A Genome Wide Association Study on Age at First Calving Using High Density Single Nucleotide Polymorphism Chips in Hanwoo (Bos taurus coreanae)

Article information

Asian-Australas J Anim Sci.. 2014;27(10):1406-1410
*Corresponding Author: Jong-Joo Kim. Tel: +82-53-810-3027, Fax: +82-53-801-3027, E-mail: kimjj@ynu.ac.kr

School of Biotechnology, Yeungnam University, Gyeongsan 712-749, Korea


The two authors equally contributed to this work

Received 2014 April 14; Revised 2014 June 25; Accepted 2014 July 14.


Age at first calving is an important trait for achieving earlier reproductive performance. To detect quantitative trait loci (QTL) for reproductive traits, a genome wide association study was conducted on the 96 Hanwoo cows that were born between 2008 and 2010 from 13 sires in a local farm (Juk-Am Hanwoo farm, Suncheon, Korea) and genotyped with the Illumina 50K bovine single nucleotide polymorphism (SNP) chips. Phenotypes were regressed on additive and dominance effects for each SNP using a simple linear regression model after the effects of birth-year-month and polygenes were considered. A forward regression procedure was applied to determine the best set of SNPs for age at first calving. A total of 15 QTL were detected at the comparison-wise 0.001 level. Two QTL with strong statistical evidence were found at 128.9 Mb and 111.1 Mb on bovine chromosomes (BTA) 2 and 7, respectively, each of which accounted for 22% of the phenotypic variance. Also, five significant SNPs were detected on BTAs 10, 16, 20, 26, and 29. Multiple QTL were found on BTAs 1, 2, 7, and 14. The significant QTLs may be applied via marker assisted selection to increase rate of genetic gain for the trait, after validation tests in other Hanwoo cow populations.


Reproductive performance is a good determinant of productivity and sustainability in farm animals (Olsen et al., 2011). The period that a cow needs to reach its maturity and to reproduce for the first time is referred to as age at first calving (van Raden and Klaaskate, 1993). When the heifer reaches its puberty, it is inseminated at the onset and, after a gestation period, the heifer calves for the first time. Age at first calving (AFC) is an index of animal’s entry into the production system for the next generation, and therefore, reducing the AFC is of crucial importance to enhance reproductive performance of animals (Sasaki et al., 2014).

Most of the agriculturally important traits are complex in genetic architecture due to their polygenic nature, i.e. many genes with small effects, which makes it difficult to identify and characterize individual quantitative trait loci (QTL) for genetic prediction (Datewyler et al., 2010; Hayes et al., 2010). However, due to rapid advancement of high throughput genotyping technologies, high density single nucleotide polymorphisms (SNPs) chips became available, which made it possible to perform whole genome association (WGA) studies in farm animals (Goddard and Hayes, 2009). The WGA analyses would significantly aid our understanding of genetic architecture for complex traits (Lee et al., 2013).

Due to low heritability of the reproduction traits, they are less responsive to genetic progress through selection. Also, the traits are sex-limited or measured only at maturity (Pryce et al., 1997; Dematawewa and Berger, 1998; Rox-strom et al., 2001). Therefore, for efficient genetic improvement of the reproductive traits, genome selection could be a good option, as it allows selection decisions at early stage (Goddard and Hayes, 2009) by performing WGA studies (Ball, 2005). Currently, several genome-wide SNP panels are available for cattle, enabling QTL mapping or prediction of animal’s genetic merit for the traits of interest (Goddard and Hayes, 2009).

In cattle, there are some reports on WGA analyses about AFC (Olsen et al., 2011; Sasaki et al., 2014). Results of WGA study on AFC, date at first service, calving interval and some fertility traits were also reported in Japanese Black cattle (Minozzi et al., 2013; Sasaki et al., 2014). Olsen et al. (2011) reported WGA results for fertility traits in Norwegian cattle. Hawken et al. (2012) reported WGA results in Brahman and some tropical breeds to identity QTL regions for age at puberty. Age at first service was also investigated with WGA analyses in Holstein by Hoglund et al. (2009) and Daetwyler et al. (2008). However, there is no report about WGA results for reproductive traits in Hanwoo.

The aim of this study was to identify QTL for AFC in Korean native cattle, Hanwoo (Bos taurus coreanae), through the WGA analysis.


Animals and phenotype data

Ninety six Hanwoo heifers were chosen in a local (Jook-Am Hanwoo) farm, Suncheon, Korea. The heifers were born between 2008 and 2010 from 13 sires. The heifers were raised in large pen with the capacity of about 100 individuals and fed ad libitum with a total mixed ration concentrates including 35% Italian ryegrass silage, 35% rice straws, 24% (yeast fermented) bran, and 6% brewers’ grain. During the experimental period, estrus cycle of each heifer was checked mainly by one veterinarian, and the first artificial insemination was applied to the heifers at approximately 15 months of age. After pregnancy, the heifers were raised in the same large pens before pregnancy, sometimes driven out of the pens and pastured during days. Age at first calving was calculated between the birth days of a heifer and the first calf of the heifer.

Molecular data

The heifers were genotyped with the Illumina 50K bovine SNP chips according to the Infinium HD Assay Ultra Protocol (Matukumalli et al., 2009). Among the SNPs in the chip, some SNPs on autosomal chromosomes were screened and removed before WGA testing, which met the following three criteria; i) the number of genotype group with one or none (e.g. only AA genotypes and no AB or BB), ii) with a minor allele frequency less than 0.05, and iii) with the proportion of genotyped individuals less than 90%.

Statistical analysis

A SAS general linear model procedure (SAS version 9.1) was used to test pre-adjust the animal phenotypes before WGA testing. Then, the significant effect of birth year-season was fitted as a fixed effect in Animal model with a genome-relationship (G) matrix that was constructed using the R (version 2.15.0) subroutine (Van Raden, 2008). After the residuals of the phenotype were obtained using ASREML program (version 3.0), the residuals were regressed on additive and dominance effects of each SNP under a linear regression model (PLINK version 1.07). In the model, SNP genotypes with AA, AB, and BB were assigned as 1, 0, and −1 for the additive, and 0, 1, and 0 for the dominance effect, respectively. For significance threshold, 0.1% point-wise p value from F distributions was applied for each SNP test.

Among the significant SNPs, the best set of SNP markers were selected using the forward regression procedure (Neter et al., 1990), because some of the significant SNPs would yield redundant information due to linkage disequilibrium between closely linked SNPs. Inclusion of each SNP into the model was determined at α = 0.05 level. For significant SNPs, phenotypic variation explained by each (jth) SNP (S2snpj) was calculated as i=13αi2fi-μ2, where i denotes each genotype, αi is the estimated additive or dominance effects, e.g.ᾱ, , and +ᾱ for BB, AB, and AA genotypes, respectively, for the SNP with both additive (ᾱ) and dominance () effects, fi is the frequency of ith genotype, μ is the population mean that can be expressed as +ᾱ (fAAfBB)+ d̄ fAB (Falconer and Mackay, 1996). Proportion of phenotypic variance due to the jth SNP was then estimated as S2snpj/S2p, in which S2p was phenotype variance that was obtained from residual values of the trait after adjusting the fixed and polygenic effects. The proportion of phenotype variance due to all of the significant SNPs was estimated as ((RSSnoSNPSSSNPs)/RSSnoSNP), where RSS was residual sum of squares for the model with or without fitting the significant SNPs of the trait.


Summary statistics of the AFC that were recorded on the 96 heifers are described in Table 1. The average AFC was 717.5±67.4, ranging from the minimum of 607 to the maximum of 898 days.

Summary statistics for age at first calving in 96 Hanwoo heifers

A total of 15 SNPs for AFC were detected (Table 2). The significant SNPs were located on ten bovine chromosomes (BTAs), among which 9 SNPs were located on BTAs 1, 2, 7, and 14. The two significant SNPs with strong statistical evidence were located at 128.9 Mb of BTA2 and at 111.1 Mb of BTA7, respectively, each of which accounted for 22% phenotypic variance (Table 2).

Identities, effects and positions of the significant SNPs for age at first calving in a population of Hanwoo heifers

Among the 15 significant SNPs, ten SNPs had (over-) dominance mode of gene action, i.e. the estimates of dominance effect were greater than of additive effect (Table 2). This result indicates that dominance effect of the QTL is highly appreciable in reproductive traits in Hanwoo, which is in agreement with the previous reports about strong effects of heterosis on reproduction in cattle (Cundiff et al., 1974).

The number of heifers for WGA analysis in this study was very limited (n = 96), which did not enable detection of many QTL (SNP) for AFC with strong statistical evidence. As sample size is one of major factor in efficient genome wide association (GWA) analysis (Goddard and Hayes, 2009), more samples are needed to detect other QTL, and the detected QTL need to be validated by adding further samples. However, our results provide a good preliminary report regarding the first genome scan on reproductive traits in Hanwoo heifers.

Schulman et al. (2011) performed a GWA study in the Finish Ayrshire breed to identify QTL for non-return rate, time from first to last insemination, and time from calving to first insemination (in days) for cows. They identified significant SNPs on BTAs 1, 2, 4, 8, 12, 13, 20, 24, and 27. In this study, we detected significant SNPs on BTAs 1, 2, 7, 9, 10, 14, 16, 20, 26, and 29 (Table 2). Hawken et al. (2012) reported three QTL for age at puberty in Brahman and tropical breeds at 43.2 Mb, 81.3 Mb and 140.8 cM of BTA1, while one significant SNP for AFC SNP was detected at 27 Mb of BTA1 (Table 2).

Minozzi et al. (2013) reported one SNP for days to first service at 136.18 Mb of BTA 2 in Italian Holstein cattle. In the similar region, we found one AFC SNP (129 Mb). Sasaki et al. (2014) reported a significant SNP at 118.24 Mb of BTA 2 that was associated with AFC in Japanese Black Cattle, which is closely located at the region (121.5 Mb) where one AFC SNP was detected in this study. Datewyler et al. (2008) reported one QTL for calving age at first service at 24 cM of BTA2 in Holstein. We also detected an AFC SNP at 25 Mb of the same chromosome (Table 2).

In this study, we detected two AFC SNPs at 22 Mb and 111 Mb of BTA 7, respectively (Table 2). Holmberg and Anderson-Eklund (2006) reported QTL for still birth at 22 Mb of BTA7 in Swedish dairy cattle and Hoglund et al. (2009) found one QTL for number of inseminations per conception in Holstein cattle at 111.6 cM of BTA 7.

Sahana et al. (2010) reported a genomic region spanning 40.4 to 52.2 Mb of BTA 9 that was associated with the interval from first to last insemination in Holstein cattle, while Olsen et al. (2011) reported a QTL for retained placenta at 52.64 Mb of BTA 9 in Norwegian Red cattle. We also detected one AFC SNP at 52 Mb of BTA 9 in Hanwoo heifers (Table 2).

One significant SNP for AFC was detected at 88 Mb of BTA10 in this study. Kuhn et al. (2013) reported a QTL for dystocia and still birth at 80 to 87 cM of the same chromosome in Holstein.

On BTA 14, two significant AFC SNPs were detected at 6 Mb and 63 Mb, respectively (Table 2). Hawken et al. (2012) reported one QTL for age at puberty at 61.9 Mb of BTA14 in Brahman and Tropical Composite breed.

McClure et al. (2010) found one SNP for percent of unassisted birth in first calf heifers at 49.1 Mb of BTA 20 in Angus cattle. In the similar region (50 Mb), an AFC SNP was detected in this study (Table 2).

On BTA 26, McClure et al. (2010) detected a significant SNP for percent of unassisted birth in first calf heifers at 18.8 Mb in Angus cattle. We also found an AFC SNP at 22 Mb of the same chromosome. In a proximal chromosomal region (17.8 Mb), Peters et al. (2013) reported one SNP for the first service conception in Brangus cattle.

On BTA29, a significant AFC SNP was detected at 12 Mb, while Hawken et al. (2012) found a significant SNP for post-partum anestrus at 11.3 Mb of the same chromosome in Brahman and tropical composite breeds.


Early AFC results in short generation interval, which would improve lifetime productivity of cows. We identified 15 significant QTL for age at first caving in Hanwoo heifers, and in similar chromosomal regions where the QTL were detected, many QTL for the related reproductive traits were reported in other cattle breeds. The identified QTL would aid in performing marker assisted selection (MAS) for genetic improvement of AFC in Hanwoo cows. However, implementation of effective and efficient MAS requires validation of the detected QTL and characterization of the QTL e.g. causal mutations responsible for AFC.


This research was supported by the Technology Development Program for Agriculture and Forestry (Project No. 311016-3), Ministry of Agriculture, Forestry, and Fisheries, Republic of Korea, 2013.


Ball RD. 2005;Experimental designs for reliable detection of linkage disequilibrium in unstructured random population association studies. Genetics 170:859–873.
Cundiff LV, Gregory KE, Koch RM. 1974;Effects of heterosis on reproduction in Hereford, Angus and Shorthorn cattle. J Anim Sci 38:711–727.
Daetwyler HD, Pong-Wong R, Villanueva B, Woolliams JA. 2010;The impact of genetic architecture on genome-wide evaluation methods. Genetics 185:1021–1031.
Daetwyler HD, Schenkel FS, Sargolzaei M, Robinson JA. 2008;A genome scan to detect quantitative trait loci for economically important traits in Holstein cattle using two methods and a dense single nucleotide polymorphism map. J Dairy Sci 91:3225–3236.
Dematawewa CM, Berger PJ. 1998;Genetic and phenotypic parameters for 305-day yield, fertility, and survival in Holsteins. J Dairy Sci 81:2700–2709.
Falconer DS, Mackay TFC. 1996. Introduction to Quantitative Genetics 4th edth ed. Pearson & Prentice Hall. London, UK:
Goddard ME, Hayes BJ. 2009;Mapping genes for complex traits in domestic animals and their use in breeding programmes. Nat Rev Genet 10:381–391.
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. 2012;Genome-wide association studies of female reproduction in tropically adapted beef cattle. J Anim Sci 90:1398–1410.
Hayes BJ, Pryce J, Chamberlain AJ, Bowman PJ, Goddard ME. 2010;Genetic architecture of complex traits and accuracy of genomic prediction: coat colour, milk-fat percentage and type in Holstein cattle as contrasting model traits. PLoS Genet 6(9):e1001139.
Hoglund JK, Guldbrandtsen B, Su G, Thomsen B, Lund MS. 2009;Genome scan detects quantitative trait loci affecting female fertilitytraits in Danish and Swedish Holstein cattle. J Dairy Sci 92:2136–2143.
Holmberg M, Andersson-Eklund L. 2006;Quantitative trait loci affecting fertility and calving traits in Swedish dairy cattle. J Dairy Sci 89:3664–3671.
Kuhn C, Bennewitz J, Reinsch N, Xu N, Thomsen H, Looft C, Brockmann GA, Schwerin M, Weimann C, Hiendleder S, Erhardt G, Medjugorac I, Forster M, Brenig B, Reinhardt F, Reents R, Russ I, Averdunk G, Blumel J, Kalm E. 2003;Quantitative trait loci mapping of functional traits in the German Holstein cattle population. J Dairy Sci 86:360–368.
Lee SH, Choi BH, Lim D, Gondro C, Cho YM, Dang CG, Sharma A, Jang GW, Lee KT, Yoon D, Lee HK, Yeon SH, Yang BS, Kang HS, Hong SK. 2013;Genome-wide association study identifies major loci for carcass weight on BTA14 in Hanwoo (Korean cattle). PLoS ONE 8(10):e74677. 10.1371/journal.pone.0074677.
Matukumalli LK, Lawley CT, Schnabel RD, Taylor JF, Allan MF, Heaton MP, Connell JO, Moore SS, Smith TPL, Sonstegard TS, Van Tassell CP. 2009;Development and characterization of high density SNP genotyping assay for cattle. PLoS One 4(4):e5350.
McClure MC, Morsci NS, Schnabel RD, Kim JW, Yao P, Rolf MM, McKay SD, Gregg SJ, Chapple RH, Northcutt SL, Taylor JF. 2010;A genome scan for quantitative trait loci influencing carcass, postnatal growth and reproductive traits in commercial Angus cattle. Anim Genet 41:597–607.
Minozzi G, Nicolazzi EL, Stella A, Biffani S, Negrini R, Lazzari1 B, Ajmone-Marsan P, Williams JL. 2013;Genome wide analysis of fertility and production traits in italian holstein cattle. PLoS ONE 8(11):e80219. 10.1371/journal.pone.0080219.
Neter J, Wasserman W, Kutner MH. 1990. Applied Linear Statistical Models 3rd edth ed. McGraw-Hill/Irwin. Boston, MA, USA:
Olsen HG, Hayes BJ, Kent MP, Nome T, Svendsen M, Larsgard AG, Lien S. 2011;Genome-wide association mapping in Norwegian Red cattle identifies quantitative trait loci for fertility and milk production onBTA12. Anim Genet 42:466–474.
Peters SO, Kizilkaya K, Garrick DJ, Fernando RL, Reecy J, Weaber RL, Silver GA, Thomas MG. 2013;Heritability and Bayesian genome wide association study of first service conception and pregnancy in Brangus heifers. J Anim Sci 91:605–612.
Pryce JE, Veerkamp RF, Thompson R, Hill WG, Simm G. 1997;Genetic aspects of common health disorders and measures of fertility in Holstein Friesian cattle. Anim Sci 65:353–360.
Rox-strom A, Strandberg E, Berglund B, Emanuelson U, Philipsson J. 2001;Genetic and environmental correlations among female fertility traits and milk production in different parities of Swedish Red and White dairy cattle. Acta Agric Scand A Anim Sci 51:7–14.
Sahana G, Guldbrandtsen B, Bendixen C, Lund MS. 2010;Genome wide association mapping for female fertility traits in Danish and Swedish Holstein cattle. Anim Genet 4:579–588.
Schnabel RD, Sonstegard TS, Taylor JF, Ashwell MS. 2005;Whole genome scan to detect QTL for milk production, conformation, fertility and functional traits in two US Holstein families. Anim Genet 36:408–416.
Sasaki S, Ibi T, Ikeda S, Sugimoto Y. 2014;A genome-wide association study reveals a quantitative trait locus for age at first calving in delta/notch-like EGF repeat containing on chromosome 2 in Japanese Black cattle. Anim Genet 45:285–287.
Schulman NF, Sahana G, Iso-Touru T, McKay SD, Schnabel RD, Lund MS, Taylor JF, Virta J, Vilkki JH. 2011;Mapping of fertility traits in Finnish Ayrshire by genome-wide association analysis. Anim Genet 42:263–269.
Van Raden PM, Klaaskate EJH. 1993;Genetic evaluation of length of productive life including predicted longevity of live cows. J Dairy Sci 76:2758–2764.
Van Raden PM. 2008;Efficient methods to compute genomic predictions. J Dairy Sci 91:4414–4423.

Article information Continued

Table 1

Summary statistics for age at first calving in 96 Hanwoo heifers

Trait Average Std. Dev Min Max CV
Age at first calving (d) 717.5 67.4 607 898 9.2

Std. Dev, standard deviation; Min, minimum; Max, maximum; CV, coefficient of variation (%).

Table 2

Identities, effects and positions of the significant SNPs for age at first calving in a population of Hanwoo heifers

SNP marker1 SNP2 BTA kb3 Additive4 SE Dominance5 SE −log10P6 %σp2,7
ARS-BFGL-NGS-98203 [C/A] 1 1,010 −7.22 6.09 32.78 8.21 3.50 0.15
ARS-BFGL-NGS-115015 [C/T] 1 27,181 −13.97 5.54 28.45 7.91 3.35 0.14
BTA-46908-no-rs [C/A] 2 24,734 43.18 10.94 −33.72 12.55 3.48 0.21
UA-IFASA-7185 [T/G] 2 121,503 −2.03 5.78 35.27 7.89 4.73 0.19
ARS-BFGL-NGS-8433 [A/G] 2 128,973 36.44 9.65 −23.33 11.15 3.31 0.22
Hapmap57628-rs29027155 [A/G] 7 21,845 34.82 10.93 −50.81 12.85 3.40 0.15
Hapmap51497-BTA-80591 [A/G] 7 111,106 28.54 6.18 −34.50 8.24 5.76 0.22
Hapmap32423-BTA-92864 [T/C] 9 52,133 20.89 8.72 5.67 10.48 3.16 0.14
Hapmap60876-rs29013997 [A/C] 10 87,684 15.94 5.34 12.63 7.80 3.02 0.12
Hapmap27709-BTC-057052 [C/T] 14 5,925 7.22 11.03 −36.42 12.91 3.01 0.13
Hapmap51719-BTA-35195 [G/A] 14 62,751 5.27 5.68 26.89 7.68 3.02 0.13
ARS-BFGL-NGS-93176 [T/C] 16 55,239 −13.47 10.82 48.42 13.46 3.77 0.16
Hapmap48608-BTA-111028 [A/G] 20 49,810 1.59 5.31 −30.18 7.65 3.39 0.14
BTB-00932019 [C/A] 26 21,870 −18.95 6.26 30.31 8.48 3.35 0.14
UA-IFASA-9371 [G/T] 29 12,185 −26.15 8.64 2.89 11.22 3.14 0.12

SNP, single nucleotide polymorphism; BTA, bovine chromosomes; SE, standard error.


SNP marker annotations and their positions.


Alternative nucleotides.


Estimates of additive and dominance effects of the SNP. The estimates are allele substitution effects of the former against the latter SNP in the SNP column.


Negative logarithm of base ten for the comparison-wise p-value of the test statistic against the null hypothesis of no SNP effect under the one SNP model.


Proportion of phenotypic variance explained by each SNP or across all of the significant SNPs for the given trait.