- Research article
- Open Access
Estimation of the genetic parameters for semen traits in Chinese Holstein bulls
BMC Genetics volume 20, Article number: 51 (2019)
Semen traits are important for the widespread use of superior bulls. Thus, the objective of this study was to estimate the heritability of five semen traits, ejaculate volume (VE), progressive sperm motility (SM), sperm concentration (SC), number of sperm (NSP), and number of progressive motile sperm (NMSP), and their genetic correlations (rg). The dataset being studied consisted of 1450 Chinese Holstein bulls with full pedigree information, born between 1996 and 2012, representing 11 AI centers. Genetic parameters were estimated in a multivariate analysis using the average information restricted maximum likelihood estimation of variance (AI-REML).
The estimates of heritability for VE, SM, SC, NSP, and NMSP were 0.15, 0.12, 0.22, 0.16 and 0.12, respectively. The genetic correlations among the five semen traits ranged from 0.02 (VE and SC) to 0.99 (NSP and NMSP).
Our findings provide useful information on the heritability of semen traits in Holstein bulls and the relationships among them, and should assist in selection for improvement of semen traits in Chinese Holstein bulls.
Female fertility has been widely studied for breeding purposes in many countries [1,2,3], whereas less focus has been placed on male fertility. In cattle, male fertility can be measured directly on semen or indirectly from females. Several studies have been conducted on bull fertility [4,5,6,7,8]. Unlike the male fertility traits measured from females, semen traits are measured directly in males and can be used as indicators of male fertility. Semen traits are complex and are influenced by genetic and environmental factors, such as season, age, bull handlers, and the interval between consecutive semen collections in days [9,10,11]. The first published estimates of the heritability of semen traits were low to medium (ranging from 0.02 to 0.31) , and the heritability of semen traits has been studied in many cattle populations [7, 10, 13,14,15,16]. However, there were big differences in the estimates from these studies.
Conception rate in the Holstein cattle has declined dramatically recently in some countries [3, 17]. Because conception rate is affected by semen characteristics such as sperm volume and sperm motility, the estimation of genetic parameters for semen traits can offer important information that could be used to improve the conception rate in Holstein dairy cattle. However, since the first AI center was built in China (1973), there has been little focus on semen traits, and no studies have been done on semen quality or the genetic parameters of semen traits in Chinese Holstein bulls. Therefore, the objectives of this study were: 1) to estimate the heritability of semen traits in Chinese Holstein cattle, and 2) to estimate the genetic and phenotypic correlations among semen traits. The results will help us better understand the genetic basis of semen traits, contributing to their genetic improvement in Holstein cattle.
Data from a total of 1450 Chinese Holstein bulls were used in this study. The bulls were from 11 artificial insemination (AI) centers in China, and the number of bulls in each center ranged from 69 to 202. All the bulls were born between 1996 and 2012, and five semen traits were recorded across years (Fig. 1), ejaculate volume (VE) (mL) was read directly from a graduated collection tube (graduation: 1/10 mL); progressive sperm motility (SM) was the proportion of forward moving sperm, which was evaluated under microscope by experienced technicians; sperm concentration (SC) was measured using a spectrophotometer; number of sperm per ejaculate (NSP) was calculated by multiplying VE and SC, and the number of motile sperm (NMSP) was calculated by multiplying NSP and SM. In each AI center, ejaculations were conducted once or twice per day, often by different semen handlers, and the interval between the consecutive collections for each bull was 30 min to 1 h. The approaches used for semen collection were standardized across AI centers.
The pedigree information was extracted from the Dairy Data Centre of China (http://www.holstein.org.cn/) by tracing back as many generations as possible for the bulls. The number of animals in the pedigree file was 12,985.
For each of the five semen traits, the model for the genetic parameters was:
where yghijklmno is the record of the semen phenotype; μ is the overall intercept; Recyeag is the fixed effect of the record year from 2003 to 2013; BulAgeh is the fixed effect of age of bull at collection, where the five age classes (h = 1 to 5) were defined as 10–24 months, 25–36 months, 37–60 months, 61–96 months, and older than 96 months; Seai is the fixed effect of the season of collection, where the four seasons (i = 1 to 4) were defined as December to February, March to May, June to August, and September to November; ColIntj is the fixed effect of the interval of between consecutive collections for each bull (j = 1 to 5), which was defined as 0 d, 1–4 d, 5–7 d, 7–11 d, and longer than 11 d; Numk is the fixed effect of the number of collections on the respective collection day (k = 1 to 2); BulStal is the fixed effect of AI station (l = 1 to 11); po is the random permanent environment effect with var.(p) ~ N(0, Iσ2pe), where σ2pe is the permanent environmental variance and I is the identity matrix; an is the random additive genetic effect with var.(a) ~ N(0, Aσ2a), where σ2a is the additive genetic variance and A is the relationship matrix; BulHanm (BulStal) is the effect of BulHanm (m = 1 to 54), as a random effect of bull handlers, nested within AI station; eghijklmno is the random residual effect with var.(e) ~ N(0, Iσ2e), where σ2e is the residual variance. (Co) variance components were estimated using the DMU software package  with the average information restricted maximum likelihood (AI-REML) .
Results and discussions
The total number of ejaculates was 332,632 and the average number of recorded ejaculates per bull was 229. The bulls had an average semen volume of 7.22 mL, progressive sperm motility of 67%, sperm concentration of 11.37 × 10 9/mL, total sperm of 83.52 × 10 9, and total motility sperm of 59.21 × 10 9 (Table 1). There was a significant difference in the five semen traits among the 11 AI centers (P < 0.01).
The estimated genetic parameters for semen traits are presented in Table 2. The heritability estimates for VE, SM, SC, NSP, and NMSP were 0.15, 0.12, 0.22, 0.16, and 0.12, respectively. All the estimates were significantly greater than zero as shown by the standard error of these estimates (Table 2). A low genetic correlation (0.02) was found between VE and SC. SM had a moderate genetic correlation with all the other traits (ranged from 0.18 to 0.36). The genetic correlations among all other semen traits were high (from 0.65 to 0.99). Similarly, a low phenotypic correlation (0.09) was also found between VE and SC. SM had a moderate phenotypic correlation with all the other traits (ranged from 0.26 to 0.56). The phenotypic correlations among the other semen traits were high (from 0.63 to 0.96). All phenotypic correlations were in the same direction as the corresponding genetic correlations.
In general, our heritability estimates for VE, SC, SM, NSP, and NMSP ranged from 0.12 to 0.22, which implies that there was a substantial environmental influence on these traits . Our heritability estimate for VE was consistent with previous studies [7, 10, 14, 18]. However, our value was lower than the estimate reported by Ducrocq et al. (1995) (0.65), both those estimates were based on mean values from more than 10 records per individual in young Normande bulls. The heritability for SM (0.12) was also in accordance with the results of some previous studies [12, 14, 19], but was lower than estimates from others [7, 10, 15, 16]. These differences may be due to different breeds, genetic backgrounds, and different models used for parameter estimation. Fewer heritability estimates have been published for NSP and NMSP. For NSP, heritability was estimated to be 0.09 by Druet et al. (2009), 0.22 by Gredler et al. (2007), and 0.24 by Knights et al. (1984), which are close to our estimate of 0.16. For NSP and NMSP, our heritability estimates are 0.16 and 0.12, respectively,which are much lower than the 0.38 and 0.49 for young bulls, and 0.54 and 0.64 for older bulls reported by Mathevon et al. (1998), which was based on a small sample size of198 bulls separated into young bulls (up to 30 months old) and older bulls (between 4 and 6 years old) . The heritability of NMSP (0.12) was higher than the 0.04 value reported by Gredler et al. (2007). This difference may be attributed to a different breed—their breed was Austrian dual-purpose Simmental (Fleckvieh)—and the sample size in their study (12,746 ejaculates from 301 bulls) was smaller than in this study . The heritability of 0.22 for SC was in accordance with values from other studies on Holstein bulls [7, 18].
In general, the genetic correlations among all five semen traits were positive. The genetic correlation between VE and SM (0.36) was similar to the value of 0.31 reported by Gredler et al. (2007), whereas the negative correlations of − 0.20 and − 0.38 were reported by Druet et al. (2009) and Kealey et al. (2006), respectively. The genetic correlation between SC and NSP (0.46, s.e. 0.18) reported by Druet et al.  was lower than the estimate of 0.71 (s.e. 0.16) found in the present study; the estimates for the genetic correlation between SC and SM differed only slightly between these studies (0.12 v 0.18, respectively) . The results of genetic correlations in our study suggest that selection for one semen trait should have no unfavourable effects on the other semen traits measured.
The low heritability of semen traits (0.12 to 0.22) indicates a substantial environmental influence; this suggests that potential environmental factors that may affect semen quality, such as bull handler, season, and the interval between semen collections need to be controlled carefully so that genetic differences are more accurately assessed. Therefore, to improve semen quality and increase the quantity of superior bulls, environmental influences such as bull handlers, season, and the interval of two continual semen collections for each bull, should be optimized. In addition, other environmental influences that are not routinely recorded in AI stations, such as temperature, humidity, and feed, should be recorded, and studies of candidate genes for semen traits and GWAS can be performed with the aim of selecting bulls with superior semen quality and quantity.
In this study, we estimated the genetic parameters of five semen traits in Chinese Holstein bulls for the first time. The heritability of semen traits was varied from 0.12 to 0.22. The correlations among the semen traits were all positive. NSP had a strong genetic correlation with all other semen traits and its heritability was 0.16, which suggests the focus of genetic improvement ought to be around NSP. Our results provide valuable insights into the genetic basis of semen traits in the production period and should assist in designing selection programs to improve semen traits in Chinese Holstein bulls.
Availability of data and materials
The datasets analyzed during the current study are available from the corresponding author on a reasonable request.
The number of motile sperm () was calculated by multiplying NSP and SM
Number of sperm per ejaculate was calculated by multiplying VE and SC
Sperm concentration was measured using a spectrophotometer
Progressive sperm motility was the proportion of forward moving sperm, which was evaluated under microscope by experienced technicians
Ejaculate volume (mL) was read directly from a graduated collection tube (graduation: 1/10 mL)
Weigel KA. Prospects for improving reproductive performance through genetic selection. Anim Reprod Sci. 2006;96(3):323–30.
Miglior F, Muir BL, Van Doormaal BJ. Selection indices in Holstein cattle of various countries. J Dairy Sci. 2005;88(3):1255–63.
Jamrozik J, Fatehi J, Kistemaker GJ, Schaeffer LR. Estimates of genetic parameters for Canadian Holstein female reproduction traits. J Dairy Sci. 2005;88(6):2199–208.
Lyons RE, Loan NT, Dierens L, Fortes MRS, Kelly M, McWilliam SS, Li Y, Bunch RJ, Harrison BE, Barendse W, et al. Evidence for positive selection of taurine genes within a QTL region on chromosome X associated with testicular size in Australian Brahman cattle. BMC Genet. 2014;15(1):6.
Fortes MRS, Reverter A, Kelly M, McCulloch R, Lehnert SA. Genome-wide association study for inhibin, luteinizing hormone, insulin-like growth factor 1, testicular size and semen traits in bovine species. Andrology. 2013;1(4):644–50.
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. Biol Reprod. 2012;87(3):58.
Druet T, Fritz S, Sellem E, Basso B, Gérard O, Salas-Cortes L, Humblot P, Druart X, Eggen A. Estimation of genetic parameters and genome scan for 15 semen characteristics traits of Holstein bulls. J Anim Breed Genet. 2009;126(4):269–77.
Blaschek M, Kaya A, Zwald N, Memili E, Kirkpatrick BW. A whole-genome association analysis of noncompensatory fertility in Holstein bulls. J Dairy Sci. 2011;94(9):4695–9.
Stälhammar E-M, Janson L, Philipsson J. Genetic studies on fertility in a.I. bulls. I. Age, season and genetic effects on semen characteristics in young bulls. Anim Reprod Sci. 1989;19(1):1–17.
Mathevon M, Buhr MM, Dekkers JCM. Environmental, management, and genetic factors affecting semen production in Holstein bulls. J Dairy Sci. 1998;81(12):3321–30.
Fuerst-Waltl B, Schwarzenbacher H, Perner C, Sölkner J. Effects of age and environmental factors on semen production and semen quality of Austrian Simmental bulls. Anim Reprod Sci. 2006;95(1–2):27–37.
Smith BA, Brinks JS, Richardson GV. Estimation of genetic parameters among breeding soundness examination components and growth traits in yearling Bulls1. J Anim Sci. 1989;67(11):2892–6.
Suchocki T, Szyda J. Genome-wide association study for semen production traits in Holstein-Friesian bulls. J Dairy Sci. 2015;98(8):5774–80.
Kealey CG, MacNeil MD, Tess MW, Geary TW, Bellows RA. Genetic parameter estimates for scrotal circumference and semen characteristics of line 1 Hereford bulls1. J Anim Sci. 2006;84(2).
Gredler B, Fuerst C, Fuerst-Waltl B, Schwarzenbacher H, Sölkner J. Genetic parameters for semen production traits in Austrian dual-purpose Simmental bulls. Reprod Domest Anim. 2007;42(3):326–8.
Ducrocq V, Humblot P. Genetic characteristics and evolution of semen production of young Normande bulls. Livest Prod Sci. 1995;41(1):1–10.
Lucy M. Reproductive loss in high-producing dairy cattle: where will it end? J Dairy Sci. 2001;84(6):1277–93.
Makulska J, Hagger C, Künzi N, Kupferschmied HU. Genetic and environmental influences on semen traits in a.I. bulls. Reprod Domest Anim. 1993;28(4):279–84.
Knights SA, Baker RL, Gianola D, Gibb JB. Estimates of Heritabilities and of genetic and phenotypic correlations among growth and reproductive traits in yearling Angus Bulls1. J Anim Sci. 1984;58(4):887–93.
Madsen P, Jensen J. A User’s guide to DMU: a package for analyzing multivariate mixed models, version 6. Danish Institute of Agricultural Sciences, Dept of Animal Breeding and Genetics, Research Centre Foulum. 2006;8830:27.
Jensen J, Mäntysaari EA, Madsen P, Thompson R. Residual maximum likelihood estimation of (co) variance components in multivariate mixed linear models using average information. Jour Ind Soc Ag Statistics. 1997;49:215–36.
We appreciate eleven AI centers in China for providing the data and pedigree for the Chinese Holstein bulls.
This work is supported by the 863 project (2013AA102504), the National Science and Technology Programs of China (2011BAD28B02), National Key Technologies R&D Program (2012BAD12B01), Beijing Dairy Industry Innovation Team, China Agricultural Research System (CARS-37), and Xinjiang Province Key Technology Integration and Demonstration Program (201230116). We are deeply grateful to all donors who participated in this program. The funding body did not participate in the design of the study and collection, analysis, interpretation of data or in writing the manuscript.
Ethics approval and consent to participate
The phenotypic data were collected from eleven AI centers in China with the owners’ consents approved by the Institutional Animal Care and Use Committee (IACUC) at China Agricultural University.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.