Emphasis on resilience in dairy cattle breeding: Possibilities and consequences

A first characterization of the microbiota-resilience link in swine

BACKGROUND

The gut microbiome plays a crucial role in understanding complex biological mechanisms, including host resilience to stressors. Investigating the microbiota-resilience link in animals and plants holds relevance in addressing challenges like adaptation of agricultural species to a warming environment. This study aims to characterize the microbiota-resilience connection in swine. As resilience is not directly observable, we estimated it using four distinct indicators based on daily feed consumption variability, assuming animals with greater intake variation may face challenges in maintaining stable physiological status. These indicators were analyzed both as linear and categorical variables. In our first set of analyses, we explored the microbiota-resilience link using PERMANOVA, α-diversity analysis, and discriminant analysis. Additionally, we quantified the ratio of estimated microbiota variance to total phenotypic variance (microbiability). Finally, we conducted a Partial Least Squares-Discriminant Analysis (PLS-DA) to assess the classification performance of the microbiota with indicators expressed in classes.

RESULTS

This study offers four key insights. Firstly, among all indicators, two effectively captured resilience. Secondly, our analyses revealed robust relationship between microbial composition and resilience in terms of both composition and richness. We found decreased α-diversity in less-resilient animals, while specific amplicon sequence variants (ASVs) and KEGG pathways associated with inflammatory responses were negatively linked to resilience. Thirdly, considering resilience indicators in classes, we observed significant differences in microbial composition primarily in animals with lower resilience. Lastly, our study indicates that gut microbial composition can serve as a reliable biomarker for distinguishing individuals with lower resilience.

CONCLUSION

Our comprehensive analyses have highlighted the host-microbiota and resilience connection, contributing valuable insights to the existing scientific knowledge. The practical implications of PLS-DA and microbiability results are noteworthy. PLS-DA suggests that host-microbiota interactions could be utilized as biomarkers for monitoring resilience. Furthermore, the microbiability findings show that leveraging host-microbiota insights may improve the identification of resilient animals, supporting their adaptive capacity in response to changing environmental conditions. These practical implications offer promising avenues for enhancing animal well-being and adaptation strategies in the context of environmental challenges faced by livestock populations. Video Abstract.

Unraveling the genomic background of resilience based on variability in milk yield and milk production levels in North American Holstein cattle through genome-wide association study and Mendelian randomization analyses

Breeding more resilient animals will benefit the dairy cattle industry in the long term, especially as global climate changes become more severe. Previous studies have reported genetic parameters for various milk yield-based resilience indicators, but the underlying genomic background of these traits remain unknown. In this study, we conducted GWAS of 62,029 SNPs with 4 milk yield-based resilience indicators, including the weighted occurrence frequency (wfPert) and accumulated milk losses (dPert) of milk yield perturbations, and log-transformed variance (LnVar) and lag-1 autocorrelation (rauto) of daily yield residuals. These variables were previously derived from 5.6 million daily milk yield records from 21,350 lactations (parities 1-3) of 11,787 North American Holstein cows. The average daily milk yield (ADMY) throughout lactation was also included to compare the shared genetic background of resilience indicators with milk yield. The differential genetic background of these indicators was first revealed by the significant genomic regions identified and significantly enriched biological pathways of positional candidate genes, which confirmed the genetic difference among resilience indicators. Interestingly, the functional analyses of candidate genes suggested that the regulation of intestinal homeostasis is most likely affecting resilience derived based on variability in milk yield. Based on Mendelian randomization analyses of multiple instrumental SNPs, we further found an unfavorable causal association of ADMY with LnVar. In conclusion, the resilience indicators evaluated are genetically different traits, and there are causal associations of milk yield with some of the resilience indicators evaluated. In addition to providing biological insights into the molecular regulation mechanisms of resilience derived based on variability in milk yield, this study also indicates the need for developing selection indexes combining multiple indicator traits and taking into account their genetic relationship for breeding more resilient dairy cattle.

Combined mechanistic and genetic modelling to benchmark body reserve traits as proxies of dairy cows' lifetime efficiency in grass-based production systems

Improving lifetime lactation efficiency of dairy cows by selection is difficult due to the complexity of this trait and the existence of genotype-by-environment interactions. This study aimed at assessing the relevance of traits derived from body reserves as lifetime efficiency indicators under contrasting nutritional environments. Given the absence of large-scale datasets covering a panel of feeding regimes, phenotypes were simulated for populations of 20 000 dairy cows using a mechanistic bioenergetic model. Ten phenotypes were computed for third-lactation cows. Analysed phenotypes comprised total milk production, lactation efficiency, BW at calving (BWcalv), DM intake (DMI) and interval between first insemination and conception. Five traits described levels and changes of body reserves at different periods during lactation. Lifetime lactation efficiency was computed for all cows (Life_Eff). Three nutritional environments were defined considering a grass-based production system with seasonal calving: a high non-limiting scenario (HS) mimicking ad libitum access to feed and two limiting environments with moderate (MS) and low (LS) feed offer. Variance components were estimated for all traits within and between environments using REML. Heritabilities estimated for milk production, lactation efficiency, BWcalv and DMI were moderate in the different environments (0.27-0.35 ± 0.04). The heritability of body reserve levels and dynamics were moderate in the HS and MS scenarios (0.23-0.30 ± 0.03) and lower in the LS scenario (0.14-0.25 ± 0.03). The heritability of Life_Eff was low in the HS environment (0.07 ± 0.01) and slightly increased in the limiting environments. All genetic correlations estimated between environments were moderate to high (≥0.66 ± 0.07), suggesting low to moderate genotype-by-environment interactions. Estimated genetic correlations were moderate between Life_Eff and body reserve levels (from 0.39 to 0.51 ± 0.08) and moderate but negative between Life_Eff and change in body reserves traits (-0.27 to -0.37 ± 0.09) in the HS environment. The genetic correlations between Life_Eff and body reserve levels increased to higher values in the limiting environments. In contrast, genetic correlations between Life_Eff and the changes in body reserves were closer to zero. In conclusion, this study showed that body reserve levels were relevant proxies of lifetime irrespective of the environment. In contrast, changes in body reserves that reflected energy mobilisation in early lactation were less informative about lifetime efficiency in environments with severe feed restrictions.

Carbon footprint based on lifetime productivity for future cows selected for resilience to climate-related disturbances

Droughts, which can affect feed production, are projected to become more common under future climate conditions. In light of this, breeding cattle resilient to changes in feeding regimens is increasingly becoming an important topic. Body reserves can play a crucial role when feed resources are limited. We simulated populations of dairy cows selected with 2 different breeding goals: one reflecting the current breeding goal and the other placing weight on minimum level of body reserves in early lactation and change in body reserves during lactation. We considered this latter as a breeding goal for resilience. We used the 2 dynamic simulation programs ADAM and AQAL to predict performance of the cows after selection. In AQAL, we modeled moderate and severe drought by decreasing feed quality and quantity offered to cows during one year. We compared cows selected with the 2 breeding goals under 3 environments: without disturbances related to climate and with moderate and severe drought. In the environments without disturbances and the moderate drought, the cows selected with the current breeding goal had higher lifetime lactation efficiency (energy invested in milk/energy acquired from feed) and lower carbon footprint per kilogram of protein in milk and meat than cows selected for resilience. However, with severe drought, cows selected for resilience had higher lifetime lactation efficiency and lower carbon footprint per kilogram of protein in milk and meat than those selected with the current breeding goal. This suggests that cows selected for high productive performance do not perform well under very limiting conditions, leading to increased climate impact. The importance of inclusion of body reserves as a resilience trait in dairy cattle breeding depends on the future environment in which the cows will be used.

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.

Relationship between feed efficiency and resilience in dairy ewes subjected to acute underfeeding

Selection of dairy sheep based on production levels has caused a loss of rusticity, which might compromise their future resilience to nutritional challenges. Although refocusing breeding programs toward improved feed efficiency (FE) is expected, more-efficient ewes also seem to be more productive. As a first step to examine the relationship between FE and resilience in dairy sheep, in this study we explored the variation in the response to and the recovery from an acute nutritional challenge in high-yielding Assaf ewes phenotypically divergent for FE. First, feed intake, milk yield and composition, and body weight changes were recorded individually over a 3-wk period in a total of 40 sheep fed a total mixed ration (TMR) ad libitum. Data were used to calculate their FE index (FEI, defined as the difference between the actual and predicted intake estimated through net energy requirements for maintenance, production, and weight change). The highest and lowest FE ewes (H-FE and L-FE groups, respectively; 10 animals/group) were selected and then subjected to the nutritional challenge (i.e., withdrawing the TMR and limiting their diet only to the consumption of straw for 3 d). Afterward, sheep were fed again the TMR ad libitum. Temporal patterns of variation in performance traits, and ruminal fermentation and blood parameters were examined. A good consistency between FEI, residual feed intake, and feed conversion ratio was observed. Results supported that H-FE were more productive than L-FE sheep at similar intake level. Average time trends of milk yield generated by a piecewise model suggest that temporal patterns of variation in this trait would be related to prechallenge production level (i.e., H-FE presented quicker response and recovery than L-FE). Considering all studied traits, the overall response to and recovery from underfeeding was apparently similar or even better in H-FE than in L-FE. This would refute the initial hypothesis of a poorer resilience of more-efficient sheep to an acute underfeeding. However, the question remains whether a longer term feed restriction might impair the ability of H-FE ewes to maintain or revert to a high-production status, which would require further research.

Genomic-based genetic parameters for resilience across lactations in North American Holstein cattle based on variability in daily milk yield records

Considering the increasing challenges imposed by climate change and the need to improve animal welfare, breeding more resilient animals capable of better coping with environmental disturbances is of paramount importance. In dairy cattle, resilience can be evaluated by measuring the longitudinal occurrences of abnormal daily milk yield throughout lactation. Aiming to estimate genetic parameters for dairy cattle resilience, we collected 5,643,193 daily milk yield records on automatic milking systems (milking robots) and milking parlors across 21,350 lactations 1 to 3 of 11,787 North American Holstein cows. All cows were genotyped with 62,029 SNPs. After determining the best fitting models for each of the 3 lactations, daily milk yield residuals were used to derive 4 resilience indicators: weighted occurrence frequency of yield perturbations (wfPert), accumulated milk losses of yield perturbations (dPert), and log-transformed variance (LnVar) and lag-1 autocorrelation (r_auto) of daily yield residuals. The indicator LnVar presented the highest heritability estimates (±standard error), ranging from 0.13 ± 0.01 in lactation 1 to 0.15 ± 0.02 in lactation 2; the other 3 indicators had relatively lower heritabilities across the 3 lactations (0.01-0.06). Based on bivariate analyses of each resilience indicator across lactations, stronger genetic correlations were observed between lactations 2 and 3 (0.88-0.96) than between lactations 1 and 2 or 3 (0.34-0.88) for dPert, LnVar, and r_auto. For the pairwise comparisons of different resilience indicators within each lactation, dPert had the strongest genetic correlations with wfPert (0.64) and r_auto (0.53) in lactation 1, whereas the correlations in lactations 2 and 3 were more variable and showed relatively high standard errors. The genetic correlation results indicated that different resilience indicators across lactations might capture additional biological mechanisms and should be considered as different traits in genetic evaluations. We also observed favorable genetic correlations of these resilience indicators with longevity and Net Merit index, but further biological validation of these resilience indicators is needed. In conclusion, this study provided genetic parameter estimates for different resilience indicators derived from daily milk yields across the first 3 lactations in Holstein cattle, which will be useful when potentially incorporating these traits in dairy cattle breeding schemes.