Selection for environmental variance of litter size in rabbits

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.

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.

Differential patterns in runs of homozygosity in two mice lines under divergent selection for environmental variability for birth weight

Runs of homozygosity (ROH) are defined as long continuous homozygous stretches in the genome which are assumed to originate from a common ancestor. It has been demonstrated that divergent selection for variability in mice is possible and that low variability in birth weight is associated with robustness. To analyse ROH patterns and ROH-based genomic inbreeding, two mouse lines that were divergently selected for birth weight variability for 26 generations were used, with: 752 individuals for the high variability line (H-Line), 766 individuals for the low variability line (L-Line) and 74 individuals as a reference population. Individuals were genotyped using the high density Affymetrix Mouse Diversity Genotyping Array. ROH were identified using both the sliding windows (SW) and the consecutive runs (CR) methods. Inbreeding coefficients were calculated based on pedigree (FPED ) information, on ROH identified using the SW method (FROHSW ) and on ROH identified using the CR method (FROHCR ). Differences in genomic inbreeding were not consistent across generations and these parameters did not show clear differences between lines. Correlations between FPED and FROH were high, particularly for FROHSW . Moreover, correlations between FROHSW and FPED were even higher when ROH were identified with no restrictions in the number of heterozygotes per ROH. The comparison of FROH estimates between either of the selected lines were based on significant differences at the chromosome level, mainly in chromosomes 3, 4, 6, 8, 11, 15 and 19. ROH-based inbreeding estimates that were computed using longer homozygous segments had a higher relationship with FPED . Differences in robustness between lines were not attributable to a higher homozygosis in the L-Line, but maybe to the different distribution of ROH at the chromosome level between lines. The analysis identified a set of genomic regions for future research to establish the genomic basis of robustness.

Macro- and micro-genetic environmental sensitivity of yearling weight in Angus beef cattle

Australian beef cattle experience variable conditions, which may give rise to genotype-by-environment interactions depending on the genotypes' macro- and/or micro-genetic environmental sensitivity (GES). Macro-GES gives rise to genotype-by-environment interactions across definable and shared environments, while micro-GES causes heritable variation of phenotypes, e.g., the performance of progeny from one sire may be more variable than other sires. Yearling weight (YW) is a key trait in Australian Angus cattle that may be impacted by both macro- and micro-GES. Current models for genetic evaluation of YW attempt to account for macro-GES by fitting sire-by-herd interactions (S × H). Variation in micro-GES had not yet been estimated for YW in Australian Angus. The aim of this study was to estimate genetic variation due to macro- and micro-GES in YW of Australian Angus cattle. A reaction norm with contemporary group effects as the environmental covariate was fitted either as an alternative to or in combination with a random S × H effect to account for macro-GES. Double hierarchical generalised linear models (DHGLM), fitted as sire models, were used to estimate the genetic variance of the dispersion as a measure of micro-GES. Variation due to both macro- and micro-GES were found in YW. The variance of the slope of the reaction norm was 0.02-0.03 (SEs 0.00), while the S × H variance accounted for 7% of the phenotypic variance in all models. Results showed that both a random S × H effect and a reaction norm should be included to account for both macro-GES and the additional variation captured by an S × H effect. The heritability of the dispersion on the measurement scale ranged from 0.06 to 0.10 (SEs 0.00) depending on which model was used. It should therefore be possible to alter both macro- and micro-GES of YW in Australian Angus through selection. However, care should be taken to ensure an appropriate data structure when including sire-by-herd interactions in the mean part of a DHGLM; otherwise, it can cause biased estimates of micro-GES.

LTP2 hypomorphs show genotype-by-environment interaction in early seedling traits in Arabidopsis thaliana

Isogenic individuals can display seemingly stochastic phenotypic differences, limiting the accuracy of genotype-to-phenotype predictions. The extent of this phenotypic variation depends in part on genetic background, raising questions about the genes involved in controlling stochastic phenotypic variation. Focusing on early seedling traits in Arabidopsis thaliana, we found that hypomorphs of the cuticle-related gene LIPID TRANSFER PROTEIN 2 (LTP2) greatly increased variation in seedling phenotypes, including hypocotyl length, gravitropism and cuticle permeability. Many ltp2 hypocotyls were significantly shorter than wild-type hypocotyls while others resembled the wild-type. Differences in epidermal properties and gene expression between ltp2 seedlings with long and short hypocotyls suggest a loss of cuticle integrity as the primary determinant of the observed phenotypic variation. We identified environmental conditions that reveal or mask the increased variation in ltp2 hypomorphs and found that increased expression of its closest paralog LTP1 is necessary for ltp2 phenotypes. Our results illustrate how decreased expression of a single gene can generate starkly increased phenotypic variation in isogenic individuals in response to an environmental challenge.

Micro-genetic environmental sensitivity across macro-environments of chickens reared in Burkina Faso and France

BACKGROUND

Commercial poultry production systems follow a pyramidal structure with a nucleus of purebred animals under controlled conditions at the top and crossbred animals under commercial production conditions at the bottom. Genetic correlations between the same phenotypes on nucleus and production animals can therefore be influenced by differences both in purebred-crossbred genotypes and in genotype-by-environment interactions across the two environments, known as macro-genetic environmental sensitivity (GES). Within each environment, genotype-by-environment interactions can also occur due to so-called micro-GES. Micro-GES causes heritable variation in phenotypes and decreases uniformity. In this study, genetic variances of body weight (BW) and of micro-GES of BW and the impacts of purebred-crossbred differences and macro-environmental differences on micro-GES of BW were estimated. The dataset contained three subpopulations of slow-growing broiler chickens: purebred chickens (PB) reared in France, and crossbred chickens reared in France (FR) under the same conditions as PB or reared in Burkina Faso (BF) under local conditions. The crossbred chickens were offspring of the same dam line and had PB as their sire line.

RESULTS

Estimates of heritability of BW and micro-GES of BW were 0.54 (SE of 0.02) and 0.06 (0.01), 0.67 (0.03) and 0.03 (0.01), and 0.68 (0.04) and 0.02 (0.01) for the BF, FR, and PB subpopulations, respectively. Estimates of the genetic correlations for BW between the three subpopulations were moderately positive (0.37 to 0.53) and those for micro-GES were weakly to moderately positive (0.01 to 0.44).

CONCLUSIONS

The results show that the heritability of the micro-GES of BW varies with macro-environment, which indicates that responses to selection are expected to differ between macro-environments. The weak to moderate positive genetic correlations between subpopulations indicate that both macro-environmental differences and purebred-crossbred differences can cause re-ranking of sires based on their estimated breeding values for micro-GES of BW. Thus, the sire that produces the most variable progeny in one macro-environment may not be the one that produces the most variable offspring in another. Similarly, the sire that produces the most variable purebred progeny may not produce the most variable crossbred progeny. The results highlight the need for investigating micro-GES for all subpopulations included in the selection scheme, to ensure optimal genetic gain in all subpopulations.

LTP2 hypomorphs show genotype-by-environment interaction in early seedling traits in Arabidopsis thaliana

Isogenic individuals can display seemingly stochastic phenotypic differences, limiting the accuracy of genotype-to-phenotype predictions. The extent of this phenotypic variation depends in part on genetic background, raising questions about the genes involved in controlling stochastic phenotypic variation. Focusing on early seedling traits in Arabidopsis thaliana, we found that hypomorphs of the cuticle-related gene LTP2 greatly increased variation in seedling phenotypes, including hypocotyl length, gravitropism and cuticle permeability. Many ltp2 hypocotyls were significantly shorter than wild-type hypocotyls while others resembled the wild type. Differences in epidermal properties and gene expression between ltp2 seedlings with long and short hypocotyls suggest a loss of cuticle integrity as the primary determinant of the observed phenotypic variation. We identified environmental conditions that reveal or mask the increased variation in ltp2 hypomorphs, and found that increased expression of its closest paralog LTP1 is necessary for ltp2 phenotypes. Our results illustrate how decreased expression of a single gene can generate starkly increased phenotypic variation in isogenic individuals in response to an environmental challenge.

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.

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.

Genetic parameters for fleece uniformity in alpacas

Fiber diameter is the main selection objective and criterion in alpaca breeding programs, but it can vary across anatomic regions of the animal. As fiber diameter is usually registered from a unique sample from the mid side of the body, fiber diameter variability within fleece is never addressed and phenotypic and genetic differences may exist for fleece uniformity in alpaca populations. The objective of this work was to estimate the genetic parameters of fleece uniformity in an alpaca population. Fiber diameters measured in three different locations were used as repeated records of the same animal and studied for fitting a model that considers heterogeneous the residual variance of the model. Also, the logarithm of the standard deviation of the three measures was used as a measure of the fleece variability. Estimate of the additive genetic variance of the environmental variability was 0.43±0.14, enough high to suggest the existence of wide room to select for fleece uniformity. Genetic correlation of the trait with its environmental variability was 0.76±0.13 showing that fleece uniformity will be indirectly selected when aiming to reduce the fiber diameter. In the light of these parameters, and due to the cost of registering and the cost of opportunity, it looks no worthy to include uniformity as a selection criterion in alpaca breeding programs.

Differential sensitivity of climate conditions on birth weight genetic values in mice divergently selected for birth weight residual variance

After 32 generations of a divergent selection experiment for residual variance of birth weight in mice, two divergent lines were thus obtained: the heterogeneous line (H-line) and the homogeneous line (L-line). Throughout the generations, differences were observed between the two lines in traits such as litter size, survival at weaning, and birth weight variability caused by unidentified environmental conditions. The L-line exhibited advantages in terms of higher survival rates, larger litter sizes, and less sensitivity to changes in food intake. The study is an examination of the effects of climate as an environmental factor on the performance of these animals. Climate factors including maximum, minimum, and mean temperature (T), humidity (H), and TH index; at three stages (the fecundation, a week before the parturition and the parturition), were linked to a birth weight dataset consisting of 22,614 records distributed as follows: 8,853 corresponding to the H-line, 12,649 to the L-line, and 1,112 to the initial population. Out of the 27 analyzed climatic variables, the maximum temperature 1 wk before parturition (MXTW) was identified as the most influential when comparing heteroscedastic models with the deviance information criterion. The order of Legendre polynomial to apply in the following random regression model was tested by a cross-validation using homoscedastic models. Finally, MXTW was compared on how it affected the two divergent lines by analyzing predicted breeding values (PBV) obtained from a random regression heteroscedastic model. The mean PBV of the H-line in the first generation showed a range of 0.070 g with a negative slope, which was 35 times higher than the range obtained for the L-line, which varied within 0.002 g. In the last generation of selection, the H-line exhibited greater instability of PBV across temperatures, with a difference of 0.101 g between the maximum and minimum mean PBV, compared to 0.017 g for the L-line. The standard deviations of the slopes in the H-line were more dispersed than in the L-line. Unlike the H-line, the L-line had slopes that were not significantly different from 0 throughout the generations of selection, indicating greater stability in response to MXTW variations. The H-line exhibited a higher sensitivity to changes in MXTW, particularly in birth weight, with the L-line being more stable. The selection for uniformity of birth weight could lead to less sensitive animals under environmental changes.

Partial diallel cross for assessing genetic merit of local rabbit breed

This study was carried out to estimate general combining ability (GCA) and specific combining ability (SCA) of pre-and post-weaning traits from a partial diallel cross in three rabbit breeds: Local rabbit (L), Flemish Giant (F), and Rex (R). Body weight at 0 (BW0), 30 (BW30), 42 (BW42), and 63 (BW63) days of age, average daily gains from 0 to 30 d of age (ADG0-30), from 30 to 42 d of age (ADG30-42), from 42 to 63 d of age (ADG42-63), and from 30 to 63 d of age (ADG30-63), litter size at birth (LS0) and at weaning (LS42), and mortality at weaning (MR42) were studied in crossing LL, FF, RR, LF, LR, and FR. Local breed had the highest GCA for BW0, BW30, BW42, and average daily gain before weaning (ADG0-30, and ADG30-42) compared to Flemish Giant and Rex, while GCA of Local breed for average daily gain after weaning (ADG42-63), litter size (LS0, and LS42), and mortality (MR42) was higher than for Rex and similar to that of Flemish Giant. Crossing LF and LR showed higher SCA for BW30, BW42, BW63, ADG0-30, ADG42-63 and ADG30-63 than FR. In conclusion, based on GCA and SCA, the Indonesian Local breed has a high genetic potential in the crossing with Flemish Giant and Rex breeds.

Genetic factors of functional traits

Selection of functional traits is a challenge for researchers, but an increasingly necessary objective due to the growing concern regarding animal welfare and overcoming the problems of reducing antibiotic use in rabbit production without undermining the animals’ productivity. The aim of this review is to discuss the genetic control of resistance to diseases, longevity and variability of birth weight within a litter, or litter size variability at birth within doe, describing the selection programmes and the first results from a multi-omics analysis of resistance/susceptibility to diseases. The heritability is around 0.13 for longevity, 0.01 for uniformity in birth weight, 0.09 for litter size variability and around 0.11 for disease resistance. Genetic correlations between functional traits and production traits are mostly no different from zero, or are moderately favourable in some cases. Six selection programmes developed in three countries are reviewed. Line foundation with high pressure for selection or divergent selection experiments are different methodologies used, and favourable responses to selection have been achieved. Genomics studies have revealed associations in regions related to immune system functionality and stress in lines selected for litter size variability. Knowledge of the role of gut microbiota in the rabbit’s immune response is very limited. A multi-omics approach can help determine the microbial mechanisms in regulation immunity genes of the host.

Selection for environmental variance of litter size in rabbits involves genes in pathways controlling animal resilience

Background

Environmental variance (V_E) is partially under genetic control, which means that the V_E of individuals that share the same environment can differ because they have different genotypes. Previously, a divergent selection experiment for V_E of litter size (LS) during 13 generations in rabbit yielded a successful response and revealed differences in resilience between the divergent lines. The aim of the current study was to identify signatures of selection in these divergent lines to better understand the molecular mechanisms and pathways that control V_E of LS and animal resilience. Three methods (F_ST, ROH and varLD) were used to identify signatures of selection in a set of 473 genotypes from these rabbit lines (377) and a base population (96). A whole-genome sequencing (WGS) analysis was performed on 54 animals to detect genes with functional mutations.

Results

By combining signatures of selection and WGS data, we detected 373 genes with functional mutations in their transcription units, among which 111 had functions related to the immune system, stress response, reproduction and embryo development, and/or carbohydrate and lipid metabolism. The genes TTC23L, FBXL20, GHDC, ENSOCUG00000031631, SLC18A1, CD300LG, MC2R, and ENSOCUG00000006264 were particularly relevant, since each one carried a functional mutation that was fixed in one of the rabbit lines and absent in the other line. In the 3ʹUTR region of the MC2R and ENSOCUG00000006264 genes, we detected a novel insertion/deletion (INDEL) variant.

Conclusions

Our findings provide further evidence in favour of V_E as a measure of animal resilience. Signatures of selection were identified for V_E of LS in genes that have a functional mutation in their transcription units and are mostly implicated in the immune response and stress response pathways. However, the real implications of these genes for V_E and animal resilience will need to be assessed through functional analyses.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00653-y.

Evaluation of direct and maternal responses in reproduction traits based on different selection strategies for postnatal piglet survival in a selection experiment

BACKGROUND

Postnatal piglet survival is important both in economic and animal welfare terms. It is influenced by the piglet's own direct genetic effects and by maternal genetic effects of the dam, associated with milk production and mothering abilities. These genetic effects might be correlated, affected by other non-genetic factors and unfavourably associated with other reproduction traits such as litter size, which makes the development of optimal breeding strategies a challenge. To identify the optimum selection strategy for piglet survival, a selection experiment was carried out to compare responses in survival and reproduction traits to selection on only direct, only maternal, or both genetic effects of postnatal survival. The data of the experiment were recorded from outdoor reared pigs, with first- and second-generation sires selected based on their estimated breeding values for maternal and direct effects of postnatal survival of indoor reared offspring, respectively, with the opportunity to identify potential genotype-by-environment interaction.

RESULTS

A Bayesian multivariate threshold-linear model that was fitted to data on 22,483 piglets resulted in significant (Pr(h2 > 0) = 1.00) estimates of maternal and direct heritabilities between 0.12 and 0.18 for survival traits and between 0.29 and 0.36 for birth weight, respectively. Selection for direct genetic effects resulted in direct and maternal responses in postnatal survival of 1.11% ± 0.17 and - 0.49% ± 0.10, respectively, while selection for maternal genetic effects led to greater direct and maternal responses, of 5.20% ± 0.34 and 1.29% ± 0.20, respectively, in part due to unintentional within-litter selection. Selection for both direct and maternal effects revealed a significant lower direct response (- 1.04% ± 0.12) in comparison to its expected response from single-effect selection, caused by interactions between direct and maternal effects.

CONCLUSIONS

Selection successfully improved post- and perinatal survival and birth weight, which indicates that they are genetically determined and that genotype-by-environment interactions between outdoor (experimental data) and indoor (selection data) housed pigs were not important for these traits. A substantially increased overall (direct plus maternal) response was obtained using selection for maternal versus direct or both direct and maternal effects, suggesting that the maternal genetic effects are the main limiting factor for improving piglet survival on which selection pressure should be emphasized.

Effects of genetics and early-life mild hypoxia on size variation in farmed gilthead sea bream (Sparus aurata)

The present study evaluated, in an 18-month gilthead sea bream trial, the time course effects of genetics on individual size variation and growth compensation processes in families selected by heritable growth in the PROGENSA^® breeding program. Families categorized as fast, intermediate, and slow growing had different growth trajectories with a more continuous growth in fast growth families. This feature was coincident with a reduced size variation at the beginning of the trial that clustered together the half-sib families sharing the same father. Regression analysis evidenced that the magnitude of compensatory growth was proportional to the initial size variation with no rescaling of families at this stage. By contrast, the finishing growth depensation process can mask, at least partially, the previous size convergence. This reflects the different contribution across the production cycle of genetics in growth. How early-life experiences affect growth compensation at juvenile stages was also evaluated in a separate cohort, and intriguingly, a first mild-hypoxia pulse at 60-81 days post-hatching (dph) increased survival rates by 10%, preventing growth impairment when fish were exposed to a second hypoxia episode (112-127 dph). The early hypoxia experience did not have a negative impact on growth compensatory processes at juvenile stages. By contrast, a diminished capacity for growth compensation was found with repeated or late hypoxia experiences. All this reinforces the use of size variation as a main criterion for improving intensive fish farming and selective breeding.

Correlated Response to Selection for Litter Size Residual Variability in Rabbits’ Body Condition

A divergent selection experiment for residual variance of litter size at birth was carried out in rabbits during twelve generations. Residual variance of litter size was estimated as the within-doe variance of litter size after pre-correction for year and season as well as parity and lactation status effects. The aim of this work was to study the correlated response to selection for litter size residual variability in body condition from mating to weaning. Body condition is related directly to an animal’s fat deposits. Perirenal fat is the main fat deposit in rabbits. Individual body weight (IBW) and perirenal fat thickness (PFT) were used to measure body condition at second mating, delivery, 10 days after delivery, and weaning. Litter size of the first three parities was analyzed. Both lines decreased body condition between mating to delivery; however, the decrease in body condition at delivery was lower in the low line, despite this line having higher litter size at birth (+0.54 kits, p = 0.93). The increment of body condition between delivery and early lactation was slightly higher in the low line. On the other hand, body condition affected success of females’ receptivity and fertility at the third mating, e.g., receptive females showed a higher IBW and PFT than unreceptive ones (+129 g and +0.28 mm, respectively), and fertile females had a higher IBW and PFT than unfertile ones (+82 g and +0.28 mm, respectively). In conclusion, the does selected for reducing litter size variability showed a better deal with situations of high-energy demand, such as delivery and lactation, than those selected for increasing litter size variability, which would agree with the better health and welfare condition in the low line.

Inflammatory Correlated Response in Two Lines of Rabbit Selected Divergently for Litter Size Environmental Variability

A divergent selection experiment for environmental variance of litter size variance was carried out in rabbits over thirteen generations. The aim of this study was to evaluate the inflammatory response in the two lines of the experiment, in order to analyse the effect of selection on susceptibility to diseases after challenging to stressful situations, such as 24 h after the first delivery. A total of 78 females were used in this study, 39 from each line. The line selected for litter size heterogeneity (the high line) showed lower white blood leukocyte count (WBC; -0.87 × 103/µL), lower percentage of basophils (-0.11%), higher concentration of TNF-α (+13.8 pg/mL), and greater concentration of CRP (+38.1 µg/mL) than the line selected for litter size homogeneity (the low line). The high line had also higher concentrations of bilirubin, cholesterol, gamma-glutamyl transferase (GGT) and alkaline phosphatase (ALP) compared to the low line (difference between lines were +0.08 µmol/L, +0.14 µmol/L, +0.35 U/L and +2.4 U/L, respectively). The high line showed higher inflammatory response than the low line, in accordance with a larger susceptibility to infectious disorders. In conclusion, the line selected to increase litter size environmental variability seems to have poor capacity coping with environmental stressors. Therefore, selection for litter size environmental variability can be a useful way to improve animal welfare.

Selection Response in a Divergent Selection Experiment for Birth Weight Variability in Mice Compared with a Control Line

Birth weight (BW) in animal production is an economically important trait in prolific species. The laboratory mouse (Mus musculus) is used as an experimental animal because it is considered a suitable model for prolific species such as rabbits and pigs. Two mouse lines were divergently selected for birth weight variability with a third line of non-selected control population of the same origin as the animals starting the experiment. The objective of this study was, therefore, to compare and evaluate the differential response of each line. The animals were from the 17th generation of both low and high BW variability lines of the divergent selection experiment, including in addition animals from the control line. The dataset contained 389 records from 48 litters of the high line, 734 records from 73 litters of the low line, and 574 records from 71 litters of the control line. The studied traits were as follows: the BW, the BW variance, the BW standard deviation, the BW coefficient of variation within-litter, the weaning weight (WW), the litter size at birth and at weaning, the weight gain, and the preweaning survival. The model included the line effect jointly with the parturition number and its interaction, the linear and quadratic LS as covariates except for the LS trait itself when analyzing litter traits, as well as the pup sex when analyzing individual traits. The low line had a lower BW and WW, but a higher litter size, and greater robustness owing to a higher survival at weaning. As a model of livestock animals, the findings from this experiment led to a proposal of selection for pig production that would combine an increase in litter size with higher survival and welfare. Compared with the control line, a much higher response was observed in the low variability line than in the high line, making it extremely satisfactory given that homogeneity provides advantages in terms of animal welfare and robustness.

Identification of functional mutations associated with environmental variance of litter size in rabbits

Background

Environmental variance (V_E) is partly under genetic control and has recently been proposed as a measure of resilience. Unravelling the genetic background of the V_E of complex traits could help to improve resilience of livestock and stabilize their production across farming systems. The objective of this study was to identify genes and functional mutations associated with variation in V_E of litter size (LS) in rabbits. To achieve this, we combined the results of a genome-wide association study (GWAS) and a whole-genome sequencing (WGS) analysis using data from two divergently selected rabbit lines for high and low V_E of LS. These lines differ in terms of biomarkers of immune response and mortality. Moreover, rabbits with a lower V_E of LS were found to be more resilient to infections than animals with a higher V_E of LS.

Results

By using two GWAS approaches (single-marker regression and Bayesian multiple-marker regression), we identified four genomic regions associated with V_E of LS, on chromosomes 3, 7, 10, and 14. We detected 38 genes in the associated genomic regions and, using WGS, we identified 129 variants in the splicing, UTR, and coding (missense and frameshift effects) regions of 16 of these 38 genes. These genes were related to the immune system, the development of sensory structures, and stress responses. All of these variants (except one) segregated in one of the rabbit lines and were absent (n = 91) or fixed in the other one (n = 37). The fixed variants were in the HDAC9, ITGB8, MIS18A, ENSOCUG00000021276 and URB1 genes. We also identified a 1-bp deletion in the 3′UTR region of the HUNK gene that was fixed in the low V_E line and absent in the high V_E line.

Conclusions

This is the first study that combines GWAS and WGS analyses to study the genetic basis of V_E. The new candidate genes and functional mutations identified in this study suggest that the V_E of LS is under the control of functions related to the immune system, stress response, and the nervous system. These findings could also explain differences in resilience between rabbits with homogeneous and heterogeneous V_E of litter size.

Patterns of environmental variance across environments and traits in domestic cattle

The variance in phenotypic trait values is a product of environmental and genetic variation. The sensitivity of traits to environmental variation has a genetic component and is likely to be under selection. However, there are few studies investigating the evolution of this sensitivity, in part due to the challenges of estimating the environmental variance. The livestock literature provides a wealth of studies that accurately partition components of phenotypic variance, including the environmental variance, in well-defined environments. These studies involve breeds that have been under strong selection on mean phenotype in optimal environments for many generations, and therefore represent an opportunity to study the potential evolution of trait sensitivity to environmental conditions. Here, we use literature on domestic cattle to examine the evolution of micro-environmental variance (CVR-the coefficient of residual variance) by testing for differences in expression of CVR in animals from the same breed reared in different environments. Traits that have been under strong selection did not follow a null expectation of an increase in CVR in heterogenous environments (e.g., grazing), a pattern that may reflect evolution of increased uniformity in heterogeneous environments. When comparing CVR across environments of different levels of optimality, here measured by trait mean, we found a reduction in CVR in the more optimal environments for both life history and growth traits. Selection aimed at increasing trait means in livestock breeds typically occurs in the more optimal environments, and we therefore suspect that the decreased CVR is a consequence of evolution of the expression of micro-environmental variance in this environment. Our results highlight the heterogeneity in micro-environmental variance across environments and point to possible connections to the intensity of selection on trait means.

Genetic and Environmental Factors Influencing Gestation Length in Lusitano Horses

Breeding and parturition records collected over a period of 35 years in the Alter Real stud of Lusitano horses were used to calculate gestation length (GL). The 1027 gestations by 209 mares mated to 60 stallions had a mean GL of 338.1 ± 9.26 days. The mixed model analysis of variance indicated that the sex of the foal and inbreeding of the dam and foal had no significant effect on GL (P > .05). On the other hand, GL increased linearly with mare age, with an estimated regression coefficient of 0.155 ± 0.069 days/year (P < .05). Year and conception month affected GL (P < .05), with longer gestations observed when breeding occurred between January and April, followed by a decline of about 5 and 10 days in GL of mares bred in May and in June-July, respectively. Important differences were observed between stallions (P < .05), with most stallions resulting in a distribution of GL in an interval of ±5 days. The inclusion of the mare as a random effect in the mixed model resulted in an estimated repeatability of GL of 0.427, indicating that mares tend to be regular in having long or short gestations across their lifetime. Variance components estimated in an Animal Model resulted in heritability estimates of 0.39 for maternal genetic effects and 0.19 for direct genetic effects, with no association between the two components. Overall, the mare seems to have the major genetic influence on GL in Lusitano horses, but environmental factors, such as month of conception and also the age of the mare, have a noticeable effect.

Genetics and genomics of uniformity and resilience in livestock and aquaculture species: A review

Genetic control of residual variance offers opportunities to increase uniformity and resilience of livestock and aquaculture species. Improving uniformity and resilience of animals will improve health and welfare of animals and lead to more homogenous products. Our aims in this review were to summarize the current models and methods to study genetic control of residual variance, genetic parameters and genomic results for residual variance and discuss future research directions. Typically, the genetic coefficient of variation is high (median = 0.27; range 0-0.86) and the heritability of residual variance is low (median = 0.01; range 0-0.10). Higher heritabilities can be achieved when increasing the number of records per animal. Divergent selection experiments have supported the feasibility of selecting for high or low residual variance. Genomic studies have revealed associations in regions related to stress, including those from the heat shock protein family. Although the number of studies is growing, genetic control of residual variance is still poorly understood, but big data and genomics offer great opportunities.

The effect of environmental enrichment on behavioral variability depends on genotype, behavior, and type of enrichment

Non-genetic individuality in behavior, also termed intragenotypic variability, has been observed across many different organisms. A potential cause of intragenotypic variability is sensitivity to minute environmental differences during development, which are present even when major environmental parameters are kept constant. Animal enrichment paradigms often include the addition of environmental diversity, whether in the form of social interaction, novel objects or exploratory opportunities. Enrichment could plausibly affect intragenotypic variability in opposing ways: it could cause an increase in variability due to the increase in microenvironmental variation, or a decrease in variability due to elimination of aberrant behavior as animals are taken out of impoverished laboratory conditions. In order to test these hypothesis, we assayed five isogenic Drosophila melanogaster lines raised in control and mild enrichment conditions, and one isogenic line under both mild and intense enrichment conditions. We compared the mean and variability of six behavioral metrics between our enriched fly populations and the laboratory housing control. We found that enrichment often caused a small increase in variability across most of our behaviors, but that the ultimate effect of enrichment on both behavioral means and variabilities was highly dependent on genotype and its interaction with the particular enrichment treatment. Our results support previous work on enrichment that presents a highly variable picture of its effects on both behavior and physiology.

The Statistical Scale Effect as a Source of Positive Genetic Correlation Between Mean and Variability: A Simulation Study

The selection objective for animal production is the highest income with the lowest production cost, while ensuring the highest animal welfare. A selection experiment for environmental variability of birth weight in mice showed a correlated response in the mean after 20 generations starting from a crossed panmictic population. The relationship between the birth weight and its environmental variability explained the correlated response. The scale effect represents a potential cause of this correlation. The relationship between the mean and the variability implies: the higher the mean, the higher the variability. The study was to quantify by simulation the genetic correlation between a trait and its environmental variability. This can be attributable to the scale effect in a range of coefficients of variation and heritabilities between 0.05 and 0.50. The resulting genetic correlation ranged from 0.1335 to 0.7021 being the highest for the highest heritability and the lowest CV. The scale effect for a trait with heritability between 0.25 and 0.35 and CV between 0.15 and 0.25 generated a genetic correlation between 0.43 and 0.57. The genetic coefficient of variation (GCV) affecting residual variability was modulated by the strength reducing the impact of the scale effect. GCV ranged from 0.0050 to 1.4984. The strength of the scale effect might be in the range between 0 and 1. The scale effect would explain many reported genetic correlation and the additive genetic variance for the variability. This is relevant when increasing the mean of a trait jointly with the reduction of its variability.

Litter Survival Differences between Divergently Selected Lines for Environmental Sensitivity in Rabbits

Simple Summary

Two rabbit lines are divergently selected for increasing or decreasing the variability of litter size at birth. Decreasing the litter size variability produces more resilient females with less sensitivity to diseases, being an indirect selection way to improve environmental sensitivity. The kits’ survival rate at weaning was higher in the homogeneous line. Moreover, this line led to a greater uniformity of the kits’ weight at weaning, although the weight variability at birth was higher, which could be due to a higher lactation capacity of the homogeneous line.

Abstract

A divergent selection experiment on environmental sensitivity was performed in rabbits. The aim was to estimate the correlated response in kit weight and survival, litter weight, and weight distance from birth to weaning. The weight distance was calculated as the absolute value of the differences between the individual value and the mean value of its litter. The relationship between the probability of survival at 4 d of age, and the weight at birth, was studied. Environmental sensitivity was measured as litter size variability. A total of 2484 kits from 127 does from the low line, and 1916 kits of 114 does from the high line of the 12th generation were weighed. Both of the lines showed similar individual and litter weights at birth and weaning, and a similar survival rate at birth, and at 4 d of age. The survival rate at weaning was higher in the low line (0.67 and 0.62; P = 0.93). The weight distance was higher at birth, but lower at weaning in the low line (47.8 g and 54.1 g; P = 0.98). When the weight at birth was high, the kits had a higher survival rate. In conclusion, selection for environmental sensitivity showed a correlated response in the kits’ survival, and in the homogeneity of litter weight at weaning.

Effect of feed restriction on the environmental variability of birth weight in divergently selected lines of mice

BACKGROUND

Selection of mice for decreased environmental variability of birth weight has achieved higher survivability and larger litter size as a correlated response to canalized selection, which suggests higher welfare and robustness, and animals that are more homogeneous. However, in these studies, animals were not exposed to an environmental challenge. To demonstrate the advantages of this mouse line with a low environmental variability of birth weight, animals from two divergent lines (high and low variability of birth weight) were subjected to feed restriction. The objective of this study was to use these divergent lines to compare their response in terms of robustness against an environmental challenge. At weaning, 120 females, i.e. four full-sib females from 10 random litters of three consecutive generations of selection, were chosen from these divergent lines. The total number of females was divided into four groups, which were subjected to a feeding regimen by imposing different levels of feed restriction (i.e. 75, 90 and 85% of full ad libitum feed across three generations, respectively) in different combinations during the growth and reproduction periods.

RESULTS

Animals from the "low" line were less sensitive to a change in feed level than those from the "high" line. Regarding reproduction, the "low" line performed better in terms of number of females having parturitions, number of parturitions, and litter size. Imposing a feed restriction on female mice during their growth period did not affect the birth weight of their pups. The "low" line was preferred because of its higher reproductive efficiency and survival under an environmental challenge.

CONCLUSIONS

Selection for decreased environmental variability of birth weight produces animals that are less sensitive to environmental conditions, which can be interpreted as having greater robustness.

Stage- and thermal-specific genetic architecture for preadult viability in natural populations of Drosophila melanogaster

Studying the processes affecting variation for preadult viability is essential to understand the evolutionary trajectories followed by natural populations. This task requires focusing on the complex nature of the phenotype-genotype relationship by taking into account usually neglected aspects of the phenotype and recognizing the modularity between different ontogenetic stages. Here, we describe phenotypic variability for viability during the larval and pupal stages in lines derived from three natural populations of Drosophila melanogaster, as well as the variability for phenotypic plasticity and canalization at two different rearing temperatures. The results indicate that the three populations present significant phenotypic differences for preadult viability. Furthermore, distinct aspects of the phenotype (means, plasticity, canalization, plasticity of canalization) are affected by different genetic bases underlying changes in viability in a stage- and environment-specific manner. These findings explain the generalized maintenance of genetic variability for this fitness trait.

Invited review: Piglet survival: benefits of the immunocompetence

Piglet mortality has a negative impact on animal welfare and public acceptance. Moreover, the number of weaned piglets per sow mainly determines the profitability of piglet production. Increased litter sizes are associated with lower birth weights and piglet survival. Decreased survival rates and performance of piglets make the control of diseases and infections within pig production even more crucial. Consequently, selection for immunocompetence becomes an important key aspect within modern breeding programmes. However, the phenotypic recording of immune traits is difficult and expensive to realize within farm routines. Even though immune traits show genetic variability, only few examples exist on their respective suitability within a breeding programme and their relationships to economically important production traits. The analysis of immune traits for an evaluation of immunocompetence to gain a generally improved immune response is promising. Generally, in-depth knowledge of the genetic background of the immune system is needed to gain helpful insights about its possible incorporation into breeding programmes. Possible physiological drawbacks for enhanced immunocompetence must be considered with regards to the allocation theory and possible trade-offs between the immune system and performance. This review aims to discuss the relationships between the immunocompetence of the pig, piglet survival as well as the potential of these traits to be included into a breeding strategy for improved robustness.

Analysis of lambing distribution in the Ripollesa sheep breed. I. Development and comparison of circular von Mises models

Circular data originates in a wide range of scientific fields and can be analyzed on the basis of directional statistics and special distributions wrapped around the circumference. However, both propensity to transform non-linear to linear data and complexity of directional statistics limited the generalization of the circular paradigm in the animal breeding framework, among others. Here, we generalized a circular mixed (CM) model within the context of Bayesian inference. Three different parametrizations with different hierarchical structures were developed on basis of the von Mises distribution; moreover, both goodness of fit and predictive ability from each parametrization were compared through the analyses of 110 116 lambing distribution records collected from Ripollesa sheep herds between 1976 and 2017. The naive circular (NC) model only accounted for population mean and homogeneous circular variance, and reached the lowest goodness-of-fit and predictive ability. The CM model assumed a hierarchical structure for the population mean by accounting for systematic (ewe age and lambing interval) and permanent environmental sources of variation (flock-year-season and ewe). This improved goodness of fit by reducing both the deviance information criterion (DIC; -2520 units) and the mean square error (MSE; -12.4%) between simulated and predicted lambing data when compared against the NC model. Finally, the last parametrization expanded CM model by also assuming a hierarchical structure with systematic and permanent environmental factors for the variance parameter of the von Mises distribution (i.e. circular canalization (CC) model). This last model reached the best goodness of fit to lambing distribution data with a DIC estimate 5425 units lower than the one for NC model (MSE reduced 13.2%). The same pattern revealed when models were compared in terms of predictive ability. The superiority revealed by CC model emphasized the relevance of heteroskedasticity for the analysis of lambing distribution in the Ripollesa breed, and suggested potential applications for the sheep industry, even genetic selection for canalization. The development of CM models on the basis of the von Mises distribution has allowed to integrate flexible hierarchical structures accounting for different sources of variation and affecting both mean and dispersion terms. This must be viewed as a useful statistical tool with multiple applications in a wide range of research fields, as well as the livestock industry. The next mandatory step should be the inclusion of genetic terms in the hierarchical structure of the models in order to evaluate their potential contribution to current selection programs.

Correlated response to selection for litter size environmental variability in rabbits' resilience

Resilience is the ability of an animal to return soon to its initial productivity after facing diverse environmental challenges. This trait is directly related to animal welfare and it plays a key role in fluctuations of livestock productivity. A divergent selection experiment for environmental variance of litter size has been performed successfully in rabbits over ten generations. The objective of this study was to analyse resilience indicators of stress and disease in the divergent lines of this experiment. The high line showed a lower survival rate at birth than the low line (-4.1%). After correcting by litter size, the difference was -3.2%. Involuntary culling rate was higher in the high than in the low line (+12.4%). Before vaccination against viral haemorrhagic disease or myxomatosis, concentration of lymphocytes, C-reactive protein (CRP), complement C3, serum bilirubin, triglycerides and cholesterol were higher in the high line than in the low line (difference between lines +4.5%, +5.6 µg/ml, +4.6 mg/ml, +7.9 mmol/l, +0.3 mmol/l and +0.4 mmol/l). Immunological and biochemical responses to the two vaccines were similar. After vaccination, the percentage of lymphocytes and CRP concentration were higher in the low line than in the high one (difference between lines +4.0% and +13.1 µg/ml). The low line also showed a higher increment in bilirubin and triglycerides than the high line (+14.2 v. +8.7 mmol/l for bilirubin and +0.11 v. +0.01 mmol/l for triglycerides); these results would agree with the protective role of bilirubin and triglycerides against the larger inflammatory response found in this line. In relation to stress, the high line had higher basal concentration of cortisol than the low line (+0.2ng/ml); the difference between lines increased more than threefold after the injection of ACTH 1 to 24, the increase being greater in the high line (+0.9 ng/ml) than in the low line (+0.4 ng/ml). Selection for divergent environmental variability of litter size leads to dams with different culling rate for reproductive causes and different kits' neonatal survival. These associations suggest that the observed fitness differences are related to differences in the inflammatory response and the corticotrope response to stress, which are two important components of physiological adaptation to environmental aggressions.

Adaptation to environmental stress at different timescales

Environments are changing rapidly, and to cope with these changes, organisms have to adapt. Adaptation can take many shapes and occur at different speeds, depending on the type of response, the trait, the population, and the environmental conditions. The biodiversity crisis that we are currently facing illustrates that numerous species and populations are not capable of adapting with sufficient speed to ongoing environmental changes. Here, we discuss current knowledge on the ability of animals and plants to adapt to environmental stress on different timescales, mainly focusing on thermal stress and ectotherms. We discuss within-generation responses that can be fast and induced within minutes or hours, evolutionary adaptations that are often slow and take several generations, and mechanisms that lay somewhere in between and that include epigenetic transgenerational effects. To understand and predict the impacts of environmental change and stress on biodiversity, we suggest that future studies should (1) have an increased focus on understanding the type and speed of responses to fast environmental changes; (2) focus on the importance of environmental fluctuations and the predictability of environmental conditions on adaptive capabilities, preferably in field studies encompassing several fitness components; and (3) look at ecosystem responses to environmental stress and their resilience when disturbed.

Correlated response in body condition and energy mobilisation in rabbits selected for litter size variability

A divergent selection experiment on litter size variability (high and low lines) was performed in rabbits over seven generations. The aim of this study was to evaluate the correlated responses to selection in body condition and fat reserves mobilisation. Litter size variability was estimated as phenotypic variance of litter size within female after correcting for the year-season and the parity-lactation status effects. A total of 226 females were used in this study, of which 158 females were used to measure body condition and energy mobilisation. Body condition was measured as BW and perirenal fat thickness. Females were stimulated with the adrenergic isoproterenol. Mobilisation capacity of fat reserves was measured by the lipolytic potential, defined as the increment in non-esterified fatty acids (NEFA) levels from basal concentration until adrenergic stimulation at mating, delivery and 10 days after delivery of the second reproductive cycle. Females were classified as survivor or non-survivor when they were culled for sanitary reasons or died before the third kindling. Data were analysed using Bayesian methodology. Survivor females presented higher BW than the non-survivor females at delivery (238 g, P=1.00) and 10 days after delivery (276 g, P=1.00). They also showed higher perirenal fat thickness at 10 days after delivery (0.62 mm, P=1.00). At delivery, basal NEFA levels was lower in survivor than non-survivor females (-0.18 mmol/l, P=1.00), but their lipolytic potential (∆NEFA) was higher (0.08 mmol/l, P=0.94). Body weight was similar between lines in survivor females. Perirenal fat thickness was lower in the high line than in the low line at delivery (-0.23 mm, P=0.90) and 10 days after delivery (-0.28 mm, P=0.92). The high line exhibited higher NEFA (0.10 mmol/l, P=0.93) and lower ∆NEFA (-0.08 mmol/l, P=0.92) than the low line at delivery. The low line showed a favourable correlated response to selection on body condition and fat reserves mobilisation. In conclusion, the low line selected for litter size variability seems to adapt better to adverse conditions, as it has a greater capacity to mobilise energy reserves at delivery than the high line. Females that adequately manage their body reserves and perform energy mobilisation correctly have a lower risk of dying or being culled.

Feed and reproductive efficiency differences between divergently selected lines for birthweight environmental variability in mice

Sustainability has come to play an important role in agricultural production. A way to combine efficiency with sustainability might be by searching for robust animals that can be selected for the homogeneity of certain traits. Furthermore, the optimization of feed efficiency is one of the challenges to improve livestock genetics programmes, but this might compromise reproductive efficiency. Animals from two divergent mouse lines, regarding variability of birthweight, were used to check whether homogeneity was also related to both feed and reproductive efficiency. The objective of this study was to use these divergent lines of mice to compare them with their feed efficiency and the reproductive capacity. Animal weight, weight gain, feed intake, relative intake and cumulated transformation index were considered as feed efficiency traits. Animals from the low line had both lower weight and feed intake from 21 to 56 days. They had a worse transforming index in the three last weeks when litter size was fitted as an effect of the model, but the lines become similar if the higher litter size of the low line was not included. Reproductively, the low line performed better considering the number of females having parturitions, the number of parturitions, and with higher litter size and survival in both parturitions. Hence, the low variability line was preferred because of reproductive efficiency without seriously affecting its feed efficiency. Homogeneity seemed to be related to robustness with similar feed efficiency but higher reproductive efficiency.

Correlated response in early embryonic development in rabbits selected for litter size variability

<p>A divergent selection experiment for litter size variability was carried out in rabbits. The litter size variability was estimated as the phenotypic variance of litter size within female. The aim of this study was to assess the effect of selecting for litter size variability on early embryonic development and survival after 7 generations of divergent selection (high and low variability lines). A total of 30 non-lactating multiparous does per line were used. The ovulation rate and early embryonic development were analysed using Bayesian methodology. Ovulation rate was not affected by the selection process. At 28 h of gestation, embryonic development and survival were similar in both lines. At 48 h of gestation, the majority of embryos in the high line were in the early morulae stage. The high line had a higher proportion of early morulae (79.54 vs. 53.43%; P=0.94) and a lower proportion of compacted morulae (20.46 vs. 46.57%; P=0.93%) than the low line. At 72 h of gestation, the high line had 1.59 fewer embryos than the more homogeneous line (P=0.85), due to reduced embryonic survival (0.60 vs. 0.74; P=0.93). The high line continued to show a higher proportion of early morulae (21.01 vs. 3.69%; P=0.93) and a lower proportion of compacted morulae and blastocysts (78.99 vs. 96.31%; P=0.94) than the low line at 72 h of gestation, indicative of reduced embryonic development. In conclusion, selection for homogeneity in litter size had a positive impact on embryonic traits.</p>

Correlated response in litter size components in rabbits selected for litter size variability

A divergent selection experiment for the environmental variability of litter size (Ve) over seven generations was carried out in rabbits at the University Miguel Hernández of Elche. The Ve was estimated as the phenotypic variance within the female, after correcting for year-season and parity-lactation status. The aim of this study was to analyse the correlated responses to selection in litter size components. The ovulation rate (OR) and number of implanted embryos (IE) in females were measured by laparoscopy at 12 day of the second gestation. At the end of the second gestation, the total number of kits born was measured (TB). Embryonic (ES), foetal (FS) and prenatal (PS) survival were computed as IE/OR, TB/IE and TB/OR, respectively. A total of 405 laparoscopies were performed. Data were analysed using Bayesian methodology. The correlated response to selection for litter size environmental variability in terms of the litter size components was estimated as either genetic trends, estimated by computing the average estimated breeding values for each generation and each line, or the phenotypic differences between lines. The OR was similar in both lines. However, after seven generations of selection, the homogenous line showed more IE (1.09 embryos for genetic means and 1.23 embryos for phenotypic means) and higher ES than the heterogeneous one (0.07 for genetic means and 0.08 for phenotypic means). The probability of the phenotypic differences between lines being higher than zero (p) was 1.00 and .99, respectively. A higher uterine overcrowding of embryos in the homogeneous line did not penalize FS; as a result, this line continued to show a greater TB (1.01 kits for genetic means and 1.30 kits for phenotypic means, p = .99, in the seventh generation). In conclusion, a decrease in litter size variability showed a favourable effect on ES and led to a higher litter size at birth.