Correlated genetic trends for production and welfare traits in a mouse population divergently selected for birth weight environmental variability

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.

Genomic inbreeding measures applied to a population of mice divergently selected for birth weight environmental variance

This study aimed to compare different inbreeding measures estimated from pedigree and molecular data from two divergent mouse lines selected for environmental birth weight during 26 generations. Furthermore, the performance of different approaches and both molecular and pedigree data sources for estimating Ne were tested in this population. A total of 1,699 individuals were genotyped using a high-density genotyping array. Genomic relationship matrices were used to calculate molecular inbreeding: Nejati-Javaremi (F NEJ), Li and Horvitz (F L&H), Van Raden method 1 (F VR1) and method 2 (F VR2), and Yang (F YAN). Inbreeding based on runs of homozygosity (F ROH) and pedigree inbreeding (F PED) were also computed. F ROH, F NEJ, and F L&H were also adjusted for their average values in the first generation of selection and named F ROH0, F NEJ0, and F L&H0. ∆F was calculated from pedigrees as the individual inbreeding rate between the individual and his parents (∆F PEDt) and individual increases in inbreeding (∆F PEDi). Moreover, individual ∆F was calculated from the different molecular inbreeding coefficients (∆F NEJ0, ∆F L&H, ∆F L&H0, ∆F VR1, ∆F VR2, ∆F YAN, and ∆F ROH0). The Ne was obtained from different ∆F, such as Ne PEDt, Ne PEDi, Ne NEJ0, Ne L&H, Ne L&H0, Ne VR1, Ne VR2, Ne YAN, and Ne ROH0. Comparing with F PED , F ROH , F NEJ and F VR2 overestimated inbreeding while F NEJ0 , F L&H , F L&H0 , F VR1 and F YAN underestimated inbreeding. Correlations between inbreeding coefficients and ∆F were calculated. F ROH had the highest correlation with F PED (0.89); F YAN had correlations >0.95 with all the other molecular inbreeding coefficients. Ne PEDi was more reliable than Ne PEDt and presented similar behaviour to Ne L&H0 and Ne NEJ0. Stable trends in Ne were not observed until the 10th generation. In the 10th generation Ne PEDi was 42.20, Ne L&H0 was 45.04 and Ne NEJ0 was 45.05 and in the last generation these Ne were 35.65, 35.94 and 35.93, respectively F ROH presented the highest correlation with F PED, which addresses the identity by descent probability (IBD). The evolution of Ne L&H0 and Ne NEJ0 was the most similar to that of Ne PEDi. Data from several generations was necessary to reach a stable trend for Ne, both with pedigree and molecular data. This population was useful to test different approaches to computing inbreeding coefficients and Ne using molecular and pedigree data.

Embryo survival and fertility differ in lines divergently selected for birth weight homogeneity in mice

The selection of animals for lower environment sensitivity around the optimum trait value can also provide benefits in productivity and welfare. A divergent selection experiment for birth weight environmental variability in mice was successfully conducted over 17 generations. Animals from low variability selected line (L-line) were more robust by having a higher litter size and survival at weaning in a common breeding environment, than those from high variability line (H-line). The objective of this study was to analyze the differences between those divergently selected lines for embryo and fetal survival and for fertility and prolificacy rate. To study embryo survival and ovulation rate, a total of 98 females (34 H-line and 50 L-line) were studied in four generations of the divergent experiment. To analyze fetal survival and fertility rate, 378 female mice (138 H-line and 240 L-line) in 10 generations or the divergent experiment were studied. Ultrasound scans were performed at day 14 of gestation to establish the number of total fetal and the embryo absorptions. Mortality was addressed as the difference between litter size at birth and the number of fetuses at 14 days of gestation. The number of pregnant females in the first 3 days after mating was used to measure fertility. A linear model was also fitted to analyze embryo mortality, litter size, and the number of embryos at 14 days of gestation. A categorical model was then used to study fertility, including line, generation, and its interaction as effects. Despite the fact that there were no significant differences in the ovulation rate, litter size at birth was significantly higher in the L-line than in the H-line (9.82 vs. 8.36 pups, p < 0.001). Moreover, embryo mortality was significantly lower in the L-line than in the H-line (1.39 vs. 2.87 fetuses, p < 0.001). L-line females were more fertile (53.49% vs. 23.26% for the H-line). According to these results, the line selected for low environmental variance would be preferable for robustness and animal welfare.

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.

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.

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

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

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.

Relationship between litter size at birth and within-litter birth weight characteristics in laboratory mice as pilot animal for pig

We investigated the relationship between litter size at birth and within-litter birth weight (BW) characteristics of laboratory mice as a pilot mammal for pig. We obtained records of number born alive (NBA) and total and mean litter BW (LWB, MWB), and maximum and minimum values of within-litter BW (MaxIWB, MinIWB), range and standard deviation (Range, SDIWB), skewness (Skew), and kurtosis (Kurt) of within-litter BW for 656 litters at first parity. Pearson's correlations of NBA were highly positive with LWB (0.92), weakly negative with MWB (-0.31), MaxIWB (-0.19), and MinIBW (-0.33), and those of MWB were negligible with Range, SDIWB, Skew, and Kurt (-0.10 to 0.06). Estimated heritabilities, treated as dam traits, were 0.32 for NBA, 0.39 for LWB, 0.24 for MWB, 0.28 for MaxIWB, 0.05 for MinIWB, 0.16 for Range, 0.17 for SDIWB, and 0.00 for Skew and Kurt. Estimated genetic correlation between NBA and LWB was high (0.95). Therefore, LWB could be promising for efficiently improving NBA. The estimated genetic correlation of NBA was negligible with MWB (0.00), positive with MaxIWB (0.10), Range (0.48), and SDIWB (0.36), and negative with MinIWB (-0.36), suggesting that selection for increased NBA brings larger SDIWB and lighter MinIWB.

Genetic parameter and genetic trend estimates for reproductive traits in Wistar rats

ABSTRACT The aim of this study was to estimate genetic parameters and genetic trends for reproductive traits in Wistar rats. A total of 1,167 data records from 283 females over six generations of monogamous mating pairs was used. Heritability and genetic correlation were estimated through Bayesian inference and genetic trends were calculated by linear regression of breeding values over generations. Heritability estimates for litter size at birth (LS), calving interval (CI), pup mortality (PM) and maternal cannibalism (CAN) presented low magnitude, ranging from 0.01 to 0.13. CAN presented high and positive genetic correlation with LS and PM (0.77 and 0.78, respectively). On the other hand, all the other estimated genetic correlations were not significant. Genetic trend was positive for LS (+0.0900 pups per generation), and negative for PM and CAN (-1.0085 and -0.5217 pups per generation, respectively). For CI the genetic trend was not significant. It is recommended to increase selection intensity on dams in this Wistar rat population in order to accelerate the genetic progress.

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.

Challenging the selection for consistency in the rank of endurance competitions

Background

Control of the environmental variability by genetic selection offers possibilities for new selection objectives for productive traits. This methodology aims at reducing heterogeneity in productive traits and has been applied to several traits and species for which animal homogeneity is profitable. In horse breeding programmes, rank in competitions is a common selection objective but has been challenging to model. In this study, the parameters of environmental variability for the rank of a horse were computed to analyse the capability of a horse to maintain the best ranking across competitions that consist of long-distance races in which the adapted physical condition of the horse is essential. The genetic component of the environmental variance for the rank in endurance competitions was evaluated, which resulted in proposing a new transformation of horse scores in competitions.

Results

Homogeneous and heterogeneous variance models were compared by assaying three random effects that affect both the rank and its variability, using endurance ride data consisting of 2863 records. The pedigree relationship matrix contained 5931 animals. The rank trait was transformed into a normalized variable to prevent false estimates of the genetic correlation by inappropriate artificial skewness. The models included the number of participants in the race, sex, and age as systematic effects. The rider, the rider-horse interaction, or an environmental permanent effect were tested as random effects, in addition to additive genetic and residual effects. The models were analysed using the GSEVM program. Estimates of heritability for rank ranged from 0.12 to 0.15. The heterogeneous variance model that fitted the rider was assessed as the best model based on the deviance information criterion. Estimates of genetic variance for rank variability ranged from 0.12 to 0.13. The genetic correlation between the rank and its environmental variability was low and did not differ from 0.

Conclusions

These results offer an opportunity to select animals for canalization by reducing the variability of race results and achieving the best positions, which could be a new selection objective by weighting estimated breeding values for rank and its variability in a selection index.

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 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.

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.

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.

Modulating birth weight heritability in mice

Expected genetic response is proportional to the heritability of the trait, and this parameter is considered inherent of a specific trait in a particular population. However, models assuming heterogeneity in residual variance lead to different estimates of heritability across combinations of systematic (environmental) effects. Modifying the residual variance of the birth weight by artificial selection was shown to be feasible in a divergent selection experiment in mice. The objectives of this work were to 1) estimate the evolution of the heritability of birth weight in mice in the mentioned experiment, and 2) estimate different heritability regarding systematic effects. Data came from eleven generations of a divergent selection experiment to modify the residual variability of birth weight in mice. A total of 15,431 birth weight records from 959 females and 1,641 litters in combination with 14,786 pedigree records were used. The model used for analysis included generation, litter size, sex, and parity number as systematic effects. Each record of birth weight was assigned to the mother of the pup in the model which assumes that the residual variance is heterogeneous and partially under genetic control. Differences in heritability between lines reached values of 0.06 in the last generations. Choosing the most extreme values of systematic effects, the birth weight heritability ranged from 0.04 to 0.22. From these results, the possibility of modulating the heritability for this trait could be explored in 1 of 2 ways: selecting to decrease the residual variability, or choosing the specific levels of the systematic effects.

Selection for environmental variance of litter size in rabbits

Background

In recent years, there has been an increasing interest in the genetic determination of environmental variance. In the case of litter size, environmental variance can be related to the capacity of animals to adapt to new environmental conditions, which can improve animal welfare.

Results

We developed a ten-generation divergent selection experiment on environmental variance. We selected one line of rabbits for litter size homogeneity and one line for litter size heterogeneity by measuring intra-doe phenotypic variance. We proved that environmental variance of litter size is genetically determined and can be modified by selection. Response to selection was 4.5% of the original environmental variance per generation. Litter size was consistently higher in the Low line than in the High line during the entire experiment.

Conclusions

We conclude that environmental variance of litter size is genetically determined based on the results of our divergent selection experiment. This has implications for animal welfare, since animals that cope better with their environment have better welfare than more sensitive animals. We also conclude that selection for reduced environmental variance of litter size does not depress litter size.