Body Weight Deviations as Indicator for Resilience in Layer Chickens

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.

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.

Sequence variants associated with resilient responses in growing pigs

The current work aimed to identify genomic regions and candidate genes associated with resilience in pigs. In previous work, we proposed the body weight deviation from the expected growth curve (ΔBW) and the increase of the positive acute-phase protein haptoglobin (ΔHP) after a vaccine challenge as resilience indicators which may be improved through selective breeding in pigs. Individuals with steady growth rate and minor activation of haptoglobin (high ΔBW and low ΔHP values) were considered resilient. In contrast, pigs with perturbed growth rate and high activation of haptoglobin (low ΔBW and high ΔHP values) were considered susceptible. Both ∆BW and ∆HP were simultaneously considered to select the most resilient (N = 40) and susceptible (N = 40) pigs. A genome-wide association study was carried out for the pigs' response classification to the challenge test using whole-genome sequence data (7,760,720 variants). Eleven associated genomic regions were identified, harbouring relevant candidate genes related to the immune response (such as pro- and anti-inflammatory responses) and growth pathways. These associated genomic regions harboured 41 potential functional mutations (frameshift, splice donor, splice acceptor, start loss and stop loss/gain) in candidate genes. Overall, this study advances our knowledge about the genetic determinism of resilience, highlighting its polygenic nature and strong relationship with immunity and growth.

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.

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.

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.

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.

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.

Disentangling the dynamics of energy allocation to develop a proxy for robustness of fattening pigs

BACKGROUND

There is a growing need to improve robustness of fattening pigs, but this trait is difficult to phenotype. Our first objective was to develop a proxy for robustness of fattening pigs by modelling the longitudinal energy allocation coefficient to growth, with the resulting environmental variance of this allocation coefficient considered as a proxy for robustness. The second objective was to estimate its genetic parameters and correlations with traits under selection and with phenotypes that are routinely collected. In total, 5848 pigs from a Pietrain NN paternal line were tested at the AXIOM boar testing station (Azay-sur-Indre, France) from 2015 to 2022. This farm is equipped with an automatic feeding system that records individual weight and feed intake at each visit. We used a dynamic linear regression model to characterize the evolution of the allocation coefficient between the available cumulative net energy, which was estimated from feed intake, and cumulative weight gain during the fattening period. Longitudinal energy allocation coefficients were analysed using a two-step approach to estimate both the genetic variance of the coefficients and the genetic variance in their residual variance, which will be referred to as the log-transformed squared residual (LSR).

RESULTS

The LSR trait, which could be interpreted as an indicator of the response of the animal to perturbations/stress, showed a low heritability (0.05 ± 0.01), a high favourable genetic correlation with average daily growth (- 0.71 ± 0.06), and unfavourable genetic correlations with feed conversion ratio (- 0.76 ± 0.06) and residual feed intake (- 0.83 ± 0.06). Segmentation of the population in four classes using estimated breeding values for LSR showed that animals with the lowest estimated breeding values were those with the worst values for phenotypic proxies of robustness, which were assessed using records routinely collected on farm.

CONCLUSIONS

Results of this study show that selection for robustness, based on estimated breeding values for environmental variance of the allocation coefficients to growth, can be considered in breeding programs for fattening pigs.

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.

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.

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.

Uganda chicken genetic resources: I. phenotypic and production characteristics

The extent of diversity in the indigenous chicken breeds of Uganda was assessed for their potential utilisation in breeding programmes. A total of 293 indigenous-chicken-keeping households in villages across 35 districts forming 12 sub-regions of the four regions were randomly sampled for 586 mature chickens of both sexes. While only 20.8% of households were headed by women, 40.6% of indigenous chicken keepers were women. The production objectives mainly focused on chickens as sources of income from sales and household food. The chickens were predominantly managed in a scavenging (94.2%) feeding system in a mainly (96.9%) mixed crop-livestock system. The average flock size was 19.80 ± 1.21 chickens with 3.83 ± 0.29 laying hens, each producing an average of 13.41 ± 0.20 eggs/clutch and 40.3 ± 0.6 eggs/annum of 5.20 ± 0.03-month egg production age. Normal-feather strains predominated (>90%), with scattered incidences of naked neck, frizzles, polydactyl, and ptilopody traits in both sexes. Dark (49.0% hen; 43.8% cock) and white (38.3% hen; 42.4% cock) skin colours were most common among the chickens compared to yellow skin. However, yellow-coloured shanks were proportionally the most observed (41% cock; 29% hen). Orange and brown iris (eye) colours were the most common in both sexes. The hens commonly had small round earlobes with varying colours, while cocks had large oval-shaped, mainly red (70%) earlobes. The single-comb type was dominant in both sexes, with wattles almost universally present. Frizzle and polydactyl allele frequencies were significantly lower (p < 0.05) than the expected Mendelian proportions, indicating a possible state of endangerment. Meanwhile, the estimated allele frequencies of ptilopody, tufted-crest, and rose comb alleles in the population were similar (p > 0.05) to the expected Mendelian frequencies. However, these strains did not show any significant (p > 0.05) influence on the body weight or the linear morphometric estimates except for being marginally higher than the normal strains. The phenotypic correlations of body weight and morphometric traits ranged from 0.457 to 0.668 and 0.292 to 0.454 in cocks and hens, respectively. These findings provide hints about the prospects for improved performance with modifications in the production environment. The wide phenotypic diversity would support management efforts for their sustainable utilisation and preservation.

Growth performance, survivability and profitability of improved smallholder chicken genetics in Nigeria: A COVID-19 intervention study

The impact of COVID-19 pandemic on smallholder farming households (SFH) includes increased poverty, and loss of livelihoods. Provision of livestock to SFH is a helpful intervention to mitigate this impact. This study provided a total of 150 smallholder poultry farmers, randomly selected from three states (Kebbi, Nasarawa, and Imo) in Nigeria, with ten 5-week-old chickens (mixed sexes) each, of either FUNAAB Alpha or Noiler chicken genetics. The improved, dual-purpose chickens were evaluated for growth performance (GP), survivability and profitability. The birds were managed under semi-scavenging production system. Body weight, mortality, and cost of production (COP) were recorded every 4 weeks until 21 weeks of age. Profitability was a function of the COP, and the selling price for live-birds (cocks). Body weight of Noiler (1,927 g) birds was not significantly (p > .05) higher than FUNAAB Alpha (1,792 g) at 21 weeks. Agroecology and genetics had significant (p < .05) effects on GP and survivability. Survivability of FUNAAB Alpha was higher (p < .05) than Noiler, with Nasarawa (81%-96%), having the highest (p < .0001) survival rate compared to Imo (62%-81%), and Kebbi (58%-75%). At 21 weeks, the number of cocks and hens differed significantly (p < .05) within the states (Imo: 2.4 ± .2 and 5.4 ± .3; Kebbi: 2.6 ± .2 and 5.5 ± .3; and Nasarawa: 2.9 ± .2 and 5.8 ± .3). Nasarawa (NGN 7,808; USD 19) ranked best for profitability, followed by Kebbi (NGN 6,545; USD 16) and Imo (NGN 5,875; USD 14). Overall, this study demonstrates that provision of improved chickens to vulnerable SFH in Nigeria holds great potential for economic growth, and resilience during emergencies, such as the COVID-19 pandemic.

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.

Genetic and phenotypic correlations between Aleutian disease tests with body weight, growth, and feed efficiency traits in mink

The ineffectiveness of vaccination, medicine, and culling strategy leads mink farmers to control Aleutian disease (AD) by selecting AD-resilient mink based on AD tests. However, the genetic background of AD tests and their correlations with economically important or AD-resilient traits are limited. This study estimated the genetic and phenotypic correlations between four AD tests and seven body weight (BW) traits, six growth parameters from the Richards growth model, and eight feed-related traits. Univariate models were used to test the significance (P < 0.05) of fixed effects (sex, color type, AD test year, birth year, and row-by-year), random effects (additive genetic, maternal genetic, and permanent environmental), and a covariate of age using ASReml 4.1. Likewise, pairwise bivariate analyses were conducted to estimate the phenotypic and genetic correlations among the studied traits. Both antigen- and virus capsid protein-based enzyme-linked immunosorbent assay tests (ELISA-G and ELISA-P) showed significant (P < 0.05) moderate positive genetic correlations (±SE) with maturation rate (from 0.36 ± 0.18 to 0.38 ± 0.19). ELISA-G showed a significant negative genetic correlation (±SE) with average daily gain (ADG, -0.37 ± 0.16). ELISA-P showed a significant positive moderate genetic correlation (±SE) with off-feed days (DOF, 0.42 ± 0.17). These findings indicated that selection for low ELISA scores would reduce the maturation rate, increase ADG (by ELISA-G), and minimize DOF (by ELISA-P). The iodine agglutination test (IAT) showed significant genetic correlations with DOF (0.73 ± 0.16), BW at 16 weeks of age (BW16, 0.45 ± 0.23), and BW at harvest (HW, -0.47 ± 0.20), indicating that selection for lower IAT scores would lead to lower DOF and BW16, and higher HW. These estimated genetic correlations suggested that the selection of AD tests would not cause adverse effects on the growth, feed efficiency, and feed intake of mink. The estimates from this study might strengthen the previous finding that ELISA-G could be applied as a reliable and practical indicator trait in the genetic selection of AD-resilient mink in AD-positive farms.

Assessing various environmental descriptors with respect to genotype x environment interaction for milk production traits in Bavarian Fleckvieh cattle

Genotype x environment (GxE) interaction for production traits in Fleckvieh cattle was assessed by means of various environmental descriptors (EDs). It was also of particular interest to search for EDs useful for studying differing robustness or resilience of individuals which implies reasonable GxE interaction. The set of studied EDs included farm/herd environment (e.g. herd size, housing/feeding regimes, herd production level), geographical location (e.g. height above sea level), temperature humidity index and fat-to-protein ratio. Milk, fat and protein yield deviations for the first test day, the first 100 days and the 305-day-yield of the first lactation were used as phenotypes. Genetic correlations were estimated with bi- and multivariate sire models in case of categorical EDs and reaction norm sire models in case of continuous EDs. Further, rank correlations of bulls' estimated breeding values were calculated for different environmental levels/classes. Almost all estimated genetic correlations were significantly different from 1, but not <0.93. There were some exceptions for extremely different classes/levels of some EDs (e.g. average herd-year production, organic vs conventional farm systems, way of concentrated feed supply). Rank correlations were rarely below 0.95. In summary, no substantial GxE interactions for milk production traits were found with the studied EDs.

The immune-neuroendocrine system, a key aspect of poultry welfare and resilience

There is increasing societal concern regarding the negative impact of intensive poultry production on animal welfare, human health, and on the environment. This is leading to the inclusion of animal welfare as an imperative aspect for sustainable production. Certain environmental factors may challenge domesticated birds, resulting in poor health and welfare status. Resilience is the capacity to rapidly return to prechallenge status after coping with environmental stressors, thus resilient individuals have better chances to maintain good health and welfare. Immune-neuroendocrine system, thoroughly characterized in the domestic bird species, is the physiological scaffold for stress coping and health maintenance, influencing resilience and linking animal welfare status to these vital responses. Modern domestic bird lines have undergone specific genetic selective pressures for fast-growing, or high egg-production, leading to a diversity of birds that differ in their coping capacities and resilience. Deepening the knowledge on pro/anti-inflammatory milieus, humoral/cell-mediated immune responses, hormonal regulations, intestinal microbial communities and mediators that define particular immune and neuroendocrine configurations will shed light on coping strategies at the individual and population level. The understanding of the profiles leading to differential coping and resilience potential will be highly relevant for improving bird health and welfare in a wider range of challenging scenarios and, therefore, crucial to scientifically tackle long term sustainability.

Broiler resilience to colibacillosis is affected by incubation temperature and post-hatch feeding strategy

Colibacillosis is a poultry disease that negatively affects welfare and causes economic losses. Treatment with antibiotics raises concerns on antimicrobial resistance. Consequently, alternative approaches to enhance poultry resilience are needed. Access to feed and water directly after hatch (early feeding) may enhance resilience at later ages. Additionally, a high eggshell temperature (EST) during mid incubation may improve chick quality at hatch, supporting potential positive effects of early feeding. Effects of EST [37.8°C (control) or 38.9°C (higher)] during mid-incubation (embryo days 7–14) and feeding strategy (early feeding or 48 h delayed feeding) were tested in a 2 × 2 factorial arrangement. At hatch, ~ 1,800 broilers were divided over 36 pens and grown for 6 wk. At d 8 post hatch, avian pathogenic E. coli (APEC) was inoculated intratracheally as model to investigate broiler resilience against respiratory diseases. Incidence and severity of colibacillosis, local infection, and systemic infection were assessed at 6 moments between 3 h and 7 d postinoculation. Broilers were weighed daily during 13 d postinoculation and weekly thereafter. At higher EST, early feeding resulted in higher incidence of systemic infection compared to delayed feeding whereas at control EST, systemic infection was not different between feeding strategies. Regardless of EST, early compared to delayed feeding resulted in lower incidence of local infection, fewer BW deviations, and higher growth until d 35. In conclusion, early feeding could be considered as a strategy to enhance broiler resilience, but only when EST is not too high.

A genome-wide screen for resilient responses in growing pigs

Background

There is a growing interest to decipher the genetic background of resilience and its possible improvement through selective breeding. The objective of the present study was to provide new insights into the genetic make-up of resilience in growing pigs by identifying genomic regions and candidate genes associated with resilience indicators. Commercial Duroc pigs were challenged with an attenuated Aujeszky vaccine at 12 weeks of age. Two resilience indicators were used: deviation from the expected body weight at 16 weeks of age given the growth curve of non-vaccinated pigs (∆BW) and the increase in acute-phase protein haptoglobin at four days post-vaccination (∆HP). Genome-wide association analyses were carried out on 445 pigs, using genotypes at 41,165 single nucleotide polymorphisms (SNPs) and single-marker and Bayesian multiple-marker regression approaches.

Results

Genomic regions on pig chromosomes 2, 8, 9, 11 (∆BW) and 8, 9, 13 (∆HP) were found to be associated with the resilience indicators and explained high proportions of their genetic variance. The genomic regions that were associated explained 27 and 5% of the genetic variance of ∆BW and ∆HP, respectively. These genomic regions harbour promising candidate genes that are involved in pathways related to immune response, response to stress, or signal transduction (CD6, PTGDR2, IKZF1, RNASEL and MYD88), and growth (GRB10 and LCORL).

Conclusions

Our study identified novel genomic regions that are associated with two resilience indicators (∆BW and ∆HP) in pigs. These associated genomic regions harbour potential candidate genes involved in immune response and growth pathways, which emphasise the strong relationship between resilience and immune response.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00739-1.

A multi-suckling system combined with an enriched housing environment during the growing period promotes resilience to various challenges in pigs

Little is known about the impact of social and environmental enrichment on improving livestock resilience, i.e. the ability to quickly recover from perturbations. We evaluated the effect of an alternative housing system (AHS) on resilience of pigs, as compared to conventional housing (CONV). The AHS consisted of multi-litter housing during lactation, delayed weaning, extra space allowance and environmental enrichment at all times. We assessed recovery to a 2 h-transport challenge, an LPS injection, 2 h-heat stress and a biopsy wound in 96 pigs. Additionally, indicators of long-term “wear and tear” on the body were determined. AHS pigs had better physiological recoveries with quicker returns to baseline in the transport and LPS challenges, showed lower cortisol accumulation in hairs and lower variance in weight gain over the experimental period compared to conventionally-housed (CONV) pigs. They also had higher levels of natural antibodies binding KLH than CONV pigs. Their response to heat stress revealed a different strategy compared to CONV pigs. Taken together, AHS pigs appear to be more resilient and experience less chronic stress. Enhancing welfare by provision of social and environmental enrichment that better meets the behavioural needs of pigs seems to be a promising approach to improve their resilience.

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

Background

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

Results

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

Conclusions

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

Supplementary Information

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

Divergent selection for natural antibodies in poultry in the presence of a major gene

Background

Natural antibodies (NAb) are antibodies that are present in a healthy individual without requiring previous exposure to an exogenous antigen. Selection for high NAb levels might contribute to improved general disease resistance. Our aim was to analyse the genetic background of NAb based on a divergent selection experiment in poultry, and in particular the effect of a polymorphism in the TLR1A gene.

Methods

The study population consisted of a base population from a commercial pure-bred elite white leghorn layer line and seven generations of birds from a High and Low selection line. Birds were selected for total KLH-binding NAb titer (IgTotal). An enzyme-linked immunosorbent assay was performed to determine NAb titers in blood plasma for IgTotal and the antibody isotypes IgM and IgG. NAb titers were available for 10,878 birds. Genotypes for a polymorphism in TLR1A were determined for chickens in generations 5, 6 and 7. The data were analysed using mixed linear animal models.

Results

The heritability estimate for IgM was 0.30 and higher than that for IgG and IgTotal (0.12). Maternal environmental effects explained 2 to 3% of the phenotypic variation in NAb. Selection for IgTotal resulted in a genetic difference between the High and Low line of 2.4 titer points (5.1 genetic standard deviation) in generation 7. For IgM, the selection response was asymmetrical and higher in the Low than the High line. The frequency of the TLR1A C allele was 0.45 in the base population and 0.66 and 0.04 in generation 7 of the High and Low line, respectively. The TLR1A polymorphism had large and significant effects on IgTotal and IgM. Estimated genotypic effects suggest full dominance of the TLR1A C allele. Significant TLR1A by generation interactions were detected for IgM and IgTotal.

Conclusions

The effect of a polymorphism in the TLR1A gene on IgTotal and IgM NAb was confirmed. Furthermore, we provide experimental verification of changes in allele frequencies at a major gene with dominant gene action on a quantitative trait that is subjected to mass selection. TLR1A by generation interactions indicate sensitivity to environmental factors.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00715-9.

Impact of selection for birth weight variability on reproductive longevity: A mice model

Uniformity, understood as a similar performance in relevant livestock traits, such as birth weight within the litter, is being included as one of the selection objectives in breeding programmes, especially for polytocous livestock species. A divergent selection experiment for birth weight within‐litter variability in mice during 23 generations showed that homogeneous animals were better for litter size, survival and feed efficiency but less heavy than heterogeneous animals. The aim of this study was to compare the reproductive longevity in both divergent lines as time to the end of the reproductive period. Two generations from both lines with an initial number of 43 females and 43 males were mated one to one and stayed together to have consecutive parturitions until the end of the reproductive life. Females were discarded when the time elapsed from the last parturition was longer than 63 days. The time to the end of the reproductive period between both lines was compared by fitting a Cox proportional hazard regression model adjusting for line, generation and its interaction. The rate of parturitions in both lines was also compared using a Prentice–Williams–Peterson model adjusted for the same effects. The low variability line was associated with a higher parturition rate, e.g., adjusted hazard ratio was 2.93 (95% CI 2.17–3.94). The Cox model showed that the low variability females also presented benefits of time to the end of the reproductive period, with an adjusted hazard ratio of 0.26 (95% CI 0.16–0.41). The median of reproductive days was 55.50 in the high variability line whilst the median was 252.50 days for the low variability females. The homogeneous line presented important reproductive advantages suggesting higher robustness and animal welfare. Further research should confirm whether the findings presented here of a better performance in the low variability line could be properly applied to some livestock species.

Development of resilience indicator traits based on daily step count data for dairy cattle breeding

Background

Resilient animals are minimally affected by disturbances, such as diseases and heat stress, and quickly recover. Daily activity data can potentially indicate resilience, because resilient animals likely keep variations due to disturbances that threat animal homeostasis at a low magnitude. We used daily step count of cows to define resilience indicators based on theory, exploratory analysis and literature, and then investigated if they can be used to genetically improve resilience by estimating heritability and repeatability, and genetic associations with other resilience-related traits, i.e. health traits, longevity, fertility, and body condition score (BCS).

Results

Two groups of resilience indicators were defined: indicators describing (1) mean step count at different lactation stages for individual cows, and (2) fluctuations in step count from individual step count curves. Heritability estimates were highest for resilience indicators describing mean step count, from 0.22 for the 2-week period pre-partum to 0.45 for the whole lactation. High mean step count was consistently, but weakly, genetically correlated with good health, fertility, and longevity, and high BCS. Heritability estimates of resilience indicators describing fluctuations ranged from 0.01 for number of step count drops to 0.15 for the mean of negative residuals from individual curves. Genetic correlations with health traits, longevity, fertility, and BCS were mostly weak, but were moderate and favorable for autocorrelation of residuals (− 0.33 to − 0.44) and number of step count drops (− 0.44 to − 0.56) with hoof health, fertility, and BCS. Resilience indicators describing variability of residuals and mean of negative residuals showed strong genetic correlations with mean step count (0.86 to 0.95, absolute), which suggests that adjustment for step count level is needed. After adjustment, ‘mean of negative residuals’ was highly genetically correlated with hoof health, fertility, and BCS.

Conclusions

Mean step count, autocorrelation and mean of negative residuals showed most potential as resilience indicators based on resilience theory, heritability, and genetic associations with health, fertility, and body condition score. Other resilience indicators were heritable, but had unfavorable genetic correlations with several health traits. This study is an important first step in the exploration of the use of activity data to breed more resilient livestock.

Genome-wide association study of disease resilience traits from a natural polymicrobial disease challenge model in pigs identifies the importance of the major histocompatibility complex region

Infectious diseases cause tremendous financial losses in the pork industry, emphasizing the importance of disease resilience, which is the ability of an animal to maintain performance under disease. Previously, a natural polymicrobial disease challenge model was established, in which pigs were challenged in the late nursery phase by multiple pathogens to maximize expression of genetic differences in disease resilience. Genetic analysis found that performance traits in this model, including growth rate, feed and water intake, and carcass traits, as well as clinical disease phenotypes, were heritable and could be selected for to increase disease resilience of pigs. The objectives of the current study were to identify genomic regions that are associated with disease resilience in this model, using genome-wide association studies and fine-mapping methods, and to use gene set enrichment analyses to determine whether genomic regions associated with disease resilience are enriched for previously published quantitative trait loci, functional pathways, and differentially expressed genes subject to physiological states. Multiple quantitative trait loci were detected for all recorded performance and clinical disease traits. The major histocompatibility complex region was found to explain substantial genetic variance for multiple traits, including for growth rate in the late nursery (12.8%) and finisher (2.7%), for several clinical disease traits (up to 2.7%), and for several feeding and drinking traits (up to 4%). Further fine mapping identified 4 quantitative trait loci in the major histocompatibility complex region for growth rate in the late nursery that spanned the subregions for class I, II, and III, with 1 single-nucleotide polymorphism in the major histocompatibility complex class I subregion capturing the largest effects, explaining 0.8–27.1% of genetic variance for growth rate and for multiple clinical disease traits. This single-nucleotide polymorphism was located in the enhancer of TRIM39 gene, which is involved in innate immune response. The major histocompatibility complex region was pleiotropic for growth rate in the late nursery and finisher, and for treatment and mortality rates. Growth rate in the late nursery showed strong negative genetic correlations in the major histocompatibility complex region with treatment or mortality rates (−0.62 to −0.85) and a strong positive genetic correlation with growth rate in the finisher (0.79). Gene set enrichment analyses found genomic regions associated with resilience phenotypes to be enriched for previously identified disease susceptibility and immune capacity quantitative trait loci, for genes that were differentially expressed following bacterial or virus infection and immune response, and for gene ontology terms related to immune and inflammatory response. In conclusion, the major histocompatibility complex and other quantitative trait loci that harbor immune-related genes were identified to be associated with disease resilience traits in a large-scale natural polymicrobial disease challenge. The major histocompatibility complex region was pleiotropic for growth rate under challenge and for clinical disease traits. Four quantitative trait loci were identified across the class I, II, and III subregions of the major histocompatibility complex for nursery growth rate under challenge, with 1 single-nucleotide polymorphism in the major histocompatibility complex class I subregion capturing the largest effects. The major histocompatibility complex and other quantitative trait loci identified play an important role in host response to infectious diseases and can be incorporated in selection to improve disease resilience, in particular the identified single-nucleotide polymorphism in the major histocompatibility complex class I subregion.

Breeding for disease resilience: opportunities to manage polymicrobial challenge and improve commercial performance in the pig industry

Disease resilience, defined as an animal's ability to maintain productive performance in the face of infection, provides opportunities to manage the polymicrobial challenge common in pig production. Disease resilience can deliver a number of benefits, including more sustainable production as well as improved animal health and the potential for reduced antimicrobial use. However, little progress has been made to date in the application of disease resilience in breeding programs due to a number of factors, including (1) confusion around definitions of disease resilience and its component traits disease resistance and tolerance, and (2) the difficulty in characterizing such a complex trait consisting of multiple biological functions and dynamic elements of rates of response and recovery from infection. Accordingly, this review refines the definitions of disease resistance, tolerance, and resilience based on previous studies to help improve the understanding and application of these breeding goals and traits under different scenarios. We also describe and summarize results from a "natural disease challenge model" designed to provide inputs for selection of disease resilience. The next steps for managing polymicrobial challenges faced by the pig industry will include the development of large-scale multi-omics data, new phenotyping technologies, and mathematical and statistical methods adapted to these data. Genome editing to produce pigs resistant to major diseases may complement selection for disease resilience along with continued efforts in the more traditional areas of biosecurity, vaccination and treatment. Altogether genomic approaches provide exciting opportunities for the pig industry to overcome the challenges provided by hard-to-manage diseases as well as new environmental challenges associated with climate change.

Evaluation of the phenotypic and genomic background of variability based on litter size of Large White pigs

BACKGROUND

The genetic background of trait variability has captured the interest of ecologists and animal breeders because the genes that control it could be involved in buffering various environmental effects. Phenotypic variability of a given trait can be assessed by studying the heterogeneity of the residual variance, and the quantitative trait loci (QTL) that are involved in the control of this variability are described as variance QTL (vQTL). This study focuses on litter size (total number born, TNB) and its variability in a Large White pig population. The variability of TNB was evaluated either using a simple method, i.e. analysis of the log-transformed variance of residuals (LnVar), or the more complex double hierarchical generalized linear model (DHGLM). We also performed a single-SNP (single nucleotide polymorphism) genome-wide association study (GWAS). To our knowledge, this is only the second study that reports vQTL for litter size in pigs and the first one that shows GWAS results when using two methods to evaluate variability of TNB: LnVar and DHGLM.

RESULTS

Based on LnVar, three candidate vQTL regions were detected, on Sus scrofa chromosomes (SSC) 1, 7, and 18, which comprised 18 SNPs. Based on the DHGLM, three candidate vQTL regions were detected, i.e. two on SSC7 and one on SSC11, which comprised 32 SNPs. Only one candidate vQTL region overlapped between the two methods, on SSC7, which also contained the most significant SNP. Within this vQTL region, two candidate genes were identified, ADGRF1, which is involved in neurodevelopment of the brain, and ADGRF5, which is involved in the function of the respiratory system and in vascularization. The correlation between estimated breeding values based on the two methods was 0.86. Three-fold cross-validation indicated that DHGLM yielded EBV that were much more accurate and had better prediction of missing observations than LnVar.

CONCLUSIONS

The results indicated that the LnVar and DHGLM methods resulted in genetically different traits. Based on their validation, we recommend the use of DHGLM over the simpler method of log-transformed variance of residuals. These conclusions can be useful for future studies on the evaluation of the variability of any trait in any species.

Fleece rot and dermatophilosis (lumpy wool) in sheep: opportunities and challenges for new vaccines

During prolonged wetting of the fleece, proliferation of bacterial flora often dominated by Pseudomonas aeruginosa or Dermatophilus congolensis can induce dermatitis and fleece damage termed fleece rot and dermatophilosis respectively, which predispose sheep to blowfly strike. A large research effort in the 1980s and 1990s on vaccines to control fleece rot and dermatophilosis met with limited success. This review examines theoretical and technological advances in microbial ecology, pathogenesis, immunology, vaccine development and the characterisation of microbial virulence factors that create new opportunities for development of vaccines against these diseases. Genomic technologies have now created new opportunities for examining microbial dynamics and pathogen virulence in dermatitis. An effective vaccine requires the combination of appropriate antigens with an adjuvant that elicits a protective immune response that ideally provides long-lasting protection in the field. A clinical goal informed by epidemiological, economic and animal welfare values is needed as a measure of vaccine efficacy. Due to dependence of fleece rot and dermatophilosis on sporadic wet conditions for their expression, vaccine development would be expedited by in vitro correlates of immune protection. The efficacy of vaccines is influenced by genetic and phenotypic characteristics of the animal. Advances in understanding vaccine responsiveness, immune defence in skin and immune competence in sheep should also inform any renewed efforts to develop new fleece rot and dermatophilosis vaccines. The commercial imperatives for new vaccines are likely to continue to increase as the animal welfare expectations of society intensify and reliance on pharmacotherapeutics decrease due to chemical resistance, market pressures and societal influences. Vaccines should be considered part of an integrated disease control strategy, in combination with genetic selection for general immune competence and resistance to specific diseases, as well as management practices that minimise stress and opportunities for disease transmission. The strategy could help preserve the efficacy of pharmacotherapeutics as tactical interventions to alleviate compromised welfare when adverse environmental conditions lead to a break down in integrated strategic disease control. P. aeruginosa and D. congolensis are formidable pathogens and development of effective vaccines remains a substantial challenge.

A Methodology to Quantify Resilience in Growing Pigs

Simple Summary

The concept of resilience can be defined as the ability of animals to maintain their productivity despite the stressors that might occur during the rearing period. Despite the growing interest in the genetic determinism of resilience and its possible inclusion in selection objectives, there is no straightforward method to measure this trait. Several quantifying methods estimate resilience based on productivity such as body weight or milk production, or non-productivity related traits such as immunity or stress responses. The objective of this study was to elaborate novel resilience indicators in growing pigs based on both productivity (body weight) and non-productivity (acute-phase proteins) related traits. We propose the deviation from the expected growth curve and the increment of the acute-phase protein haptoglobin, after applying a common vaccine, as resilience indicators in growing pigs under standard rearing conditions. We showed that the suggested resilience indicators are under a genetic control, show a substantial variability in the population, and may be improved through selection.

Abstract

There is a growing concern about the genetic determinism of resilience and its possible implementation in breeding programs. The objective of our study was to elaborate novel resilience indicators in growing pigs based on the deviation from the expected growth curve and the increment of the acute-phase protein haptoglobin (HP) after applying a common vaccine. A total of 445 pigs were vaccinated with an attenuated Aujeszky vaccine at 12 weeks of age. Deviation from the expected body weight (ΔBW) given the growth curve of unvaccinated pigs at 28 days post-vaccination (DPV) and the increment of HP at 4 DPV (ΔHP) were suggested as resilience indicators. Challenged pigs that maintained their productivity and had a minor activation of HP were deemed resilient, whereas pigs that had low ∆BW values and a high activation of HP were deemed susceptible. Pigs were also classified based on ∆BW and ∆HP relative to the expected BW at 28 DPV and to the basal level of HP, respectively. The concordance was high between both methods, indicating that ΔBW and ΔHP are not sensitive to the animal’s expected BW nor the basal level of HP. The heritability estimates were moderate for ∆BW (0.33) and low-to-moderate for ∆HP (0.16). Our study suggests ΔBW and ΔHP as novel resilience indicators in pigs. The suggested indicators capture different aspects of resilience, are easy to measure, and are genetically controlled. Thus, they may be improved through selective breeding. Further analyses are needed to validate our findings.

Validation of resilience indicators by estimating genetic correlations among daughter groups and with yield responses to a heat wave and disturbances at herd level

Resilient cows are minimally affected in their functioning by disturbances, and if affected, they quickly recover. Previously, the variance and autocorrelation of daily deviations from a lactation curve were proposed as resilience indicators. These traits were heritable and genetically associated with good health and longevity. However, it was unknown if selection for these indicators would lead to desired changes in the phenotype. The first aim of this study was to investigate if forward prediction of the resilience indicators in another environment was possible. Therefore, the resilience indicator records were split into 2 subsets, each containing half of the daughters of each sire, split within sire into cows that calved in early year-seasons and cows that calved in more recent year-seasons. Genetic correlations between the subsets were then estimated for each resilience indicator. The second aim was to estimate genetic correlations between the resilience indicators and traits describing production responses to actual disturbances. The disturbances were a heat wave in July 2015 and yield disturbances at herd level. The latter were selected by decreases in mean yield of all primiparous cows in a herd, indicating that a disturbance occurred. The data set used for calculation of the resilience indicators and the traits describing yield responses contained 62,932,794 daily milk yield records on 199,104 primiparous cows. Genetic correlations (r_g) between recent and earlier daughter groups were 1 for both resilience indicators, which suggests that selection will result in changes in the phenotype in the next generation. Furthermore, low variance was genetically correlated with weak response in milk yield to both the heat wave and herd disturbances (r_g 0.47 to 0.97). Low autocorrelation was genetically correlated with reduced perturbation length and quick recovery after the heat wave and herd disturbances (0.28 to 0.97). These results suggest that variance and autocorrelation cover different aspects of resilience, and should be combined in a resilience index. In conclusion, genetic selection for the resilience indicators will likely result in favorable changes in the traits themselves, and in response and recovery to actual disturbances, which confirms that they are useful resilience indicators.

Detection of unrecorded environmental challenges in high-frequency recorded traits, and genetic determinism of resilience to challenge, with an application on feed intake in lambs

BACKGROUND

Resilient animals can remain productive under different environmental conditions. Rearing in increasingly heterogeneous environmental conditions increases the need of selecting resilient animals. Detection of environmental challenges that affect an entire population can provide a unique opportunity to select animals that are more resilient to these events. The objective of this study was two-fold: (1) to present a simple and practical data-driven approach to estimate the probability that, at a given date, an unrecorded environmental challenge occurred; and (2) to evaluate the genetic determinism of resilience to such events.

METHODS

Our method consists of inferring the existence of highly variable days (indicator of environmental challenges) via mixture models applied to frequently recorded phenotypic measures and then using the inferred probabilities of the occurrence of an environmental challenge in a reaction norm model to evaluate the genetic determinism of resilience to these events. These probabilities are estimated for each day (or other time frame). We illustrate the method by using an ovine dataset with daily feed intake (DFI) records.

RESULTS

Using the proposed method, we estimated the probability of the occurrence of an unrecorded environmental challenge, which proved to be informative and useful for inclusion as a covariate in a reaction norm animal model. We estimated the breeding values for sensitivity of the genetic potential for DFI of animals to environmental challenges. The level and slope of the reaction norm were negatively correlated (- 0.46 ± 0.21).

CONCLUSIONS

Our method is promising and appears to be viable to identify unrecorded events of environmental challenges, which is useful when selecting resilient animals and only productive data are available. It can be generalized to a wide variety of phenotypic records from different species and used with large datasets. The negative correlation between level and slope indicates that a hypothetical selection for increased DFI may not be optimal depending on the presence or absence of stress. We observed a reranking of individuals along the environmental gradient and low genetic correlations between extreme environmental conditions. These results confirm the existence of a G [Formula: see text] E interaction and show that the best animals in one environmental condition are not the best in another one.

Stress, health and the welfare of laying hens

It is essential to understand responses to stress and the impact of stress on physiological and behavioural functioning of hens, so as to assess their welfare. The current understanding of stress in laying hens is comprehensively reviewed here. Most research on stress in hens has focussed on the activity of the adrenal glands, with the most common approach being to measure corticosterone, which is the predominant glucocorticoid produced by birds in response to stress. While these measures are useful, there is a need to understand how the brain regulates stress responses in hens. A greater understanding of the sympathoadrenal system and its interaction with the hypothalamo–pituitary–adrenal axis is required. There is also a lack of knowledge about the many other peptides and regulatory systems involved in stress responses in hens. The usefulness of understanding stress in hens in terms of assessing welfare depends on appreciating that different stressors elicit different responses and that there are often differences in responses to, and impacts of, acute and chronic stress. It is also important to establish the actions and fate of stress hormones within target tissues. It is the consequences of these actions that are important to welfare. A range of other measures has been used to assess stress in hens, including a ratio of heterophils to lymphocytes and haematocrit:packed cell-volume ratio and measures of corticosterone or its metabolites in eggs, excreta, feathers and the secretions of the uropygial gland. Measures in eggs have proffered varying results while measures in feathers may be useful to assess chronic stress. There are various studies in laying hens to indicate impacts of stress on the immune system, health, metabolism, appetite, and the quality of egg production, but, generally, these are limited, variable and are influenced by the management system, environment, genetic selection, type of stressor and whether or not the birds are subjected to acute or chronic stress. Further research to understand the regulation of stress responses and the impact of stress on normal functioning of hens will provide important advances in the assessment of stress and, in turn, the assessment of welfare of laying hens.

Genetic analysis of resilience indicators based on milk yield records in different lactations and at different lactation stages

Resilience is the ability of cows to cope with disturbances, such as pathogens or heat waves. To breed for improved resilience, it is important to know whether resilience genetically changes throughout life. Therefore, the aim was to perform a genetic analysis on 2 resilience indicators based on data from 3 periods of the first lactation (d 11-110, 111-210, and 211-340) and the first 3 full lactations, and to estimate genetic correlations with health traits. The resilience indicators were the natural log-transformed variance (LnVar) and lag-1 autocorrelation (r_auto) of daily deviations in milk yield from an expected lactation curve. Low LnVar and r_auto indicate low variability in daily milk yield and quick recovery, and were expected to indicate good resilience. Data of 200,084 first, 155,784 second, and 89,990 third lactations were used. Heritabilities were similar based on different lactation periods (0.12-0.15 for LnVar, 0.05-0.06 for r_auto). However, the heritabilities of the resilience indicators based on full first lactation were higher than those based on lactation periods (0.20 for LnVar, 0.08 for r_auto), due to lower residual variances. Heritabilities decreased from 0.20 in full lactation 1 to 0.19 in full lactation 3 for LnVar and from 0.08 to 0.06 for r_auto. For LnVar, as well as for r_auto, the strongest genetic correlation between lactation periods was between period 2 and 3 (0.97 for LnVar, 0.96 for r_auto) and the weakest between period 1 and 3 (0.81 for LnVar, 0.65 for r_auto). Similarly, for both traits the genetic correlation between full lactations was strongest between lactations 2 and 3 (0.99 for LnVar, 0.95 for r_auto) and weakest between lactations 1 and 3 (0.91 for LnVar, 0.71 for r_auto). For LnVar, genetic correlations with resilience-related traits, such as udder health, ketosis, and longevity, adjusted for correlations with milk yield, were almost always favorable (-0.59 to 0.02). In most cases these genetic correlations were stronger based on full lactations than on lactation periods. Genetic correlations were similar across full lactations, but the correlation with udder health increased substantially from -0.31 in lactation 1 to -0.51 in lactation 3. For r_auto, genetic correlations with resilience-related traits were always favorable in lactation period 1 and in most full lactations, but not in the other lactation periods. However, correlations were weak (-0.27 to 0.15). Therefore, as a resilience indicator for breeding, LnVar is preferred over r_auto. A multitrait index based on estimated breeding values for LnVar in lactations 1, 2, and 3 is recommended to improve resilience throughout the lifetime of a cow.

Lactation curve model with explicit representation of perturbations as a phenotyping tool for dairy livestock precision farming

In the context of dairy farming, ruminant females often face challenges inducing perturbations that affect their performance and welfare. A key issue is how to assess the effect of perturbations and provide metrics to quantify how animals cope with their environment. Milk production dynamics are good candidates to address this issue: i) they are easily accessible, ii) overall dynamics throughout lactation process are well described and iii) perturbations are visible through milk losses. In this study, a perturbed lactation model (PLM) with explicit representation of perturbations was developed. The model combines two components: i) the unperturbed lactation model that describes a theoretical lactation curve, assumed to reflect female production potential and ii) the perturbation model that describes all the deviations from the unperturbed lactation model with four parameters: starting date, intensity and shape (collapse and recovery). To illustrate the use of the PLM as a phenotyping tool, it was fitted on a data set of 319 complete lactations from 181 individual dairy goats. A total of 2 354 perturbations were detected, with an average of 7.40 perturbations per lactation. Loss of milk production for the whole lactation due to perturbations varied between 2 and 19% of the milk production predicted by the unperturbed lactation model. The number of perturbations was not the major factor explaining the loss of milk production, suggesting that there are different types of animal response to challenges. By incorporating explicit representation of perturbations in a lactation model, it was possible to determine for each female the potential milk production, characteristics of each perturbation and milk losses induced by perturbations. Further, it was possible to compare animals and analyze individual variability. The indicators produced by the PLM are likely to be useful to move from raw data to decision support tools in dairy production.

Exploring simultaneous perturbation profiles in milk yield and body weight reveals a diversity of animal responses and new opportunities to identify resilience proxies

Livestock husbandry aims to manage the environment in which animals are reared to enable them to express their production potential. However, animals are often confronted with perturbations that affect their performance. Evaluating effects of these perturbations on animal performance could provide metrics to quantify and understand how animals cope with their environment, and therefore to better manage them. Body weight (BW) and milk yield (MY) dynamics over lactation may be used for this purpose. The goal of this study was to estimate an unperturbed performance trajectory using a differential smoothing approach on both MY and BW time series, and then to identify the perturbations and extract their phenotypic features. Daily MY and BW records from 490 primiparous Holstein cows from 33 commercial French herds were used. From the fitting procedure, estimated unperturbed performance trajectories of BW and MY were clustered into 3 groups. After the fitting procedure, 1,754 deviations were detected in the MY time series and 964 were detected in the BW time series across all cows. Overall, 425 of these deviations were detected during the same period (±10 d) in both MY and BW time series, 76 of which started at the same time. Results suggest that combining various individual dynamic measures and revealing the relationship that exists between them could be of great value in obtaining reliable estimates of resilience components in large populations.

Why breed disease-resilient livestock, and how?

BACKGROUND

Fighting and controlling epidemic and endemic diseases represents a considerable cost to livestock production. Much research is dedicated to breeding disease resilient livestock, but this is not yet a common objective in practical breeding programs. In this paper, we investigate how future breeding programs may benefit from recent research on disease resilience.

MAIN BODY

We define disease resilience in terms of its component traits resistance (R: the ability of a host animal to limit within-host pathogen load (PL)) and tolerance (T: the ability of an infected host to limit the damage caused by a given PL), and model the host's production performance as a reaction norm on PL, depending on R and T. Based on this, we derive equations for the economic values of resilience and its component traits. A case study on porcine respiratory and reproductive syndrome (PRRS) in pigs illustrates that the economic value of increasing production in infectious conditions through selection for R and T can be more than three times higher than by selection for production in disease-free conditions. Although this reaction norm model of resilience is helpful for quantifying its relationship to its component traits, its parameters are difficult and expensive to quantify. We consider the consequences of ignoring R and T in breeding programs that measure resilience as production in infectious conditions with unknown PL-particularly, the risk that the genetic correlation between R and T is unfavourable (antagonistic) and that a trade-off between them neutralizes the resilience improvement. We describe four approaches to avoid such antagonisms: (1) by producing sufficient PL records to estimate this correlation and check for antagonisms-if found, continue routine PL recording, and if not found, shift to cheaper proxies for PL; (2) by selection on quantitative trait loci (QTL) known to influence both R and T in favourable ways; (3) by rapidly modifying towards near-complete resistance or tolerance, (4) by re-defining resilience as the animal's capacity to resist (or recover from) the perturbation caused by an infection, measured as temporal deviations of production traits in within-host longitudinal data series.

CONCLUSIONS

All four alternatives offer promising options for genetic improvement of disease resilience, and most rely on technological and methodological developments and innovation in automated data generation.

Quantifying Individual Response to PRRSV Using Dynamic Indicators of Resilience Based on Activity

Pigs are faced with various perturbations throughout their lives, some of which are induced by management practices, others by natural causes. Resilience is described as the ability to recover from or cope with a perturbation. Using these data, activity patterns of an individual, as well as deviations from these patterns, can potentially be used to quantify resilience. Dynamic indicators of resilience (DIORs) may measure resilience on a different dimension by calculating variation, autocorrelation and skewness of activity from the absolute activity data. The aim of this study was to investigate the potential of using DIORs of activity, such as average, root mean square error (RMSE), autocorrelation or skewness as indicators of resilience to infection with the Porcine Reproductive and Respiratory Syndrome Virus (PRRSV). For this study, individual activity was obtained from 232 pigs equipped with ear tag accelerometers and inoculated with PRRSV between seven and 9 weeks of age. Clinical scores were assigned to each individual at 13 days post-challenge and used to distinguish between a resilient and non-resilient group. Mortality post-challenge was also recorded. Average, RMSE, autocorrelation and skewness of activity were calculated for the pre- and post-challenge phases, as well as the change in activity level pre- vs. post-challenge (i.e., delta). DIORs pre-challenge were expected to predict resilience to PRRSV in the absence of PRRSV infection, whereas DIORs post-challenge and delta were expected to reflect the effect of the PRRSV challenge. None of the pre-challenge DIORs predicted morbidity or mortality post-challenge. However, a higher RMSE in the 3 days post-challenge and larger change in level and RMSE of activity from pre- to post-challenge tended to increase the probability of clinical signs at day 13 post-infection (poor resilience). A higher skewness post-challenge (tendency) and a larger change in skewness from pre- to post-challenge increased the probability of mortality. A decrease in skewness post-challenge lowered the risk of mortality. The post-challenge DIOR autocorrelation was neither linked to morbidity nor to mortality. In conclusion, results from this study showed that post-challenge DIORs of activity can be used to quantify resilience to PRRSV challenge.