Genome-wide association study for birth weight in Nellore cattle points to previously described orthologous genes affecting human and bovine height
- Yuri T Utsunomiya†1,
- Adriana S do Carmo†1,
- Roberto Carvalheiro2,
- Haroldo HR Neves1,
- Márcia C Matos1,
- Ludmilla B Zavarez1,
- Ana M Pérez O’Brien3,
- Johann Sölkner3,
- John C McEwan4,
- John B Cole5,
- Curtis P Van Tassell6,
- Flávio S Schenkel7,
- Marcos VGB da Silva8,
- Laercio R Porto Neto9, 10,
- Tad S Sonstegard6 and
- José F Garcia11Email author
© Utsunomiya et al.; licensee BioMed Central Ltd. 2013
Received: 19 December 2012
Accepted: 5 June 2013
Published: 13 June 2013
Birth weight (BW) is an economically important trait in beef cattle, and is associated with growth- and stature-related traits and calving difficulty. One region of the cattle genome, located on Bos primigenius taurus chromosome 14 (BTA14), has been previously shown to be associated with stature by multiple independent studies, and contains orthologous genes affecting human height. A genome-wide association study (GWAS) for BW in Brazilian Nellore cattle (Bos primigenius indicus) was performed using estimated breeding values (EBVs) of 654 progeny-tested bulls genotyped for over 777,000 single nucleotide polymorphisms (SNPs).
The most significant SNP (rs133012258, PGC = 1.34 × 10-9), located at BTA14:25376827, explained 4.62% of the variance in BW EBVs. The surrounding 1 Mb region presented high identity with human, pig and mouse autosomes 8, 4 and 4, respectively, and contains the orthologous height genes PLAG1, CHCHD7, MOS, RPS20, LYN, RDHE2 (SDR16C5) and PENK. The region also overlapped 28 quantitative trait loci (QTLs) previously reported in literature by linkage mapping studies in cattle, including QTLs for birth weight, mature height, carcass weight, stature, pre-weaning average daily gain, calving ease, and gestation length.
This study presents the first GWAS applying a high-density SNP panel to identify putative chromosome regions affecting birth weight in Nellore cattle. These results suggest that the QTLs on BTA14 associated with body size in taurine cattle (Bos primigenius taurus) also affect birth weight and size in zebu cattle (Bos primigenius indicus).
KeywordsGWAS Birth weight Bos primigenius indicus Nellore cattle Stature
Birth weight (BW) is an economically important trait in beef cattle, and is usually the first characteristic measured in a calf. Birth weight is associated with growth-related traits , mature size  and carcass weight, thus being a valuable production indicator, as well as a selection criterion to improve calving ease [3, 4].
Despite the beef industry’s pursuit of animals with rapid growth, yielding heavier carcasses, selection for these objectives needs to be properly balanced against selection for reproductive traits, which have great economic importance in beef cattle production systems [5, 6]. While low estimated breeding values (EBVs) for BW are associated with reduced calf viability  and lower growth rates [3, 7], the use of sires with high EBVs for BW on dams with small pelvic size may result in higher rates of dystocia  and increased perinatal mortality . These antagonisms result from the strong association of birth weight with the body size of the calf, i.e., with the stature of the animal . Calving difficulties can result from a mismatch between pelvic opening and calf size . This relationship between BW with reproductive and growth/size traits highlights the importance of understanding the underlying genetic architecture of BW.
Birth weight exhibits sufficient variability and heritability in the Nellore breed (Bos primigenius indicus), with an average of 29.8 ± 2.7 kg  and estimated heritability between 0.25 and 0.33 [1, 9, 10]. Despite the low frequency of dystocia in Nellore cows, BW has been recorded and used to monitor genetic trend. One selection strategy of the breeding programs in Brazil has been to preferentially use sires with higher EBVs for weaning and yearling weights, but with low or close to average EBVs for BW . The identification of major genes and variants affecting multiple weight and carcass traits or influencing BW alone would be of help to balance these conflicting goals, because BW is positively correlated with weaning and yearling weights .
In the past two decades, linkage studies attempting to map quantitative trait loci (QTLs) affecting weight, growth, or stature in cattle have been published (e.g. [12–16]). The release of the reference bovine genome , the discovery of common single nucleotide polymorphisms (SNPs) across breeds [18, 19], and the availability of high-throughput microarrays have enhanced the process of mapping loci that affect complex traits. This has led to several population-based investigations of associations between weight/growth/height phenotypes with genome-wide variants in different cattle breeds [15, 20–23].
In particular, two regions of the bovine genome associated with stature and growth have been highlighted recently. The first, located on Bos primigenius taurus (BTA) autosome 6 [20, 23], shelters the orthologous genes NCAPG and LCORL, which have been also found to be associated with adult height in humans [24, 25]. The second, located on BTA14 [15, 21–23], contains the genes PLAG1, CHCHD7, RDHE2, MOS, RPS20, LYN, PENK and TGS1, that were previously found to affect stature in both cattle and humans [22, 24, 26–28]. Importantly, the majority of the genome-wide association studies (GWAS) reported in literature were conducted in the humpless subspecies of cattle (Bos primigenius taurus, known as taurine cattle), and GWAS in the humped bovine subspecies (Bos primigenius indicus, often referred as indicine or zebu cattle) are only now emerging, especially because the first SNP microarrays were optimized for taurine cattle .
In this paper, results from a genome-wide scan for SNPs associated with BW variation in Nellore cattle using EBVs of progeny-tested Brazilian bulls are reported. As EBVs take into account information from performance of the individual, progeny and parents, pedigree relationships, and systematic management and environmental factors, they can be used as composite phenotypes for proceeding with association analyses (see, e.g., ). The objective of this study was to identify putative SNP associated with differences in BW and to explore the genomic regions around them to unravel prospective functional relationships among weight, fertility, and growth/size traits.
Genotype informativeness and quality control
From the initial set of 777,961 SNPs, 42,669 (5.5%) were non-autosomal markers. Fifty four autosomal SNPs with redundant genomic coordinates were identified and excluded from further analyses. The identity by state (IBS) check revealed no unexpected sample duplicates. A total of 223,309 (30.4%) markers were excluded due to minor allele frequency (MAF) < 0.02. The number of SNPs excluded due to SNP Call Rate (CRSNP) < 0.98 and Fisher’s exact test P-value for Hardy-Weinberg Equilibrium (HWE) < 1 × 10-5 were 122,611 (16.7%) and 13,194 (1.8%), respectively. Five individuals were removed due to low Call Rate (CRIND). The final dataset included data for 649 individuals and 434,020 SNPs.
Descriptive statistics of dependent variables
Based on the results for the Shapiro-Wilk test, there was no evidence that EBVs deviated from normality (P = 0.415), and no outliers were observed. Average accuracy was 0.87 ± 0.11, with minimum, median and maximum accuracies of 0.51, 0.91 and 0.99, respectively. After fitting the regression in formula (2), there was no evidence against the hypotheses of normally distributed (P = 0.57) and homoskedastic residuals (studentized Breusch-Pagan test, P = 0.11). These findings suggest that the dependent variables used were reliable and did not violate possible assumptions of the statistical analyses used hereafter.
Summary of parameters and statistics estimated for the identified significant SNPs
Ensembl variant ID
Illumina probe ID
BTA14 Position (bp)
T 2 c
List of genes within the 1 Mb region surrounding the most significant SNP (rs133012258)
Distance from SNP (kb)
U6 spliceosomal RNA
Short-chain dehydrogenase/reductase family 16C member 6
Epidermal retinol dehydrogenase 2
Coiled-coil-helix-coiled-coil-helix domain-containing protein 7
Pleiomorphic adenoma gene 1
V-mosMoloneymurine sarcoma viral oncogenehomolog
U1 spliceosomal RNA
40S ribosomalprotein S20
Small nucleolar RNA U54
QTLdb hits within the 1 Mb region surrounding the most significant SNP (rs133012258)
Body weight (birth)
Pre-weaning average daily gain
Calving ease (maternal)
Longissimus muscle area
Fat thickness at the 12th rib
Abnormal flavor intensity
Milk fat yield (or percentage)
Milk protein yield (or percentage)
Somatic cell score
Five SNPs on BTA14 were identified as associated with BW in Nellore cattle (P < 1.15 × 10-7), whose surrounding region has been shown to contain many QTLs, genes and variants affecting stature-related traits in cattle by several independent studies [12–15, 21–23]. More particularly, the genes PLAG1, CHCHD7, RDHE2, MOS, RPS20, LYN and PENK have been found to influence both human and cattle height [22–28].
The BTA14 region pointed out by the present study has also been shown to be associated with reproductive traits. Cole et al.  reported a QTL on BTA14 associated with stillbirth, which also has been associated with body size in dairy cattle [8, 30], but found no effect on stature or other conformation traits on that chromosome. The region also associates with many fertility and growth-related traits in the indicine breed Brahman, for example scrotal circumference [31, 32], age at the first corpus luteum[31, 33], blood levels of insulin-like growth factor 1 (IGF1) [32, 33] and hip height .
A significant SNP was found within intron 2 of the XKR4 gene in the present study. Lindholm-Perry et al.  identified five SNPs near XKR4 associated with feed intake and gain in crossbred steers. Bolormaa et al.  found five SNPs in a narrow region of BTA14 encompassing XKR4 associated with rump fat thickness measured at the P8 position (CHILLP8) in seven breeds of cattle, including taurine, indicine and composite breeds. The authors found that four of these SNPs were also associated with CHILLP8 in a confirmatory sample of 1,338 animals, including Angus, Hereford and Brahman cattle. Furthermore, Porto Neto et al.  performed a replication study using samples of Belmont Red, Santa Gertrudis and Brahman animals genotyped for SNPs within XKR4 and found that although the SNP effect may vary depending on the breed, the variant rs42646708 (BTA14:24573257) explain around 1.3% of CHILLP8 variance in cattle. This SNP is also located within intron 2 of XKR4, only 17.6 kb apart from the intronic SNP detected in the present study, which strongly suggests XKR4 as a candidate gene for being further explored in future studies of weight and carcass traits in Nellore cattle.
The most significant SNP (rs133012258, PGC = 1.34 × 10-9) was found to explain 4.62% of the variance in sires EBVs, with a 95% CI of 2.12-8.09%. One hundred and eighty loci associated with human adult height explain only 10% of the phenotypic variance together, while individual loci account for 0.4% or less . SNPs analyzed by  within a nearby BTA14 region explain from 0.29 to 2.53% of the bovine stature variability, and the quantitative trait nucleotides (QTN) spanning MOS, CHCHD7 and PLAG1 described by  explain from 1.10 to 3.50% of height in Jersey and Holstein breeds. Furthermore, the genome-wide survey performed by  provided strong evidence for two QTL on BTA14 and BTA21 that together explain at least 10% of the variation of EBVs for calving ease in the German Fleckvieh.
Considering that multiple stature-related traits are governed by variants with small effects, and that the genomic region identified in this study has been previously found to be associated with several of these traits, the putative SNP detected in the present analysis can be considered as a marker in linkage disequilibrium (LD) with major untyped (i.e., not probed by the SNP assay used) causative variants affecting BW and other height-associated traits in Nellore cattle, and further studies would be needed to determine if the QTNs reported by  are also segregating in the Nellore population. Also, future investigations are needed to better characterize the effect of nearby SNPs on other weight and carcass traits in Nellore cattle, as it is not clear yet how the putative pleiotropic effect of these variants would be used towards balancing conflicting selection goals for birth, weaning and yearling weights. Although we cannot confirm that the allele substitution effects of these SNPs work in the same direction for all three traits, because only birth weight was analyzed here, these findings suggest that the SNPs identified would be key polymorphisms to be monitored over time. In a scenario where the SNP effects have the same direction in all three traits, one strategy could be avoiding strong positive selection or drifting of the allele that contributes to higher BW EBVs, and identify and promote positive selection of other variants that have effects on weaning and yearling weights only.
The high identity found in the alignment of this BTA14 region against other mammalian species genomes suggests that these orthologous genes are located in a conserved syntenic block which may have arisen and been maintained after speciation from a common ancestor of the mammal clade. Moreover, the evidence for variants associated with growth and stature within this BTA14 region in both taurine and zebu cattle raises two hypotheses: 1) these variants have been introgressed into Nellore via historical admixture with taurine Creole cattle in the maternal line, and was maintained in the breed in spite of several generations of backcrossing; 2) these are ancient polymorphisms, probably already segregating in the founder population of wild Aurochs (Bos primigenius) before subspecies formation.
Regarding functional meaning, the set of genes reported participate in diverse growth and tumor development mechanisms. Among these genes, PLAG1 is the most appealing functional candidate. It is an oncogene that encodes a transcription factor broadly expressed during fetal development, but is down-regulated at birth . It interacts with several growth factors controlling body size, including IGF2. In addition, PLAG1 knock-out mice have been shown to have marked growth retardation and reduced fertility . In a replication study,  confirmed the findings reported by , demonstrating association of growth rate and early life and peripubertal body weight with PLAG1 polymorphisms, supporting its status as a key regulator of mammalian growth.
The lack of significant association between BW and SNPs within other previously described weight- and height-related chromosome regions in the present study should not be interpreted as a lack of existence of true association, but rather it might be due to limitations specific to this study. Firstly, because complex trait mapping requires large sample sizes and only 649 bulls were analyzed here. Secondly, the significance level adopted was highly stringent, which may have caused inflation of type II errors. In spite of these limitations, it was possible to demonstrate that a well-characterized chromosome region affecting human and taurine cattle stature also associates with BW in a zebu breed. The release of a Bos primigenius indicus reference genome assembly, as well as the application of re-sequencing and replication studies would help improve resolution to narrow down the genomic region as close as possible to the true causative variants.
This study is believed to be the first genome-wide association study applying a high-density SNP panel to identify putative chromosome regions affecting birth weight in zebu cattle. The findings presented, which are strongly supported by the literature, point to orthologous genes already known to affect growth- and stature-related traits in both humans and cattle, which may shelter ancient polymorphisms responsible for variation in those traits since before cattle subspecies divergence.
Estimated breeding values
Estimated breeding values for BW were obtained from routine genetic evaluations using performance and pedigree data from the Aliança database , containing data from different commercial Nellore breeding programs, including more than 250 farms distributed across Brazil and Paraguay. The genetic evaluation for BW was calculated using a subset of that data that included 542,918 animals, born from 1985 to 2011, and distributed in approximately 5,000 distinct contemporary groups. These data were collected in 243 grazing-based herds in Brazil. Estimated breeding values were obtained using an animal model that included fixed effects for the age of dam at calving and contemporary group (defined as animals from the same herd, born in the same year and season, and belonging to the same management group at birth, and sex), as well as random effects that include direct additive genetic, maternal additive genetic, maternal permanent environmental and residual error effects. The variance ratios required to solve the mixed model equations were computed based on restricted maximum likelihood (REML) estimates of the variance components from previous studies in this population. Only EBVs of progeny-tested bulls whose accuracy (i.e., square root of reliability, calculated based on prediction error variance estimates) was ≥ 0.50 were used for sample collection and genotyping (described later). The majority of the bulls were used under artificial insemination service.
Genotyping, informativeness, and quality assurance
A total of 654 progeny-tested Nellore bulls were genotyped with the Illumina® BovineHD Genotyping BeadChip assay, according to the manufacturer’s protocol. Genotype calls (i.e. successfully determined genotypes) were defined as genotypes with GenCall Scores greater than 0.70, using the validated standard cluster file provided by the manufacturer. As chromosomes X, Y and mtDNA present different mode of inheritance from the rest of the genome, only autosomal markers with unique genomic coordinates were included into the analyses. After this initial screening, potential duplicated samples were determined by calculating the proportion of alleles identical by state (IBS) shared between all pairs of individuals. Any pair of samples with IBS ≥ 0.95 for 2,000 randomly sampled markers was considered unexpected duplicates, and resulted in the exclusion of both members of the pair. Individual SNPs were removed from the dataset if they did not exhibit: 1) minor allele frequency (MAF) greater than or equal to 0.02, 2) Fisher’s exact test P-value for Hardy-Weinberg Equilibrium (HWE) greater than or equal to 1 × 10-5 (i.e. extremely deviating from HWE, suggesting potential genotyping error) or 3) Call rate (CRSNP) of at least 98%. After the SNP pruning, individuals exhibiting call rate (CRIND) below 90% were also removed. These procedures and many others described later were performed in the R v2.15.0 environment , using combinations of functions from the R base, locally developed scripts, and the GenABEL v1.7-2 package .
Assessment of population substructure
Where is the estimated genomic kinship between individuals i and j, L is the total number of loci used for the calculation, p l is the reference allele frequency for locus l, and gl,i and gl,j are the locus l genotypes for individuals i and j, respectively (coded as 0, 1 or 2 reference alleles). Calculations were based on 10,000 randomly sampled markers using GenABEL . The calculated genomic kinship coefficients within the yielded n x n symmetric matrix (where n is the total number of samples) were then transformed to squared Euclidean distances, and the dissimilarities between the subjects within the matrix were captured in n-1 dimensional spaces of n observations (eigenvectors), via classical multidimensional scaling .
A clustering analysis was applied to the two eigenvectors that explained the largest proportion of the data variance using the k-means algorithm  implemented in R. Individuals were clustered into 2 groups, and the association between the prior information on breeding program and the k-means clustering results was tested using Pearson’s χ2 with Yates’ continuity correction in order to see if the algorithm could reproduce the known breeding programs subgroups. Additionally, an F-test for homogeneity of variance between subpopulations and a t-test for difference between subpopulation means, defining subpopulations either as k-means assignments or breeding program of origin, was performed to determine if there was confounding due to stratification.
The EBVs were again weighted by their respective accuracies. The conservative Bonferroni adjustment for multiple testing (α=0.05/N, where N is the number of tests, i.e., number of SNPs) was used to reject the null hypothesis (β g = 0, i.e. there is no association between the SNP and the EBVs), which resulted in an adjusted significance of α = 1.15 × 10-7.
Exploratory view of significant SNPs
Where p and q are the allele frequencies and S 2 is the sample EBVs variance. The upper and lower limits of the estimated 95% CI for were used to derive a 95% CI for % .
The genomic region containing the most significant SNP of a peak was explored by inspecting a 1 Mb window around the location of this SNP using the BioMart tool and the Ensembl genes 69 database  to interrogate 500 kb to each side of the marker using the UMD v3.1 assembly. The cattle QTLdb database  was also examined to find out if the significant SNP mapped against any previously described bovine QTL. Additionally, the alignments of the UMD v3.1 assembly sequence of the 1 Mb window against the human (Homo sapiens, GRCh37 assembly), pig (Sus scrofa, Sscrofa10.2 assembly) and mouse (Mus musculus, GRCm38 assembly) genome builds were inspected in the Ensembl Comparative genomics alignments and Comparative genomics synteny tools in order to determine if any homologous genes were present in the putative region.
Estimated breeding value
Single nucleotide polymorphism
Bos primigenius taurus
Genome-wide association study
Identity by state
Minor allele frequency
Principal coordinates analysis
- Q-Q plot:
Quantitative trait locus
Rump fat thickness measured at the P8 position
Quantitative trait nucleotide
We thank Guilherme Penteado Coelho Filho and Daniel Biluca for technical assistance in sample acquisition. This research was supported by: Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) - process 560922/2010-8 and 483590/2010-0; Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) - process 2011/16643-2 and 2010/52030-2; Next-Generation BioGreen 21 Program (No. PJ008196), Rural Development Administration, Republic of Korea.
Mention of trade name proprietary product or specified equipment in this article is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the authors or their respective institutions.
- Boligon AA, Albuquerque LG, Mercadante MEZ, Lôbo RB: Herdabilidades e correlações entre pesos do nascimento à idade adulta em rebanhos da raça Nelore. Rev. Bras. Zootec. 2009, 38 (12): 2320-2326. 10.1590/S1516-35982009001200005.View ArticleGoogle Scholar
- Meyer K: Estimates of genetic parameters for mature weight of Australian beef cows and its relationship to early growth and skeletal measures. Livest Prod Sci. 1995, 44: 125-137. 10.1016/0301-6226(95)00067-4.View ArticleGoogle Scholar
- Bourdon RM, Brinks JS: Genetic, environmental and phenotypic relationships among gestation length, birth weight, growth traits and age at first calving in beef cattle. J Anim Sci. 1982, 55 (3): 543-553.PubMedGoogle Scholar
- Eriksson S, Näsholm A, Johansson K, Philipsson J: Genetic parameters for calving difficulty, stillbirth, and birth weight for Hereford and Charolais at first and later parities. J Anim Sci. 2004, 82: 375-383.PubMedGoogle Scholar
- Phocas F, Bloch C, Chapelle P, Bécherel F, Renand G, Ménissier F: Developing a breeding objective for a French purebred beef cattle selection programme. Livest Prod Sci. 1998, 57: 49-65. 10.1016/S0301-6226(98)00157-2.View ArticleGoogle Scholar
- Gutiérrez JP, Goyache F, Fernández I, Alvarez I, Royo LJ: Genetic relationships among calving ease, calving interval, birth weight, and weaning weight in the Asturiana de los Valles beef cattle breed. J Anim Sci. 2007, 85: 69-75. 10.2527/jas.2006-168.View ArticlePubMedGoogle Scholar
- Cook BR, Tess MW, Kress DD: Effects of selection strategies using heifer pelvic area and sire birth weight expected progeny difference on dystocia in first-calf heifers. J Anim Sci. 1993, 71: 602-607.PubMedGoogle Scholar
- Johanson JM, Berger PJ: Birth weight as a predictor of calving ease and perinatal mortality in Holstein cattle. J Dairy Sci. 2003, 86 (11): 3745-3755. 10.3168/jds.S0022-0302(03)73981-2.View ArticlePubMedGoogle Scholar
- Nobre PR, Misztal I, Tsuruta S, Bertrand JK, Silva LO, Lopes OS: Analyses of growth curves of Nellore cattle by multiple-trait and random regression models. J Anim Sci. 2003, 81: 918-926.PubMedGoogle Scholar
- Albuqueque LG, Meyer K: Estimates of direct and maternal genetic effects for weights from birth to 600 days of age in Nelore cattle. J Anim Breed Genet. 2001, 118: 83-92. 10.1046/j.1439-0388.2001.00279.x.View ArticleGoogle Scholar
- Aliança: Sumário de touros 2011 Aliança Nelore. 2011, Porto Alegre, Brazil: GenSys Consultores AssociadosGoogle Scholar
- Spelman RJ, Huisman AE, Singireddy SR, Coppieters W, Arranz J, Georges M, Garrick DJ: Short communication: quantitative trait loci analysis on 17 nonproduction traits in the New Zealand dairy population. J Dairy Sci. 1999, 82 (11): 2514-2516. 10.3168/jds.S0022-0302(99)75503-7.View ArticlePubMedGoogle Scholar
- Kneeland J, Li C, Basarab J, Snelling WM, Benkel B, Murdoch B, Hansen C, Moore SS: Identification and fine mapping of quantitative trait loci for growth traits on bovine chromosomes 2, 6, 14, 19, 21, and 23 within one commercial line of Bos taurus. J Anim Sci. 2004, 82 (12): 3405-3414.PubMedGoogle Scholar
- Maltecca C, Weigel KA, Khatib H, Cowan M, Bagnato A: Whole-genome scan for quantitative trait loci associated with birth weight, gestation length and passive immune transfer in a Holstein x Jersey crossbred population. Anim Genet. 2009, 40 (1): 27-34. 10.1111/j.1365-2052.2008.01793.x.View ArticlePubMedGoogle Scholar
- McClure MC, Morsci NS, Schnabel RD, Kim JW, Yao P, Rolf MM, McKay SD, Gregg SJ, Chapple RH, Northcutt SL, Taylor JF: A genome scan for quantitative trait loci influencing carcass, post-natal growth and reproductive traits in commercial Angus cattle. Anim Genet. 2010, 41 (6): 597-607. 10.1111/j.1365-2052.2010.02063.x.View ArticlePubMedGoogle Scholar
- Cole JB, Wiggans GR, Ma L, Sonstegard TS, Lawlor TJ, Crooker BA, Van Tassell CP, Yang J, Wang S, Matukumalli LK, Da Y: Genome-wide association analysis of thirty one production, health, reproduction and body conformation traits in contemporary U.S. Holstein cows. BMC Genomics. 2011, 12: 408-10.1186/1471-2164-12-408.PubMed CentralView ArticlePubMedGoogle Scholar
- Elsik CG, Tellam RL, Worley KC, Gibbs RA, Muzny DM, Weinstock GM, Adelson DL, Eichler EE, Elnitski L, Guigó R, Hamernik DL, Kappes SM, Lewin HA, Lynn DJ, Nicholas FW, Reymond A, Rijnkels M, Skow LC, Zdobnov EM, Schook L, Womack J, Alioto T, Antonarakis SE, Astashyn A, Chapple CE, Chen HC, Chrast J, Câmara F, Ermolaeva O, Bovine Genome Sequencing and Analysis Consortium, et al: The genome sequence of taurine cattle: a window to ruminant biology and evolution. Science. 2009, 324 (5926): 522-528.PubMed CentralView ArticlePubMedGoogle Scholar
- Gibbs RA, Taylor JF, Van Tassell CP, Barendse W, Eversole KA, Gill CA, Green RD, Hamernik DL, Kappes SM, Lien S, Matukumalli LK, McEwan JC, Nazareth LV, Schnabel RD, Weinstock GM, Wheeler DA, Ajmone-Marsan P, Boettcher PJ, Caetano AR, Garcia JF, Hanotte O, Mariani P, Skow LC, Sonstegard TS, Williams JL, Diallo B, Hailemariam L, Martinez ML, Morris CA, Bovine HapMap Consortium, et al: Science. 2009, 324 (5926): 528-532.View ArticlePubMedGoogle Scholar
- Matukumalli LK, Lawley CT, Schnabel RD, Taylor JF, Allan MF, Heaton MP, O’Connell J, Moore SS, Smith TP, Sonstegard TS, Van Tassell CP: Development and characterization of a high density SNP genotyping assay for cattle. PLoS One. 2009, 4 (4): e5350-10.1371/journal.pone.0005350.PubMed CentralView ArticlePubMedGoogle Scholar
- Snelling WM, Allan MF, Keele JW, Kuehn LA, McDaneld T, Smith TP, Sonstegard TS, Thallman RM, Bennett GL: Genome-wide association study of growth in crossbred beef cattle. J Anim Sci. 2010, 88 (3): 837-848. 10.2527/jas.2009-2257.View ArticlePubMedGoogle Scholar
- Pausch H, Flisikowski K, Jung S, Emmerling R, Edel C, Götz KU, Fries R: Genome-wide association study identifies two major loci affecting calving ease and growth-related traits in cattle. Genetics. 2011, 187: 289-297. 10.1534/genetics.110.124057.PubMed CentralView ArticlePubMedGoogle Scholar
- Pryce JE, Hayes BJ, Bolormaa S, Goddard ME: Polymorphic regions affecting human height also control stature in cattle. Genetics. 2011, 187: 981-984. 10.1534/genetics.110.123943.PubMed CentralView ArticlePubMedGoogle Scholar
- Nishimura S, Watanabe T, Mizoshita K, Tatsuda K, Fujita T, Watanabe N, Sugimoto Y, Takasuga A: Genome-wide association study identified three major QTL for carcass weight including the PLAG1-CHCHD7 QTN for stature in Japanese Black cattle. BMC Genet. 2012, 13: 40-PubMed CentralView ArticlePubMedGoogle Scholar
- Gudbjartsson DF, Walters GB, Thorleifsson G, Stefansson H, Halldorsson BV, Zusmanovich P, Sulem P, Thorlacius S, Gylfason A, Steinberg S, Helgadottir A, Ingason A, Steinthorsdottir V, Olafsdottir EJ, Olafsdottir GH, Jonsson T, Borch-Johnsen K, Hansen T, Andersen G, Jorgensen T, Pedersen O, Aben KK, Witjes JA, Swinkels DW, den Heijer M, Franke B, Verbeek AL, Becker DM, Yanek LR, Becker LC, et al: Many sequence variants affecting diversity of adult human height. Nat Genet. 2008, 40: 609-615. 10.1038/ng.122.View ArticlePubMedGoogle Scholar
- Weedon MN, Lango H, Lindgren CM, Wallace C, Evans DM, Mangino M, Freathy RM, Perry JR, Stevens S, Hall AS, Samani NJ, Shields B, Prokopenko I, Farrall M, Dominiczak A, Johnson T, Bergmann S, Beckmann JS, Vollenweider P, Waterworth DM, Mooser V, Palmer CN, Morris AD, Ouwehand WH, Zhao JH, Li S, Loos RJ, Diabetes Genetics Initiative; Wellcome Trust CaseControl Consortium, et al: Genome-wide association analysis identifies 20 loci that influence adult height. Nat Genet. 2008, 40: 575-583. 10.1038/ng.121.PubMed CentralView ArticlePubMedGoogle Scholar
- Lettre G, Jackson AU, Gieger C, Schumacher FR, Berndt SI, Sanna S, Eyheramendy S, Voight BF, Butler JL, Guiducci C, Illig T, Hackett R, Heid IM, Jacobs KB, Lyssenko V, Uda M, Boehnke M, Chanock SJ, Groop LC, Hu FB, Isomaa B, Kraft P, Peltonen L, Salomaa V, Diabetes Genetics Initiative; FUSION; KORA; Prostate, Lung Colorectal and Ovarian Cancer Screening Trial; Nurses’ HealthStudy; SardiNIA, et al: Identification of ten loci associated withheight highlights new biological pathways in human growth. Nat Genet. 2008, 40: 584-591. 10.1038/ng.125.PubMed CentralView ArticlePubMedGoogle Scholar
- Karim L, Takeda H, Lin L, Druet T, Arias JA, Baurain D, Cambisano N, Davis SR, Farnir F, Grisart B, Harris BL, Keehan MD, Littlejohn MD, Spelman RJ, Georges M, Coppieters W: Variants modulating the expression of a chromosome domain encompassing PLAG1 influence bovine stature. Nat Genet. 2011, 43: 405-413. 10.1038/ng.814.View ArticlePubMedGoogle Scholar
- Littlejohn M, Grala T, Sanders K, Walker C, Waghorn G, Macdonald K, Coppieters W, Georges M, Spelman R, Hillerton E, Davisand S, Snell R: Genetic variation in PLAG1 associates with early life body weight and peripubertal weight and growth in Bos taurus. Anim Genet. 2012, 43 (5): 591-594. 10.1111/j.1365-2052.2011.02293.x.View ArticlePubMedGoogle Scholar
- Garrick DJ, Taylor JF, Fernando RL: Deregressing estimated breeding values and weighting information for genomic regression analyses. Genet Sel Evol. 2009, 41: 55-10.1186/1297-9686-41-55.PubMed CentralView ArticlePubMedGoogle Scholar
- Cole JB, VanRaden PM, O’Connell JR, Van Tassell CP, Sonstegard TS, Schnabel RD, Taylor JF, Wiggans GR: Distribution and location of genetic effects for dairy traits. J Dairy Sci. 2009, 92: 2931-2946. 10.3168/jds.2008-1762.View ArticlePubMedGoogle Scholar
- Fortes MRS, Lehnert SA, Bolormaa S, Reich C, Fordyce G, Corbet NJ, Whan V, Hawken RJ, Reverter A: Finding genes for economically important traits: Brahman cattle puberty. Animal Production Science. 2012, 52: 143-150. 10.1071/AN11165.View ArticleGoogle Scholar
- Fortes MR, 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. BiolReprod. 2012, 87 (3): 58-Google Scholar
- Hawken RJ, Zhang YD, Fortes MRS, Collis E, Barris WC, Corbet NJ, Williams PJ, Fordyce G, Holroyd RG, Walkley JRW, Barendse W, Johnston DJ, Prayaga KC, Tier B, Reverter A, Lehnert SA: Genome-wide association studies of female reproduction in tropically adapted beef cattle. J Anim Sci. 2012, 90: 1398-1410. 10.2527/jas.2011-4410.View ArticlePubMedGoogle Scholar
- Lindholm-Perry AK, Kuehn LA, Smith TP, Ferrell CL, Jenkins TG, Freetly HC, Snelling WM: A region on BTA14 that includes the positional candidate genes LYPLA1, XKR4 and TMEM68 is associated with feed intake and growth phenotypes in cattle. Anim Genet. 2012, 43 (2): 216-219. 10.1111/j.1365-2052.2011.02232.x.View ArticlePubMedGoogle Scholar
- Bolormaa S, Porto Neto LR, Zhang YD, Bunch RJ, Harrison BE, Goddard ME, Barendse W: A genome-wide association study of meat and carcass traits in Australian cattle. J Anim Sci. 2011, 89 (2): 297-309.Google Scholar
- Porto Neto LR, Bunch RJ, Harrison BE, Barendse W: Variation in the XKR4 gene was significantly associated with subcutaneous rump fat thickness in indicine and composite cattle. Anim Genet. 2012, 43 (6): 785-789. 10.1111/j.1365-2052.2012.02330.x.View ArticlePubMedGoogle Scholar
- Lango Allen H, Estrada K, Lettre G, Berndt SI, Weedon MN, Rivadeneira F, Willer CJ, Jackson AU, Vedantam S, Raychaudhuri S, Ferreira T, Wood AR, Weyant RJ, Segrè AV, Speliotes EK, Wheeler E, Soranzo N, Park JH, Yang J, Gudbjartsson D, Heard-Costa NL, Randall JC, Qi L, Vernon Smith A, Mägi R, Pastinen T, Liang L, Heid IM, Luan J, Thorleifsson G, et al: Hundreds of variants clustered in genomic loci and biological pathways affect human height. Nature. 2010, 467: 832-838. 10.1038/nature09410.PubMed CentralView ArticlePubMedGoogle Scholar
- Van Dyck F, Declercq J, Braem CV, Van de Ven WJ: PLAG1, the prototype of the PLAG gene family: versatility in tumour development. Int J Oncol. 2007, 30: 765-774.PubMedGoogle Scholar
- Hensen K, Braem C, Declercq J, Van Dyck F, Dewerchin M, Fiette L, Denef C, Van de Ven WJM: Targeted disruption of the murine PLAG1 proto-oncogene causes growth retardation and reduced fertility. Dev Growth Differ. 2004, 46: 459-470. 10.1111/j.1440-169x.2004.00762.x.View ArticlePubMedGoogle Scholar
- R Development Core Team: R: A language and environment for statistical computing. 2008, Vienna, Austria: R Foundation for Statistical Computing, http://www.R-project.org,Google Scholar
- Aulchenko YS, Ripke S, Isaacs A, Van Duijn CM: GenABEL: an R library for genome-wide association analysis. Bioinformatics. 2007, 23 (10): 1294-1296. 10.1093/bioinformatics/btm108.View ArticlePubMedGoogle Scholar
- Amin N, Van Duijn CM, Aulchenko YS: A genomic background based method for association analysis in related individuals. PLoS One. 2007, 2: e1274-10.1371/journal.pone.0001274.PubMed CentralView ArticlePubMedGoogle Scholar
- Astle W, Balding DJ: Population structure and cryptic relatedness in genetic association studies. Stat Sci. 2009, 24 (4): 451-471. 10.1214/09-STS307.View ArticleGoogle Scholar
- Mardia KV: Some properties of classical multi-dimensional scaling. Communications on Statistics-Theory and Methods. 1978, A7 (13): 1233-1241.View ArticleGoogle Scholar
- Hartigan JA, Wong MA: A K-means clustering algorithm. Applied Statistics. 1979, 28: 100-108. 10.2307/2346830.View ArticleGoogle Scholar
- Aulchenko YS, de Koning D-J, Haley C: Genomewide rapid association using mixed model and regression: a fast and simple method for genomewide pedigree-based quantitative trait loci association analysis. Genetics. 2007, 177: 577-585. 10.1534/genetics.107.075614.PubMed CentralView ArticlePubMedGoogle Scholar
- Devlin B, Roeder K: Genomic control for association studies. Biometrics. 1999, 55: 997-1004. 10.1111/j.0006-341X.1999.00997.x.View ArticlePubMedGoogle Scholar
- Kinsella RJ, Kähäri A, Haider S, Zamora J, Proctor G, Spudich G, Almeida-King J, Staines D, Derwent P, Kerhornou A, Kersey P, Flicek P: Ensembl BioMarts: a hub for data retrieval across taxonomic space. Database (Oxford). 2011, 2011: Bar030-10.1093/database/bar030.View ArticleGoogle Scholar
- Hu ZL, Park CA, Wu XL, Reecy JM: Animal QTLdb: an improved database tool for livestock animal QTL/association data dissemination in the post-genome era. Nucleic Acids Res. 2012, in pressGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.