Genetic parameters and evaluation of mortality and slaughter rate in Holstein and Jersey cows

The longevity of dairy cattle has economic, animal welfare, and health implications and is influenced by the frequency of mortality on the farm and sale for slaughter. In this study cows removed from the herd due to death or slaughter during the lactation were coded 1 and cows that were not terminated were coded 0. Genetic parameters for mortality rates (MR) and slaughter rates (SR) were estimated for Holstein (H) and Jersey (J) breeds by applying both linear (LM) and threshold (TM) sire models using about 1.2 million H and 286,000 J cows. Estimated breeding values (EBV) for MR and SR were predicted using animal models to assess the opportunity for selection and genetic trends. Cow termination data, recorded between 1990 and 2020 on a voluntary basis by Australian dairy farmers, were analyzed. Cow MR has increased from below 1% in the 1990s to 4.1% and 3.6% in recent years in H and J cows, respectively. Most dead cows (∼36%) left the herd before 120 d of lactation, while cows that were slaughtered left the herd toward the end of the lactation. Using the LM, heritability (h2) estimates for MR were lower (1%) than those for SR (2%-3.5%). When h2 were estimated using a TM, the estimates for both traits varied between 4% and 20%, suggesting that the difference in incidence level is one of the reasons for the difference in the h2 values between MR and SR. Early test-day milk yield (MY) and 305-d MY (305-d MY) have unfavorable genetic correlations (0.32-0.41) with MR in both breeds. The genetic correlations of calving interval with MR were stronger (0.54-0.68) than with SR (0.28-0.45) suggesting that poor fertility can serve as an early indicator of poor cow health that may lead to increased risk of death. High early test-day somatic cell count is genetically associated with increased likelihood of slaughter (0.24-0.46), but not with increased likelihood of death. In H, 305-d protein yield (PY) had the strongest genetic correlation (-0.34 to -0.40) with SR whereas in J, both 305-d PY and fat yield showed high genetic (-0.64 to -0.70) and moderate environmental (-0.35 to -0.37) correlations with SR. The genetic correlation of removal from the herd due to death and slaughter was negative (-0.3) in J and zero in H. Strong selection for improved fertility and survival and less selection emphasis for MY, has led to an improvement in the genetic trend for cow MR in H and the trend in J has stabilized. Although genetic evaluations for cow MR are feasible, the reliabilities of the EBV are low and the level of cow MR in Australia are relatively low compared with similar countries. Therefore, genetic evaluation for survival based on mortality and slaughter data could be sufficient in the current selection circumstances where breeding objectives are broadly defined. Nevertheless, all Australian farmers should be encouraged to continue recording mortality and slaughter data for monitoring of the trends and for future development of genetic evaluations.

Insight into selective breeding for robustness based on field survival records: New genetic evaluation of survival traits in pacific white shrimp (Penaeus vannamei) breeding line

Survival can be considered a relatively 'old' trait in animal breeding, yet commonly neglected in aquaculture breeding because of the simple binary records and generally low heritability estimates. Developing routine genetic evaluation systems for survival traits however, will be important for breeding robust strains based on valuable field survival data. In the current study, linear multivariate animal model (LMA) was used for the genetic analysis of survival records from 2-year classes (BL2019 and BL2020) of pacific white shrimp (Penaeus vannamei) breeding lines with data collection of 52, 248 individuals from 481 fullsib families. During grow-out test period, 10 days intervals of survival data were considered as separate traits. Two survival definitions, binary survivability (S) and continuous survival in days (SL), were used for the genetic analysis of survival records to investigate; 1) whether adding more survival time information could improve estimation of genetic parameters; 2) the trajectory of survival heritability across time, and 3) patterns of genetic correlations of survival traits across time. Levels of heritability estimates for both S and SL were low (0.005-0.076), while heritability for survival day number was found to be similar with that of binary records at each observation time and were highly genetically correlated (r g > 0.8). Heritability estimates of body weight (BW) for BL2019 and BL2020 were 0.486 and 0.373, respectively. Trajectories of survival heritability showed a gradual increase across the grow-out test period but slowed or reached a plateau during the later grow-out test period. Genetic correlations among survival traits in the grow-out tests were moderate to high, and the closer the times were between estimates, the higher were their genetic correlations. In contrast, genetic correlations between both survival traits and body weight were low but positive. Here we provide the first report on the trajectory of heritability estimates for survival traits across grow-out stage in aquaculture. Results will be useful for developing robust improved pacific white shrimp culture strains in selective breeding programs based on field survival data.

Analysis of culling reasons and age at culling in Australian dairy cattle

Context A thorough analysis of the reasons for culling was made to understand the phenotypic trend in herd life. In addition, identification of culling reasons could enable to develop a strategy for further evaluation of longevity in Australian dairy cows. Aims The aim of this study was to investigate the main causes of culling in Australian dairy herds and thereby to assess the trend of reason-specific culling over time. Methods Culling reasons in Australian dairy cattle were studied based on culling records from 1995 through 2016. A total of 2452124 individual cow culling observations were obtained from Datagene, Australia, of which 2140337 were Holstein and 311787 were from Jersey cows. A binary logistic regression model was used to estimate effects of breed and age and the trend of a particular culling reason over time. Key results The most important culling reasons identified over the 21-year period were infertility (17.0%), mastitis (12.9%), low production (9.3%), sold for dairy purpose (6.4%) and old age (6.2%), whereas 37.4% were ‘other reasons not reported’. The average age at culling was nearly the same for Holstein (6.75 years) and Jersey (6.73 years) cows. The estimated age at culling was slightly increased for Holstein cows (by 3.7 days) and somewhat decreased for Jersey cows (by 11 days) over the last two decades. The probability of culling cows for infertility and low production was high in early parities and consistently declined as age advanced, and culling due to mastitis was higher in older cows. The trend of main culling reasons over time was evaluated, indicating that the probability of culling due to infertility has progressively increased over the years in both breeds, and culling for mastitis in Jersey cows has also increased. Culling of cows due to low production sharply decreased from 2.5 to –8% for Holstein and from 73 to 60% for Jersey cows over the 21-year period. Conclusions Culling age has changed only little in both breeds whereas culling reasons have changed over the last two decades, with low production becoming a less important reason for culling and infertility becoming more important for Holstein and Jersey breeds. Implications Due to changes of culling reasons, there could be a change in the meaning of survival over time as well. As a result, genetic correlation with survival and other traits might be changed and accuracy and bias of genetic evaluations could be affected.

Beef trait genetic parameters based on old and recent data and its implications for genomic predictions in Italian Simmental cattle

This study aimed to evaluate the changes in variance components over time to identify a subset of data from the Italian Simmental (IS) population that would yield the most appropriate estimates of genetic parameters and breeding values for beef traits to select young bulls. Data from bulls raised between 1986 and 2017 were used to estimate genetic parameters and breeding values for four beef traits (average daily gain [ADG], body size [BS], muscularity [MUS], and feet and legs [FL]). The phenotypic mean increased during the years of the study for ADG, but it decreased for BS, MUS, and FL. The complete dataset (ALL) was divided into four generational subsets (Gen1, Gen2, Gen3, and Gen4). Additionally, ALL was divided into two larger subsets: the first one (OLD) combined data from Gen1 and Gen2 to represent the starting population, and the second one (CUR) combined data from Gen3 and Gen4 to represent a subpopulation with stronger ties to the current population. Genetic parameters were estimated with a four-trait genomic animal model using a single-step genomic average information restricted maximum likelihood algorithm. Heritability estimates from ALL were 0.26 ± 0.03 for ADG, 0.33 ± 0.04 for BS, 0.55 ± 0.03 for MUS, and 0.23 ± 0.03 for FL. Higher heritability estimates were obtained with OLD and ALL than with CUR. Considerable changes in heritability existed between Gen1 and Gen4 due to fluctuations in both additive genetic and residual variances. Genetic correlations also changed over time, with some values moving from positive to negative or even to zero. Genetic correlations from OLD were stronger than those from CUR. Changes in genetic parameters over time indicated that they should be updated regularly to avoid biases in genomic estimated breeding values (GEBV) and low selection accuracies. GEBV estimated using CUR variance components were less biased and more consistent than those estimated with OLD and ALL variance components. Validation results indicated that data from recent generations produced genetic parameters that more appropriately represent the structure of the current population, yielding accurate GEBV to select young animals and increasing the likelihood of higher genetic gains.

Changes in genetic parameters for fitness and growth traits in pigs under genomic selection

Genomic selection increases accuracy and decreases generation interval, speeding up genetic changes in the populations. However, intensive changes caused by selection can reduce the genetic variation and can strengthen undesirable genetic correlations. The purpose of this study was to investigate changes in genetic parameters for fitness traits related with prolificacy (FT1) and litter survival (FT2 and FT3), and for growth (GT1 and GT2) traits in pigs over time. The data set contained 21,269 (FT1), 23,246 (FT2), 23,246 (FT3), 150,492 (GT1), and 150,493 (GT2) phenotypic records obtained from 2009 to 2018. The pedigree file included 369,776 animals born between 2001 and 2018, of which 39,103 were genotyped. Genetic parameters were estimated with bivariate models (FT1-GT1, FT1-GT2, FT2-GT1, FT2-GT2, FT3-GT1, and FT3-GT2) using 3-yr sliding subsets. With a Bayesian implementation using the GIBBS3F90 program computations were performed as genomic analysis (GEN) or pedigree-based analysis (PED), that is, with or without genotypes, respectively. For GEN (PED), the changes in heritability from the first to the last year interval, that is, from 2009–2011 to 2015–2018 were 8.6 to 5.6 (7.9 to 8.8) for FT1, 7.8 to 7.2 (7.7 to 10.8) for FT2, 11.4 to 7.6 (10.1 to 7.5) for FT3, 35.1 to 16.5 (32.5 to 23.7) for GT1, and 35.9 to 16.5 (32.6 to 24.1) for GT2. Differences were also observed for genetic correlations as they changed from −0.31 to −0.58 (−0.28 to −0.73) for FT1-GT1, −0.32 to −0.50 (−0.29 to −0.74) for FT1-GT2, −0.27 to −0.45 (−0.30 to −0.65) for FT2-GT1, −0.28 to −0.45 (−0.32 to −0.66) for FT2-GT2, 0.14 to 0.17 (0.11 to 0.04) for FT3-GT1, and 0.14 to 0.18 (0.11 to 0.05) for FT3-GT2. Strong selection in pigs reduced heritabilities and emphasized the antagonistic genetic relationships between fitness and growth traits. With genotypes considered, heritability estimates were smaller and genetic correlations were greater than estimates with only pedigree and phenotypes. When selection is based on genomic information, genetic parameters estimated without this information can be biased because preselection is not accounted for by the model.

Predicting survival in dairy cattle by combining genomic breeding values and phenotypic information

Advances in technology and improved data collection have increased the availability of genomic estimated breeding values (gEBV) and phenotypic information on dairy farms. This information could be used for the prediction of complex traits such as survival, which can in turn be used in replacement heifer management. In this study, we investigated which gEBV and phenotypic variables are of use in the prediction of survival. Survival was defined as survival to second lactation, plus 2 wk, a binary trait. A data set was obtained of 6,847 heifers that were all genotyped at birth. Each heifer had 50 gEBV and up to 62 phenotypic variables that became gradually available over time. Stepwise variable selection on 70% of the data was used to create multiple regression models to predict survival with data available at 5 decision moments: distinct points in the life of a heifer at which new phenotypic information becomes available. The remaining 30% of the data were kept apart to investigate predictive performance of the models on independent data. A combination of gEBV and phenotypic variables always resulted in the model with the highest Akaike information criterion value. The gEBV selected were longevity, feet and leg score, exterior score, udder score, and udder health score. Phenotypic variables on fertility, age at first calving, and milk quantity were important once available. It was impossible to predict individual survival accurately, but the mean predicted probability of survival of the surviving heifers was always higher than the mean predicted probability of the nonsurviving group (difference ranged from 0.014 to 0.028). The model obtained 2.0 to 3.0% more surviving heifers when the highest scoring 50% of heifers were selected compared with randomly selected heifers. Combining phenotypic information and gEBV always resulted in the highest scoring models for the prediction of survival, and especially improved early predictive performance. By selecting the heifers with the highest predicted probability of survival, increased survival could be realized at the population level in practice.