Estimation of genetic parameters of fertility traits in Chinese Holstein cattle

Prediction of likelihood of conception in dairy cows using milk mid-infrared spectra collected before the first insemination and machine learning algorithms

Accurate and ex-ante prediction of cows' likelihood of conception (LC) based on milk composition information could improve reproduction management on dairy farms. Milk composition is already routinely measured by mid-infrared (MIR) spectra, which are known to change with advancing stages of pregnancy. For lactating cows, MIR spectra may also be used for predicting the LC. Our objectives were to classify the LC at first insemination using milk MIR spectra data collected from calving to first insemination and to identify the spectral regions that contribute the most to the prediction of LC at first insemination. After quality control, 4,866 MIR spectra, milk production, and reproduction records from 3,451 Holstein cows were used. The classification accuracy and area under the curve (AUC) of 6 models comprising different predictors and 3 machine learning methods were estimated and compared. The results showed that partial least square discriminant analysis (PLS-DA) and random forest had higher prediction accuracies than logistic regression. The classification accuracy of good and poor LC cows and AUC in herd-by-herd validation of the best model were 76.35 ± 10.60% and 0.77 ± 0.11, respectively. All wavenumbers with values of variable importance in the projection higher than 1.00 in PLS-DA belonged to 3 spectral regions, namely from 1,003 to 1,189, 1,794 to 2,260, and 2,300 to 2,660 cm-1. In conclusion, the model can predict LC in dairy cows from a high productive TMR system before insemination with a relatively good accuracy, allowing farmers to intervene in advance or adjust the insemination schedule for cows with a poor predicted LC.

Genetic parameters for milk yield and reproductive traits in Honduran Holstein cattle

The aim of this study was to estimate the (co)variance components and genetic parameters for milk yield adjusted to 305d (MY305), calving-to-conception interval (CCI), number of services per conception (NSC) and calving interval (CI) of Honduran Holstein cows, by fitting a bivariate animal model using Maximum Restricted Likelihood procedures. Model included the fixed effects of calving number, the contemporary calving group (farm-season-year of calving and the cow age as covariate). The estimated means and standard deviations for MY, CCI, NSC and CI were, 5098.60 ± 1564.32 kg, 168.27 ± 104.71 days, 2.46 ± 1.69 services, and 448.73 ± 109.16 days, respectively; and their estimated heritabilities were 0.21 ± 0.05, 0.03 ± 0.028, 0.02 ± 0.024 and 0.06 ± 0.04, respectively. The genetic correlations between MY305 and CCI, NSC and CI were positive and antagonist, with values of 0.64 ± 0.52, 0.99 ± 0.56, and 0.32 ± 0.24 respectively. Even though moderate to low heritability was estimated for MY305, systematic selection for milk yield, with a reduction in reproductive efficiency, if considered as the only selection criterion is important to be considered. By including reproductive traits and considering permanent environment effects into the breeding program, might yield a slow, but constant and permanent improvement over time.

Estimation of genetic parameters for fertility traits in Chinese Holstein of south China

Introduction: Chinese Holstein in South China suffer heat stress for a long period, which leads to evolutionary differences from Chinese Holstein in North China. The aim of this study was to estimate the genetic parameters of fertility traits for Chinese Holstein in South China. Methods: A total of 167,840 Chinese Holstein heifers and cows from Guangming Animal Husbandry Co., LTD farms were used in this study. The fertility traits analyzed were calving interval (CI), days open (DO), age of first service (AFS), age of first calving (AFC), calving to first insemination (CTFS), first insemination to conception (FSTC), gestation length (GL), non-return rate to 56 days (NRR), and number of services (NS). Results: The descriptive statistics revealed that the same trait in heifers performed better than in cows, which was consistent with the other studies. The heritabilities of fertility traits in this study ranged from close to 0 (for NS of cows) to 0.2474 (for AFC of heifers). The genetic correlation of NRR between heifers and cows was 0.9993, which indicates that the NRR for heifers and cows could be treated as one trait in this population. Conclusion: The heritabilities of fertility traits in Chinese Holstein in south China were quite different from the heritabilities of fertility traits in North China. NRR56, NS, AFC, and CI were suggested to be included into the selection index to improve fertility performance of Chinsese Holstein of south China. The results of this study could provide genetic parameters for the animal breeding program of Chinese Holstein in the south of China.

Genetic parameters for health traits and their association with fertility and milk production in Chinese Holsteins

Herd health is one of the key problems influencing the efficiency of the dairy industry. Genetic selection, with a focus on animal health, is important for herd improvement. This study aimed to estimate genetic parameters for health traits and their correlations with fertility and milk production traits in dairy cattle. Based on records from 58,549 lactating cows calved between 2015 and 2021, a total of 24 health traits (six composite health traits and 18 independent health traits), four fertility traits and five milk production traits were analysed. First, linear and threshold animal models were used to estimate the variance components and heritabilities of the health traits. Second, a bivariate linear animal model was used to estimate genetic correlations among all 24 health traits. Finally, a bivariate linear animal model based on records from the first lactation was used to estimate the correlations between health traits and fertility or milk production traits. The results showed that all health traits had low heritabilities, ranging from 0.002 (0.001) to 0.048 (0.004) in the linear model and from <0.001 (0.021) to 0.226 (0.035) in the threshold model. Genetic correlations between health traits across categories were generally low, whereas the relatively high genetic correlations were found between health traits within the same category. In this study, only a few significant and moderate genetic correlations were observed between health traits and fertility or milk production traits. Clinical mastitis showed relatively moderate correlations with fertility traits, ranging from 0.277 (0.113) (interval from first to last insemination) to 0.401 (0.104) (calving interval). Moreover, there were moderate genetic correlations between hoof health and milk production traits. The results from the current study will support balanced dairy breeding to genetically improve disease resistance in dairy cows.

Estimation of genetic parameters and genetic change of first parity reproductive traits in Alpine × Beetal goats

The study of reproductive traits is crucial for improving genetic potential of goats because of their significant utility in meat production. Hence, genetic analysis was conducted for reproductive traits on Alpine × Beetal goats using animal model for first parity data. Information on 1462 reproductive records were collected over five decades from ICAR-National Dairy Research Institute, Karnal, Haryana (1971-2021). Single-trait and multi-trait animal models were used for genetic analysis. Estimates of (co)variance components and genetic parameters were obtained using Gibbs Sampler for Animal Model due to non-normal distribution of data. Six single-trait animal models (including or excluding maternal and environmental effects) were fitted and best models were determined based on Deviance Convergence Criterion values. The prolificacy for the A × B goats for first parity data was 32%, having 68% single births, 31% twins and 1% triplets/quadruplets. The least squares mean for age at first service (AFS), age at first kidding (AFK), service period (SP), dry period (DP), gestation length (GL), kidding interval (KI), litter weight (LW), number of kids born (NKB) and number of females kids born (NFKB) in first parity were 546.15 ± 4.10 days, 679.05 ± 4.07 days, 226.51 ± 4.02 days, 67.96 ± 2.76 days, 150.74 ± 0.13 days, 362.53 ± 3.35 days, 3.99 ± 0.04 kg, 1.32 ± 0.02 and 0.64 ± 0.02, respectively. The heritability estimates obtained from best model for AFS, AFK, GL, KI, SP, and DP were 0.12 ± 0.00, 0.10 ± 0.00, 0.09 ± 0.01, 0.03 ± 0.00, 0.04 ± 0.00, and 0.05 ± 0.00, respectively. For NKB, NFKB and LW, heritability estimates were 0.16 ± 0.01, 0.03 ± 0.03, and 0.04 ± 0.00, respectively. These results imply lower heritability estimates for reproductive traits, and hence meagre scope for selection for further improvement. Maternal effects contributed significantly for traits such as GL, NKB and NFKB. Genetic correlation for number of female kids born was negative with SP and DP which is favourable. Furthermore, genetic correlation was negative for dry period and litter weight which is favourable as number of kids born and litter weight are traits of direct economic importance. Results reveal high genetic potential of this breed for meat industry owing to high prolificacy, provided consistent efforts are made for the genetic improvement of this germplasm.

Exploring milk loss and variability during environmental perturbations across lactation stages as resilience indicators in Holstein cattle

Genetic selection for resilience is essential to improve the long-term sustainability of the dairy cattle industry, especially the ability of cows to maintain their level of production when exposed to environmental disturbances. Recording of daily milk yield provides an opportunity to develop resilience indicators based on milk losses and fluctuations in daily milk yield caused by environmental disturbances. In this context, our study aimed to explore milk loss traits and measures of variability in daily milk yield, including log-transformed standard deviation of milk deviations (Lnsd), lag-1 autocorrelation (Ra), and skewness of the deviations (Ske), as indicators of general resilience in dairy cows. The unperturbed dynamics of milk yield as well as milk loss were predicted using an iterative procedure of lactation curve modeling. Milk fluctuations were defined as a period of at least 10 successive days of negative deviations in which milk yield dropped at least once below 90% of the expected values. Genetic parameters of these indicators and their genetic correlation with economically important traits were estimated using single-trait and bivariate animal models and 8,935 lactations (after quality control) from 6,816 Chinese Holstein cows. In general, cows experienced an average of 3.73 environmental disturbances with a milk loss of 267 kg of milk per lactation. Each fluctuation lasted for 19.80 ± 11.46 days. Milk loss traits are heritable with heritability estimates ranging from 0.004 to 0.061. The heritabilities differed between Lnsd (0.135-0.250), Ra (0.008-0.058), and Ske (0.001-0.075), with the highest heritability estimate of 0.250 ± 0.020 for Lnsd when removing the first and last 10 days in milk in a lactation (Lnsd2). Based on moderate to high genetic correlations, lower Lnsd2 is associated with less milk losses, better reproductive performance, and lower disease incidence. These findings indicate that among the variables evaluated, Lnsd2 is the most promising indicator for breeding for improved resilience in Holstein cattle.

Estimation of Genetic Parameters for Heifer and Cow Fertility Traits Derived from On-Farm AI Service Records of South African Holstein Cattle

This study aimed to derive additional fertility traits from service data and estimate their genetic parameters for their possible inclusion in the South African Holstein cattle breeding programs. Service records (n = 64,464) were collected from 18 South African Holstein herds using on-farm automated milk recording systems. Using a multivariate model, the data were used to estimate heritabilities and correlations among several fertility traits. The pedigree data consisted of information on 18,592 animals born between 1981 and 2013. Heritability estimates observed were low to moderate (0.02 ± 0.00 to 0.24 ± 0.00), indicating that there is some genetic basis for the explored fertility traits to warrant selection. The genetic correlations observed between fertility traits were generally favorable, with some high correlations between age at first service (AFS) and services per conception for heifers (SPCh) (0.73 ± 0.00) and between days from calving to first service (CFS) and services per conception for cows (SPC) (0.90 ± 0.01). Positive genetic correlations indicate that improvement in one trait is coupled with a correlated genetic increase in another trait. The studied fertility traits could be used in addition to AFC and CI to serve as a basis for the selection of reproduction in dairy cattle to minimize selection bias.

A Post-GWAS Functional Analysis Confirming Effects of Three BTA13 Genes CACNB2, SLC39A12, and ZEB1 on Dairy Cattle Reproduction

In our previous GWAS of Chinese and Nordic dairy cattle, genes CACNB2, SLC39A12, and ZEB1 locating on BTA 13 were suggested as candidate genes for reproduction. In this study, validation of these associations was performed in an independent population with records of nine reproductive traits. More importantly, functions of these genes in the reproductive process were verified by employing the expression data of ovarian follicles. The potential variants within the three genes were firstly detected in 68 Chinese Holstein bulls, and then screened in 1,588 Chinese Holstein cows using the KASP (Kompetitive allele-specific PCR) method. There were nine variants with polymorphisms in CACNB2, five in SLC39A12, and four in ZEB1, respectively, of which one SNP was in the upstream regulatory region, two in exon region, four in downstream regulatory region, and 11 SNPs in intronic regions. Amongst the 18 variants, g.33267056T/G in CACNB2 explained the largest phenotypic variance for age at first calving (0.011%), interval from first to last insemination (0.004%), and calving ease (0.002%), while g.32751518G/A in SLC39A12 contributed the most to stillbirth in heifers (0.038%). Two haplotype blocks were constructed for CACNB2 while one each for SLC39A12 and ZEB1, which were significantly associated with five reproductive traits, including age at the first service, age at the first calving, calving ease in heifers and cows, and the interval from calving to the first insemination. We then studied the profile of gene expression in granulosa cells isolated from four developmental stages of ovarian follicles from eight dairy cows. All three genes were differentially expressed between ovarian follicles with different sizes (p < 0.05), indicating their potential roles in the reproductive process of dairy cows. This study successfully demonstrated the associations of three BTA 13 genes CACNB2, SLC39A12, and ZEB1 with reproduction and further examined their expression levels in ovarian follicles directly. These findings can be beneficial for the ongoing genomic selection program for reproductive traits which have long been considered as traits that are difficult to achieve genetic improvement due to the lack of efficient genetic markers.

Estimation of Genetic Parameters for Female Fertility Traits in the Polish Holstein-Friesian Population

Simple Summary

Female fertility is an important functional trait in dairy cattle. The aim of this study was to estimate genetic parameters for some fertility traits of Polish Holstein-Friesian cows. The results indicate that a few analyzed traits could be of use in genetic evaluation of Polish Holstein-Friesian cows.

Abstract

The objective of this study was to estimate genetic parameters for the analyzed fertility traits of Polish Holstein-Friesian primiparous and multiparous cows, as a step toward the incorporation of new traits into routine genetic evaluation. Lactation records of 116,836 Polish Holstein-Friesian cows were studied. The records cover the first three lactations of all cows. The cows, daughters of 2884 sires, were calved from 2006 to 2020. The conception rate for heifers (CRh) and for cows in the second (CR2) and third parity (CR3), the interval from first calving to first insemination (CTFI), the days open (DO), and the interval from first to successful insemination for heifers (FSh) and for cows in second (FS2) and third (FS3) parity were analyzed. The BLUPf90 package and a Bayesian method via Gibbs sampling were used to estimate (co)variance components. In general, all heritabilities were low and ranged from 0.013 (CTFI) to 0.038 (FS2). The heritability of conception rate and interval from first to successful insemination was slightly lower for heifers than for cows. Genetic correlations were moderate to high with two exceptions: correlation of CTFI with CRh (0.033) and with FSh (−0.051). The results indicate that a few analyzed traits could be used in genetic evaluation of Polish Holstein-Friesian cows. It is suggested to complement the selection index with traits for both heifers and cows, such as the interval from first to successful insemination (i.e., FSh, FS2 and FS3), in order to enable effective improvement of female fertility scores in the Polish Holstein-Friesian population.

Genetic analyses of blood β-hydroxybutyrate predicted from milk infrared spectra and its association with longevity and female reproductive traits in Holstein cattle

Ketosis is one of the most prevalent and complex metabolic disorders in high-producing dairy cows and usually detected through analyses of β-hydroxybutyrate (BHB) concentration in blood. Our main objectives were to evaluate genetic parameters for blood BHB predicted based on Fourier-transform mid-infrared spectra from 5 to 305 d in milk, and estimate the genetic relationships of blood BHB with 7 reproduction traits and 6 longevity traits in Holstein cattle. Predicted blood BHB records of 11,609 Holstein cows (after quality control) were collected from 2016 to 2019 and used to derive 4 traits based on parity number, including predicted blood BHB in all parities (BHBp), parity 1 (BHB1), parity 2 (BHB2), and parity 3+ (BHB3). Single- and multitrait repeatability models were used for estimating genetic parameters for the 4 BHB traits. Random regression test-day models implemented via Bayesian inference were used to evaluate the daily genetic feature of BHB variability. In addition, genetic correlations were calculated for the 4 BHB traits with reproduction and longevity traits. The heritability estimates of BHBp, BHB1, BHB2, and BHB3 ranged from 0.100 ± 0.026 (± standard error) to 0.131 ± 0.023. The BHB in parities 1 to 3+ were highly genetically correlated and ranged from 0.788 (BHB1 and BHB2) to 0.911 (BHB1 and BHB3). The daily heritability of BHBp ranged from 0.069 to 0.195, higher for the early and lower for the later lactation periods. A similar trend was observed for BHB1, BHB2, and BHB3. There are low direct genetic correlations between BHBp and selected reproductive performance and longevity traits, which ranged from -0.168 ± 0.019 (BHBp and production life) to 0.157 ± 0.019 (BHBp and age at first calving) for the early lactation stage (5 to 65 d). These direct genetic correlations indicate that cows with higher BHBp (greater likelihood of having ketosis) in blood usually have shorter production life (-0.168 ± 0.019). Cows with higher fertility and postpartum recovery, such as younger age at first calving (0.157 ± 0.019) and shorter interval from calving to first insemination in heifer (0.111 ± 0.006), usually have lower BHB concentration in the blood. Furthermore, the direct genetic correlations change across parity and lactation stage. In general, our results suggest that selection for lower predicted BHB in early lactation could be an efficient strategy for reducing the incidence of ketosis as well as indirectly improving reproductive and longevity performance in Holstein cattle.

Genetic and Genomic Analyses of Service Sire Effect on Female Reproductive Traits in Holstein Cattle

Fertility and reproductive performance are key drivers of dairy farm profitability. Hence, reproduction traits have been included in a large majority of worldwide dairy cattle selection indexes. The reproductive traits are lowly heritable but can be improved through direct genetic selection. However, most scientific studies and dairy cattle breeding programs have focused solely on the genetic effects of the dam (GED) on reproductive performance and, therefore, ignored the contribution of the service sire in the phenotypic outcomes. This study aimed to investigate the service sire effects on female reproductive traits in Holstein cattle from a genomic perspective. Genetic parameter estimation and genome-wide association studies (GWAS) were performed for the genetic effect of service sire (GESS) on conception rate (CR), 56-day non-return rate (NRR56), calving ease (CE), stillbirth (SB), and gestation length (GL). Our findings indicate that the additive genetic effects of both sire and dam contribute to the phenotypic variance of reproductive traits measured in females (0.0196 vs. 0.0109, 0.0237 vs. 0.0133, 0.0040 vs. 0.0289, 0.0782 vs. 0.0083, and 0.1024 vs. 0.1020 for GESS and GED heritability estimates for CR, NRR56, CE, SB, and GL, respectively), and these two genetic effects are positively correlated for SB (0.1394) and GL (0.7871). Interestingly, the breeding values for GESS on insemination success traits (CR and NRR56) are unfavorably and significantly correlated with some production, health, and type breeding values (ranging from -0.449 to 0.274), while the GESS values on calving traits (CE, SB, and GL) are usually favorably associated with those traits (ranging from -0.493 to 0.313). One hundred sixty-two significant single-nucleotide polymorphisms (SNPs) and their surrounding protein-coding genes were identified as significantly associated with GESS and GED, respectively. Six genes overlapped between GESS and GED for calving traits and 10 genes overlapped between GESS for success traits and calving traits. Our findings indicate the importance of considering the GESS when genetically evaluating the female reproductive traits in Holstein cattle.

Genetic parameters and genetic trends for reproductive traits in Egyptian buffalo

Reproductive traits are important for farm profitability because failure to reproduce is the primary reason for culling animals. Study objectives were to estimate genetic parameters and evaluate the trends for reproductive traits. Age at first calving (AFC), gestation length (GL), postpartum interval to pregnancy (PPIP), calving interval (CI) and calving ease score (CE) were recorded. A total of 38,906 records were available from 2426 buffalo cows. There was evaluation of genetic parameters using three models. The first model was applied to the first three parities fitting animal as a random effect. There was also a repeatability model utilized where data from all parities were evaluated to estimate heritability and repeatability. There was also a bivariate model to estimate genetic correlations between reproductive traits. Heritability estimates ranged from 0.0001 to 0.12 for PPIP and CE, respectively. Repeatability estimates were low to moderate ranging from 0.13 to 0.38 for PPIP and GL, respectively. There were close genetic correlations of 0.99 and - 0.93 between PPIP-CI and GL-CE, respectively. Genetic correlations between the other reproductive traits were low to moderate. Genetic trends for all reproductive traits were negative with and of a small magnitude, and regression coefficients were significant only for AFC and PPIP. The results from the current study supported the urgent need, not only for genetic or genomic selection improvement programs, but also for improving the farm management practices for reproductive traits in Egyptian buffalo.

Genetic parameters for rectal temperature, respiration rate, and drooling score in Holstein cattle and their relationships with various fertility, production, body conformation, and health traits

Genetic selection for improved climatic resilience is paramount to increase the long-term sustainability of high-producing dairy cattle, especially in face of climate change. Various physiological indicators, such as rectal temperature (RT), respiration rate score (RR), and drooling score (DS), can be used to genetically identify animals with more effective coping mechanisms in response to heat stress events. In this study, we investigated genetic parameters for RT, RR (score from 1-3), and DS (score from 1-3). Furthermore, we assessed the genetic relationship among these indicators and other economically important traits for the dairy cattle industry. After data editing, 59,265 (RT), 30,290 (RR), and 30,421 (DS) records from 13,592 lactating Holstein cows were used for the analyses. Variance components were estimated based on a multiple-trait repeatability animal model. The heritability ± standard error estimate for RT, RR, and DS was 0.06 ± 0.01, 0.04 ± 0.01, and 0.02 ± 0.01, respectively, whereas their repeatability was 0.19, 0.14, and 0.14, respectively. Moderate genetic correlations of RR with RT and DS (0.26 ± 0.11 and 0.25 ± 0.16) and nonsignificant correlation between RT and DS (-0.11 ± 0.14) were observed. Furthermore, the approximate genetic correlations between RT, RR, and DS with 12 production, 29 conformation, 5 fertility and reproduction, 5 health, and 9 longevity-indicator traits were assessed. In general, the approximate genetic correlations calculated were low to moderate. In summary, 3 physiological indicators of heat stress response were measured in a large number of animals and shown to be lowly heritable. There is a value in developing a selection index including all the 3 indicators to improve heat tolerance in dairy cattle. All the unfavorable genetic relationships observed between heat tolerance and other economically important traits can be accounted for in a selection index to enable improved climatic resilience while also maintaining or increasing productivity in Holstein cattle.

Genetics of heifer reproductive traits in Japanese Black cattle

Objective

The objective of this study was to identify environmental factors strongly associated with and to estimate genetic parameters of reproductive traits in Japanese Black heifers.

Methods

Data included reproduction records of Japanese Black heifers born between 2004 and 2014. First service non-return rate (NRR) to 56 days from first to successful insemination (FS), number of services per conception (IN), age at first calving (AFC) and gestation length were analyzed with the use of the general linear model. Genetic parameters were estimated with the use of the univariate animal model of the residual maximum likelihood.

Results

Averages of reproductive traits over eleven years were assessed, and the effects of farm, year, month, artificial insemination technician and interaction of farm×year on the traits were determined. Estimated heritability of FS was very low and that of AFC was higher than that of the other traits. A close genetic relation was observed among NRR, IN, and FS; however, their heritabilities were very low. AFC shows favorable genetic correlation with IN and FS.

Conclusion

Low heritabilities of most reproductive traits in Japanese Black heifers are strongly influenced by farm management practices, and that large residual variances make genetic evaluation difficult. Among the reproductive traits, AFC is potentially more useful for genetic improvement of heifer reproductive traits because it has high heritability and favorable genetic correlations with IN and FS.

Variance components and correlations of female fertility traits in Chinese Holstein population

BACKGROUND

The objective of the present study was to estimate (co)variance components of female fertility traits in Chinese Holsteins, considering fertility traits in different parities as different traits. Data on 88,647 females with 215,632 records (parities) were collected during 2000 to 2014 from 32 herds in the Sanyuan Lvhe Dairy Cattle Center, Beijing, China. The analyzed female fertility traits included interval from calving to first insemination, interval from first to last insemination, days open, conception rate at first insemination, number of inseminations per conception and non-return rates within 56 days after first insemination.

RESULTS

The descriptive statistics showed that the average fertility of heifers was superior to that of cows. Moreover, the genetic correlations between the performances of a trait in heifers and in cows were all moderate to high but far from one, which suggested that the performances of a trait in heifers and cows should be considered as different but genetically correlated traits in genetic evaluations. On the other hand, genetic correlations between performances of a trait in different parities of cows were greater than 0.87, with only a few exceptions, but variances were not homogeneous across parities for some traits. The estimated heritabilities of female fertility traits were low; all were below 0.049 (except for interval from calving to first insemination). Additionally, the heritabilities of the heifer interval traits were lower than those of the corresponding cow interval traits. Moreover, the heritabilities of the interval traits were higher than those of the threshold traits when measuring similar fertility functions. In general, estimated genetic correlations between traits were highly consistent with the biological categories of the female fertility traits.

CONCLUSIONS

Interval from calving to first insemination, interval from first to last insemination and non-return rates within 56 days after first insemination are recommended to be included in the selection index of the Chinese Holstein population. The parameters estimated in the present study will facilitate the development of a genetic evaluation system for female fertility traits to improve the reproduction efficiency of Chinese Holsteins.