Genotype-by-environment interactions in beef and dairy cattle populations: A review of methodologies and perspectives on research and applications

Modern livestock production systems are characterized by a greater focus on intensification, involving managing larger numbers of animals to achieve higher productive efficiency and animal health and welfare within herds. Therefore, animal breeding programs need to be strategically designed to select animals that can effectively enhance production performance and animal welfare across a range of environmental conditions. Thus, this review summarizes the main methodologies used for assessing the levels of genotype-by-environment interaction (G × E) in cattle populations. In addition, we explored the importance of integrating genomic and phenotypic information to quantify and account for G × E in breeding programs. An overview of the structure of cattle breeding programs is provided to give insights into the potential outcomes and challenges faced when considering G × E to optimize genetic gains in breeding programs. The role of nutrigenomics and its impact on gene expression related to metabolism in cattle are also discussed, along with an examination of current research findings and their potential implications for future research and practical applications. Out of the 116 studies examined, 60 and 56 focused on beef and dairy cattle, respectively. A total of 83.62% of these studies reported genetic correlations across environmental gradients below 0.80, indicating the presence of G × E. For beef cattle, 69.33%, 24%, 2.67%, 2.67%, and 1.33% of the studies evaluated growth, reproduction, carcass and meat quality, survival, and feed efficiency traits, respectively. By contrast, G × E research in dairy cattle populations predominantly focused on milk yield and milk composition (79.36% of the studies), followed by reproduction and fertility (19.05%), and survival (1.59%) traits. The importance of G × E becomes particularly evident when considering complex traits such as heat tolerance, disease resistance, reproductive performance, and feed efficiency, as highlighted in this review. Genomic models provide a valuable avenue for studying these traits in greater depth, allowing for the identification of candidate genes and metabolic pathways associated with animal fitness, adaptation, and environmental efficiency. Nutrigenetics and nutrigenomics are emerging fields that require extensive investigation to maximize our understanding of gene-nutrient interactions. By studying various transcription factors, we can potentially improve animal metabolism, improving performance, health, and quality of products such as meat and milk.

Estimation of Genetic Parameters for Milk Production Rate and Its Stability in Holstein Population

Milk production rate (MPR) refers to the rate of milk secretion per hour (kg/h), calculated from the harvested milk yield and milking interval, and it is considered an appropriate measure to evaluate the production potential of cows. The objective of this study was to estimate the phenotypic and genetic parameters of milk production rate traits. In this study, the milking records of 4760 Holstein cows were collected, and four milk yield traits and six milk production rate traits were defined. The MIXED procedure was used to detect the impacts of non-genetic effects on milk yield and milk production rate traits, including parity, measured season and lactation stage. Variance and covariance components for milk yield and milk production rate traits were estimated using a univariate linear repeatability model. Parity, measurement season and lactation stage had significant effects (p < 0.01) on milk yield, milk production rate and its stability. Milk yield and milk production traits had high heritability, and ranged from 0.25 to 0.39. The stability of milk production rate had low heritability (0.04~0.05). Milk production rate is beneficial for the devolving novel trait in dairy breeding and provides new insights for herd management and genetic selection of production performance of dairy cattle.

Metabolic and productive adaptive response of beef cows to successive short-nutritional challenges

This study aimed to analyze the response of lactating beef cows to repeated short nutritional challenges with their performance parameters and plasma metabolites. Multiparous lactating beef cows were subjected to three repeated nutritional challenges in the fourth month of lactation. Each challenge consisted of a 4-d feed restriction (55% of their average energy and protein requirements), followed by a 3-d refeeding period (100% requirements). Cows were classified into two groups differing in their performance (milk yield) and metabolic adaptation [non esterified fatty acids (NEFA) and β-hydroxybutyrate (BHB)] to diet changes (metabolic response, MR): High and Low MR cows, where the High MR cows showed a faster and larger response to diet changes than the Low MR cows (P < 0.001). The loss in milk yield during restriction was the smallest in challenge 1 (P < 0.001). Milk urea increased during restriction in challenges 1 and 2 (P < 0.001). The High MR cows had greater NEFA concentrations than their Low MR counterparts during restrictions, and greater BHB concentrations during the restriction of challenge 2 (P < 0.001). Restriction increased NEFA, BHB (only in the High MR cows) and urea (P < 0.01). During refeeding, both milk yield and plasma metabolites recovered basal values (P > 0.05). These results highlight the ability of beef cows to respond to and recover from successive short-term nutrient restrictions, and that despite a certain degree of sensitization of milk yield may have occurred, there were only minimal changes in the metabolic strategies triggered to cope with repeated underfeeding.

Comparison of 3 mathematical models to estimate lactation performance in dairy cows

Milk yield dynamics and production performance reflect how dairy cows cope with their environment. To optimize farm management, time series of individual cow milk yield have been studied in the context of precision livestock farming, and many mathematical models have been proposed to translate raw data into useful information for the stakeholders of the dairy chain. To gain better insights on the topic, this study aimed at comparing 3 recent methods that allow one to estimate individual cow potential lactation performance, using daily data recorded by the automatic milking systems of 14 dairy farms (7 Holstein, 7 Italian Simmental) from Belgium, the Netherlands, and Italy. An iterative Wood model (IW), a perturbed lactation model (PLM), and a quantile regression (QR) were compared in terms of estimated total unperturbed (i.e., expected) milk production and estimated total milk loss (relative to unperturbed yield). The IW and PLM can also be used to identify perturbations of the lactation curve and were thus compared in this regard. The outcome of this study may help a given end-user in choosing the most appropriate method according to their specific requirements. If there is a specific interest in the post-peak lactation phase, IW can be the best option. If one wants to accurately describe the perturbations of the lactation curve, PLM can be the most suitable method. If there is need for a fast and easy approach on a very large dataset, QR can be the choice. Finally, as an example of application, PLM was used to analyze the effect of cow parity, calving season, and breed on their estimated lactation performance.

Trait development and genetic parameters of resilience indicators based on variability in milk consumption recorded by automated milk feeders in North American Holstein calves

The implementation of automated milk feeders (AMF) on precision dairy farms has enabled efficient management of large numbers of group-housed replacement calves with reduced labor requirements and improved calf welfare. In this study, we investigated the feasibility of deriving calf resilience indicators based on variability in milk consumption using data from 10,076 North American Holstein calves collected between 2015 and 2021. We modeled and evaluated deviations in observed and predicted daily milk consumption trajectories as indicators of resilience to environmental perturbations. We also analyzed average milk intake and the number of treatments for bovine respiratory disease (BRD) and their genetic correlations with the derived resilience parameters. Milk consumption was recorded using the Förster-Technik AMF. Deviations in cumulative milk intake were modeled using various methods, including quantile regression and the Gompertz function. Ten resilience indicators were derived to quantify the degree and duration of perturbations, including amplitude, perturbation time, recovery time, and deviation velocities. After data editing, genomic data from 9,273 calves and pedigree information from 10,076 calves with 321,388 phenotypic records were used to estimate genetic parameters for 12 traits, including 10 calf resilience indicators as well as average milk intake and treatments for bovine respiratory disease. Substantial phenotypic variability was observed for all calf resilience indicators derived and genetic parameters related to these novel resilience indicators were estimated. The heritability estimates for the resilience traits are as follows: amplitude of the deviation (in L) 0.047 (0.032, 0.064) (HPD interval), perturbation time of deviation (in d) 0.011 (0.0056, 0.016), recovery time of the deviation (in d) 0.025 (0.016, 0.035), maximum velocity of perturbation (L/d) 0.039 (0.024, 0.053), average velocity of perturbation (L/d) 0.038 (0.022, 0.050), area between the curves (L x d) 0.039 (0.027, 0.054), recovery ratio 0.053 (0.036, 0.072), deviation variance 0.049 (0.32, 0.068), log-deviation variance 0.027 (0.016, 0.044), deviation auto-correlation 0.010 (0.0042, 0.017) and number of deviation occurrences 0.023 (0.0094, 0.036). Some of the highlighted genetic correlations observed with average milk consumption include amplitude: 0.569 (0.474, 0.666), perturbation time: -0.534 (-0.73, -0.342), and average velocity: 0.554 (0.432, 0.672). Similarly, the genetic correlations between the number of times treated for BRD with perturbation time was 0.494 (0.251, 0.723), -0.294 (-0.52, -0.095) with number of deviations, and 0.348 (0.131, 0.578) with deviation autocorrelation. This study highlights the genetic influence on various resilience traits in calves, including amplitude, perturbation time, recovery time, and velocity measures of the perturbation. Our findings suggest the need for prioritizing genetic selection based on traits like recovery time, which exhibits higher heritability and a moderate genetic correlation with the number of times a calf is treated for BRD. The combination of AMF data, mathematical modeling, and genomic evaluation provides a comprehensive framework for assessing and breeding more resilient dairy calves in the face of environmental and health challenges.

Validating statistical properties of resilience indicators derived from simulated longitudinal performance measures of farmed animals

Resilience is commonly defined as the ability of an individual to be minimally affected or to quickly recover from a challenge. Improvement of animals' resilience is a vital component of sustainable livestock production but has so far been hampered by the lack of established quantitative resilience measures. Several studies proposed that summary statistics of the deviations of an animal's observed performance from its target performance trajectory (i.e., performance in the absence of challenge) may constitute suitable quantitative resilience indicators. However, these statistical indicators require further validation. The aim of this study was to obtain a better understanding of these resilience indicators in their ability to discriminate between different response types and their dependence on different response characteristics of animals, and data recording features. To this purpose, milk-yield trajectories of individual dairy cattle differing in resilience, without and when exposed to a short-term challenge, were simulated. Individuals were categorised into three broad response types (with individual variation within each type): Fully Resilient animals, which experience no systematic perturbation in milk yield after challenge, Non-Resilient animals whose milk yield permanently deviates from the target trajectory after challenge and Partially Resilient animals that experience temporary perturbations but recover. The following statistical resilience indicators previously suggested in the literature were validated with respect to their ability to discriminate between response types and their sensitivity to various response features and data characteristics: logarithm of mean of squares (LMS), logarithm of variance (LV), skewness (S), lag-1 autocorrelation (AC1), and area under the curve (AUC) of deviations. Furthermore, different methods for estimating unknown target trajectories were evaluated. All of the considered resilience indicators could distinguish between the Fully Resilient response type and either of the other two types when target trajectories were known or estimated using a parametric method. When the comparison was between Partially Resilient and Non-Resilient, only LMS, LV, and AUC could correctly rank the response types, provided that the observation period was at least twice as long as the perturbation period. Skewness was in general the least reliable indicator, although all indicators showed correct dependency on the amplitude and duration of the perturbations. In addition, all resilience indicators except for AC1 were robust to lower frequency of measurements. In general, parametric methods (quantile or repeated regression) combined with three resilience indicators (LMS, LV and AUC) were found the most reliable techniques for ranking animals in terms of their resilience.

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 (AFS), and leg accelerometer) to track behavioral changes in response during the time when the stressors were present. The AFS 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 calf to novelty and fear stimuli. Factors extracted from a principal component analysis on the behaviors from the personality test were utilized to represent personality traits: Factor 1 ('Fearful'), Factor 2 ('Active') and Factor 3 ('Explorative'). The magnitude change 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 utilized to determine if 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 to 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, than 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.

Dare to be resilient: the key to future pesticide-free orchards?

Considering the urgent need for more sustainable fruit tree production, it is high time to find durable alternatives to the systematic use of phytosanitary products in orchards. To this end, resilience can deliver a number of benefits. Relying on a combination of tolerance, resistance, and recovery traits, disease resilience appears as a cornerstone to cope with the multiple pest and disease challenges over an orchard's lifetime. Here, we describe resilience as the capacity of a tree to be minimally affected by external disturbances or to rapidly bounce back to normal functioning after being exposed to these disturbances. Based on a literature survey largely inspired from research on livestock, we highlight different approaches for dissecting phenotypic and genotypic components of resilience. In particular, multisite experimental designs and longitudinal measures of so-called 'resilience biomarkers' are required. We identified a list of promising biomarkers relying on ecophysiological and digital measurements. Recent advances in high-throughput phenotyping and genomics tools will likely facilitate fine scale temporal monitoring of tree health, allowing identification of resilient genotypes with the calculation of specific resilience indicators. Although resilience could be considered as a 'black box' trait, we demonstrate how it could become a realistic breeding goal.

Quantification of the effect of in utero events on lifetime resilience in dairy cows

Currently, the dairy industry is facing many challenges that could affect its sustainability, including climate change and public perception of the industry. As a result, interest is increasing in the concept of identifying resilient animals, those with a long productive lifespan, as well as good reproductive performance and milk yield. There is much evidence that events in utero, that is, the developmental origins of health and disease hypothesis, alter the life-course health of offspring and we hypothesized that these could alter resilience in calves, where resilience is identified using lifetime data. The aim of this study was to quantify lifetime resilience scores (LRS) using an existing scoring system, based on longevity with secondary corrections for age at first calving and calving interval, and to quantify the effects of in utero events on the LRS using 2 datasets. The first was a large dataset of cattle on 83 farms in Great Britain born from 2006 to 2015 and the second was a smaller, more granular dataset of cattle born between 2003 and 2015 in the Langhill research herd at Scotland's Rural College. Events during dam's pregnancy included health events (lameness, mastitis, use of an antibiotic or anti-inflammatory medication), the effect of heat stress as measured by temperature-humidity index, and perturbations in milk yield and quality (somatic cell count, percentage fat, percentage protein and fat:protein ratio). Daughters born to dams that experienced higher temperature-humidity indexes while they were in utero during the first and third trimesters of pregnancy had lower LRS. Daughter LRS were also lower where milk yields or median fat percentages in the first trimester were low, and when milk yields were high in the third trimester. Dam LRS was positively associated with LRS of their offspring; however, as parity of the dam increased, LRS of their calves decreased. Similarly, in the Langhill herd, dams of a higher parity produced calves with lower LRS. Additionally, dams that recorded a high maximum locomotion score in the third trimester of pregnancy were negatively associated with lower calf LRS in the Langhill herd. Our results suggest that events that occur during pregnancy have lifelong consequences for the calf's lifetime performance. However, experience of higher temperature-humidity indexes, higher dam LRS, and mothers in higher parities explained a relatively small proportion of variation in offspring LRS, which suggests that other factors play a substantial role in determining calf LRS. Although "big data" can contain a considerable amount of noise, similar findings between the 2 datasets indicate it is likely these findings are real.

Breeding for resilience in finishing pigs can decrease tail biting, lameness and mortality

BACKGROUND

Previous research showed that deviations in longitudinal data are heritable and can be used as a proxy for pigs' general resilience. However, only a few studies investigated the relationship between these resilience traits and other traits related to resilience and welfare. Therefore, this study investigated the relationship between resilience traits derived from deviations in longitudinal data and traits related to animal resilience, health and welfare, such as tail and ear biting wounds, lameness and mortality.

RESULTS

In our experiment, 1919 finishing pigs with known pedigree (133 Piétrain sires and 266 crossbred dams) were weighed every 2 weeks and scored for physical abnormalities, such as lameness and ear and tail biting wounds (17,066 records). Resilience was assessed via deviations in body weight, deviations in weighing order and deviations in observed activity during weighing. The association between these resilience traits and physical abnormality traits was investigated and genetic parameters were estimated. Deviations in body weight had moderate heritability estimates (h2 = 25.2 to 36.3%), whereas deviations in weighing order (h2 = 4.2%) and deviations in activity during weighing (h2 = 12.0%) had low heritability estimates. Moreover, deviations in body weight were positively associated and genetically correlated with tail biting wounds (rg = 0.22 to 0.30), lameness (rg = 0.15 to 0.31) and mortality (rg = 0.19 to 0.33). These results indicate that events of tail biting, lameness and mortality are associated with deviations in pigs' body weight evolution. This relationship was not found for deviations in weighing order and activity during weighing. Furthermore, individual body weight deviations were positively correlated with uniformity at the pen level, providing evidence that breeding for these resilience traits might increase both pigs' resilience and within-family uniformity.

CONCLUSIONS

In summary, our findings show that breeding for resilience traits based on deviations in longitudinal weight data can decrease pigs' tail biting wounds, lameness and mortality while improving uniformity at the pen level. These findings are valuable for pig breeders, as they offer evidence that these resilience traits are an indication of animals' general health, welfare and resilience. Moreover, these results will stimulate the quantification of resilience via longitudinal body weights in other species.

Genetic parameters for novel climatic resilience indicators derived from automatically-recorded vaginal temperature in lactating sows under heat stress conditions

BACKGROUND

Longitudinal records of automatically-recorded vaginal temperature (TV) could be a key source of data for deriving novel indicators of climatic resilience (CR) for breeding more resilient pigs, especially during lactation when sows are at an increased risk of suffering from heat stress (HS). Therefore, we derived 15 CR indicators based on the variability in TV in lactating sows and estimated their genetic parameters. We also investigated their genetic relationship with sows' key reproductive traits.

RESULTS

The heritability estimates of the CR traits ranged from 0.000 ± 0.000 for slope for decreased rate of TV (SlopeDe) to 0.291 ± 0.047 for sum of TV values below the HS threshold (HSUB). Moderate to high genetic correlations (from 0.508 ± 0.056 to 0.998 ± 0.137) and Spearman rank correlations (from 0.431 to 1.000) between genomic estimated breeding values (GEBV) were observed for five CR indicators, i.e. HS duration (HSD), the normalized median multiplied by normalized variance (Nor_medvar), the highest TV value of each measurement day for each individual (MaxTv), and the sum of the TV values above (HSUA) and below (HSUB) the HS threshold. These five CR indicators were lowly to moderately genetically correlated with shoulder skin surface temperature (from 0.139 ± 0.008 to 0.478 ± 0.048) and respiration rate (from 0.079 ± 0.011 to 0.502 ± 0.098). The genetic correlations between these five selected CR indicators and sow reproductive performance traits ranged from - 0.733 to - 0.175 for total number of piglets born alive, from - 0.733 to - 0.175 for total number of piglets born, and from - 0.434 to - 0.169 for number of pigs weaned. The individuals with the highest GEBV (most climate-sensitive) had higher mean skin surface temperature, respiration rate (RR), panting score (PS), and hair density, but had lower mean body condition scores compared to those with the lowest GEBV (most climate-resilient).

CONCLUSIONS

Most of the CR indicators evaluated are heritable with substantial additive genetic variance. Five of them, i.e. HSD, MaxTv, HSUA, HSUB, and Nor_medvar share similar underlying genetic mechanisms. In addition, individuals with higher CR indicators are more likely to exhibit better HS-related physiological responses, higher body condition scores, and improved reproductive performance under hot conditions. These findings highlight the potential benefits of genetically selecting more heat-tolerant individuals based on CR indicators.

Key determinants of adaptive strategies of goats to a 2-day nutritional challenge during early lactation

Little is known about the key determinants of the physiological adaptations to environmental challenges and how these determinants interact. We evaluated how the response/recovery profiles to a short-term nutritional challenge during early lactation are affected by early-life nutritional strategies in dairy goats divergently selected for functional longevity. We used 72 females, split into two cohorts, daughters of Alpine bucks divergently selected for functional longevity. The females from the two lines were fed with two divergent diets, normal vs low-energy, from weaning until the middle of first gestation, and then fed with the same standard diet. Individual BW, body condition score, morphology, and plasma samples were collected from birth to first kidding. The adaptative physiological strategy to a nutritional challenge was assessed via a 2-day feed restriction challenge, during early lactation, which consisted of a five-day control period on a standard lactation diet followed by a 2-day challenge with straw-only feeding and then a 10-day recovery period on a standard lactation diet. During the challenge, DM intake, BW, milk yield (MY), and plasma and milk metabolite composition were recorded daily. Linear mixed-effects models were used to analyze all traits, considering the individual nested in the cohort as a random effect and the 2 × 2 treatments (i.e., line and rearing diet) and litter size as fixed effects. Linear mixed-effects models using a piecewise arrangement were used to analyze the response/recovery profiles to nutritional challenge. Random parameters estimated for each individual, using the mixed-effects models without the fixed effects of rearing diet and genetic line, were used in a stepwise model selection based on R2 to identify key determinants of an individual's physiological adaptations to environmental challenges. Differences in stature and body reserves created by the two rearing diets diminished during late gestation and the 5-day control period. Genetic line did not affect body reserves during the rearing phase. Rearing diet and genetic line slightly affected the recovery profiles of evaluated traits and had no effects on prechallenge and response to challenge profiles. The prekidding energy status measures and MY before challenge were selected as strong predictors of variability in response-recovery profiles of milk metabolites that have strong links with body energy dynamics (i.e., isoCitrate, ß-hydroxybutyrate, choline, cholesterol, and triacylglycerols; R2 = 35%). Our results suggested that prekidding energy status and MY are key determinants of adult resilience and that rearing diet and genetic line may affect adult resilience insofar as they affect the animals' energy status.

Resilience indicator traits in 3 dairy cattle breeds in Baden-Württemberg

In recent years, research in animal breeding has increasingly focused on the topic of resilience, which is expected to continue in the future due to the need for high-yielding, healthy, and robust animals. In this context, an established approach is the calculation of resilience indicator traits with time series analyses. Examples are the variance and autocorrelation of daily milk yield in dairy cows. We applied this methodology to the German dairy cow population. Data from the 3 breeds (German Holstein, German Fleckvieh, and German Brown Swiss) were obtained, which included 13,949 lactations from 36 farms from the state Baden-Württemberg in Germany working with automatic milking systems. Using the milk yield data, the daily absolute milk yields, deviations between observed and expected daily milk yields, and relative proportions of daily milk yields in relation to lactation performance were calculated. We used the variance and autocorrelation of these data as phenotypes in our statistical analyses. We estimated a heritability of 0.047 for autocorrelation and heritabilities between 0.026 and 0.183 for variance-based indicator traits. Furthermore, significant breed differences could be observed, with a tendency of better resilience in Brown Swiss. The breed differences can be due to both genetic and environmental factors. A high value of a variance-based indicator trait indicates a low resilience. Performance traits were positively correlated with variance-based indicator traits calculated from absolute daily milk yields, but they were negatively correlated with variance-based indicators calculated from relative daily milk yields. Thus, they can be considered as different traits. Although variance-based indicators based on absolute daily milk yields were affected by the performance level, variance-based indicators based on relative daily milk yields were corrected for the performance level and also showed higher heritabilities. Thus, they seem to be more suitable for practical use. Further studies need to be conducted to calculate the correlations between resilience indicators, functional traits, and health traits.

Review: Environmental enrichment builds functional capacity and improves resilience as an aspect of positive welfare in production animals

The success of the animal in coping with challenges, and in harnessing opportunities to thrive, is central to its welfare. Functional capacity describes the capacity of molecules, cells, organs, body systems, the whole animal, and its community to buffer against the impacts of environmental perturbations. This buffering capacity determines the ability of the animal to maintain or regain functions in the face of environmental perturbations, which is recognised as resilience. The accuracy of physiological regulation and the maintenance of homeostatic balance underwrite the dynamic stability of outcomes such as biorhythms, feed intake, growth, milk yield, and egg production justifying their assessment as indicators of resilience. This narrative review examines the influence of environmental enrichments, especially during developmental stages in young animals, in building functional capacity and in its subsequent expression as resilience. Experience of enriched environments can build skills and competencies across multiple functional domains including but not limited to behaviour, immunity, and metabolism thereby increasing functional capacity and facilitating resilience within the context of challenges such as husbandry practices, social change, and infection. A quantitative method for measuring the distributed property of functional capacity may improve its assessment. Methods for analysing embedded energy (emergy) in ecosystems may have utility for this goal. We suggest functional capacity provides the common thread that links environmental enrichments with an ability to express resilience and may provide a novel and useful framework for measuring and reporting resilience. We conclude that the development of functional capacity and its subsequent expression as resilience is an aspect of positive animal welfare. The emergence of resilience from system dynamics highlights a need to shift from the study of physical and mental states to the study of physical and mental dynamics to describe the positive dimension of animal welfare.

Genetic analysis of phenotypic indicators for heat tolerance in crossbred dairy cattle

Climate change-induced rise in global temperatures has intensified heat stress on dairy cattle and is contributing to the generally observed low milk productivity. Selective breeding aimed at enhancing animals' ability to withstand rising temperatures while maintaining optimal performance is crucial for ensuring future access to dairy products. However, phenotypic indicators of heat tolerance are yet to be effectively factored into the objectives of most selective breeding programs. This study investigated the response of milk production to changing heat load as an indication of heat tolerance and the influence of calving season on this response in multibreed dairy cattle performing in three agroecological zones Kenya. First-parity 7-day average milk yield (65 261 milk records) of 1 739 cows were analyzed. Based on routinely recorded weather data that were accessible online, the Temperature-Humidity Index (THI) was calculated and used as a measure of heat load. THI measurements used represented averages of the same 7-day periods corresponding to each 7-day average milk record. Random regression models, including reaction norm functions, were fitted to derive two resilience indicators: slope of the reaction norm (Slope) and its absolute value (Absolute), reflecting changes in milk yield in response to the varying heat loads (THI 50 and THI 80). The genetic parameters of these indicators were estimated, and their associations with average test-day milk yield were examined. There were no substantial differences in the pattern of milk yield response to heat load between cows calving in dry and wet seasons. Animals with ≤50% Bos taurus genes were the most thermotolerant at extremely high heat load levels. Animals performing in semi-arid environments exhibited the highest heat tolerance capacity. Heritability estimates for these indicators ranged from 0.06 to 0.33 and were mostly significantly different from zero (P < 0.05). Slope at THI 80 had high (0.64-0.71) negative correlations with average daily milk yield, revealing that high-producing cows are more vulnerable to heat stress and vice versa. A high (0.63-0.74) positive correlation was observed between Absolute and average milk yield at THI 80. This implied that low milk-producing cows have a more stable milk production under heat-stress conditions and vice versa. The study demonstrated that the slope of the reaction norms and its absolute value can effectively measure the resilience of crossbred dairy cattle to varying heat load conditions. The implications of these findings are valuable in improving the heat tolerance of livestock species through genetic selection.

Clustering of countries based on dairy productivity characteristics of Holstein cattle for breeding material selection

Background and Aim

The aim of any breeding process is to create a herd based on certain parameters that reflect an ideal animal vision. Targeted herding involves selecting the source of breeding material to be imported from another country. Therefore, there is a problem in selecting a breeding material importer to rapidly form a uterine canopy with the required properties. The purpose of this study was to evaluate a set of predictive milk productivity traits in Holstein cattle across countries.

Materials and Methods

This research was based on records of 819,358 recorded animals from 28 countries born after January 1, 2018, from open databases. We used the Euclidean metric to construct dendrograms characterizing the similarity of countries according to the complex milk productivity traits of the daughters of bulls. The Ward method was used to minimize intracluster variance when forming clusters and constructing the corresponding diagrams. Principal component analysis was used to reduce dimensionality and eliminate the effect of multicollinearity. The principal components were selected using the Kaiser-Harris criteria.

Results

A ranking of multidimensional complex milk productivity traits in different countries over the past 5 years was performed. A group of leading countries led by the USA was established according to the studied indicators, and the possible reasons for such a division into groups were described.

Conclusion

The pressure of purposeful artificial selection prevails in comparison with the pressure of natural selection concerning milk productivity traits in a certain group of countries, which allows specialists to choose suppliers when buying breeding animals and materials. The findings are based solely on data from recorded animals, which may not represent the entire breed population within each country, especially in regions where record-keeping may be inconsistent. It is expected that further studies will include regional data from large enterprises not part of Interbull, with mandatory verification and validation. An important element of such work is seen as the ability to compare the milk productivity of populations from different countries using a different scale, as well as studying the differentiation of countries by other selection traits of dairy.

A comprehensive meta-analysis of genetic parameters for resilience and productivity indicator traits in Holstein cattle

Selection for resilience indicator (RIND) traits in Holstein cattle is becoming an important breeding objective as the worldwide population is expected to be exposed to increased environmental stressors due to both climate change and changing industry standards. However, genetic correlations between RIND and productivity indicator (PIND) traits, which are already being selected for and have the most economic value, are often unfavorable. As a result, it is necessary to fully understand these genetic relationships when incorporating novel traits into selection indices, so that informed decisions can be made to fully optimize selection for both groups of traits. In the past 2 decades, there have been many estimates of RIND traits published in the literature, albeit in small populations. To provide valuable pooled summary estimates, a random-effects meta-analysis was conducted for heritability and genetic correlation estimates for PIND and RIND traits in worldwide Holstein cattle. In total, 926 heritability estimates for 9 PIND and 27 RIND traits, along with 362 estimates of genetic correlation (PIND × RIND traits) were collected. Resilience indicator traits were grouped into the following subgroups: Metabolic Diseases, Hoof Health, Udder Health, Fertility, Heat Tolerance, Longevity, and Other. Pooled estimates of heritability for PIND traits ranged from 0.201 ± 0.05 (energy-corrected milk) to 0.377 ± 0.06 (protein content), while pooled estimates of heritability for RIND traits ranged from 0.032 ± 0.02 (incidence of lameness, incidence of milk fever) to 0.497 ± 0.05 (measures of body weight). Pooled estimates of genetic correlations ranged from -0.360 ± 0.25 (protein content vs. milk acetone concentration) to 0.535 ± 0.72 (measures of fat-to-protein ratio vs. milk acetone concentration). Additionally, out of 243 potential genetic correlations between PIND and RIND traits that could have been reported, only 40 had enough published estimates to implement the meta-analysis model. Our results confirmed that the interactions between PIND and RIND traits are complex, and all relationships should be evaluated when incorporating novel traits into selection indices. This study provides a valuable reference for breeders looking to incorporate RIND traits for Holstein cattle into selection indices.

Responses of selected plasma metabolites to a two-day nutritional challenge of goats divergently selected for functional longevity

Understanding the extent to which genetics × environment plays a role in shaping individual strategies to environmental challenges is of considerable interest for future selection of more resilient animals. Accordingly, the objective of this study was to evaluate the metabolic responses to a nutritional challenge of goats divergently selected for functional longevity based on plasma metabolites and the repeatability of these responses across 2 experimental farms and years. We carried out 6 different experimental trials from years 2018 to 2022 (4 trials on site Bourges (2018-21) and 2 trials (2021-22) on site Grignon) in which 267 first kidding goats, daughters of Alpine bucks divergently selected for functional longevity, longevity plus (n = 137), and longevity minus (n = 130), were exposed to a 2-d nutritional challenge in early lactation. The experiments consisted of a 5 or 7-d control period (pre-challenge) on a standard lactation diet followed by a 2-d nutritional challenge with straw-only feeding and then a 7 or 10-d recovery period on a standard lactation diet, for site Bourges and Grignon, respectively. During the challenge plasma metabolite composition was recorded daily. Linear mixed-effects models were used to analyze all traits, considering the individual as a random effect and the 2x2 treatments (i.e., genetic line and year nested in site) and litter size as fixed effects. The linear mixed-effects model using a piecewise arrangement was used to analyze the response/recovery profiles to the nutritional challenge. Random parameters estimated for each individual, using the mixed-effects models without the fixed effects of genetic line, were used in a Sparse Partial Least Square Discriminant Analysis (sPLS-DA) to compare the goat metabolism response to the challenge on a multivariate scale. The plasma metabolites, glucose, β-hydroxybutyrate (BHB), and nonesterified fatty acids (NEFA), and urea concentrations responded to the 2-d nutritional challenge. Selection for functional longevity did not affect plasma glucose, NEFA, BHB, and urea response/recoveries to a 2-d nutritional challenge. However, site, trial, and litter size affected these responses. Moreover, the plasma metabolites seem not to fully recover to prechallenge levels after the recovery phase. The sPLS-DA analysis did not discriminate between the 2 longevity lines. We observed meaningful between-individuals' variability in plasma BHB, especially on the prechallenge and rate of response and rate of recovery from the 2-d nutritional challenge (CV = 26.2%, 36.1%, and 41.2%, repeatability = 0.749, 0.322, and 0.741, respectively). Plasma NEFA recovery from challenge also demonstrated high between-individuals' variability (CV = 16.4%, repeatability = 0.323). Selection for functional longevity did not affect plasma metabolites responses to a 2-d nutritional challenge in dairy goats. Plasma NEFA and BHB response/recovery presented high between-individuals' variability, indicating individual adaptative characteristics to nutritional challenges not related to the environmental conditions but to inherent individual characteristics.

Microbiome composition as a potential predictor of longevity in rabbits

BACKGROUND

Longevity and resilience are two fundamental traits for more sustainable livestock production. These traits are closely related, as resilient animals tend to have longer lifespans. An interesting criterion for increasing longevity in rabbit could be based on the information provided by its gut microbiome. The gut microbiome is essential for regulating health and plays crucial roles in the development of the immune system. The aim of this research was to investigate if animals with different longevities have different microbial profiles. We sequenced the 16S rRNA gene from soft faeces from 95 does. First, we compared two maternal rabbit lines with different longevities; a standard longevity maternal line (A) and a maternal line (LP) that was founded based on longevity criteria: females with a minimum of 25 parities with an average prolificacy per parity of 9 or more. Second, we compared the gut microbiota of two groups of animals from line LP with different longevities: females that died/were culled with two parities or less (LLP) and females with more than 15 parities (HLP).

RESULTS

Differences in alpha and beta diversity were observed between lines A and LP, and a partial least square discriminant analysis (PLS-DA) showed a high prediction accuracy (> 91%) of classification of animals to line A versus LP (146 amplicon sequence variants (ASV)). The PLS-DA also showed a high prediction accuracy (> 94%) to classify animals to the LLP and HLP groups (53 ASV). Interestingly, some of the most important taxa identified in the PLS-DA were common to both comparisons (Akkermansia, Christensenellaceae R-7, Uncultured Eubacteriaceae, among others) and have been reported to be related to resilience and longevity.

CONCLUSIONS

Our results indicate that the first parity gut microbiome profile differs between the two rabbit maternal lines (A and LP) and, to a lesser extent, between animals of line LP with different longevities (LLP and HLP). Several genera were able to discriminate animals from the two lines and animals with different longevities, which shows that the gut microbiome could be used as a predictive factor for longevity, or as a selection criterion for these traits.

Genetic analysis of lactation consistency in US Holsteins using temporal variation in daily milk weights

The ability of a dairy cow to perform reliably over time is an interesting trait to include in dairy cattle breeding programs aimed at improving dairy cow resilience. Consistency, defined as the quality of performing as expected each day of the lactation, could be highly associated with resilience, defined as animal's ability to maintain health and performance in the presence of environmental challenges, including pathogens, heat waves, and nutritional changes. A total of 51,415,022 daily milk weights collected from 2018 to 2023 were provided for 255,191 multiparous Holstein cows milked 3 times daily in conventional parlor systems on farms in 32 states. The temporal variance (TEMPVAR) of milk yield from 5 to 305 d postpartum was computed as the log-transformed variance of daily deviations between observed and expected individual milk weights. Lower values of TEMPVAR imply smaller day-to-day deviations from expectations, indicating consistent performance, whereas larger values indicate inconsistent performance. Expected daily milk weights were computed using 3 nonparametric and parametric regression models: (1) loceally estimated scatterplot smoothing regression with a 0.75 span; (2) polynomial quantile regression using the median (0.5 quantile), and (3) polynomial quantile regression using a 0.7 quantile. The univariate statistical model included age at first calving and herd-year-season as fixed effects and cow as a random effect. Heritability estimates (standard errors) of TEMPVAR phenotypes calculated over the entire lactation ranged between 0.227 (0.011) and 0.237 (0.011), demonstrating that cows are genetically predisposed to display consistent or inconsistent performance. Estimated genetic correlations calculated using a multiple-trait model between TEMPVAR traits and between lactations were high (>0.95), indicating TEMPVAR is repeatable across lactations and robust to the model used to compute expected daily milk yield. Higher TEMPVAR phenotypes reflect more variation in performance, hence greater inconsistency, which is undesirable. Therefore, correlations between predicted transmitting abilities (PTA) for TEMPVAR and milk yield of 0.57 indicate that high-producing cows exhibit more day-to-day variation in performance. Correlations with productive life and livability were -0.38 and -0.48, respectively. Correlations between PTA for TEMPVAR and those of postpartum health traits were also negative, ranging from -0.41 to -0.08. Given that health traits are derived from disease resistance measurements, and higher health trait PTA are preferred, our results indicate that more consistent cows tend to have fewer health problems and greater longevity. Overall, our findings suggest that temporal variation in daily milk weights can be used to identify consistent animals that maintain expected performance throughout the lactation, which will enable selection for greater resilience to management and environmental perturbations.

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.

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.

Investigating the relationship between fluctuations in daily milk yield as resilience indicators and health traits in Holstein cattle

Disease-related milk losses directly affect dairy herds' profitability and the production efficiency of the dairy industry. Therefore, this study aimed to quantify phenotypic variability in milk fluctuation periods related to diseases and to explore milk fluctuation traits as indicators of disease resilience. By combining high-frequency daily milk yield data with disease records of cows that were treated and recovered from the disease, we estimated milk variability trends within a fixed period around the treatment day of each record for 5 diseases: udder health, reproductive disorders, metabolic disorders, digestive disorders, and hoof health. The average milk yield decreased rapidly from 6 to 8 d before the treatment day for all diseases, with the largest milk reduction observed on the treatment day. Additionally, we assessed the significance of milk fluctuation periods highly related to diseases by defining milk fluctuations as a period of at least 10 consecutive days in which milk yield fell below 90% of the expected milk production values at least once. We defined the development and recovery phases of milk fluctuations using 3,847 milk fluctuation periods related to disease incidences, and estimated genetic parameters of milk fluctuation traits, including milk losses, duration of the fluctuation, variation rate in daily milk yield, and standard deviation of milk deviations for each phase and their genetic correlation with several important traits. In general, the disease-related milk fluctuation periods lasted 21.19 ± 10.36 d with a milk loss of 115.54 ± 92.49 kg per lactation. Compared with the development phase, the recovery phase lasted an average of 3.3 d longer, in which cows produced 11.04 kg less milk and exhibited a slower variation rate in daily milk yield of 0.35 kg/d. There were notable differences in milk fluctuation traits depending on the disease, and greater milk losses were observed when multiple diseases occurred simultaneously. All milk fluctuation traits evaluated were heritable with heritability estimates ranging from 0.01 to 0.10, and moderate to high genetic correlations with milk yield (0.34 to 0.64), milk loss throughout the lactation (0.22 to 0.97), and resilience indicator (0.39 to 0.95). These results indicate that cows with lower milk losses and higher resilience tend to have more stable milk fluctuations, which supports the potential for breeding for more disease-resilient cows based on milk fluctuation traits. Overall, this study confirms the high effect of diseases on milk yield variability and provides insightful information about their relationship with relevant traits in Holstein cattle. Furthermore, this study shows the potential of using high-frequency automatic monitoring of milk yield to assist on breeding practices and health management in dairy cows.

Exploring genetic variation in potential indicators of resilience in sheep using fibre diameter measured along the wool staple

Production animals are increasingly exposed to a wide variety of disturbances that can compromise their productivity, health and well-being. As a result, there is a growing need to be able to select animals that are more resilient to environmental disturbances. Fibre diameter variation measured along a wool staple is expected to contain information about how resilient sheep are to the disturbances of their internal and external environment. This study aimed to develop potential resilience indicators from fibre diameter variation, estimate their genetic parameters and assess whether these traits are genetically correlated across three age stages. The study used 6 140 Merino sheep from the Sheep Cooperative Research Centre Information Nucleus Flocks recorded at yearling, 2 years old, and adult ages. Eight potential traits were defined based on theory, literature and exploratory analysis, which were suggested to capture the animal's ability to resist, respond and recover from potential disturbances. Genetic evaluation of the traits was conducted using pedigree-based animal models. The traits were shown to be low to moderately heritable (0.01-0.33) when examined at each of the three age stages. The potential indicators were generally well correlated with one another within age stages. Further, the genetic correlation between the same trait measured at different age stages was moderate to high between yearling and 2 years old (0.35-0.94) and between 2 years old and adults (0.18-0.70), while slightly lower between yearling and adult estimates (0.09-0.62). These results suggest that selection for resilience indicators from fibre diameter is possible; however, further studies are warranted to refine the trait definitions and validate these indicators against other measures of health, fitness and productive performance.

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.

Consistency of dry matter intake in Holstein cows: Heritability estimates and associations with feed efficiency

Resilience can be defined as the capacity to maintain performance or bounce back to normal functioning after a perturbation, and studying fluctuations in daily feed intake may be an effective way to identify resilient dairy cows. Our goal was to develop new phenotypes based on daily dry matter intake (DMI) consistency in Holstein cows, estimate genetic parameters and genetic correlations with feed efficiency and milk yield consistency, and evaluate their relationships with production, longevity, health, and reproduction traits. Data consisted of 397,334 daily DMI records of 6,238 lactating Holstein cows collected from 2007 to 2022 at 6 research stations across the United States. Consistency phenotypes were calculated based on the deviations from expected daily DMI for individual cows during their respective feeding trials, which ranged from 27 to 151 d in duration. Expected values were derived from different models, including simple average, quadratic and cubic quantile regression with a 0.5 quantile, and locally estimated scatterplot smoothing (LOESS) regression with span parameters 0.5 and 0.7. We then calculated the log of variance (log-Var-DMI) of daily deviations for each model as the consistency phenotype. Consistency of milk yield was also calculated, as a reference, using the same methods (log-Var-Milk). Genetic parameters were estimated using an animal model, including lactation, days in milk and cohort as fixed effects, and animal as random effect. Relationships between log-Var-DMI and traits currently considered in the US national genetic evaluation were evaluated using Spearman's rank correlations between sires' breeding values. Heritability estimates for log-Var-DMI ranged from 0.11 ± 0.02 to 0.14 ± 0.02 across models. Different methods (simple average, quantile regressions, and LOESS regressions) used to calculate log-Var-DMI yielded very similar results, with genetic correlations ranging from 0.94 to 0.99. Estimated genetic correlations between log-Var-DMI and log-Var-Milk ranged from 0.51 to 0.62. Estimated genetic correlations between log-Var-DMI and feed efficiency ranged from 0.55 to 0.60 with secreted milk energy, from 0.59 to 0.63 with metabolic body weight, and from 0.26 to 0.31 with residual feed intake (RFI). Relationships between log-Var-DMI and the traits in the national genetic evaluation were moderate and positive correlations with milk yield (0.20 to 0.21), moderate and negative correlations with female fertility (-0.07 to -0.20), no significant correlations with health and longevity, and favorable correlations with feed efficiency (-0.23 to -0.25 with feed saved and 0.21 to 0.26 with RFI). We concluded that DMI consistency is heritable and may be an indicator of resilience. Cows with lower variation in the difference between actual and expected daily DMI (more consistency) may be more effective in maintaining performance in the face of challenges or perturbations, whereas cows with greater variation in observed versus expected daily DMI (less consistency) are less feed efficient and may be less resilient.

Stress-induced loss of social resilience in honeybee colonies and its implications on fitness

Stressors may lead to a shift in the timing of life-history events of species, causing a mismatch with optimal environmental conditions, potentially reducing fitness. In honeybees, the timing of brood rearing and nest emergence in late winter/early spring is critical as colonies need to grow fast after winter to prepare for reproduction. However, the effects of stress on these life-history events in late winter/early spring and the possible consequences are not well understood. Therefore, we tested whether (i) honeybee colonies shift timing of brood rearing and nest emergence as response to stressors, and (ii) if there is a consequent loss of social resilience, reflected in colony fitness (survival, growth and reproduction). We monitored stressed (high load of the parasitic mite Varroa destructor or nutrition restricted) colonies and presumably non-stressed colonies from the beginning of 2020 till spring of 2021. We found that honeybee colonies do not shift the timing of brood rearing and nest emergence in spring as a coping mechanism to stressors. However, we show that there is loss of social resilience in stressed colonies, leading to reduced growth and reproduction. Our study contributes to better understanding the effects of stressors on social resilience in eusocial organisms.

Identification and characterization of unknown disturbances in a structured population using high-throughput phenotyping data and measurement of robustness: application to growing pigs

Improving the robustness of animals has become a priority in breeding due to climate change, new societal demands, and the agroecological transition. Components of animal robustness can be extracted from the analysis of the adaptive response of an animal to disturbance using longitudinal data. Nonetheless, this response is a function of animal robustness as well as of disturbance characteristics (intensity and duration). To correctly assess an animal's robustness potential, it is therefore useful to know the characteristics of the disturbances it faces. The UpDown method, which detects and characterizes unknown disturbances at different levels of organization of the population (e.g., individual, pen, and batch disturbances), has been proposed for this purpose. Furthermore, using the outputs of the method, it is possible to extract proxies of the robustness of animals. In this context, the objective of the study was to evaluate the performances of the UpDown method to detect and characterize disturbances and quantify the robustness of animals in a genetic framework using different sets of simulations, and to apply this method to real pig longitudinal data recorded during the fattening period (body weight, cumulative feed intake, and feeding rate). Based on the simulations, the specificity of the UpDown method was high (>0.95). Its sensitivity increased with the level of organization exposed (from 0.23 to 0.32 for individual disturbances, from 0.45 to 0.59 for pen disturbances, and from 0.77 to 0.88 for batch disturbances). The UpDown method also showed a good ability to characterize detected disturbances. The average time interval between the estimated and true start date or duration of the disturbance was lower than 3 d. The correlation between the true and estimated intensity of the disturbance increased with the hierarchical level of organization (on average, 0.41, 0.78, and 0.83 for individual, pen, and batch disturbance, respectively). The accuracy of the estimated breeding values of the proxies for robustness extracted from the analysis of individual trajectories over time were moderate (lower than 0.33). Applied to real data, the UpDown method detected different disturbances depending on the phenotype analyzed. The heritability of the proxies of robustness were low to moderate (ranging from 0.11 to 0.20).

Milk yield residuals and their link with the metabolic status of dairy cows in the transition period

The transition period is one of the most challenging periods in the lactation cycle of high-yielding dairy cows. It is commonly known to be associated with diminished animal welfare and economic performance of dairy farms. The development of data-driven health monitoring tools based on on-farm available milk yield development has shown potential in identifying health-perturbing events. As proof of principle, we explored the association of these milk yield residuals with the metabolic status of cows during the transition period. Over 2 yr, 117 transition periods from 99 multiparous Holstein-Friesian cows were monitored intensively. Pre- and postpartum dry matter intake was measured and blood samples were taken at regular intervals to determine β-hydroxybutyrate, nonesterified fatty acids (NEFA), insulin, glucose, fructosamine, and IGF1 concentrations. The expected milk yield in the current transition period was predicted with 2 previously developed models (nextMILK and SLMYP) using low-frequency test-day (TD) data and high-frequency milk meter (MM) data from the animal's previous lactation, respectively. The expected milk yield was subtracted from the actual production to calculate the milk yield residuals in the transition period (MRT) for both TD and MM data, yielding MRTTD and MRTMM. When the MRT is negative, the realized milk yield is lower than the predicted milk yield, in contrast, when positive, the realized milk yield exceeded the predicted milk yield. First, blood plasma analytes, dry matter intake, and MRT were compared between clinically diseased and nonclinically diseased transitions. MRTTD and MRTMM, postpartum dry matter intake and IGF1 were significantly lower for clinically diseased versus nonclinically diseased transitions, whereas β-hydroxybutyrate and NEFA concentrations were significantly higher. Next, linear models were used to link the MRTTD and MRTMM of the nonclinically diseased cows with the dry matter intake measurements and blood plasma analytes. After variable selection, a final model was constructed for MRTTD and MRTMM, resulting in an adjusted R2 of 0.47 and 0.73, respectively. While both final models were not identical the retained variables were similar and yielded comparable importance and direction. In summary, the most informative variables in these linear models were the dry matter intake postpartum and the lactation number. Moreover, in both models, lower and thus also more negative MRT were linked with lower dry matter intake and increasing lactation number. In the case of an increasing dry matter intake, MRTTD was positively associated with NEFA concentrations. Furthermore, IGF1, glucose, and insulin explained a significant part of the MRT. Results of the present study suggest that milk yield residuals at the start of a new lactation are indicative of the health and metabolic status of transitioning dairy cows in support of the development of a health monitoring tool. Future field studies including a higher number of cows from multiple herds are needed to validate these findings.

Potential for quantifying general environmental resilience of dairy cattle in sub-Saharan Africa using deviations in milk yield

Introduction: Genetic improvement of general resilience of dairy cattle is deemed as a part of the solution to low dairy productivity and poor cattle adaptability in sub-Saharan Africa (SSA). While indicators of general resilience have been proposed and evaluated in other regions, their applicability in SSA remains unexplored. This study sought to test the viability of utilizing log-transformed variance (LnVar), autocorrelation (rauto), and skewness (Skew) of deviations in milk yield as indicators of general resilience of dairy cows performing in the tropical environment of Kenya. Methods: Test-day milk yield records of 2,670 first-parity cows performing in three distinct agroecological zones of Kenya were used. To predict expected milk yield, quantile regression was used to model lactation curve for each cow. Subsequently, resilience indicators were defined based on actual and standardized deviations of observed milk yield from the expected milk yield. The genetic parameters of these indicators were estimated, and their associations with longevity and average test-day milk yield were examined. Results: All indicators were heritable except skewness of actual and standardized deviation. The log-transformed variance of actual (LnVar1) and standardized (LnVar2) deviations had the highest heritabilities of 0.19 ± 0.04 and 0.17 ± 0.04, respectively. Auto-correlation of actual (rauto1) and standardized (rauto2) deviations had heritabilities of 0.05 ± 0.03 and 0.07 ± 0.03, respectively. Weak to moderate genetic correlations were observed among resilience indicators. Both rauto and Skew indicators had negligible genetic correlations with both longevity and average test-day milk yield. LnVar1 and LnVar2 were genetically associated with better longevity (rg = -0.47 ± 0.26 and -0.49 ± 0.26, respectively). Whereas LnVar1 suggested that resilient animals produce lower average test-day milk yield, LnVar2 revealed a genetic association between resilience and higher average test-day milk yield. Discussion: Log transformed variance of deviations in milk yield holds a significant potential as a robust resilience indicator for dairy animals performing in SSA. Moreover, standardized as opposed to actual deviations should be employed in defining resilience indicators because the resultant indicator does not inaccurately infer that low-producing animals are inherently resilient. This study offers an opportunity for enhancing the productivity of dairy cattle performing in SSA through selective breeding for resilience to environmental stressors.

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.

Association between management practices and estimated mastitis incidence and milk losses on robotic dairy farms

This study aims to describe the relation between farm-level management factors and estimated farm-level mastitis incidence and milk loss traits (MIMLT) at dairy farms with automated milking systems. In this observational study, 43 commercial dairy farms in Belgium and the Netherlands were included and 148 'management and udder health related variables' were obtained during a farm visit through a farm audit and survey. The MIMLT were estimated from milk yield data. Quarter-level milk yield perturbations that were caused by presumable mastitis cases (PMC) were selected based on quarter-level milk yield and electrical conductivity. On average, 57.6 ± 5.4% of the identified milk yield perturbations complied with our criteria. From these PMC, 3 farm-level MIMLT were calculated over a one-year period around the farm visit date: (1) the 'average number of PMC per cow per year', (2) the 'absolute milk loss per cow per day', calculated as the farm-level sum of all milk losses during PMC in one year, divided by the average number of lactating cows and the number of days, and (3) the 'relative milk loss', calculated as the farm-level sum of milk losses during PMC in one year, divided by the estimated total production in the absence of PMC. The 'average number of PMC per cow per year' was on average 1.81 ± 0.47. The PMC caused an average milk loss of 0.77 ± 0.26 kg per lactating cow per day, which corresponded to an average production loss of 2.38 ± 0.82% of the expected production in the absence of PMC. We performed a principal component regression (PCR) analysis to link the 3 MIMLT to the 'management and udder health related variables', whilst reducing the multicollinearity and the number of dimensions. The first principal component was mainly related to 'milking system brand, maintenance and settings'. The second component mainly linked to average productivity and somatic cell counts, whereas the third component mainly contained variables linked with mastitis management, treatment, and biosecurity. The 3 PCR models had R² ranging from 0.46 (for absolute milk loss per cow per day) to 0.57 (for relative milk loss). For all models, the second PC had the largest effect size. This analysis raises awareness of the impact of management factors on a factual basis and provides handles to take management actions to improve udder health.

Multivariate analysis of milk metabolite measures shows potential for deriving new resilience phenotypes

In a context of growing interest in breeding more resilient animals, a noninvasive indicator of resilience would be very valuable. We hypothesized that the time-course of concentrations of several milk metabolites through a short-term underfeeding challenge could reflect the variation of resilience mechanisms to such a challenge. We submitted 138 one-year-old primiparous goats, selected for extreme functional longevity (i.e., productive longevity corrected for milk yield [60 low longevity line goats and 78 high longevity line goats]), to a 2-d underfeeding challenge during early lactation. We measured the concentration of 13 milk metabolites and the activity of 1 enzyme during prechallenge, challenge, and recovery periods. Functional principal component analysis summarized the trends of milk metabolite concentration over time efficiently without preliminary assumptions concerning the shapes of the curves. We first ran a supervised prediction of the longevity line of the goats based on the milk metabolite curves. The partial least square analysis could not predict the longevity line accurately. We thus decided to explore the large overall variability of milk metabolite curves with an unsupervised clustering. The large year × facility effect on the metabolite concentrations was precorrected for. This resulted in 3 clusters of goats defined by different metabolic responses to underfeeding. The cluster that showed higher β-hydroxybutyrate, cholesterol, and triacylglycerols increase during the underfeeding challenge was associated with poorer survival compared with the other 2 clusters. These results suggest that multivariate analysis of noninvasive milk measures show potential for deriving new resilience phenotypes.

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.

Candidate genes associated with heat stress and breeding strategies to relieve its effects in dairy cattle: a deeper insight into the genetic architecture and immune response to heat stress

The need for food products of animal origin is increasing worldwide. Satisfying these needs in a way that has minimal impact on the environment requires cutting-edge technologies and techniques to enhance the genetic quality of cattle. Heat stress (HS), in particular, is affecting dairy cattle with increasing frequency and severity. As future climatic challenges become more evident, identifying dairy cows that are more tolerant to HS will be important for breeding dairy herds that are better adapted to future environmental conditions and for supporting the sustainability of dairy farming. While research into the genetics of HS in the context of the effect of global warming on dairy cattle is gaining momentum, the specific genomic regions involved in heat tolerance are still not well documented. Advances in omics information, QTL mapping, transcriptome profiling and genome-wide association studies (GWAS) have identified genomic regions and variants associated with tolerance to HS. Such studies could provide deeper insights into the genetic basis for response to HS and make an important contribution to future breeding for heat tolerance, which will help to offset the adverse effects of HS in dairy cattle. Overall, there is a great interest in identifying candidate genes and the proportion of genetic variation associated with heat tolerance in dairy cattle, and this area of research is currently very active worldwide. This review provides comprehensive information pertaining to some of the notable recent studies on the genetic architecture of HS in dairy cattle, with particular emphasis on the identified candidate genes associated with heat tolerance in dairy cattle. Since effective breeding programs require optimal knowledge of the impaired immunity and associated health complications caused by HS, the underlying mechanisms by which HS modulates the immune response and renders animals susceptible to various health disorders are explained. In addition, future breeding strategies to relieve HS in dairy cattle and improve their welfare while maintaining milk production are discussed.

Animal board invited review: Quantification of resilience in farm animals

Resilience, when defined as the capacity of an animal to respond to short-term environmental challenges and to return to the prechallenge status, is a dynamic and complex trait. Resilient animals can reinforce the capacity of the herd to cope with often fluctuating and unpredictable environmental conditions. The ability of modern technologies to simultaneously record multiple performance measures of individual animals over time is a huge step forward to evaluate the resilience of farm animals. However, resilience is not directly measurable and requires mathematical models with biologically meaningful parameters to obtain quantitative resilience indicators. Furthermore, interpretive models may also be needed to determine the periods of perturbation as perceived by the animal. These applications do not require explicit knowledge of the origin of the perturbations and are developed based on real-time information obtained in the data during and outside the perturbation period. The main objective of this paper was to review and illustrate with examples, different modelling approaches applied to this new generation of data (i.e., with high-frequency recording) to detect and quantify animal responses to perturbations. Case studies were developed to illustrate alternative approaches to real-time and post-treatment of data. In addition, perspectives on the use of hybrid models for better understanding and predicting animal resilience are presented. Quantification of resilience at the individual level makes possible the inclusion of this trait into future breeding programmes. This would allow improvement of the capacity of animals to adapt to a changing environment, and therefore potentially reduce the impact of disease and other environmental stressors on animal welfare. Moreover, such quantification allows the farmer to tailor the management strategy to help individual animals to cope with the perturbation, hence reducing the use of pharmaceuticals, and decreasing the level of pain of the animal.

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.

A promising resilience parameter for breeding: the use of weight and feed trajectories in growing pigs

BACKGROUND

Increasing resilience is a priority in modern pig breeding. Recent research shows that general resilience can be quantified via variability in longitudinal data. The collection of such longitudinal data on weight, feed intake and feeding behaviour in pigs has been facilitated by the development of technologies such as automated feeding stations. The goal of this study was to investigate resilience traits, which were estimated as deviations from longitudinal weight, feed intake and feeding behaviour data during the finishing phase. A dataset with 324,207 records between the age of 95 and 155 days on 5,939 Piétrain pigs with known pedigree and genomic information was used. We provided guidelines for a rigid quality control of longitudinal body weight data, as we found that outliers can significantly affect results. Gompertz growth curve analysis, linear modelling and trajectory analyses were used for quantifying resilience traits.

RESULTS

To our knowledge, this is the first study comparing resilience traits from longitudinal body weight, feed intake and feeding behaviour data in pigs. We demonstrated that the resilience traits are lowly to moderately heritable for deviations in body weight (h² = 2.9%-20.2%), in feed intake (9.4%-23.3%) and in feeding behaviour (16.2%-28.3%). Additionally, these traits have good predictive abilities in cross-validation analyses. Deviations in individual body weight and feed intake trajectories are highly correlated (r_g = 0.78) with low to moderate favourable genetic correlations with feed conversion ratio (r_g = 0.39-0.49). Lastly, we showed that some resilience traits, such as the natural logarithm of variances of observed versus predicted body weights (lnvar_weight), are more robust to lower observation frequencies and are repeatable over three different time periods of the finishing phase.

CONCLUSIONS

Our results will help future studies investigating resilience traits and resilience-related traits. Moreover, our study provides first results on standardization of quality control and efficient data sampling from automated feeding station data. Our findings will be valuable for breeding organizations as they offer evidence that pigs' general resilience can be selected on with good accuracy. Moreover, this methodology might be extended to other species to quantify resilience based on longitudinal data.

Differences across herds with different dairy breeds in daily milk yield based proxies for resilience

Global sustainability issues such as climate change, biodiversity loss and food security require food systems to become more resource efficient and better embedded in the local environment. This needs a transition towards more diverse, circular and low-input dairy farming systems with animals best suited to the specific environmental conditions. When varying environmental challenges are posed to animals, cows need to become resilient to disturbances they face. This resilience of dairy cows for disturbances can be quantified using sensor features and resilience indicators derived from daily milk yield records. The aim of this study was to explore milk yield based sensor features and resilience indicators for different cattle groups according to their breeds and herds. To this end, we calculated 40 different features to describe the dynamics and variability in milk production of first parity dairy cows. After correction for milk production level, we found that various aspects of the milk yield dynamics, milk yield variability and perturbation characteristics indeed differed across herds and breeds. On farms with a lower breed proportion of Holstein Friesian across cows, there was more variability in the milk yield, but perturbations were less severe upon critical disturbances. Non-Holstein Friesian breeds had a more stable milk production with less (severe) perturbations. These differences can be attributed to differences in genetics, environments, or both. This study demonstrates the potential to use milk yield sensor features and resilience indicators as a tool to quantify how cows cope with more dynamic production conditions and select animals for features that best suit a farms' breeding goal and specific environment.

Exploring the value of genomic predictions to simultaneously improve production potential and resilience of farmed animals

Sustainable livestock production requires that animals have a high production potential but are also highly resilient to environmental challenges. The first step to simultaneously improve these traits through genetic selection is to accurately predict their genetic merit. In this paper, we used simulations of sheep populations to assess the effect of genomic data, different genetic evaluation models and phenotyping strategies on prediction accuracies and bias for production potential and resilience. In addition, we also assessed the effect of different selection strategies on the improvement of these traits. Results show that estimation of both traits greatly benefits from taking repeated measurements and from using genomic information. However, the prediction accuracy for production potential is compromised, and resilience estimates tends to be upwards biased, when families are clustered in groups even when genomic information is used. The prediction accuracy was also found to be lower for both traits, resilience and production potential, when the environment challenge levels are unknown. Nevertheless, we observe that genetic gain in both traits can be achieved even in the case of unknown environmental challenge, when families are distributed across a large range of environments. Simultaneous genetic improvement in both traits however greatly benefits from the use of genomic evaluation, reaction norm models and phenotyping in a wide range of environments. Using models without the reaction norm in scenarios where there is a trade-off between resilience and production potential, and phenotypes are collected from a narrow range of environments may result in a loss for one trait. The study demonstrates that genomic selection coupled with reaction-norm models offers great opportunities to simultaneously improve productivity and resilience of farmed animals even in the case of a trade-off.

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.

Milk metabolites can characterise individual differences in animal resilience to a nutritional challenge in lactating dairy goats

The aim of this study is built in two phases: to quantify the ability of novel milk metabolites to measure between-animal variability in response and recovery profiles to a short-term nutritional challenge, then to derive a resilience index from the relationship between these individual variations. At two different stages of lactation, sixteen lactating dairy goats were exposed to a 2-d underfeeding challenge. The first challenge was in late lactation, and the second was carried out on the same goats early in the following lactation. During the entire experiment period, samples were taken at each milking for milk metabolite measures. For each metabolite, the response profile of each goat was characterised using a piecewise model for describing the dynamic pattern of response and recovery profiles after the challenge relative to the start of the nutritional challenge. Cluster Analysis identified three types of response/recovery profiles per metabolite. Using cluster membership, multiple correspondence analyses (MCAs) were performed to further characterise response profile types across animals and metabolites. This MCA analysis identified three groups of animals. Further, discriminant path analysis was able to separate these groups of multivariate response/recovery profile type based on threshold levels of three milk metabolites: β-hydroxybutyrate, free glucose and uric acid. Further analyses were done to explore the possibility of developing an index of resilience from milk metabolite measures. Different types of performance response to short-term nutritional challenge can be distinguished using multivariate analyses of a panel of milk metabolites.

Reproductive performance of rabbit females from three paternal lines with a different potential for growth rate and resilience

A total of 197 nulliparous rabbits (from three paternal lines) were used to test potential strategies to overcome the consequences on reproduction associated with the selection for high growth rate. The R line was selected for growth rate during the growing period for 37 generations, the RF line was founded through a high selection intensity of elite animals of the R line, and the RFLP line, which was obtained by backcrossing RF animals with the LP line (a long-lived productive maternal line, characterised by high resilience). BW, perirenal fat thickness, fertility, daily feed intake, milk yield and blood metabolites of females were controlled from 1st artificial insemination to 3rd parturition. Litter size, litter weight, individual weight and feed ingestion of kits were controlled from birth to weaning. Our results show that RF females were significantly lighter than R and RFLP females throughout the trial (-5.0%; P < 0.05). Furthermore, RF animals had a higher fertility rate than RFLP females, at first cycle (+10.5 percentage points; P < 0.05). However, RFLP had a higher fertility rate than RF females at second cycle (+21.5 percentage points; P < 0.01). On average, RFLP females had higher perirenal fat thickness than R females at parturition (+3.0%; P < 0.05) and higher daily feed intake than of R and RF females during gestation and late lactation (+9.7 and +8.7%, respectively; P < 0.05). RFLP females produced more milk than R and RF females in the two first lactations (+18.5%; P < 0.001). In the first three parturitions, R females delivered fewer kits born alive (-1.7 kits than RF and RFLP; P < 0.05). In addition, R females' blood had a higher concentration of glutamine and glutamate than RFLP (+24 and +22.7%, respectively; P < 0.05). RFLP litters were heavier than both R and RF litters throughout lactation. However, R kits were heavier at birth than RF and RFLP (+7.9%). Results suggest that the foundation of a paternal line using elite animals could generate females with better early reproductive performance. In addition, backcrossing the RF line with a maternal LP line resulted in a genetic line whose females had a different resource allocation strategy to foster reproduction during the studied period.

Multiphasic mixed growth models for turkeys

Growth models are important for optimization of feed formulation and breeding programs in turkeys. The objectives of this study were 1) to compare sex and line differences in turkeys in parameter estimates of mono- and di-phasic Gompertz growth models, and 2) to evaluate mono and diphasic mixed Gompertz growth models to determine the variation in parameter estimates in a group of female line turkey toms. A total of 1,056 manually recorded weekly average body weight (BW) observations from male and female turkeys of a male and female line from weeks 1 to 20 were used for objective 1. Daily median values of automatically collected individual BW of female line turkey toms were used for objective 2 and random components associated with individual subject animals related to mature weight and/or timing of maximum gain during each phase were introduced in the Gompertz model. Growth curve shapes were different between male line toms, male line hens, female line toms, and female line hens (P < 0.001). However, inflection points were similar between male and female line toms and between male and female line hens (14.06 vs. 13.72 wk and 11.22 and 10.71 wk, respectively), while mature BW differed between lines by 6.49 and 3.81 kg for toms and hens, respectively. The normalized growth rate constant (growth rate constant corrected for mature weight) was around the same magnitude between male and female line toms (0.0031 vs. 0.0038, respectively), but slightly lower in male line hens compared to female line hens (0.0072 vs. 0.0091, respectively). Diphasic Gompertz models described growth better in all line × sex combinations compared to the monophasic models (P < 0.001) and mixed diphasic Gompertz models showed improved fit over mixed monophasic Gompertz models. The correlation structure of the random components identified that individuals with a higher mature weight had a later inflection point and lower growth rate coefficients. These results provide tools for improved breeding practices and a structure to evaluate the effects of dietary or environmental factors on growth trajectories.

Joint Models to Predict Dairy Cow Survival from Sensor Data Recorded during the First Lactation

Early predictions of cows' probability of survival to different lactations would help farmers in making successful management and breeding decisions. For this purpose, this research explored the adoption of joint models for longitudinal and survival data in the dairy field. An algorithm jointly modelled daily first-lactation sensor data (milk yield, body weight, rumination time) and survival data (i.e., time to culling) from 6 Holstein dairy farms. The algorithm was set to predict survival to the beginning of the second and third lactations (i.e., second and third calving) from sensor observations of the first 60, 150, and 240 days in milk of cows' first lactation. Using 3-time-repeated 3-fold cross-validation, the performance was evaluated in terms of Area Under the Curve and expected error of prediction. Across the different scenarios and farms, the former varied between 45% and 76%, while the latter was between 3.5% and 26%. Significant results were obtained in terms of expected error of prediction, meaning that the method provided survival probabilities in line with the observed events in the datasets (i.e., culling). Furthermore, the performances were stable among farms. These features may justify further research on the use of joint models to predict the survival of dairy cattle.

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.

Research Note: Genome-wide association study for natural antibodies and resilience in a purebred layer chicken line

Resilience is the capacity of an animal to be minimally affected by disturbances or rapidly return to the state pertained before exposure to a disturbance. Resilience indicators can be estimated from longitudinal production data, using deviations of observed from expected production levels. One component of resilience is disease resilience, which includes general disease resistance. Natural antibodies (NAbs) are an indicator trait for general disease resistance. The aim of this study was to perform a genome-wide association study (GWAS) for resilience indicators and NAbs in a Rhode Island purebred layer line and study potential overlap in genomic regions detected for these traits. For 2,494 hens, deviations (i.e., differences) between observed weekly egg production and expected weekly egg production were calculated. Resilience indicators were then defined as the natural logarithm of the variance of deviations, skewness of deviations, and lag-one autocorrelation of deviations. For a subset of 1,221 hens genotyped with the 60 K Illumina SNP BeadChip, NAbs binding keyhole-limpet hemocyanin were available (isotypes IgM and IgG). Heritabilities, estimated with a linear mixed animal model, were 0.39 for IgM and 0.20 for IgG, and ranged from 0.03 to 0.18 for the resilience indicators. No significant associations were found in the GWAS, except for a single chromosomal region for the skewness of egg deviations in wk 25 to 83 of the laying period. The absence of significant peaks for NAbs and resilience indicators suggests that there are no genes with major effect and that the traits are likely under polygenic control in this line.

Immune mechanisms, resistance genes, and their roles in the prevention of mastitis in dairy cows

Mastitis is one of the most important diseases of the mammary gland. The increased incidence of this disease in cows is due to the breeding of dairy cattle for higher yields, which is accompanied by an increased susceptibility to mastitis. Therefore, the difficulty involved with preventing this disease has increased. An integral part of current research is the elimination of mastitis in order to reduce the consumption of antibiotic drugs, thereby reducing the resistance of microorganisms and decreasing companies' economic losses due to mastitis (i.e. decreased milk yield, increased drug costs, and reduced milk supply). Susceptibility to mastitis is based on dairy cows' immunity, health, nutrition, and welfare. Thus, it is important to understand the immune processes in the body in order to increase the resistance of animals. Recently, various studies have focused on the selection of mastitis resistance genes. An important point is also the prevention of mastitis. This publication aims to describe the physiology of the mammary gland along with its immune mechanisms and to approximate their connection with potential mastitis resistance genes. This work describes various options for mastitis elimination and focuses on genetic selection and a closer specification of resistance genes to mastitis. Among the most promising resistance genes for mastitis, we consider CD14, CXCR1, lactoferrin, and lactoglobulin.