Abstract

The elite bull plays an extremely important role in the genetic progression of the dairy cow population. The previous results indicated the potential positive relationship of large scrotal circumference (SC) with improved semen volume, concentration, and motility. In order to improve bull’s semen quantity and quality by selection, it is necessary to estimate the genetic parameters of semen traits and their correlations with other conformation traits such as SC that could be used for an indirect selection. In this study, the genetic parameters of seven semen traits (n = 66,260) and nine conformation traits (n = 3,642) of Holstein bulls (n = 453) were estimated by using the bivariate repeatability animal model with the average information-restricted maximum likelihood (AI-REML) approach. The results showed that the estimated heritabilities of semen traits ranged from 0.06 (total number of motile sperm, TNMS) to 0.37 (percentage of abnormal sperm, PAS) and conformation traits ranged from 0.23 (pin width, PW) to 0.69 (hip height, HH). The highest genetic correlations were found between semen volume per ejaculation (SVPE), semen concentration per ejaculation (SCPE), total number of sperm (TNS), and TNMS traits that were 0.97, 0.98, 1.00, and 0.99, respectively. Phenotypic correlations between SC and SVPE, SCPE, TNS, and TNMS were 0.35, 0.35, 0.48, and 0.42, respectively. In summary, the moderate or high heritability of semen traits indicates that genetic improvement of semen quality by selection is feasible, where SC could be a useful trait for indirect selection or as correlated information to improve semen quantity and production in the practical bull breeding programs.

1. Introduction

The widespread use of frozen semen from bulls for artificial insemination contributes 70% to the genetic progression in dairy cattle [1]. As semen qualities directly determine the conception rates, the reproductive performances of bulls require to be continually improved to supply high-quality semen for efficient production. Therefore, selection is necessary for high-quality and large-quantity semen from genetically superior sires [2, 3]. Meyer et al. [4] estimated the genetic correlations between reproduction and growth traits in Australian beef cattle and suggested the indirect selection of larger scrotal circumference (SC) sizes for the shorter calving intervals [4]. In addition, other conformation traits such as body weight have positive correlations with viable sperm cell production that can also be used for the selection of reproduction traits [5]. Bulls’ breeding fitness may perform superiorly if their morphological conformation is better as well, so bull conformation traits have potential influences on semen production and quality [6]. Thus, understanding the heritability of these traits is crucial for effective selection programs.

In the previous studies, Nickolov and Metodiev [7] found the heritabilities of semen volume, motility, and concentration traits equal to 0.43, 0.54, and 0.61, respectively, in Holstein bulls. However, Druet et al. [8] presented a relatively low heritability of semen volume (0.15) but a moderate heritability of concentration (0.30) and a relatively high heritability of motility (0.59). Of note, a significant positive effect of SC on semen volume, concentration, and motility can be identified, for example, the significant genetic correlations with volume (0.20), concentration (0.77), and motility (0.34) in Hereford bulls [9]. In young Nellore bulls, a low negative genetic correlation (−0.24) of SC with the percentage of total sperm defects was observed [10].

The inconsistent genetic parameters of these previous studies could be caused by heterogeneous datasets, different genetic and environmental backgrounds, and various statistical models among different populations [11]. It becomes apparent that further estimations of genetic parameters in one population are needed to guarantee accurate breeding programs in the dairy industry. Until now, correlation studies between semen traits and other conformation traits such as body length and chest girth are also scarce [12]. Therefore, our study aims to estimate the genetic parameters of semen traits and body conformation traits in Chinese Holstein bulls to investigate the possible selection methods to improve semen quality. The obtained estimation results could be used for genetic improvement of semen production and quality in practical selection programs of the Holstein population.

2. Materials and Methods

2.1. Semen and Conformation Traits

The phenotype records of semen (n = 66,260) and conformation (n = 3,642) traits were collected from 453 Chinese Holstein bulls over the last 10 years from Shandong OX Livestock Breeding Co., Ltd., in Shandong province, China. After several quality control procedures, 56,980 semen records and 2,635 conformation records were retained for further analysis (Table 1). The bulls were raised in their own management systems and fed with customized pellet feed and high-quality green hay (Leymus chinensis, alfalfa, and oats) with the mineral supplement. All bulls were vaccinated and treated based on the guidelines and recommendations of the veterinarian. The distributions of seven semen traits and nine conformation traits are shown in Figures 1 and 2, respectively. The detailed information on them is listed in Supplementary File 1.

The seven semen traits included semen volume per ejaculation (SVPE), semen concentration per ejaculation (SCPE), sperm motility per ejaculation (SMPE), total number of sperm (TNS), total number of motile sperm (TNMS), percentage of abnormal sperm (PAS), and postthaw motility (PTM). Semen was collected twice a week from each bull by using an artificial vagina and then immediately stored at 37°C in a water bath. The fresh semen traits of SVPE (mL), SCPE (108/mL), SMPE (%), and PAS (%) were evaluated by the skilled technicians, where TNS = SVPE (mL) × SCPE (108/mL) and TNMS = TNS (108) × SMPE (%). After 5–7 days of cryopreservation, two straws of each bull were randomly taken and thawed at 38°C for 20 seconds. Immediately, PTM (%) was evaluated under a light microscope.

The nine conformation traits were generally measured in the middle of the month which included body weight (BW), withers height (WH), body length (BL), chest girth (CG), abdominal circumference (AC), tube circumference (TC), pin width (PW), hip height (HH), and scrotal circumference (SC). WH (cm) was measured from the highest point of the armor to the ground as the vertical distance. BL (cm) was measured from the shoulder to the end of the sciatic bone. CG (cm) was measured as the vertical circumference of the body at the posterior scapula. TC (cm) was measured as the upper third of the tube in the left lead limb. HH (cm) was measured from the center of the two waist angles to the ground. SC (cm) was measured with a special measuring ruler as the circumference of the most prominent part of the middle part of the scrotum.

2.2. Data Preprocessing for Genetic Analysis

The screening criteria of qualified semen were as follows: SVPE ≥ 2 mL, SCPE ≥ 2 × 108/mL, and SMPE ≥ 0.4. Three levels were involved in the season effects, where June, July, and August were the first level, January, February, and December were the second level, and the rest of the months were grouped into the third level. According to pedigree information obtained from the Dairy Data Center of China (https://www.holstein.org.cn) and the Canadian Dairy Network (https://www.cdn.ca), a total of 453 bulls were traced back to three generations, with 2,030 individuals included in the pedigree.

2.3. Two-Trait Genetic Analysis

A bivariate repeatability animal model was used to estimate variance and covariance components of semen and conformation traits using average information-restricted maximum likelihood (AI-REML) in the DMU package [13]. The model following the previous study [14, 15] is described as follows:where and are the vectors of phenotypes (i.e., semen or conformation traits), and are the design matrices for fixed effects, and are the vectors of fixed effects (i.e., year-season and age effects), and are the design matrices for additive genetic effects, and are the vectors of additive genetic effects that satisfy and , where and are additive genetic variances and is the relationship matrix calculated using pedigree information, and are the design matrices for permanent environmental effects, and are vectors of permanent environmental effects that satisfy and , where and are permanent environmental variances, and and are the vectors of random environmental effects that satisfy and , where and are residual variances and is the identity matrix.

The variance and covariance structure of the additive genetic, permanent environmental, and random environmental effects for two traits (i.e., semen or conformation traits) arewhere heritability () and repeatability () of the two traits are estimated as , , , and . Genetic correlation and phenotypic correlation are estimated as and , where and are covariances of the additive genetic and permanent environmental effects, respectively. The covariances of random environmental effects were not considered between semen and conformation traits, but were considered among semen traits or conformation traits.

3. Results

3.1. Statistical Description and Heritability of Semen and Conformation Traits

The statistical description of seven semen traits and nine conformation traits are presented in Table 2. The heritability, repeatability, and variance component estimates are presented in Table 3. It showed that semen traits had low to medium heritabilities ranging from 0.06 to 0.37, while conformation traits had relatively high heritabilities ranging from 0.23 to 0.69. TNMS, PTM, and SCPE traits exhibited the lower heritabilities of 0.06, 0.09, and 0.12, respectively. In contrast, AC, WH, SC, and HH traits exhibited higher heritabilities of 0.54, 0.59, 0.64, and 0.69, respectively. The repeatability ranges from 0.24 to 0.49 among semen traits and 0.27 to 0.76 among conformation traits (Table 3).

3.2. Genetic and Phenotypic Correlations of Semen and Conformation Traits

The genetic correlations (upper diagonal) and phenotypic correlations (lower diagonal) of semen and conformation traits are presented in Table 4. The phenotypic correlations among semen traits range from −0.20 to 0.72, with the highest observed correlation of SCPE with TNS and TNMS traits (0.72). For conformation traits, the phenotypic correlations range from 0.20 to 0.83, all showing positive correlations. The correlations of SM and PTM traits with conformation traits are generally weak or close to zero, while SVPE, SCPE, TNS, and TNMS traits exhibit positive correlations ranging from 0.07 to 0.48 (Table 4).

Genetic correlations among semen traits ranged from −0.34 to 0.80, where SVPE with TNS (0.79) and SCPE with TNMS (0.80) showed high correlations. Genetic correlations among conformation traits were all positive ranging from 0.18 to 0.88, where BW with CG (0.88), BW with AC (0.84), BW with BL (0.82), and BW with TC (0.82) showed high correlations. General high genetic correlations between semen traits and conformation traits were found largely from −1.00 to 1.00, where TNS with SC showed the highest correlation (1.00). Notably, genetic correlations for SC with SVPE, SCPE, TNS, and TNMS traits were 0.97, 0.98, 1.00, and 0.99, respectively, which indicates that the larger scrotal circumference can lead to a greater volume of semen, higher semen concentration, total number of sperm, and total number of motile sperm per ejaculation. Such high genetic correlations were also observed in WH and HH traits. Furthermore, BW also had a moderate genetic correlation with SVPE (0.94), TNS (0.93), and TNMS (0.89) traits (Table 4).

4. Discussion

Based on 56,980 semen records and 2,635 conformation records of 453 Chinese Holstein bulls over the last 10 years, the estimated genetic parameters of our study were favorable and comparable. For semen traits, the SVPE (0.36), SMPE (0.27), SCPE (0.12), PAS (0.37), and TNMS (0.06) heritabilities estimated in our study (Table 3) were similar to those of other studies (0.04∼0.43) [3, 8, 9, 1619]. Carvalho Filho et al. [20] reported a quite low heritability of SMPE (only 0.05) suggesting the significant influence of it by the unstable environmental factors and the high phenotypic variability. However, the PAS (0.37) and TNS (0.32) heritabilities (Table 3) were quite higher than those of the beef cattle (0.03 and 0.13) as reported by Butler et al. [19] and Berry et al. [3]. The higher estimates of heritability in our study than in other studies were likely influenced by the limited repeated phenotypic records that made it hard to estimate the larger variance component for the permanent environmental effect than for the additive genetic effect [2124].

For conformation traits, the HH heritability (0.69) was consistent with the results of Arango et al. [25]; however, the BW heritability (0.33) was clearly lower than the adult weight (0.75) of Hereford cattle [26]. Obviously, different data sources, sample sizes, and applied models may result in different estimated heritabilities of the same conformation traits [27]. Several previous studies also indicated that SC heritability from 0.40 to 0.81 [20, 28, 29] was consistent with our study (0.64) (Table 3). The high estimated heritabilities of conformation traits indicate that genetic improvement of semen production and quality by selection could be feasible if their genetic correlations with conformation traits are high enough.

The findings of genetic correlations among conformation traits underscore a robust relationship in the growth patterns across various body dimensions, showing the notable correlations of BW with BL at 0.71, CG at 0.83, and AC at 0.79 (Table 4). It highlights the significant interdependence in the development of distinct body parts and emphasizes the impact of specific dimensions on the overall conformation of each individual [27]. Moreover, our study reveals the strong genetic correlations ranging from 0.79 to 0.97 between semen and conformation traits. It indicates the significant genetic interconnection between seminal production and quality and various aspects of conformation in Chinese Holstein bulls.

Many studies have also reported the positive correlations between these two types of traits, such as 0.36 for SC with sperm motility in Angus cattle [30] and 0.20 and 0.77 for SC with semen volume and contraction, respectively, in Hereford cattle [9]. Berry et al. [3] also showed the high genetic correlations of SC with ejaculation volume, concentration, and motility (0.20, 0.77, and 0.76) in dairy and beef cattle. It was worth mentioning the quite high genetic correlations for SC with SVPE (0.97), SCPE (0.98), TNS (1.00), and TNMS (0.99) traits (Table 4) in the current study. The elevated genetic correlations suggest that the selection for favorable conformation traits could positively contribute to the enhancement of seminal production and quality, particularly for the SVPE trait. Therefore, the selection of those conformation traits could help improve semen production and quality in Holstein bulls.

Some other conformation traits were related to semen traits in strong positive correlation, specifically for WH with SVPE (0.88), SCPE (0.87), TNS (0.93), and TNMS (0.91) traits and HH with SVPE (0.87), SCPE (0.85), TNS (0.92), and TNMS (0.90) traits. It indicates that superior morphological conformation may underlie much of the variation in the semen volume, concentration, and number properties. Similarly, we identified the robust genetic correlation between BW and TNS and TNMS at 0.93 and 0.89, respectively (Table 4). These results highlight the close linkage between specific body dimensions and reproductive characteristics. Understanding these relationships is pivotal for comprehensive breeding strategies and substantial insights into the particular selection for conformation traits positively influencing both semen production and quality in Chinese Holstein bulls. Undoubtedly, the potential unstable and unfavorable estimates still exist, so further study for more accurate estimations is needed by using an enlarged population size when more phenotypic records and pedigree information are enrolled.

5. Conclusion

In summary, our study provides the estimated heritabilities and genetic correlations of semen and conformation traits in Chinese Holstein bulls. Of note, positive genetic correlations of conformation traits were found with SVPE, SCPE, TNS, and TNMS, especially for SC, BW, WH, and HH traits (0.87∼1.00). These results indicate a robust relationship of the growth traits with semen traits, so that the genetic improvement of semen production and quality traits by selection could be feasible. Especially, SC could be used as an indicator trait for indirect selection, or as an assistant trait to improve selection accuracy for semen production and quality. Such insights shed light on potential synergies between reproductive and conformation traits that emphasize the importance of considering both parts in genetic improvement programs. However, the enlarged population size is still needed to guarantee more accurate estimations in the further practical genetic improvement programs in Holstein bulls.

Abbreviations

AC:Abdominal circumference
AI-REML:Average information-restricted maximum likelihood
BL:Body length
BW:Body weight
CG:Chest girth
HH:Hip height
PAS:Percentage of abnormal sperm
PTM:Postthaw motility
PW:Pin width
SC:Scrotal circumference
SCPE:Semen concentration per ejaculation
SD:Standard deviation
SE:Standard error
SMPE:Sperm motility per ejaculation
SVPE:Semen volume per ejaculation
TC:Tube circumference
TNMS:Total number of motile sperms
TNS:Total number of sperms
WH:Withers height.

Data Availability

All data from the current study are consolidated in tables and Supplementary Materials presented in this article. If further information is required, please contact the corresponding author: Jianbin Li, [email protected].

Ethical Approval

All experimental protocols were approved by the Animal Ethics Committee of Shandong Academy of Agricultural Sciences, China. All methods were carried out in accordance with relevant guidelines and regulations.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Xiao Wang and Jian Yang contributed equally to this work.

Acknowledgments

This research was supported by the National Key R&D Program of China (2021YFF1000701-06 and 2023ZD0404901-05), the Earmarked Fund for CARS-36, and the Shandong Provincial Natural Science Foundation (ZR2021MC070 and ZR2023QC252). The authors thank Shandong OX Livestock Breeding Co., Ltd., in Shandong province, China, for providing phenotype records of semen and conformation traits.

Supplementary Materials

Supplementary File 1: the detailed information of seven semen traits and nine conformation traits. (Supplementary Materials)