Opportunities to Improve Resilience in Animal Breeding Programs

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

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

Unraveling metabolic stress response in dairy cows: genetic control of plasma biomarkers throughout lactation and the transition period

Breeding animals able to effectively respond to stress could be a long-term, sustainable, and affordable strategy to improve resilience and welfare in livestock systems. In the present study, the concentrations of 29 plasma biomarkers were used as candidate endophenotypes for metabolic stress response in single-SNP, gene- and haplotype-based GWAS using 739 healthy lactating Italian Holstein cows and 88,271 variants. Significant genetic associations were found in all the 3 GWAS approaches for plasma γ-glutamyl transferase concentration on BTA17, for paraoxonase on BTA4, and for alkaline phosphatase and zinc on BTA2. On these chromosomes, single-SNP and gene-based chromosome-wide association studies were performed, confirming GWAS findings. The signals identified for paraoxonase, γ-glutamyl transferase, and alkaline phosphatase were in proximity of the genes coding for them. The heritability of these 4 biomarkers ranged from moderate to high (from 0.39 to 0.54). Plasma biomarkers are known to undergo large changes in concentration during metabolic stress in the transition period, with an inter-individual variability in the rate of change and recovery time. Genetics may account in part for these differences. To assess this, we studied a subset of 139 periparturient cows homozygous at 3 SNPs known to be respectively associated with concentration of plasma ceruloplasmin, paraoxonase and γ-glutamyl transferase. We compared the immune-metabolic profile measured in plasma at -7, +5 and +30 d relative to calving between groups of opposite homozygotes. A significant effect of the genotype was found on paraoxonase and γ-glutamyl transferase plasma concentration at all the 3 time points. No evidence for genotype effect was detected for ceruloplasmin. Understanding the genetic control underlying metabolic stress response may suggest new approaches to foster resilience in dairy cows.

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.

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

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

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

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

Integration of selective sweeps across the sheep genome: understanding the relationship between production and adaptation traits

BACKGROUND

Livestock populations are under constant selective pressure for higher productivity levels for different selective purposes. This pressure results in the selection of animals with unique adaptive and production traits. The study of genomic regions associated with these unique characteristics has the potential to improve biological knowledge regarding the adaptive process and how it is connected to production levels and resilience, which is the ability of an animal to adapt to stress or an imbalance in homeostasis. Sheep is a species that has been subjected to several natural and artificial selective pressures during its history, resulting in a highly specialized species for production and adaptation to challenging environments. Here, the data from multiple studies that aim at mapping selective sweeps across the sheep genome associated with production and adaptation traits were integrated to identify confirmed selective sweeps (CSS).

RESULTS

In total, 37 studies were used to identify 518 CSS across the sheep genome, which were classified as production (147 prodCSS) and adaptation (219 adapCSS) CSS based on the frequency of each type of associated study. The genes within the CSS were associated with relevant biological processes for adaptation and production. For example, for adapCSS, the associated genes were related to the control of seasonality, circadian rhythm, and thermoregulation. On the other hand, genes associated with prodCSS were related to the control of feeding behaviour, reproduction, and cellular differentiation. In addition, genes harbouring both prodCSS and adapCSS showed an interesting association with lipid metabolism, suggesting a potential role of this process in the regulation of pleiotropic effects between these classes of traits.

CONCLUSIONS

The findings of this study contribute to a deeper understanding of the genetic link between productivity and adaptability in sheep breeds. This information may provide insights into the genetic mechanisms that underlie undesirable genetic correlations between these two groups of traits and pave the way for a better understanding of resilience as a positive ability to respond to environmental stressors, where the negative effects on production level are minimized.

Resistance to natural tick infestation varies with age and coat and hair traits in a tropically adapted beef cattle breed

Rhipicephalus (Boophilus) microplus causes considerable livestock production losses. Knowledge of the traits that influence tick resistance contributes to the development of breeding strategies designed to improve herd productivity. Within this context, this study evaluated the resistance of Caracu, a tropically adapted cattle breed, to R. microplus. Tick count, hair length, coat thickness, and coat color were evaluated in 202 naturally tick-infested females (cows and heifers) over a period of 18 months. Blood samples were collected from all animals during the winter season for hematological analysis. Data were analyzed using Pearson correlations, generalized linear models, and principal component analysis. Correlation coefficients of tick count with coat color, coat thickness, and hair length were estimated within each season. Hematological parameters were only included in the winter season analysis and were analyzed by the restricted maximum likelihood method using log-transformed data. No differences in blood parameters were observed between animals with and without ticks. However, tick count was negatively correlated with erythrocytes (-0.29) and hematocrit (-0.24) and positively correlated with mean corpuscular hemoglobin (0.21) and mean corpuscular hemoglobin concentration (0.25). These findings suggest that higher tick counts lead to a decrease in erythrocytes but also to an increase in the amount of hemoglobin per erythrocyte, which could reduce the damage caused by low erythrocyte levels due to tick hematophagy, delaying or preventing anemia. Although tick infestation on pasture was demonstrated by the infestation of all staff members during herd management, none of the animals exhibited high tick counts, providing evidence of resistance of Caracu animals to R. microplus. Tick infestation was influenced by age class (cows > heifers), season (spring and summer > fall and winter), coat thickness (>1.5 mm > <1.5 mm), and hair length (>6 mm > <6 mm). Three components were extracted by principal component analysis, which accounted for 69.46% of data variance. The findings of this study will contribute to the development of efficient strategies aimed at reducing economic losses due to tick infestation and could be applied in animal breeding to select for tick resistance traits, reducing chemical control strategies and consequently improving sustainable livestock production.

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.

The genetic architecture of complete blood counts in lactating Holstein dairy cows

Complete blood counts (CBCs) measure the abundance of individual immune cells, red blood cells, and related measures such as platelets in circulating blood. These measures can indicate the health status of an animal; thus, baseline circulating levels in a healthy animal may be related to the productive life, resilience, and production efficiency of cattle. The objective of this study is to determine the heritability of CBC traits and identify genomic regions that are associated with CBC measurements in lactating Holstein dairy cattle. The heritability of CBCs was estimated using a Bayes C0 model. The study population consisted of 388 cows with genotypes at roughly 75,000 markers and 16 different CBC phenotypes taken at one to three time points (n = 33, 131, and 224 for 1, 2, and 3 time points, respectively). Heritabilities ranged from 0.00 ± 0.00 (red cell distribution width) to 0.68 ± 0.06 (lymphocytes). A total of 96 different 1-Mb windows were identified that explained more than 1% of the genetic variance for at least one CBC trait, with 10 windows explaining more than 1% of the genetic variance for two or more traits. Multiple genes in the identified regions have functions related to immune response, cell differentiation, anemia, and disease. Positional candidate genes include RAD52 motif-containing protein 1 (RDM1), which is correlated with the degree of immune infiltration of immune cells, and C-X-C motif chemokine ligand 12 (CXCL12), which is critically involved in neutrophil bone marrow storage and release regulation and enhances neutrophil migration. Since animal health directly impacts feed intake, understanding the genetics of CBCs may be useful in identifying more disease-resilient and feed-efficient dairy cattle. Identification of genes responsible for variation in CBCs will also help identify the variability in how dairy cattle defend against illness and injury.

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.

Landscape genomics reveals regions associated with adaptive phenotypic and genetic variation in Ethiopian indigenous chickens

Climate change is a threat to sustainable livestock production and livelihoods in the tropics. It has adverse impacts on feed and water availability, disease prevalence, production, environmental temperature, and biodiversity. Unravelling the drivers of local adaptation and understanding the underlying genetic variation in random mating indigenous livestock populations informs the design of genetic improvement programmes that aim to increase productivity and resilience. In the present study, we combined environmental, genomic, and phenotypic information of Ethiopian indigenous chickens to investigate their environmental adaptability. Through a hybrid sampling strategy, we captured wide biological and ecological variabilities across the country. Our environmental dataset comprised mean values of 34 climatic, vegetation and soil variables collected over a thirty-year period for 260 geolocations. Our biological dataset included whole genome sequences and quantitative measurements (on eight traits) from 513 individuals, representing 26 chicken populations spread along 4 elevational gradients (6-7 populations per gradient). We performed signatures of selection analyses ([Formula: see text] and XP-EHH) to detect footprints of natural selection, and redundancy analyses (RDA) to determine genotype-environment and genotype-phenotype-associations. RDA identified 1909 outlier SNPs linked with six environmental predictors, which have the highest contributions as ecological drivers of adaptive phenotypic variation. The same method detected 2430 outlier SNPs that are associated with five traits. A large overlap has been observed between signatures of selection identified by[Formula: see text]and XP-EHH showing that both methods target similar selective sweep regions. Average genetic differences measured by [Formula: see text] are low between gradients, but XP-EHH signals are the strongest between agroecologies. Genes in the calcium signalling pathway, those associated with the hypoxia-inducible factor (HIF) transcription factors, and sports performance (GALNTL6) are under selection in high-altitude populations. Our study underscores the relevance of landscape genomics as a powerful interdisciplinary approach to dissect adaptive phenotypic and genetic variation in random mating indigenous livestock populations.

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.

The Resilient Dairy Genome Project-A general overview of methods and objectives related to feed efficiency and methane emissions

The Resilient Dairy Genome Project (RDGP) is an international large-scale applied research project that aims to generate genomic tools to breed more resilient dairy cows. In this context, improving feed efficiency and reducing greenhouse gases from dairy is a high priority. The inclusion of traits related to feed efficiency (e.g., dry matter intake [DMI]) or greenhouse gases (e.g., methane emissions [CH4]) relies on available genotypes as well as high quality phenotypes. Currently, 7 countries (i.e., Australia, Canada, Denmark, Germany, Spain, Switzerland, and United States) contribute with genotypes and phenotypes including DMI and CH4. However, combining data are challenging due to differences in recording protocols, measurement technology, genotyping, and animal management across sources. In this study, we provide an overview of how the RDGP partners address these issues to advance international collaboration to generate genomic tools for resilient dairy. Specifically, we describe the current state of the RDGP database, data collection protocols in each country, and the strategies used for managing the shared data. As of February 2022, the database contains 1,289,593 DMI records from 12,687 cows and 17,403 CH4 records from 3,093 cows and continues to grow as countries upload new data over the coming years. No strong genomic differentiation between the populations was identified in this study, which may be beneficial for eventual across-country genomic predictions. Moreover, our results reinforce the need to account for the heterogeneity in the DMI and CH4 phenotypes in genomic analysis.

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

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

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

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

Genetics and Genomics of Infectious Diseases in Key Aquaculture Species

Diseases pose a significant and pressing concern for the sustainable development of the aquaculture sector, particularly as their impact continues to grow due to climatic shifts such as rising water temperatures. While various approaches, ranging from biosecurity measures to vaccines, have been devised to combat infectious diseases, their efficacy is disease and species specific and contingent upon a multitude of factors. The fields of genetics and genomics offer effective tools to control and prevent disease outbreaks in aquatic animal species. In this study, we present the key findings from our recent research, focusing on the genetic resistance to three specific diseases: White Spot Syndrome Virus (WSSV) in white shrimp, Bacterial Necrotic Pancreatitis (BNP) in striped catfish, and skin fluke (a parasitic ailment) in yellowtail kingfish. Our investigations reveal that all three species possess substantial heritable genetic components for disease-resistant traits, indicating their potential responsiveness to artificial selection in genetic improvement programs tailored to combat these diseases. Also, we observed a high genetic association between disease traits and survival rates. Through selective breeding aimed at enhancing resistance to these pathogens, we achieved substantial genetic gains, averaging 10% per generation. These selection programs also contributed positively to the overall production performance and productivity of these species. Although the effects of selection on immunological traits or immune responses were not significant in white shrimp, they yielded favorable results in striped catfish. Furthermore, our genomic analyses, including shallow genome sequencing of pedigreed populations, enriched our understanding of the genomic architecture underlying disease resistance traits. These traits are primarily governed by a polygenic nature, with numerous genes or genetic variants, each with small effects. Leveraging a range of advanced statistical methods, from mixed models to machine and deep learning, we developed prediction models that demonstrated moderate-to-high levels of accuracy in forecasting these disease-related traits. In addition to genomics, our RNA-seq experiments identified several genes that undergo upregulation in response to infection or viral loads within the populations. Preliminary microbiome data, while offering limited predictive accuracy for disease traits in one of our studied species, underscore the potential for combining such data with genome sequence information to enhance predictive power for disease traits in our populations. Lastly, this paper briefly discusses the roles of precision agriculture systems and AI algorithms and outlines the path for future research to expedite the development of disease-resistant genetic lines tailored to our target species. In conclusion, our study underscores the critical role of genetics and genomics in fortifying the aquaculture sector against the threats posed by diseases, paving the way for more sustainable and resilient aquaculture development.

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

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

Inheritance of feed intake-based resilience traits and their correlation with production traits in Finnish pig breeds

High resilience against diseases, changing environmental conditions, and other stress factors and the ability to efficiently recover to normal status, is becoming increasingly important in pig production. Finding new phenotypes that relate to resilience is a crucial step for improving the resilience of pigs through selection. The objective of this study was to extract resilience-related phenotypes based on fluctuations in daily feed intake (DFI, g) and time spent in feeding per day (TPD, min) and to estimate the heritability of these traits and genetic correlations with production traits (PT). Resilience-related traits with high enough heritability and with either favorable or neutral genetic correlation with PT could be used in the selection program to improve the productivity and welfare of pigs. In this study, we used data from 7,347 Finnish Yorkshire, Landrace, and crossbred pigs raised at the test station. Six pig-specific resilience-related phenotypes were extracted from the individual DFI and TPD: root mean square error (RMSE), quantile regression (QR), and coefficient of variation (CV). RMSE was calculated from the differences between the actual DFI (or TPD) and the pig-specific predicted values. QR was based on the number of days that a pig belonged to the group with the lowest 5% of pigs based on DFI (or TPD), and CV was calculated over the daily observations of DFI (or TPD). PT included average daily gain (ADG, g), backfat thickness (BF, mm), and feed conversion rate (FCR, g/g). The heritability estimates for resilience-related traits varied between 0.07 ± 0.02 (QRDFI) and 0.20 ± 0.03 (RMSETPD). The genetic correlations between resilience-related traits and PT were mostly neutral, but for example, RMSEDFI had a favorable genetic correlation with FCR and BF but an unfavorable correlation with ADG. Lastly, we observed that pigs belonging to the lowest 10% group based on their breeding value (BV) for QRTPD had a lower proportion (10% incidence) of sick days compared to the highest 10% BV group (30% incidence). Therefore, pigs exhibiting small TPD variation (related to high resilience) tend to be less susceptible to sickness than pigs with large TPD variation (related to low resilience). Given its moderate heritability, neutral genetic correlation with PT, and positive effect on health, QRTPD can be considered the most promising resilience-related trait in the Finnish production system.

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.

Genotype-by-environment interactions for feed efficiency traits in Nellore cattle based on bi-trait reaction norm models

BACKGROUND

Selecting animals for feed efficiency directly impacts the profitability of the beef cattle industry, which contributes to minimizing the environmental footprint of beef production. Genetic and environmental factors influence animal feed efficiency, leading to phenotypic variability when exposed to different environmental conditions (i.e., temperature and nutritional level). Thus, our aim was to assess potential genotype-by-environment (G × E) interactions for dry matter intake (DMI) and residual feed intake (RFI) in Nellore cattle (Bos taurus indicus) based on bi-trait reaction norm models (RN) and evaluate the genetic association between RFI and DMI across different environmental gradient (EG) levels. For this, we used phenotypic information on 12,958 animals (young bulls and heifers) for DMI and RFI recorded during 158 feed efficiency trials.

RESULTS

The heritability estimates for DMI and RFI across EG ranged from 0.26 to 0.54 and from 0.07 to 0.41, respectively. The average genetic correlations (± standard deviation) across EG for DMI and RFI were 0.83 ± 0.19 and 0.81 ± 0.21, respectively, with the lowest genetic correlation estimates observed between extreme EG levels (low vs. high) i.e. 0.22 for RFI and 0.26 for DMI, indicating the presence of G × E interactions. The genetic correlation between RFI and DMI across EG levels decreased as the EG became more favorable and ranged from 0.79 (lowest EG) to 0.52 (highest EG). Based on the estimated breeding values from extreme EG levels (low vs. high), we observed a moderate Spearman correlation of 0.61 (RFI) and 0.55 (DMI) and a selection coincidence of 53.3% and 40.0% for RFI and DMI, respectively.

CONCLUSIONS

Our results show evidence of G × E interactions on feed efficiency traits in Nellore cattle, especially in feeding trials with an average daily gain (ADG) that is far from the expected of 1 kg/day, thus increasing reranking of animals.

Identification of Genomic Instability in Cows Infected with BVD Virus

An important factor for dairy cattle farmers is the profitability of cattle rearing, which is influenced by the animals' health and reproductive parameters, as well as their genomic stability and integrity. Bovine viral diarrhea (BVD) negatively affects the health of dairy cattle and causes reproductive problems. The aim of the study was to identify genomic instability in cows with reproductive disorders following infection with the BVD virus. The material for analysis was peripheral blood from Holstein-Friesian cows with reproductive problems, which had tested positive for BVD, and from healthy cows with no reproductive problems, which had tested negative for BVD. Three cytogenetic tests were used: the sister chromatid exchange assay, fragile sites assay, and comet assay. Statistically significant differences were noted between the groups and between the individual cows in the average frequency of damage. The assays were good biomarkers of genomic stability and enabled the identification of individuals with an increased frequency of damage to genetic material that posed a negative impact on their health. The assays can be used to prevent disease during its course and evaluate the genetic resistance of animals. This is especially important for the breeder, both for economic and breeding reasons. Of the three assays, the comet assay proved to be the most sensitive for identifying DNA damage in the animals.

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

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

Unravelling the genetic variability of host resilience to endo- and ectoparasites in Nellore commercial herds

BACKGROUND

Host resilience (HR) to parasites can affect the performance of animals. Therefore, the aim of this study was to present a detailed investigation of the genetic mechanisms of HR to ticks (TICK), gastrointestinal nematodes (GIN), and Eimeria spp. (EIM) in Nellore cattle that were raised under natural infestation and a prophylactic parasite control strategy. In our study, HR was defined as the slope coefficient of body weight (BW) when TICK, GIN, and EIM burdens were used as environmental gradients in random regression models. In total, 1712 animals were evaluated at five measurement events (ME) at an average age of 331, 385, 443, 498, and 555 days, which generated 7307 body weight (BW) records. Of the 1712 animals, 1075 genotyped animals were used in genome-wide association studies to identify genomic regions associated with HR.

RESULTS

Posterior means of the heritability estimates for BW ranged from 0.09 to 0.54 across parasites and ME. The single nucleotide polymorphism (SNP)-derived heritability for BW at each ME ranged from a low (0.09 at ME.331) to a moderate value (0.23 at ME.555). Those estimates show that genetic progress can be achieved for BW through selection. Both genetic and genomic associations between BW and HR to TICK, GIN, and EIM confirmed that parasite infestation impacted the performance of animals. Selection for BW under an environment with a controlled parasite burden is an alternative to improve both, BW and HR. There was no impact of age of measurement on the estimates of genetic variance for HR. Five quantitative trait loci (QTL) were associated with HR to EIM but none with HR to TICK and to GIN. These QTL contain genes that were previously shown to be associated with the production of antibody modulators and chemokines that are released in the intestinal epithelium.

CONCLUSIONS

Selection for BW under natural infestation and controlled parasite burden, via prophylactic parasite control, contributes to the identification of animals that are resilient to nematodes and Eimeria ssp. Although we verified that sufficient genetic variation existed for HR, we did not find any genes associated with mechanisms that could justify the expression of HR to TICK and GIN.

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.

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

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

Behavioral patterns as indicators of resilience after parturition in dairy cows

During the transition phase, dairy cows are susceptible to develop postpartum diseases. Cows that stay healthy or recover rapidly can be considered to be more resilient in comparison to those that develop postpartum diseases. An indication of loss of resilience will allow for early intervention with preventive and supportive measures before the onset of disease. We investigated which quantitative behavioral characteristics during the dry period could be used as indicators of reduced resilience after calving, using noninvasive Smart Tag neck and Smart Tag leg sensors in dairy cows (Nedap N.V.). We followed 180 cows during 2 wk before until 6 wk after parturition at 4 farms in the Netherlands. Serving as proxy for loss of resilience, as defined by the duration and severity of disease, a clinical assessment was performed twice weekly and blood samples were taken in the first and fifth week after parturition. For each cow, clinical and serum value deviations were aggregated into a total deficit score (TDS total). We also calculated TDS values relating to inflammation, locomotion, or metabolic problems, which were further divided into macro-mineral and liver-related deviations. Smart Tag neck and leg sensors provided continuous behavioral activity signals of which we calculated the average, variance, and autocorrelation during the dry period. Diurnal patterns in the behavioral activity signals were derived by fast Fourier transformation and the calculation of the nonperiodicity. To select significant predictors of resilience, we first performed a univariate analysis with TDS as dependent variable and the behavioral characteristics that were measured during the dry period, as potential predictors with cow as experimental unit. We included parity group as fixed effect and farm as random effect. Next, we performed multivariable analysis with only significant predictors, followed by a variable selection procedure to obtain a final linear mixed model with an optimal subset of predictors with parity group as fixed effect and farm as random effect. The TDS total was best predicted by average inactive time, nonperiodicity ruminating, nonperiodicity of bouts standing up and fast Fourier transformation stand still. Average inactive time was negatively correlated with average eating time, and these 2 predictors could be exchanged with only little difference in model performance. Our best performing model predicted TDS total at a cutoff level of 60 points, with a sensitivity of 79.5% and a specificity of 73.2% with a positive predicted value of 0.69 and a negative predicted value of 0.83. The models to predict the other TDS categories showed a lower predictive performance as compared with the TDS total model, which could be related to the limited sample size and therefore, low occurrence of problems within a specific TDS category. Furthermore, more resilient dairy cows are characterized by high averages of eating time with high regularity in rumination and low averages of inactive time. They reveal high regularity in standing time and transitions from lying to standing, in the dry period. These behaviors can be used as indicators of resilience and allow for preventive intervention during the dry period in vulnerable dairy cattle. However, further examination is still required to find clues for adequate intervention strategies.

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.

Selection for environmental variance shifted the gut microbiome composition driving animal resilience

BACKGROUND

Understanding how the host's microbiome shapes phenotypes and participates in the host response to selection is fundamental for evolutionists and animal and plant breeders. Currently, selection for resilience is considered a critical step in improving the sustainability of livestock systems. Environmental variance (V _E), the within-individual variance of a trait, has been successfully used as a proxy for animal resilience. Selection for reduced V _E could effectively shift gut microbiome composition; reshape the inflammatory response, triglyceride, and cholesterol levels; and drive animal resilience. This study aimed to determine the gut microbiome composition underlying the V _E of litter size (LS), for which we performed a metagenomic analysis in two rabbit populations divergently selected for low (n = 36) and high (n = 34) V _E of LS. Partial least square-discriminant analysis and alpha- and beta-diversity were computed to determine the differences in gut microbiome composition among the rabbit populations.

RESULTS

We identified 116 KEGG IDs, 164 COG IDs, and 32 species with differences in abundance between the two rabbit populations studied. These variables achieved a classification performance of the V _E rabbit populations of over than 80%. Compared to the high V _E population, the low V _E (resilient) population was characterized by an underrepresentation of Megasphaera sp., Acetatifactor muris, Bacteroidetes rodentium, Ruminococcus bromii, Bacteroidetes togonis, and Eggerthella sp. and greater abundances of Alistipes shahii, Alistipes putredinis, Odoribacter splanchnicus, Limosilactobacillus fermentum, and Sutterella, among others. Differences in abundance were also found in pathways related to biofilm formation, quorum sensing, glutamate, and amino acid aromatic metabolism. All these results suggest differences in gut immunity modulation, closely related to resilience.

CONCLUSIONS

This is the first study to show that selection for V _E of LS can shift the gut microbiome composition. The results revealed differences in microbiome composition related to gut immunity modulation, which could contribute to the differences in resilience among rabbit populations. The selection-driven shifts in gut microbiome composition should make a substantial contribution to the remarkable genetic response observed in the V _E rabbit populations. Video Abstract.

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

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

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

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

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

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

Going to scale—From community-based to population-wide genetic improvement and commercialized sheep meat supply in Ethiopia

Community-based breeding programs (CBBPs) have shown, at pilot scale, to be effective and beneficial in achieving genetic progress and in improving livelihoods of smallholder communities. In Ethiopia 134 sheep and goat CBBPs were operational producing their own improved rams and bucks. Based on experience the implementation of further programs is possible with appropriate private and public support. A different challenge is the efficient dissemination of the improved genetics produced in current CBBPs to create population-wide economic impact. We present a framework applied to the Ethiopian Washera sheep breed to meet this challenge. We propose the establishment of a genetic improvement structure that supports a meat commercialization model based on the integration of community-based breeding program cooperatives, client communities and complementary services such as fattening enterprises. We calculated that the recently established 28 community-based breeding programs in the Washera breeding tract can provide genetically improved rams to 22% of the four million head. To reach the whole population 152 additional CBBPs are needed. We simulated the genetic improvements obtainable in the current 28 CBBPs assuming realized genetic progress in CBBPs of a similar breed and calculated the expected additional lamb carcass meat production after 10 years of selection to be 7 tons and the accumulated discounted benefit 327 thousand USD. These benefits could be increased if the CBBPs are linked to client communities by providing them with improved rams: additional meat production would be 138 tons with a value of 3,088 thousand USD. The total meat production of the existing Washera CBBPs was calculated at 152 tons and the joint meat production of CBBPs if integrated with client communities would be 3,495 tons. A full integration model, which includes enterprises purchasing lambs for fattening, can produce up to 4,255 tons of meat. We conclude that Washera CBBPs cooperatives can benefit from a higher level of organization to produce population-wide genetic improvement and economic benefits. Unlike in the dairy and chicken industries, for low input sheep and goat smallholder systems the proposed commercialization model puts breeder cooperatives at the center of the operation. Cooperatives need to be capacitated and supported to become fully functional business ventures.

Effect of environmental variance-based resilience selection on the gut metabolome of rabbits

BACKGROUND

Gut metabolites are key actors in host-microbiota crosstalk with effect on health. The study of the gut metabolome is an emerging topic in livestock, which can help understand its effect on key traits such as animal resilience and welfare. Animal resilience has now become a major trait of interest because of the high demand for more sustainable production. Composition of the gut microbiome can reveal mechanisms that underlie animal resilience because of its influence on host immunity. Environmental variance (V_E), specifically the residual variance, is one measure of resilience. The aim of this study was to identify gut metabolites that underlie differences in the resilience potential of animals originating from a divergent selection for V_E of litter size (LS). We performed an untargeted gut metabolome analysis in two divergent rabbit populations for low (n = 13) and high (n = 13) V_E of LS. Partial least square-discriminant analysis was undertaken, and Bayesian statistics were computed to determine dissimilarities in the gut metabolites between these two rabbit populations.

RESULTS

We identified 15 metabolites that discriminate rabbits from the divergent populations with a prediction performance of 99.2% and 90.4% for the resilient and non-resilient populations, respectively. These metabolites were suggested to be biomarkers of animal resilience as they were the most reliable. Among these, five that derived from the microbiota metabolism (3-(4-hydroxyphenyl)lactate, 5-aminovalerate, and equol, N6-acetyllysine, and serine), were suggested to be indicators of dissimilarities in the microbiome composition between the rabbit populations. The abundances of acylcarnitines and metabolites derived from the phenylalanine, tyrosine, and tryptophan metabolism were low in the resilient population and these pathways can, therefore impact the inflammatory response and health status of animals.

CONCLUSIONS

This is the first study to identify gut metabolites that could act as potential resilience biomarkers. The results support differences in resilience between the two studied rabbit populations that were generated by selection for V_E of LS. Furthermore, selection for V_E of LS modified the gut metabolome, which could be another factor that modulates animal resilience. Further studies are needed to determine the causal role of these metabolites in health and disease.

Genotype by environment interaction due to heat stress in Brown Swiss cattle

Due to its geographical position and a highly variable orography, Italy is characterized by several climatic areas and thus, by many different dairy cow farming systems. Brown Swiss cattle, in this context, are a very appreciated genetic resource for their adaptability and low metabolic requirement. The significant heterogeneity in farming systems may consist of genotype by environment (G × E) interactions with neglected changes in animals' rank position. The objective of this study was to investigate G × E for heat tolerance in Brown Swiss cattle for several production traits (milk, fat, and protein yield in kilograms; fat, protein, and cheese yield in percentage) and 2 derivate traits (fat-corrected milk and energy-corrected milk). We used the daily maximum temperature-humidity index (THI) range, calculated according to weather stations' data from 2008 to 2018 in Italy, and 202,776 test-day records from 23,396 Brown Swiss cows from 639 herds. Two different methodologies were applied to estimate the effect of the environmental variable (THI) on genetic parameters: (1) the reaction norm model, which uses a continuous random covariate to estimate the animal additive effect, and (2) the multitrait model, which splits each production pattern as a distinct and correlated trait according to the first (a thermal comfort condition), third (a moderate heat stress condition), and fifth (a severe heat stress condition) mean THI value quintile. The results from the reaction norm model showed a descending trend of the additive genetic effect until THI reached the value of 80. Then we recorded an increase with high extreme THI values (THI 90). Permanent environmental variance at increasing THI values revealed an opposite trend: The plot of heritability and the ratio of animal permanent environmental variance to phenotypic variance showed that when the environmental condition worsens, the additive genetic and permanent environmental component for production traits play a growing role. The negative additive genetic correlation between slope and linear random coefficient indicates no linear relationship between the production traits or under heat stress conditions, except for milk yield and protein yield. In tridimensional wireframe plots, the extreme margin decreases until a minimum of ∼0.90 of genetic correlation in the ECM trait, showing that the magnitude of G × E interaction is greater than the other traits. Genetic correlation values in Brown Swiss suggest the possibility of moderate changes in animals' estimated breeding value in heat stress conditions. Results indicated a moderate G × E interaction but significant variability in sire response related to their production level.

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

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

Candidate genes for mastitis resistance in dairy cattle: a data integration approach

BACKGROUND

Inflammation of the mammary tissue (mastitis) is one of the most detrimental health conditions in dairy ruminants and is considered the most economically important infectious disease of the dairy sector. Improving mastitis resistance is becoming an important goal in dairy ruminant breeding programmes. However, mastitis resistance is a complex trait and identification of mastitis-associated alleles in livestock is difficult. Currently, the only applicable approach to identify candidate loci for complex traits in large farm animals is to combine different information that supports the functionality of the identified genomic regions with respect to a complex trait.

METHODS

To identify the most promising candidate loci for mastitis resistance we integrated heterogeneous data from multiple sources and compiled the information into a comprehensive database of mastitis-associated candidate loci. Mastitis-associated candidate genes reported in association, expression, and mouse model studies were collected by searching the relevant literature and databases. The collected data were integrated into a single database, screened for overlaps, and used for gene set enrichment analysis.

RESULTS

The database contains candidate genes from association and expression studies and relevant transgenic mouse models. The 2448 collected candidate loci are evenly distributed across bovine chromosomes. Data integration and analysis revealed overlaps between different studies and/or with mastitis-associated QTL, revealing promising candidate genes for mastitis resistance.

CONCLUSION

Mastitis resistance is a complex trait influenced by numerous alleles. Based on the number of independent studies, we were able to prioritise candidate genes and propose a list of the 22 most promising. To our knowledge this is the most comprehensive database of mastitis associated candidate genes and could be helpful in selecting genes for functional validation studies.

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

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

Exploration of plasma metabolite levels in healthy nursery pigs in response to environmental enrichment and disease resilience

The purpose of this study was to explore plasma metabolite levels in young healthy pigs and their potential association with disease resilience and estimate genetic and phenotypic correlation with the change in lymphocyte concentration following disease challenge. Plasma samples were collected from 968 healthy nursery pigs over 15 batches at an average of 28 ± 3.23 d of age. Forty-four metabolites were identified and quantified by nuclear magnetic resonance. Pigs were then introduced into a natural disease challenge barn, and were classified into four groups based on the growth rate of each animal in the grow-to-finish phase (GFGR) and treatment rate (TR): resilient (RES), average (MID), susceptible (SUS), and dead (pigs that died before harvest). Blood samples were collected from all pigs before and 2 wk after disease challenge and complete blood count was determined. Environmental enrichment (inedible point source objects) was provided for half of the pigs in seven batches (N = 205) to evaluate its impact on resilience and metabolite concentrations. Concentration of all metabolites was affected by batch, while entry age affected the concentration of 16 metabolites. The concentration of creatinine was significantly lower for pigs classified as "dead" and "susceptible" when compared to "average" (P < 0.05). Pigs that received enrichment had significantly lower concentrations of six metabolites compared with pigs that did not receive enrichment (P ≤ 0.05). Both, group classification and enrichment affected metabolites that are involved in the same pathways of valine, leucine, and isoleucine biosynthesis and degradation. Resilient pigs had higher increase in lymphocyte concentration after disease challenge. The concentration of plasma l-α-aminobutyric acid was significantly negatively genetically correlated with the change in lymphocyte concentration following challenge. In conclusion, creatinine concentration in healthy nursery pigs was lower in pigs classified as susceptible or dead after disease challenge, whilst l-α-aminobutyric may be a genetic biomarker of lymphocyte response after pathogen exposure, and both deserve further investigation. Batch, entry age, and environmental enrichment were important factors affecting the concentration of metabolites and should be taken into consideration in future studies.

Exploiting morphobiometric and genomic variability of African indigenous camel populations-A review

Camels (Camelus dromedarius) in Africa are adapted to arid and the semi-arid environmental conditions, and are valuable for meat, milk and fiber production. On account of the growing demand for camels in this continent, there is a need for knowledge on their phenotypic and genetic diversity. This is fundamental to sustainable herd management and utilization including the design of appropriate breeding and conservation strategies. We reviewed studies on the phenotypic and genetic characterization, breeding objectives, systems of production, productive and reproductive performances, and pathways for the sustainable rearing and use of camels in Africa. The morphological and genetic diversity, productive and reproductive abilities of African camels suggest the existence of genetic variations that can be utilized for breeds/ecotypes' genetic improvement and conservation. Possible areas of intervention include the establishment of open nucleus and community-based breeding schemes and utilization of modern reproductive technologies for the genetic improvement of milk and meat yields, sustainable management of rangelands, capacity building of the pastoralists and agro-pastoralists, institutional supports, formation of centralized conservation centres and efficient and effective marketing systems.

Exploring milk loss and variability during environmental perturbations across lactation stages as resilience indicators in Holstein cattle

Genetic selection for resilience is essential to improve the long-term sustainability of the dairy cattle industry, especially the ability of cows to maintain their level of production when exposed to environmental disturbances. Recording of daily milk yield provides an opportunity to develop resilience indicators based on milk losses and fluctuations in daily milk yield caused by environmental disturbances. In this context, our study aimed to explore milk loss traits and measures of variability in daily milk yield, including log-transformed standard deviation of milk deviations (Lnsd), lag-1 autocorrelation (Ra), and skewness of the deviations (Ske), as indicators of general resilience in dairy cows. The unperturbed dynamics of milk yield as well as milk loss were predicted using an iterative procedure of lactation curve modeling. Milk fluctuations were defined as a period of at least 10 successive days of negative deviations in which milk yield dropped at least once below 90% of the expected values. Genetic parameters of these indicators and their genetic correlation with economically important traits were estimated using single-trait and bivariate animal models and 8,935 lactations (after quality control) from 6,816 Chinese Holstein cows. In general, cows experienced an average of 3.73 environmental disturbances with a milk loss of 267 kg of milk per lactation. Each fluctuation lasted for 19.80 ± 11.46 days. Milk loss traits are heritable with heritability estimates ranging from 0.004 to 0.061. The heritabilities differed between Lnsd (0.135-0.250), Ra (0.008-0.058), and Ske (0.001-0.075), with the highest heritability estimate of 0.250 ± 0.020 for Lnsd when removing the first and last 10 days in milk in a lactation (Lnsd2). Based on moderate to high genetic correlations, lower Lnsd2 is associated with less milk losses, better reproductive performance, and lower disease incidence. These findings indicate that among the variables evaluated, Lnsd2 is the most promising indicator for breeding for improved resilience in Holstein cattle.

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.

Variations in Growth and Photosynthetic Traits of Polyploid Poplar Hybrids and Clones in Northeast China

To evaluate differences among 19 different ploidy hybrid poplar clones grown in northeast China, 21 traits related to growth traits and photosynthetic characteristics were detected and analyzed. Abundant phenotypic variations exist among and within populations, and these variations are the basis of forest tree genetic improvements. In this research, variance analysis showed that the traits except the net photosynthesis rate among the different ploidies and all the other traits exhibited significant differences among the ploidies or clones (p < 0.01). Estimation of phenotypic coefficients of variation, genotypic coefficients of variation, and repeatability is important for selecting superior materials. The larger the value, the greater the potential for material selection improvement. The repeatability of the different traits ranged from 0.88 to 0.99. The phenotypic and genotypic coefficients of variation of all the investigated traits ranged from 6.88% to 57.40% and from 4.85% to 42.89%, respectively. Correlation analysis showed that there were significant positive correlations between tree height, diameter, and volume. Transpiration rate, intercellular carbon dioxide concentration, and stomatal conductance were significantly positively correlated with each other but negatively correlated with instantaneous water use efficiency. Growth traits were weakly correlated with photosynthetic indexes. The rank correlation coefficient showed that most of the growth indicators reached a significant correlation level among different years (0.40-0.98), except 1-year-old tree height with 4-year-old tree height and 1-year-old ground diameter with 3-year-old tree height, which indicated the potential possibility for early selection of elite clones. Principal analysis results showed that the contribution rate of the first principal component was 46.606%, and 2-year-old tree height, 2-year-old ground diameter, 3-year-old tree height, 3-year-old ground diameter, 3-year-old diameter at breast height, 3-year-old volume, 4-year-old tree height, 4-year-old ground diameter, 4-year-old diameter at breast height, and 4-year-old volume showed higher vector values than other traits. With the method of multiple-trait comprehensive evaluation to evaluate clones, SX3.1, SY3.1, and XY4.2 were selected as elite clones, and the genetic gains of height, basal diameter, diameter at breast height, and volume of selected clones ranged from 12.85% to 64.87% in the fourth growth year. The results showed fundamental information for selecting superior poplar clones, which might provide new materials for the regeneration and improvement of forests in Northeast China.

Positive Welfare Indicators in Dairy Animals

Nowadays, there is growing interest in positive animal welfare not only from the view of scientists but also from that of society. The consumer demands more sustainable livestock production, and animal welfare is an essential part of sustainability, so there is interest in incorporating positive welfare indicators into welfare assessment schemes and legislation. The aim of this review is to cite all the positive welfare indicators that have been proposed for dairy animals in theory or practice. In total, twenty-four indicators were retrieved. The most promising are exploration, access to pasture, comfort and resting, feeding, and behavioral synchronicity. Qualitative behavioral assessment (QBA), social affiliative behaviors, play, maternal care, ear postures, vocalizations, visible eye white, nasal temperature, anticipation, cognitive bias, laterality, and oxytocin have been also studied in dairy ruminants. QBA is the indicator that is most often used for the on-farm welfare assessment. Among all dairy animals, studies have been performed mostly on cattle, followed by sheep and goats, and finally buffaloes. The research on camel welfare is limited. Therefore, there is a need for further research and official assessment protocols for buffaloes and especially camels.

Can breeders prevent pathogen adaptation when selecting for increased resistance to infectious diseases?

Background

Recent research shows that genetic selection has high potential to reduce the prevalence of infectious diseases in livestock. However, like all interventions that target infectious diseases, genetic selection of livestock can exert selection pressure on pathogen populations. Such selection on the pathogen may lead to escape strategies and reduce the effect of selection of livestock for disease resistance. Thus, to successfully breed livestock for lower disease prevalence, it is essential to develop strategies that prevent the invasion of pathogen mutants that escape host resistance. Here we investigate the conditions under which such “escape mutants” can replace wild-type pathogens in a closed livestock population using a mathematical model of disease transmission.

Results

Assuming a single gene that confers sufficient resistance, results show that genetic selection for resistance in livestock typically leads to an “invasion window” within which an escape mutant of the pathogen can invade. The bounds of the invasion window are determined by the frequency of resistant hosts in the population. The lower bound occurs when the escape mutant has an advantage over the wild-type pathogen in the population. The upper bound occurs when local eradication of the pathogen is expected. The invasion window is smallest when host resistance is strong and when infection with the wild-type pathogen provides cross immunity to infection with the escape mutant.

Conclusions

To minimise opportunities for pathogens to adapt, under the assumptions of our model, the aim of disease control through genetic selection should be to achieve herd-level eradication of the infection faster than the rate of emergence of escape mutants of the pathogen. Especially for microparasitic infections, this could be achieved by placing animals into herds according to their genetic resistance, such that these herds stay completely out of the invasion window. In contrast to classical breeding theory, our model suggests that multi-trait selection with gradual improvement of each trait of the breeding goal might not be the best strategy when resistance to infectious disease is part of the breeding goal. Temporally, combining genetic selection with other interventions helps to make the invasion window smaller, and thereby reduces the risk of invasion of escape mutants.

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

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

Relationship among Milk Conductivity, Production Traits, and Somatic Cell Score in the Italian Mediterranean Buffalo

The measurement of milk electrical conductivity (EC) is a relatively simple and inexpensive technique that has been evaluated as a routine method for the diagnosis of mastitis in dairy farms. The aim of this study was to obtain further knowledge on relationships between EC, production traits and somatic cell count (SCC) in Italian Mediterranean Buffalo. The original dataset included 5411 records collected from 808 buffalo cows. Two mixed models were used to evaluate both the effect of EC on MY, PP and FP and EC at test-day, and the effect of EC on somatic cell score (SCS) by using five different parameters (EC_param), namely: EC collected at the official milk recording test day (EC_day0), EC collected 3 days before official milk recording (EC_day3), and three statistics calculated from EC collected 1, 3 and 5 days before each test-day, respectively. All effects included in the model were significant for all traits, with the only exception of the effect of EC nested within parity for FP. The relationship between EC and SCS was always positive, but of different magnitude according to the parity. The regression of EC on SCS at test-day using different EC parameters was always significant except when the regression parameter was the slope obtained from a linear regression of EC collected over the 5-day period. Moreover, in order to evaluate how well the different models fit the data, three parameters were used: the Average Information Criteria (AIC), the marginal R2 and the conditional R2. According to AIC and to both the Marginal and Conditional R2, the best results were obtained when the regression parameter was the mean EC estimated over the 5-day period.