Between-herd variation in resilience and relations to herd performance

Are you ready for a challenge? Personality traits influence dairy calves' responses to disease, pain, and nutritional challenges

Dairy calves routinely experience disease, pain, and nutritional stressors such as diarrhea, dehorning, and weaning early in life. These stressors lead to changes in behavioral expression that varies in magnitude between individuals, where a greater magnitude change would suggest lower resilience in individuals to a stressor. Thus, this study first aimed to quantify the individual variation in magnitude change in feeding behaviors and activity in response to a bout of diarrhea, dehorning, and weaning. The next objective was to then investigate if personality traits were related to this magnitude of behavioral response in dairy calves, and thus their resilience toward these stressors. Calves were followed with 2 precision livestock technologies (e.g., an automatic feeding system, and leg accelerometer) to track behavioral changes in response during the time when the stressors were present. The automatic feeding system provided daily measures of milk intake, drinking speed, rewarded and unrewarded visits to the milk feeding station, and calf starter intake. The leg accelerometer provided daily measures of steps, activity index, lying time, and lying bouts. At 23 ± 3 d of age, Holstein dairy calves (n = 49) were subjected to a series of standardized personality tests that exposed the calf to novelty and fear stimuli. Factors extracted from a principal component analysis on the behaviors from the personality test were interpreted as personality traits: Factor 1 (fearful), Factor 2 (active) and Factor 3 (explorative). The magnitude changes in behaviors from the precision livestock technologies were calculated relative to the behavior performed on the day the stressor occurred (i.e., day of diagnosis, day of dehorning, day weaned). Linear regression models were used to determine whether calf scores on each factor were associated with magnitude change in behavior for each of the stressor periods with day relative to the stressor included as a repeated measure. Models were run independently for the period leading up to and following each stressor. We found that calves varied in their behavioral responses to diarrhea, dehorning, and weaning stressors, despite being reared in the same environment and experiencing consistent management procedures. Additionally, personality traits measured from standardized tests were associated with both the direction and magnitude of change in behaviors around each stressor. For instance, with diarrhea, calves that were highly fearful had a greater magnitude change in milk intake and drinking speed following diagnosis than the least fearful calves. With dehorning, calves that were highly explorative had a greater magnitude change in lying time when dehorned, but a smaller magnitude change in lying bouts and drinking speed following dehorning, compared with the least explorative calves. With weaning, calves that were highly active had a smaller magnitude change in unrewarded visits leading up to and following weaning than calves that were the least active. Each of the personality traits had a significant association with change in behavior surrounding each of the stressors evaluated, although these associations depended on the type of stressor. These results have implications for how individual calves experience each stressor and therefore individual animal welfare.

The genetics of resilience and its relationships with egg production traits and antibody traits in chickens

BACKGROUND

Resilience is the capacity of an animal to be minimally affected by disturbances or to rapidly return to its initial state before exposure to a disturbance. Resilient livestock are desired because of their improved health and increased economic profit. Genetic improvement of resilience may also lead to trade-offs with production traits. Recently, resilience indicators based on longitudinal data have been suggested, but they need further evaluation to determine whether they are indeed predictive of improved resilience, such as disease resilience. This study investigated different resilience indicators based on deviations between expected and observed egg production (EP) by exploring their genetic parameters, their possible trade-offs with production traits, and their relationships with antibody traits in chickens.

METHODS

Egg production in a nucleus breeding herd environment based on 1-week-, 2-week-, or 3-week-intervals of two purebred chicken lines, a white egg-laying (33,825 chickens) and a brown egg-laying line (34,397 chickens), were used to determine deviations between observed EP and expected average batch EP, and between observed EP and expected individual EP. These deviations were used to calculate three types of resilience indicators for two life periods of each individual: natural logarithm-transformed variance (ln(variance)), skewness, and lag-one autocorrelation (autocorrelation) of deviations from 25 to 83 weeks of age and from 83 weeks of age to end of life. Then, we estimated their genetic correlations with EP traits and with two antibody traits.

RESULTS

The most promising resilience indicators were those based on 1-week-intervals, as they had the highest heritability estimates (0.02-0.12) and high genetic correlations (above 0.60) with the same resilience indicators based on longer intervals. The three types of resilience indicators differed genetically from each other, which indicates that they possibly capture different aspects of resilience. Genetic correlations of the resilience indicator traits based on 1-week-intervals with EP traits were favorable or zero, which means that trade-off effects were marginal. The resilience indicator traits based on 1-week-intervals also showed no genetic correlations with the antibody traits, which suggests that they are not informative for improved immunity or vice versa in the nucleus environment.

CONCLUSIONS

This paper gives direction towards the evaluation and implementation of resilience indicators, i.e. to further investigate resilience indicator traits based on 1-week-intervals, in breeding programs for selecting genetically more resilient layer chickens.

Predicting dairy herd resilience on farms with conventional milking systems

This research paper addresses the problem that, thus far, there is no method available to predict herd resilience for farms that do not use automated milking systems (AMS). Recently, a methodology was developed to estimate both individual cow as well as herd resilience using daily milk yield observations at individual cow level from farms with AMS. This AMS-based method, however, is not suitable on farms that use conventional milking systems (CMS) where such individual cow milk yield observations are lacking. Therefore, this research aimed at predicting herd resilience using herd performance data that is commonly available on CMS farms. To do so, data consisting of 585 Dutch AMS farms where herd resilience estimates using the AMS-based method were available was examined. To predict herd resilience with herd performance data, only those data that are also commonly available on CMS farms were used in a 5-fold cross validation Random Forest model. These herd resilience estimates were subsequently compared with the AMS-based herd resilience estimates. Results showed that it is possible to predict with a 69.9% probability whether a herd performs with above or below average herd resilience using only variables available on CMS farms. Especially, the proportion of cows with an indication of rumen acidosis, proportion of cows with an elevated somatic cell count and the fluctuation in herd size over the years are good predictors of herd resilience. Since herd management decisions appear to affect herd resilience, a lower predicted herd resilience could be taken as a general indication that tactical or strategic management changes could be taken to improve the herd resilience.

Estimates of dairy herd health indicators of mastitis, ketosis, inter-calving interval, and fresh cow replacement in the Piedmont region, Italy

Test-day milk analysis has largely been used to study health and performance parameters in dairy cows. In this study, we estimated four health indicators of dairy cows using test-day data. Our purpose was to estimate (1) mastitis incidence rate, prevalence, and the probability of recovery; (2) the incidence proportion of ketosis; (3) the duration of inter-calving interval; and (4) the risk of a fresh cow being replaced, in a large cohort of dairy herds in the Piedmont region (Italy). We retrospectively analysed test day records of 261,121 lactating cows and 1315 herds during five years (2015-2020). Mastitis was defined by somatic cell count and ketosis by fat-to-protein ratio. Calving dates were used to calculate ICI and to estimate the removal of a fresh cow from the herd. Mixed-effect generalized linear models were used to adjust for unmeasured herd-level risk factors. The risk of mastitis increased by 120% with parity (Odds ratio [OR] = 2.20, confidence interval [CI]: 2.17 - 2.23), by 7% by months in milking (OR = 1.07, CI: 1.07 - 1.07), and even more if the cow was already affected during the same lactation (OR = 8.74, CI: 8.67 - 8.82). Lactose concentration on the previous test day was the best positive prognostic factor for mastitis recovery (OR = 1.12, CI: 1.08 - 1.17). Ketosis risk was the highest between 3rd and 4th lactations and itself increased the risk of having ICI longer than 440 days (OR = 1.12, CI: 1.02 - 1.22), and fresh-cow removal (OR = 1.75, CI: 1.58 - 1.93). Also, the removal of fresh cows was more likely when mastitis (OR = 1.31, CI: 1.19 - 1.45) or long ICI (OR = 1.34, CI: 1.22 - 1.48) occurred. For each health indicator, herd-level risk factors had an important role (18-56% of within-herd covariance). Our results indicate that milk analysis could be also useful for predicting mastitis, its cure rate, and ketosis. Cow-level risk factors are not enough to explain the risk of these issues. By studying a large population over a long period, this study provides an updated estimate of dairy cow health indicators in Piedmont (north-western Italy), useful for benchmarking dairy herds.

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.

Observational study on associations between resilience indicators based on daily milk yield in first lactation and lifetime profitability

Resilience is the ability of cows to be minimally affected by disturbances, such as pathogens, heat waves, and changes in feed quality, or to quickly recover. Obvious advantages of resilience are good animal welfare and easy and pleasant management for farmers. Furthermore, economic effects are also expected, but these remain to be determined. The goal of this study was to investigate the association between resilience and lifetime gross margin, using indicators of resilience calculated from fluctuations in daily milk yield using an observational study. Resilience indicators and lifetime gross margin were calculated for 1,325 cows from 21 herds. These cows were not alive anymore and, therefore, had complete lifetime data available for many traits. The resilience indicators were the natural log-transformed variance (LnVar) and the lag-1 autocorrelation (r_auto) of daily milk yield deviations from cow-specific lactation curves in parity 1. Good resilience is indicated by low LnVar (small yield response to disturbances) and low r_auto (quick yield recovery to baseline). Lifetime gross margin was calculated as the sum of all revenues minus the sum of all costs throughout life. Included revenues were from milk, calf value, and slaughter of the cow. Included costs were from feed, rearing, insemination, management around calving, disease treatments, and destruction in case of death on farm. Feed intake was unknown and, therefore, lifetime feed costs had to be estimated based on milk yield records. The association of each resilience indicator with lifetime gross margin, and also with the underlying revenues and costs, was investigated using analysis of covariance (ANCOVA) models. Mean daily milk yield in first lactation, herd, and year of birth were included as covariates and factors. Natural log-transformed variance had a significantly negative association with lifetime gross margin, which means that cows with stable milk yield (low LnVar, good resilience) in parity 1 generated on average a higher lifetime gross margin than cows that had the same milk yield level but with more fluctuations. The association with lifetime gross margin could be mainly attributed to higher lifetime milk revenues for cows with low LnVar, due to a longer lifespan. Unlike LnVar, r_auto was not significantly associated with lifetime gross margin or any of the underlying lifetime costs and revenues. However, it was significantly associated with yearly treatment costs, which is important for ease of management. In conclusion, the importance of resilience for total profit generated by a cow at the end of life was confirmed by the significant association of LnVar with lifetime gross margin, although effects of differences in feed efficiency between resilient and less resilient cows remain to be studied. The economic advantage can be mainly ascribed to benefits of long lifespan.

Using egg production longitudinal recording to study the genetic background of resilience in purebred and crossbred laying hens

Background

There is growing interest in using genetic selection to obtain more resilient farm animals (i.e. that are minimally affected by disturbances or rapidly recover from them). The aims of this study were to: (i) estimate the genetic parameters of resilience indicator traits based on egg production data, (ii) assess whether these traits are genetically correlated in purebreds and crossbreds, and (iii) assess the genetic correlations of these traits with egg production (EP) as total number of eggs between 25 and 83 weeks. Purebred hens (33,825 from a White Leghorn (WA) line and 34,397 from a Rhode Island (BD) line were housed in individual cages, while crossbred hens were housed in collective cages of 6 to 8 paternal half-sibs (12,852 WA and 3898 BD crossbred groups, where the name of the group refers to the line used as the sire). Deviations of a hen’s weekly egg production from the average of the corresponding batch were calculated. Resilience indicator traits investigated were the natural logarithm of the variance (LNVAR), the skewness (SKEW), and the lag-one autocorrelation (AUTO-R) of these deviations.

Results

In both purebred lines, EP was estimated to be lowly heritable (WA: 0.11 and BD: 0.12). Resilience indicators were also estimated to be lowly heritable in both lines (LNVAR: 0.10 and 0.12, SKEW: 0.04 and 0.02, AUTO-R: 0.06 and 0.08 in WA and BD, respectively). In both crossbred groups, EP, AUTO-R, and SKEW were estimated to be less heritable than in purebreds (EP: \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h^{2}$$\end{document}h2 ≤ 0.07; and resilience indicator traits: \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h^{2}$$\end{document}h2 ≤ 0.03), while LNVAR had an \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h^{2}$$\end{document}h2 estimate that was similar to or higher in crossbreds (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h^{2}$$\end{document}h2 ranged from 0.13 to 0.21) than in purebreds. In both purebreds and crossbreds, resilience indicator traits were estimated to have favorable genetic correlations with EP and between each other. For all traits and in both lines, estimates of genetic correlations between purebreds and crossbreds (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r_{pc}$$\end{document}rpc) differed from 1 and ranged from 0.16 to 0.63.

Conclusions

These results show that selection for resilience based on EP data can be considered in breeding programs for layers. Genetic improvement of resilience in crossbreds can be achieved by using information on purebreds, but would be greatly enhanced by the integration of information on crossbreds in breeding programs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00716-8.