Genome-wide scan reveals population stratification and footprints of recent selection in Nelore cattle
- Diercles F. Cardoso1Email author,
- Lucia Galvão de Albuquerque1, 2,
- Christian Reimer3,
- Saber Qanbari3,
- Malena Erbe3, 4,
- André V. do Nascimento1,
- Guilherme C. Venturini1,
- Daiane C. Becker Scalez1,
- Fernando Baldi1, 2,
- Gregório M. Ferreira de Camargo1,
- Maria E. Zerlotti Mercadante2, 5,
- Joslaine N. do Santos Gonçalves Cyrillo5,
- Henner Simianer3 and
- Humberto Tonhati1, 2
© The Author(s) 2018
Received: 30 May 2017
Accepted: 20 February 2018
Published: 2 May 2018
This study aimed at (1) assessing the genomic stratification of experimental lines of Nelore cattle that have experienced different selection regimes for growth traits, and (2) identifying genomic regions that have undergone recent selection. We used a sample of 763 animals genotyped with the Illumina BovineHD BeadChip, among which 674 animals originated from two lines that are maintained under directional selection for increased yearling body weight and 89 animals from a control line that is maintained under stabilizing selection.
Multidimensional analysis of the genomic dissimilarity matrix and admixture analysis revealed a substantial level of population stratification between the directional selection lines and the stabilizing selection control line. Two of the three tests used to detect selection signatures (FST, XP-EHH and iHS) revealed six candidate regions with indications of selection, which strongly indicates truly positive signals. The set of identified candidate genes included several genes with roles that are functionally related to growth metabolism, such as COL14A1, CPT1C, CRH, TBC1D1, and XKR4.
The current study identified genetic stratification that resulted from almost four decades of divergent selection in an experimental Nelore population, and highlighted autosomal genomic regions that present patterns of recent selection. Our findings provide a basis for a better understanding of the metabolic mechanism that underlies the growth traits, which are modified by selection for yearling body weight.
Human/animal interactions in livestock production have led to major phenotypic changes in animal traits, which are mainly due to the effect of selection . Anthropogenic selection, along with natural selection and adaptation, have increased the frequency of favorable alleles and resulted in changes in the frequencies of their surrounding loci , thereby establishing genomic signatures of positive selection. In turn, the methods used to identify this type of signature of selection represent powerful tools to detect evidence of selection during the process of animal domestication or breed formation, and ultimately they aid in elucidating the mechanisms that underlie the morphological, productive and functional traits of a species .
Cattle are domestic animals used worldwide, with two divergent groups: taurine (Bos primigenious taurus) and indicine (Bos primigenious indicus) cattle, each of which includes several breeds, which have been established for different purposes by humans. One common strategy for investigating past selection in cattle involves the systematic comparison of the genomes of breeds that have evolved under different selection regimes, e.g., dairy versus beef breeds [4, 5]. The motivation behind such comparisons is to identify strongly differentiated genomic regions between different groups of animals. However, analysis of within-breed stratification may allow the identification of signatures of selection within a single breed for which genetically distinct lines have been established. One example of such a study in cattle involved Holstein lines that were divergently selected for milk yield, in which the identified signatures of selection co-localized with known quantitative trait loci (QTL) that affect milk traits [6, 7].
Nelore cattle (originally named Ongole in India) are an indicine breed, which is predominantly raised in Brazil for meat production. Apart from the national Nelore breeding programs, there is an experimental program with three selection lines that has been running since 1978 at the APTA Beef Cattle Center—Institute of Animal Science (IZ), Sertãozinho—Brazil. This program was initiated to demonstrate to producers the benefits of selection for growth traits and its economic impact on the bovine industry . Briefly, three selection lines were established in 1980 by dividing an experimental herd representative of the main Nelore lineages that currently exist in Brazil. Since then, a line referred to as Nelore Control (NeC) has been maintained under stabilizing selection, in which animals with a yearling body weight (YBW) close to the average of the contemporary group are annually selected for reproduction. The two other lines, Nelore Selection (NeS) and Nelore Traditional (NeT), have independently undergone selection for increased YBW. While the NeS and NeC lines are maintained as closed lines, the NeT line has received sires from the two other lines (7 from NeC and 47 from NeT) and commercial herds.
After almost 40 years of selection, distinct phenotypic differences are observed between the lines undergoing stabilizing and directional selection, including differences in average body weight at various ages, body measurements, scrotal circumference and carcass quality [8–10]. Inbreeding level in these lines is controlled by adopting breeding strategies such as exclusion of siblings as sires, short mating time for each selected sire (only two seasons), and avoidance of mating between closely-related animals. In 2003, the observed coefficients of inbreeding in the three lines ranged from 4 to 4.3% , while in animals born in 2012 they ranged from 4.2 to 6.4%.
The existence of these experimental bovine lines has instigated studies on candidate genes, with the goal of detecting polymorphisms that show extreme differences in allele frequencies between lines, which may be linked to phenotypic differences [12–14]. However, genome-wide screening of high-density single nucleotide polymorphism (SNP) panels to detect selection is probably more suitable for revealing genomic regions that have undergone selection in this population and should help to understand the genetic mechanisms that underlie response to selection. Thus, the aims of this study were (1) to assess the effect of selection within this experimental Nelore population and (2) to scan for signatures of selection left by recent directional selection.
Description of the experimental selection program
Since the experimental program was established in 1980, the NeC, NeS and NeT lines are composed of 60, 120 and 170 cows with four, six and up to eight bulls, respectively, with 50% of the bulls at 2 years of age (first mating) and 50% at 3 years (second and last mating). The females that were initially assigned to each line were randomly chosen from an experimental herd of 350 cows, and since then, the annual culling rate of cows is approximately 20%, with replacement females for each line being selected based on their YBW at 550 days of age and raised on pasture. The bulls that were initially assigned to each line were chosen based on their YBW at 378 days of age, after 168-day feedlot performance tests. The bulls that were assigned to the NeC line had a YBW close to the average value for bulls with the same birth year (1977 or 1978), while those assigned to the NeS and NeT lines had a higher YBW. This bull selection scheme was maintained to the present day, by replacing annually the 50% older bulls within each line. Average annual selection intensities of 0.05 ± 0.24, 0.42 ± 0.29 and 0.56 ± 0.28 were applied to NeC, NeS and NeT females, respectively, while the corresponding values for males are equal to − 0.03 ± 0.20, 1.46 ± 0.30 and 1.60 ± 0.54.
Sampling of animals and quality control procedures
Average additive genetic relationships in the genotyped subsample
0.180 ± 0.072
0.007 ± 0.001
0.007 ± 0.001
0.115 ± 0.060
0.053 ± 0.035
0.091 ± 0.059
Each animal was genotyped using the Illumina® BovineHD BeadChip (San Diego, CA). Quality control (QC) procedures were performed using Plink (version 1.07) : samples with call rates lower than 0.90 were removed, and SNPs with call rates lower than 0.95, non-autosomal, unmapped, and duplicated SNPs or those that were monomorphic in all lines were also discarded. In addition, SNPs that might be misplaced were removed based on the decay of r2 as a function of physical distance. Following an algorithm adapted from Bohmanova et al. , this QC procedure consisted of: (1) recording the distance between each SNP and the top 10 syntenic SNPs with the highest r2 values; (2) flagging a particular SNP once for each distance greater than 10 Mb within the list of its top 10 r2 and compiling a list of candidate misplaced SNPs for each chromosome; (3) excluding from further analysis the SNPs with the highest flag counts on each chromosome; and (4) repeating steps (1), (2) and (3) with the list of candidate misplaced SNPs, until the list is empty (an example of the removal of a misplaced SNP is in Figure S1 [see Additional file 2: Figure S1]).
Finally, 529,172 SNPs remained for further analyses [see Additional file 1: Table S2]. Next, imputation of missing genotypes and haplotype phasing were performed by combining all the samples together, regardless of their respective line, using Beagle software (version 3.3.2) .
Population structure analyses
Two different methods were used to assess population stratification. First, a multidimensional scaling (MDS) analysis was applied using the genomic distance matrix and the “cmdscale” function from the “stats” package in R software . The genomic distance between two individuals was estimated as 1 minus the proportion of identical by state (IBS) alleles that they share. Second, an admixture analysis was performed to measure the proportion of individual ancestry from different numbers of hypothetical ancestral populations, using the ADMIXTURE software (version 1.04) . Thus, an unsupervised model-based approach that relies on maximum likelihood was applied to the genotype matrix, using hypothetical ancestral population (K) numbers of 2 and 3, which reflect the different selection regimes and the number of lines in the subsample used for this study, respectively.
Wright’s fixation index (F ST)
Wright’s fixation index (FST) was computed for each pairwise combination of lines based on the estimator corrected for unequal sample sizes  implemented in the “fsthet” package  in R software , using all SNPs that had passed the QC procedures. Then, FST and heterozygosity averages were estimated for genomic windows of 100 kb across the autosomes, with an overlap of 75 kb. Putative signals of divergence were detected with two empirical strategies that are typically used in studies involving FST tests, one that identifies the top-ranking windows [25, 26] and another method that identifies outliers from the background level of FST-heterozygosity [24, 27]. Thus, windows were considered to show putative signatures of selection if they contained more than five SNPs, presented an average FST within the top 1%, and outside the 95% confidence interval of their heterozygosity bin in both cross-line comparisons of NeC with one of the lines under directional selection (“NeC vs. NeS” and “NeC vs. NeT”) [see Additional file 2: Figure S2].
Haplotype-based methods for assigning signatures of selection
The extended haplotype homozygosity (EHH) metric is defined as the probability of randomly choosing two identical extended haplotypes that surround a specific locus , for which higher values are expected for selected alleles with a frequency that increased rapidly due to selection. In our study, we applied two somewhat complementary EHH-derived methods to detect signatures of selection, the cross population EHH (XP-EHH)  and the integrated haplotype homozygosity score (iHS) , both of which were estimated using the SELSCAN software (version 1.1.0) . Whereas XP-EHH is a cross population method that is useful to detect signals of divergent selection between lines, iHS is a within-population method that allows the identification of signals of ongoing positive selection.
The XP-EHH test was applied to the comparisons between NeC and each of the lines under directional selection (“NeC vs. NeS” and “NeC vs. NeT”). XP-EHH values were averaged for 100-kb windows with an overlap of 75 kb, and a cut-off threshold of 2.58 (p < 0.005) was applied to determine candidate signatures of selection. Only windows for which the average XP-EHH exceeded the upper cut-off in both cross-line comparisons and contained more than five SNPs were considered to harbor a signature of selection [see Additional file 2: Figure S3], which was characterized by the presence of some high-frequency haplotypes in the NeS and NeT lines but not in the NeC line.
The iHS test was performed by considering the lines under directional selection (NeS and NeT) together after removing SNPs with a minor allele frequency (MAF) lower than 0.01, since iHS is a within-population methodology that is not appropriate for analyzing SNPs close to fixation. One allele of each SNP was arbitrarily assigned as the “reference” allele, and absolute values of iHS (|iHS|) were then used to identify signatures of selection, as in Cohen-Zinder et al. . The “reference” allele was randomly determined based on the high similarity observed between iHS values obtained from a subset of SNPs, with or without ancestral allele information retrieved from two cattle reference allele lists [33, 34] [see Additional file 2: Figure S4]. The |iHS| values of 494,372 SNPs were divided into 20 bins based on the frequency of the assigned ancestral allele, and the values were then standardized within each bin. The |iHS| values were averaged across 100-kb windows, with an overlap of 75 kb. Windows with |iHS| scores higher than 2.58 (p < 0.005) were considered as putative signatures of selection.
QTL and gene annotation
Known QTL that overlapped with each detected signature of selection were retrieved from the CattleQTL database . Genes were identified within each putative selected region using the Ensembl cow gene set 84 . Enrichment of genes in specific biological processes and molecular pathways was then analyzed using the DAVID database (version 6.7) , however since no significant pathway was detected, the results are not shown.
Results and discussion
Polymorphisms and population structure
Summary statistics for SNPs in the three selection lines
Monomorphic SNPs (N)
Heterozygosity ± SD
0.26 ± 0.01
0.28 ± 0.01
0.28 ± 0.01
We performed an admixture analysis to estimate the proportions of ancestry in NeC, NeS and NeT with various numbers of clusters (K), which represent hypothetical ancestral populations. When K = 2, the clustering algorithm clearly separated NeC from NeS and NeT (Fig. 2b), and the gene flow into NeT became evident with K = 3 (Fig. 2c). Both MDS and admixture analyses identified NeT animals that showed a high degree of shared ancestry with the NeS cluster (Fig. 2a, c), which is likely due to the possible use of NeS sires for reproduction in the NeT line.
FST is a widely used metric for identifying population differentiation with values ranging from 0 to 1, where values close to 0 indicate more homogeneous populations and values close to 1 indicate more genetically diverged populations. The estimated average FST values in cross-line comparisons were equal to 0.05 ± 0.04 for “NeC vs. NeS”, 0.06 ± 0.04 for “NeC vs. NeT” and 0.02 ± 0.02 for “NeS vs. NeT”. These results reflect only subtle differentiation on a genome-wide scale, even after 40 years of intensive directional selection, which is consistent with the results of Flori et al. , who found minor genomic divergence between evolving lines of a single bovine breed. In spite of the low FST averages, their ranking was consistent with the population structure, with higher values being observed between NeC and one of the lines undergoing directional selection than between the two selection lines.
Signatures of selection
Reported QTL and genes overlapping the signatures selection that were detected by the FST, XP-EHH and iHS tests in Nelore cattle
Chromosome position (kb)
FST and iHS
Birth, weaning and yearling body weight 
FST and iHS
XP-EHH and iHS
FST and XP-EHH
Maintenance efficiency and partial efficiency of growth 
XP-EHH and iHS
Calving ease 
XP-EHH and iHS
Calving ease 
Some interesting genomic regions were detected exclusively based on one of the three metrics. For instance, the signature of selection on BTA6:58,700,000–58,925,000 that was revealed by the FST test overlaps with a previously reported signature of selection in Limousin cattle . The candidate gene located in this region, TBC1D1 (TBC1 domain family member 1), encodes a well-known body weight regulator in several species, including humans, mice, chickens and pigs [59–63], that is directly related to insulin-stimulated glucose uptake in muscle cells . Another candidate signature of selection identified via FST, BTA6:118,650,000–118,800,000, is linked to a QTL for average daily gain that was previously detected in the same population as that used here . This signature of selection could reflect the influence of the selection for YBW on the daily rate of weight gain. This region includes the PSAPL1 (prosaposin-like 1) gene and coincides with a copy number variant (CNV) region reported in Nelore cattle . Previous studies have indicated a high level of coincidence between signatures of selection and CNV regions in cattle, which suggests that CNV may have a role in recent selection [66, 67].
Three of the signatures of selection detected via XP-EHH in our study (BTA14:77,975,000–78,875,000, BTA14:80,525,000–84,675,000 and BTA16:69,100,000–70,900,000) correspond to signals of recent positive selection that were previously reported in Nelore cattle . In addition, the XKR4 (XK related 4) gene resides in the signature of selection on BTA14:24,150,000–24,300,000, which is a well-characterized gene associated with growth traits in cattle, including Nelore birth weight and meat tenderness [43, 57].
In addition to the signatures of selection that could be associated with directional selection for growth, the iHS test revealed signatures of selection that represent historical selection for attributes that have become fixed in the Nelore breed. For instance, the signature of selection on BTA23:25,450,000–25,625,000 overlaps with the Bola DR-alpha gene, which is involved in the bovine immune system and may be related to adaptation and high resistance to parasites. In addition, the signature of selection on BTA14:30,850,000–31,650,000 coincides with a QTL that was reported in the Nelore breed  and affects alignment and appropriate angles of feet and legs, which are phenotypes that have been strongly selected for in animals registered with the Brazilian breed association ever since the breed was first established in the country.
We describe a genomic scan in an indicine breed that has been selected for growth traits and is crucial for beef production in tropical environments. Our findings demonstrate that genetic stratification has occurred in this experimental Nelore population due to almost four decades of divergent selection. Three complementary approaches for the detection of signatures of selection were used to identify candidate regions on the autosomes, most of which overlapped with or were near known QTL and candidate genes that affect essential economic traits of beef cattle. The panel of identified loci includes well-known genes related to growth traits, such as TBC1D1, CRH, and XKR4, as well as novel candidates, such as COL14A1, PSAPL1 and CPT1C1. The candidate signatures of selection and the list of candidate genes presented here provide a basis to better understand the metabolic mechanisms that underlie the growth traits, which have been modified by selection for YBW. Moreover, these results provide candidate genes for further fine-mapping studies.
Conceived and designed the study: HT, FB, LGA, JNSGC and MEZM. Funding acquisition: LGA. Formal analysis: DFC. Support for data analysis: GCV, CR, DCBS, SQ, ME, GMFC, AVN. Guided the research: HT and HS. Original draft: DFC. All authors contributed to the interpretation of the results. All authors read and approved the final manuscript.
The authors are grateful to the APTA Beef Cattle Center - Institute of Animal Science (IZ) for providing the data and biological material used in this study.
The authors declare that they have no competing interests.
This study was approved by the Ethics Committee on Animal Use (CEUA) of the Faculty of Agrarian and Veterinary Sciences, Sao Paulo State University (UNESP), under Protocol Number 18340/16.
This work was supported by the Sao Paulo Research Foundation (FAPESP), with Grant 2009/16118-5 for financial support; and Grants 2012/24600-4 and 2014/13445-3 for the first author.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Lush JL. Animal breeding plans. 3rd ed. Ames: Iowa State College Press; 1945.Google Scholar
- Smith JM, Haigh J. The hitch-hiking effect of a favourable gene. Genet Res. 1974;23:23–35.View ArticlePubMedGoogle Scholar
- Qanbari S, Simianer H. Mapping signatures of positive selection in the genome of livestock. Livest Sci. 2014;166:133–43.View ArticleGoogle Scholar
- Flori L, Gonzatti MI, Thevenon S, Chantal I, Pinto J, Berthier D, et al. A quasi-exclusive European ancestry in the Senepol tropical cattle breed highlights the importance of the slick locus in tropical adaptation. PLoS One. 2012;7:e36133.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhao F, McParland S, Kearney F, Du L, Berry DP. Detection of selection signatures in dairy and beef cattle using high-density genomic information. Genet Sel Evol. 2015;47:49.View ArticlePubMedPubMed CentralGoogle Scholar
- Sonstegard TS, Ma L, Van Tassell CP, Kim ES, Cole JB, Wiggans GR, et al. Forty years of artificial selection in US Holstein cattle had genome-wide signatures. In Proceedings of the 9th World Congress on Genetics Applied to Livestock Production: 1–6 August 2010; Leipzig. 2010.Google Scholar
- Kim ES, Cole JB, Huson H, Wiggans GR, Van Tassell CP, Crooker BA, et al. Effect of artificial selection on runs of homozygosity in US Holstein cattle. PLoS One. 2013;8:e80813.View ArticlePubMedPubMed CentralGoogle Scholar
- Mercadante MEZ, Packer IU, Razook AG, Cyrillo JNSG, Figueiredo LA. Direct and correlated responses to selection for yearling weight on reproductive performance of Nelore cows. J Anim Sci. 2003;81:376–84.View ArticlePubMedGoogle Scholar
- Cyrillo JNSG, Razook AG, Figueiredo LA, Bonilha Neto LM, Ruggieri AC, Tonhati H. Efeitos da Seleção para Peso Pós-desmame sobre Medidas Corporais e Perímetro Escrotal de Machos Nelore de Sertãozinho (SP). Braz J Anim Sci. 2000;29:403–12.Google Scholar
- Bonilha SFM, Tedeschi LO, Packer IU, Razook AG, Alleoni GF, Nardon RF, et al. Evaluation of carcass characteristics of Bos indicus and tropically adapted Bos taurus breeds selected for postweaning weight. J Anim Sci. 2008;86:1770–80.View ArticlePubMedGoogle Scholar
- Mercadante MEZ, Razook AG. Programa de melhoramento genético do Instituto de Zootecnia. In: Pires AV. Bovinocultura de Corte. 1st ed. Piracicaba: FEALQ. 2010. pp. 1508.Google Scholar
- Ayres DR, Souza FRP, Mercadante MEZ, Fonseca LFS, Tonhati H, Cyrillo JNSG, et al. Evaluation of TFAM and FABP4 gene polymorphisms in three lines of Nellore cattle selected for growth. Genet Mol Res. 2010;9:2050–9.View ArticlePubMedGoogle Scholar
- Cardoso DF, de Souza FRP, de Camargo GMF, da Fonseca PD, Fonseca LFS, Braz CU, et al. Polymorphism analysis in genes of the somatotropic axis in Nellore cattle selected for growth. Gene. 2014;545:215–9.View ArticlePubMedGoogle Scholar
- de Souza FRP, Chiquitelli MG, da Fonseca LFS, Cardoso DF, da Silva Fonseca PD, de Camargo GMF, et al. Associations of FASN gene polymorphisms with economical traits in Nellore cattle (Bos primigenius indicus). Mol Biol Rep. 2012;39:10097–104.View ArticlePubMedGoogle Scholar
- Brinks JS, Rice FJ, Clark RT. Estimation of genetic trends in beef cattle. J Anim Sci. 1961;20:903.View ArticleGoogle Scholar
- Olivieri BF, Mercadante MEZ, Cyrillo JN, Branco RH, Bonilha SFM, de Albuquerque LG, et al. Genomic regions associated with feed efficiency indicator traits in an experimental Nellore cattle. PLoS One. 2016;11:e0164390.View ArticlePubMedPubMed CentralGoogle Scholar
- Silva RMO, Fragomeni BO, Lourenco DAL, Magalhães AFB, Irano N, Carvalheiro R, et al. Accuracies of genomic prediction of feed efficiency traits using different prediction and validation methods in an experimental Nelore cattle population. J Anim Sci. 2016;94:3613–23.View ArticlePubMedGoogle Scholar
- Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–75.View ArticlePubMedPubMed CentralGoogle Scholar
- Bohmanova J, Sargolzaei M, Schenkel FS. Characteristics of linkage disequilibrium in North American Holsteins. BMC Genomics. 2010;11:421.View ArticlePubMedPubMed CentralGoogle Scholar
- Browning SR, Browning BL. Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. Am J Hum Genet. 2007;81:1084–97.View ArticlePubMedPubMed CentralGoogle Scholar
- R Core Team R. A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing. http://www.r-project.org/ (2013). Accessed 1 Feb 2018.
- Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009;19:1655–64.View ArticlePubMedPubMed CentralGoogle Scholar
- Cockerham CC, Weir BS. Estimation of gene flow from F-statistics. Evolution. 1993;47:855–63.PubMedGoogle Scholar
- Flanagan SP, Jones AG. Constraints on the FST–heterozygosity outlier approach. J Hered. 2017;108:561–73.View ArticlePubMedGoogle Scholar
- Gholami M, Erbe M, Gärke C, Preisinger R, Weigend A, Weigend S, et al. Population genomic analyses based on 1 million SNPs in commercial egg layers. PLoS One. 2014;9:e94509.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang B, Cui L, Perez-Enciso M, Traspov A, Crooijmans RPMA, Zinovieva N, et al. Genome-wide SNP data unveils the globalization of domesticated pigs. Genet Sel Evol. 2017;49:71.View ArticlePubMedPubMed CentralGoogle Scholar
- Beaumont MA, Nichols RA. Evaluating loci for use in the genetic analysis of population structure. Proc R Soc Lond B Biol Sci. 1996;263:1619–26.View ArticleGoogle Scholar
- Sabeti PC, Reich DE, Higgins JM, Levine HZP, Richter DJ, Schaffner SF, et al. Detecting recent positive selection in the human genome from haplotype structure. Nature. 2002;419:832–7.View ArticlePubMedGoogle Scholar
- Sabeti PC, Varilly P, Fry B, Lohmueller J, Hostetter E, Cotsapas C, et al. Genome-wide detection and characterization of positive selection in human populations. Nature. 2007;449:913–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Voight BF, Kudaravalli S, Wen X, Pritchard JK. A map of recent positive selection in the human genome. PLoS Biol. 2006;4:e72.View ArticlePubMedPubMed CentralGoogle Scholar
- Szpiech ZA, Hernandez RD. selscan: an efficient multithreaded program to perform EHH-based scans for positive selection. Mol Biol Evol. 2014;31:2824–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Cohen-Zinder M, Donthu R, Larkin DM, Kumar CG, Rodriguez-Zas SL, Andropolis KE, et al. Multisite haplotype on cattle chromosome 3 is associated with quantitative trait locus effects on lactation traits. Physiol Genomics. 2011;43:1185–97.View ArticlePubMedGoogle Scholar
- Rocha D, Billerey C, Samson F, Boichard D, Boussaha M. Identification of the putative ancestral allele of bovine single-nucleotide polymorphisms. J Anim Breed Genet. 2014;131:483–6.View ArticlePubMedGoogle Scholar
- Utsunomiya YT, Pérez O’Brien AM, Sonstegard TS, Van Tassell CP, do Carmo AS, Mészáros G, et al. Detecting loci under recent positive selection in dairy and beef cattle by combining different genome-wide scan methods. PLoS One. 2013;8:e0064280.View ArticleGoogle Scholar
- Hu ZL, Park CA, Reecy JM. Developmental progress and current status of the Animal QTLdb. Nucleic Acids Res. 2015;44:D827–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Yates A, Akanni W, Amode MR, Barrell D, Billis K, Carvalho-Silva D, et al. Ensembl 2016. Nucleic Acids Res. 2016;44:D710–6.View ArticlePubMedGoogle Scholar
- Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.View ArticleGoogle Scholar
- Gutiérrez JP, Cervantes I, Goyache F. Improving the estimation of realized effective population sizes in farm animals. J Anim Breed Genet. 2009;126:327–32.View ArticlePubMedGoogle Scholar
- Vitti JJ, Grossman SR, Sabeti PC. Detecting natural selection in genomic data. Annu Rev Genet. 2013;47:97–120.View ArticlePubMedGoogle Scholar
- Evans LM, Slavov GT, Rodgers-Melnick E, Martin J, Ranjan P, Muchero W, et al. Population genomics of Populus trichocarpa identifies signatures of selection and adaptive trait associations. Nat Genet. 2014;46:1089–96.View ArticlePubMedGoogle Scholar
- He Q, Kim KW, Park YJ. Population genomics identifies the origin and signatures of selection of Korean weedy rice. Plant Biotechnol J. 2017;15:357–66.View ArticlePubMedGoogle Scholar
- Saatchi M, Schnabel RD, Taylor JF, Garrick DJ. Large-effect pleiotropic or closely linked QTL segregate within and across ten US cattle breeds. BMC Genomics. 2014;15:442.View ArticlePubMedPubMed CentralGoogle Scholar
- Magalhães AFB, de Camargo GMF, Fernandes GA, Gordo DGM, Tonussi RL, Costa RB, et al. Genome-wide association study of meat quality traits in Nellore cattle. PLoS One. 2016;11:e0157845.View ArticlePubMedPubMed CentralGoogle Scholar
- Fortes MRS, Reverter A, Hawken RJ, Bolormaa S, Lehnert SA. Candidate genes associated with testicular development, sperm quality, and hormone levels of inhibin, luteinizing hormone, and insulin-like growth factor 1 in Brahman bulls. Biol Reprod. 2012;87:58.View ArticlePubMedGoogle Scholar
- Ramayo-Caldas Y, Fortes MRS, Hudson NJ, Porto-Neto LR, Bolormaa S, Barendse W, et al. A marker-derived gene network reveals the regulatory role of PPARGC1A, HNF4G, and FOXP3 in intramuscular fat deposition of beef cattle. J Anim Sci. 2014;92:2832–45.View ArticlePubMedGoogle Scholar
- Veneroni-Gouveia G, Meirelles SL, Grossi DA, Santiago AC, Sonstegard TS, Yamagishi MEB, et al. Whole-genome analysis for backfat thickness in a tropically adapted, composite cattle breed from Brazil. Anim Genet. 2012;43:518–24.View ArticlePubMedGoogle Scholar
- Lee SH, Choi BH, Lim D, Gondro C, Cho YM, Dang CG, et al. Genome-wide association study identifies major loci for carcass weight on BTA14 in Hanwoo (Korean cattle). PLoS One. 2013;8:e74677.View ArticlePubMedPubMed CentralGoogle Scholar
- McClure MC, Ramey HR, Rolf MM, McKay SD, Decker JE, Chapple RH, et al. Genome-wide association analysis for quantitative trait loci influencing Warner-Bratzler shear force in five taurine cattle breeds. Anim Genet. 2012;43:662–73.View ArticlePubMedPubMed CentralGoogle Scholar
- de Oliveira PSN, Cesar ASM, do Nascimento ML, Chaves AS, Tizioto PC, Tullio RR, et al. Identification of genomic regions associated with feed efficiency in Nelore cattle. BMC Genet. 2014;15:100.View ArticlePubMedPubMed CentralGoogle Scholar
- Tono-Oka S, Tanase S, Miike T, Tanaka H. Transient expression of collagen type XIV during muscle development and its reappearance after denervation and degeneration. J Histochem Cytochem. 1996;44:907–18.View ArticlePubMedGoogle Scholar
- Buchanan FC, Thue TD, Yu P, Winkelman-Sim DC. Single nucleotide polymorphisms in the corticotrophin-releasing hormone and pro-opiomelancortin genes are associated with growth and carcass yield in beef cattle. Anim Genet. 2005;36:127–31.View ArticlePubMedGoogle Scholar
- Sevane N, Armstrong E, Wiener P, Pong-Wong R, Dunner S. GemQual Consortium. Polymorphisms in twelve candidate genes are associated with growth, muscle lipid profile and meat quality traits in eleven European cattle breeds. Mol Biol Rep. 2014;41:4721–31.View ArticlePubMedGoogle Scholar
- Kemper KE, Saxton SJ, Bolormaa S, Hayes BJ, Goddard ME. Selection for complex traits leaves little or no classic signatures of selection. BMC Genomics. 2014;15:246.View ArticlePubMedPubMed CentralGoogle Scholar
- Urbinati I, Stafuzza NB, Oliveira MT, Chud TCS, Higa RH, Regitano LC, et al. Selection signatures in Canchim beef cattle. J Anim Sci Biotechnol. 2016;7:29.View ArticlePubMedPubMed CentralGoogle Scholar
- Karim L, Takeda H, Lin L, Druet T, Arias JAC, Baurain D, et al. Variants modulating the expression of a chromosome domain encompassing PLAG1 influence bovine stature. Nat Genet. 2011;43:405–13.View ArticlePubMedGoogle Scholar
- Fernandes Júnior GA, Costa RB, de Camargo GMF, Carvalheiro R, Rosa GJM, Baldi F, et al. Genome scan for postmortem carcass traits in Nellore cattle. J Anim Sci. 2016;94:4087–95.View ArticleGoogle Scholar
- Utsunomiya YT, do Carmo AS, Carvalheiro R, Neves HH, Matos MC, Zavarez LB, et al. Genome-wide association study for birth weight in Nellore cattle points to previously described orthologous genes affecting human and bovine height. BMC Genet. 2013;14:52.View ArticlePubMedPubMed CentralGoogle Scholar
- Costa RB, Camargo GMF, Diaz IDPS, Irano N, Dias MM, Carvalheiro R, et al. Genome-wide association study of reproductive traits in Nellore heifers using Bayesian inference. Genet Sel Evol. 2015;47:67.View ArticlePubMedPubMed CentralGoogle Scholar
- Stone S, Abkevich V, Russell DL, Riley R, Timms K, Tran T, et al. TBC1D1 is a candidate for a severe obesity gene and evidence for a gene/gene interaction in obesity predisposition. Hum Mol Genet. 2006;15:2709–20.View ArticlePubMedGoogle Scholar
- Meyre D, Farge M, Lecoeur C, Proenca C, Durand E, Allegaert F, et al. R125 W coding variant in TBC1D1 confers risk for familial obesity and contributes to linkage on chromosome 4p14 in the French population. Hum Mol Genet. 2008;17:1798–802.View ArticlePubMedGoogle Scholar
- Chadt A, Leicht K, Deshmukh A, Jiang LQ, Scherneck S, Bernhardt U, et al. Tbc1d1 mutation in lean mouse strain confers leanness and protects from diet-induced obesity. Nat Genet. 2008;40:1354–9.View ArticlePubMedGoogle Scholar
- Wang Y, Xu HY, Gilbert ER, Peng X, Zhao XL, Liu YP, et al. Detection of SNPs in the TBC1D1 gene and their association with carcass traits in chicken. Gene. 2014;547:288–94.View ArticlePubMedGoogle Scholar
- Fontanesi L, Galimberti G, Calò DG, Fronza R, Martelli PL, Scotti E, et al. Identification and association analysis of several hundred single nucleotide polymorphisms within candidate genes for back fat thickness in Italian Large White pigs using a selective genotyping approach. J Anim Sci. 2012;90:2450–64.View ArticlePubMedGoogle Scholar
- Dokas J, Chadt A, Nolden T, Himmelbauer H, Zierath JR, Joost HG, et al. Conventional knockout of Tbc1d1 in mice impairs insulin- and AICAR-stimulated glucose uptake in skeletal muscle. Endocrinology. 2013;154:3502–14.View ArticlePubMedGoogle Scholar
- da Silva JM, Giachetto PF, da Silva LO, Cintra LC, Paiva SR, Yamagishi MEB, et al. Genome-wide copy number variation (CNV) detection in Nelore cattle reveals highly frequent variants in genome regions harboring QTLs affecting production traits. BMC Genomics. 2016;17:454.View ArticlePubMedPubMed CentralGoogle Scholar
- Pérez O’Brien AM, Utsunomiya YT, Mészáros G, Bickhart DM, Liu GE, Van Tassell CP, et al. Assessing signatures of selection through variation in linkage disequilibrium between taurine and indicine cattle. Genet Sel Evol. 2014;46:19.View ArticlePubMedPubMed CentralGoogle Scholar
- Teo YY, Fry AE, Bhattacharya K, Small KS, Kwiatkowski DP, Clark TG. Genome-wide comparisons of variation in linkage disequilibrium. Genome Res. 2009;19:1849–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Somavilla AL, Sonstegard TS, Higa RH, Rosa AN, Siqueira F, Silva LOC, et al. A genome-wide scan for selection signatures in Nellore cattle. Anim Genet. 2014;45:771–81.View ArticlePubMedGoogle Scholar
- Vargas G, Neves HHR, Camargo GMF, Cardoso V, Munari DP, Carvalheiro R. Genome-wide association study and functional analysis of feet and leg conformation traits in Nellore cattle. J Anim Sci. 2018 (accept for publication).Google Scholar