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.

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.

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.

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.