Skip to main content

ORIGINAL RESEARCH article

Front. Genet., 14 December 2023
Sec. Livestock Genomics

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

  • 1Departamento de Producción Animal, Facultad de Veterinaria, Universidad Complutense de Madrid, Madrid, Spain
  • 2Departamento de Producción Agraria, E.T.S. Ingeniería Agronómica, Alimentaria y de Biosistemas, Universidad Politécnica de Madrid, Madrid, Spain

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 (FNEJ), Li and Horvitz (FL&H), Van Raden method 1 (FVR1) and method 2 (FVR2), and Yang (FYAN). Inbreeding based on runs of homozygosity (FROH) and pedigree inbreeding (FPED) were also computed. FROH, FNEJ, and FL&H were also adjusted for their average values in the first generation of selection and named FROH0, FNEJ0, and FL&H0. ∆F was calculated from pedigrees as the individual inbreeding rate between the individual and his parents (∆FPEDt) and individual increases in inbreeding (∆FPEDi). Moreover, individual ∆F was calculated from the different molecular inbreeding coefficients (∆FNEJ0, ∆FL&H, ∆FL&H0, ∆FVR1, ∆FVR2, ∆FYAN, and ∆FROH0). The Ne was obtained from different ∆F, such as NePEDt, NePEDi, NeNEJ0, NeL&H, NeL&H0, NeVR1, NeVR2, NeYAN, and NeROH0. Comparing with FPED, FROH, FNEJ and FVR2 overestimated inbreeding while FNEJ0, FL&H, FL&H0, FVR1 and FYAN underestimated inbreeding. Correlations between inbreeding coefficients and ∆F were calculated. FROH had the highest correlation with FPED (0.89); FYAN had correlations >0.95 with all the other molecular inbreeding coefficients. NePEDi was more reliable than NePEDt and presented similar behaviour to NeL&H0 and NeNEJ0. Stable trends in Ne were not observed until the 10th generation. In the 10th generation NePEDi was 42.20, NeL&H0 was 45.04 and NeNEJ0 was 45.05 and in the last generation these Ne were 35.65, 35.94 and 35.93, respectively FROH presented the highest correlation with FPED, which addresses the identity by descent probability (IBD). The evolution of NeL&H0 and NeNEJ0 was the most similar to that of NePEDi. 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.

1 Introduction

Inbreeding appears due to mating between related individuals and is related to negative consequences because of an increase in homozygosity, such as a reduction in fitness, namely, inbreeding depression (Hedrick, 2012; Alemu et al., 2021). Therefore, the study of inbreeding is essential in many areas, e.g., animal and plant breeding (Villanueva et al., 2021), human genetics (McQuillan et al., 2012), and evolutionary (Roff, 1997) and conservation biology (Frankham et al., 2010).

Traditionally, pedigree data has been used to measure inbreeding (Howard et al., 2017; Villanueva et al., 2021); however, this only provides the expected proportions of the genome that are identical by descent (IBD) and does not capture variation due to Mendelian sampling and linkage during gamete formation (Howard et al., 2017).

The implementation of molecular technologies has led to the development of numerous genomic estimators of inbreeding coefficients (Howard et al., 2017; Alemu et al., 2021; Villanueva et al., 2021), which can distinguish between individuals with the same ascendants and differentiate inbreeding at specific regions in chromosomes (Howard et al., 2017; Villanueva et al., 2021). Different approaches have been developed to measure inbreeding coefficients based on genomic data, including methods based on homozygous stretches of DNA sequences (runs of homozygosity—ROH) (Ceballos et al., 2018). However, ROH detection is highly dependent on the parameters set by the user (Peripolli et al., 2017; Rodríguez-Ramilo et al., 2019; Villanueva et al., 2021); therefore, comparisons between studies are not straightforward. In this context, other methodologies appeared that detected homozygosity by descent (HBD) segments under a probabilistic framework based on hidden Markov models. The probability of a single nucleotide polymorphism (SNP) being in an HBD segment was modelled as a function of genotyping error, intermarker distances and allele frequencies (Druet and Gautier, 2017). In addition, other methods are based on analysing each SNP and are derived from genomic relationship matrices (GRMs) (Li and Horvitz, 1953; Nejati-Javaremi et al., 1997; VanRaden, 2008; Yang et al., 2010; VanRaden et al., 2011). GRMs are used to obtain genomic predictions in animal breeding and can be used to estimate inbreeding given (Villanueva et al., 2021). These measures, however, are designed with different constraints and can lead to very different results with several implications (Villanueva et al., 2021).

One of the indicators most used to assess the genetic diversity of a population is the effective population size (Ne) (Wright, 1931). There are some approaches based on different indicators that have been used to measure Ne: the change in allele frequencies across generations (variance Ne), increase in homozygosity (inbreeding Ne), increase in kinship (coancestry Ne) or the rate of loss of less frequent alleles (eigenvalue Ne). Ne can also be estimated using census parameters, pedigree data, individual genotypes or demographic information (Sjödin et al., 2005; Leroy et al., 2013). The relationship between the increase in inbreeding (∆F) and Ne is defined by the classic formula: Ne =12F. Traditionally, ∆F has been calculated as the increase in the inbreeding rate between two successive generations (Wright, 1931). However, other methods have been proposed to estimate the individual increase in inbreeding based on pedigree knowledge, which has demonstrated less susceptibility to mating method changes (Gutiérrez et al., 2008; Gutiérrez et al., 2009; Cervantes et al., 2011). Moreover, the availability of molecular information has led to the estimation of molecular Ne. Some estimators have been proposed based in heterozygosity excess (Pudovkin et al., 1996), linkage disequilibrium (Hill, 1981), temporal changes in allele frequency (Krimbas and Tsakas, 1971; Nei and Tajima, 1981; Pollak, 1983), half-sib and full-sib (Wang, 2009), or in approximate Bayesian computation that take into account, among other parameters, the number of alleles per locus, the linkage disequilibrium or the mean and variance of multilocus homozygosity (Tallmon et al., 2008). Wang et al., 2016 discussed some concerns about these different Ne estimators as the reliability, the interpretation of each method, or the implicit assumption in some of these methods.

A divergent selection experiment for birth weight variability in mice has been successfully carried out, creating two lines, one homogeneous and the other heterogeneous over 26 generation (Formoso-Rafferty et al., 2016a). The experimental mating design avoided sharing grandparents and was optimised as a function of selection criteria in both selected lines (Formoso-Rafferty et al., 2016b). The homogeneous line presented advantages in production, animal welfare, response to selection, and traits related to robustness, such as longevity, survival and feed efficiency (Formoso-Rafferty et al., 2017; Formoso-Rafferty et al., 2018; Formoso-Rafferty et al., 2019; Formoso-Rafferty et al., 2020; Formoso-Rafferty et al., 2022; Formoso-Rafferty et al., 2023).

This study aimed to compare different inbreeding measures estimated from pedigree and molecular data from divergent mouse lines for environmental variability in birth weight. Furthermore, the performance of Ne estimated using different approaches and both molecular and pedigree data sources was tested in this divergent population. The applications of different methodologies in this particular population structure (strong divergent selection) are intended to evaluate the suitability in populations under selection.

2 Methods

Pedigree data were obtained from 26 discrete generations of a successful divergent selection experiment for environmental variability of birth weight in mice, including an additional five previous discrete generations randomly mated. In the experiment, in each generation, 43 females were mated to one male to give a maximum of two parturitions, avoiding mating animals sharing grandparents, with an approximately 30% proportion selected within the line. In the five preceding generations, 32 males were mated to 64 females, each male with 2 females. More details of the selection process can be found in Formoso-Rafferty et al. (2016a).

A total of 1824 individuals from the 26 selection generations were genotyped using the Affymetrix Mouse Diversity Genotyping Array with 616,137 SNPs. The individuals genotyped were those selected according to the selection criteria in each generation. During quality control (QC) animals with a call rate lower than 97% were removed, leaving 1,669 for analysis. Of these, 1,587 were females and 112 were males. The first generation of selected mice was used as the reference population (a total of 70 individuals, all females). The QC was also applied to the SNPs: 3% of missing genotypes were allowed, SNPs mapped to sex chromosomes were removed, and 545,656 SNPs were retained. This set was used to detect ROH to ensure maximum genome coverage and, therefore, avoid loss of information: no minor allele frequency (MAF) was applied as recommended other authors recommended (Meyermans et al., 2020).

Additional QC was applied to determine the other genomic inbreeding coefficients estimators: SNPs presenting MAFs lower than 0.05 in the reference population were removed. Additionally, remnant SNPs with MAFs less than 0.05 among the whole population were also removed. Finally, 173,546 SNPs were kept, which were used to obtain inbreeding coefficients from different GRMs. The number of animals genotyped per generation in the selected population is presented in Table 1.

TABLE 1
www.frontiersin.org

TABLE 1. Number of animals sampled per line and generation. H-Line: high variability line. L-Line: low variability line.

2.1 Inbreeding coefficients

Inbreeding coefficients from genomic data were obtained from the diagonal elements of different GRMs and ROH.

Inbreeding based on GRMs was named using the nomenclature chosen by Villanueva et al., (2021); therefore, inbreeding coefficients were calculated as follows:

1) FNEJ uses the diagonal elements of the allelic relationship matrix of Nejati-Javaremi et al. (1997). FNEJ is calculated using all individuals (FNEJT), those from the high variability line (FNEJH) and those from the low variability line (FNEJL).

2) FNEJ0 of an individual is calculated as FNEJ0=FNEJFNEJ1st¯1FNEJ1st, where FNEJ is the Nejati-Javaremi inbreeding coefficient and FNEJ1st¯ is the average Nejati-Javaremi inbreeding coefficient of the reference population.

3) FL&H estimates the deviation of the observed frequency of homozygotes from the expected frequency of homozygotes in the reference population under Hardy-Weinberg proportions (Li and Horvitz, 1953).

4) FL&H0 of an individual k is calculated as FL&H0=FL&HFL&H1st¯1FL&H1st, where FL&H is the Li and Horvitz inbreeding coefficient of the individual and FL&H1st¯ is the average FL&H of the reference population. Therefore, FL&H0 is expected to perform similarly to that of FNEJ0.

5) FVR1 is computed from the diagonal elements of the GRM obtained according to Van Raden’s method 1 (VanRaden, 2008).

6) FVR2 is computed from the diagonal elements of the GRM obtained according to Van Raden’s method 2 (Leutenegger et al., 2003; VanRaden, 2008).

7) FYAN uses the diagonal elements from the GRM of Yang to compute inbreeding coefficients (Yang et al., 2010).

8) Sliding window algorithms were used to detect ROH in each population. The parameters set were as follows: 50 SNPs per window; one heterozygote allowed in a window; no limit of the number of heterozygotes per ROH; five missing SNPs allowed in a window; the minimum length of an ROH was 1,000 kb; the minimum number of homozygous SNPs in an ROH was set at 100; the required minimum density was set at one SNP/50 kb; the window threshold was set at 0.5; and the minimum distance between two ROHs was 1,000 kb. Inbreeding based on ROH was calculated as FROHi (McQuillan et al., 2012).

9) FROH0 is calculated as FROH0=FROHFROH1st¯1FROH1st¯, where FROH is the ROH inbreeding coefficient and FROH1st¯ is the average ROH inbreeding coefficient of the reference population.

10) Pedigree inbreeding coefficients (FPED) are defined as the probability of an individual having two identical alleles by descent and are computed following Meuwissen and Luo (1992).

In summary, we tested several methods; some were intrinsically adjusted for the allele frequencies of the reference population (FL&H, FVR1, FVR2, and FYAN), others were adjusted for mean inbreeding of the reference population (FNEJ0, FL&H0, and FROH0), and, finally, some were not adjusted (FPED, FNEJ, and FROH). The genomic inbreeding estimators described above were chosen because they have been already tested in other populations (Alemu et al., 2020; Villanueva et al., 2021; Caballero et al., 2022), except for the adjusted inbreeding estimators 2, 4 and 9, which, not been previously described before. The application of these estimators to this population with a high number of generations, a particular and very structure linked to a specific mating design, represented an opportunity to evaluate the performance of these inbreeding coefficient estimators.

Expected homozygosity (FEXP) was calculated from the frequency of the reference allele (p) per generation as FEXP=12p1p. This coefficient was calculated considering both lines together (FEXP) and within the high variability line (FEXPH) and the low variability line (FEXPL).

Our own R code was used to calculate coefficients based on GRMs. PLINK v 1.9 (Chang et al., 2015) was used to detect ROH, and ENDOG v 4.8 (Gutierrez and Goyache, 2005) was used to calculate pedigree inbreeding coefficients.

2.2 Individual increase in inbreeding and realised effective population size

Effective population size was computed for each generation from individual increases in inbreeding (Gutiérrez et al., 2008; Gutiérrez et al., 2009) based on the different inbreeding coefficients described above. When pedigree inbreeding was used, two approaches were applied to calculate the individual increase in inbreeding. One is based on the classic formula: ΔFPEDt=FPEDtFPEDt11FPEDt1, where FPEDt is the inbreeding coefficient of the individual of generation t and FPEDt1 is the average coefficient of inbreeding of their parents (Falconer and Mackay, 1996). Then, the effective population size was computed as NePEDt=12FPED¯, where FPEDt¯ is the average ΔFPEDt of n individuals included in each generation (Falconer and Mackay, 1996).

Moreover, the individual increase in inbreeding coefficient using FPED was also calculated as ΔFPEDi=11FPEDit1 (Gutiérrez et al., 2009), where FPEDi is the individual coefficient of inbreeding of i and t are the number of generations. The individual increase in inbreeding has been proposed as a measure of standardized inbreeding rate per generation (Gonzále-Recio et al., 2007). The realised effective population was computed as NePEDi=12FPEDi¯, where FPEDi¯ is the average ΔFPEDi of n individuals included in each generation (Cervantes et al., 2008).

The molecular Ne based on FNEJ0 (NeNEJ0,FL&H0 NeL&H0), FVR1 (NeVR1, FVR2 (NeVR2), FYAN (NeYAN), and FROH0 (NeROH0 was calculated as Ne=12F, where F is the average ΔFi calculated from FNEJ0, FL&H0, FVR1, FVR2, FYAN, or FROH0 of n individuals included in each generation. ΔFi was calculated as: Fi=11Fit, where Fi is the individual FNEJ0, FL&H0, FROH0, FVR1, FVR2 or FYAN coefficient of inbreeding of individual i and t is the number of generations that passed since the first generation of selection where the reference population has been defined.

The expected Ne per generation was computed as NeEXP=12FEXP, where FEXP is the expected increase in homozygosity per generation calculated as FEXP=FEXPtFEXPt11FEXPt1, where FEXPt is the expected homozygosity of the current generation and FEXPt1 is the average coefficient of expected homozygosity of the previous generation. This expected Ne was obtained considering both lines together NeEXP, also within line for the high variability line NeEXPH and for the low variability line NeEXPH.

In addition, the number of generations that passed since the founder generation to the reference population could be derived from the classic formula Ft=11Ft, where Ft is the average inbreeding of the reference population and t is the generation of the reference population. Thus, the estimated generations that had passed from the founder population to the reference population are t=ln1Ftln1F. Different t values were calculated considering FNEJ (tNEJ) and FROH (tROH); tNEJ and tROH were calculated using the average inbreeding coefficients of the reference population, and the increase in inbreeding (F=12Ne) was derived from Ne based on the sex ratio (Nes), which was calculated as Nes=4MFM+F, where M is the number of breeding males and F is the number of breeding females (Wright, 1931). Nes was calculated using the permanent number of males (32) and females (64) in the randomly mated population, which was the origin of the selected divergent lines.

Moreover, t was also calculated using the average FPED of the reference population (tPED) and the average ∆FPEDi of whole selection generations. In this case, tPED represented the number of generations elapsed if the breeding method of the selection experiment would be applied from the founder population of the recorded pedigree to the first generation of the experiment (reference population).

Pearson’s correlation coefficients were calculated between all computed inbreeding coefficients and increases in inbreeding. Correlations between the different inbreeding coefficients were calculated for the whole population and different generation groups: initial (1, 2, 3, and 4), intermediate (16, 17, 18, and 19), and most recent generations (23, 24, 25, and 26). R software was used to calculate correlations using the function “cor·” (Rizopoulos, 2007).

3 Results

Figure 1 shows the trend of average inbreeding across generations calculated using molecular and pedigree approaches, which was positive in both cases. Both approaches presented similar performances across generations, except FROH and FNEJ, which showed the greatest amount of inbreeding in the first and last generations. In the first generation, FROH was 0.54 and FNEJ was 0.61, and in the last generation, FROH was 0.68 and FNEJ was 0.73. The lowest individual inbreeding coefficient was identified for FL&H in the first (−0.02) and last generation (0.28) (Supplementary Table S1). There were not differences in the evolution of the different inbreeding coefficients in the two selection lines when analysed separately (results are not shown). We therefore decided to analyse the data together.

FIGURE 1
www.frontiersin.org

FIGURE 1. Evolution of the different inbreeding coefficients computed from pedigree and molecular information across the selection generation. FPED: inbreeding coefficient computed from pedigree information. FNEJ: inbreeding coefficient computed from the diagonal of the Nejati Javaremi Matrix. FL&H: inbreeding coefficient computed from the diagonal of the Li and Horvitz Matrix. FVR1: inbreeding coefficient computed from the diagonal of the Van Raden type 1 matrix. FVR2: inbreeding coefficient computed from the diagonal of the Van Raden type 2 matrix. FYAN: inbreeding coefficient computed from the diagonal of the Yang matrix. FROH: inbreeding coefficient computed using runs of homozygosity. FPED0, FNEJ0, FL&H0, and FROH0 were adjusted for the mean inbreeding coefficients of the first selection generation.

The evolution of expected homozygosity and FNEJ for the whole population and selected lines for divergent birth weight variability across generations are represented in Figure 2. Expected homozygosity presented a positive trend when data from all individuals was used and when selected lines were analysed separately. Moreover, FEXPH and FEXPL were similar to FNEJH and FNEJL, with 0.75 FEXPH, 0.74 FEXPL, 0.73 FNEJH, and 0.72 FNEJL in the last generation. FEXP was greater than FNEJ until the seventh generation. However, after these generations, FEXP was markedly lower than FNEJ and lower than expected homozygosity within line because of the structure created by divergent selection. The highest values of FEXP and FNEJ were reached in the last selection generation, at 0.67 and 0.72, respectively.

FIGURE 2
www.frontiersin.org

FIGURE 2. Evolution of expected homozygosity (FEXP) and mean inbreeding coefficient obtained from the Nejati-Javaremi genomic relationship matrix (FNEJ) across the selection generations. FEXP-T was calculated using all individuals. FEXP-H was calculated using individuals from the high variability line. FEXP-L was calculated using individuals from the low variability line. FNEJ-T was calculated using all genotyped individuals. FNEJ-H was calculated using individuals from the high variability line. FNEJ-L was calculated using individuals from the low variability line.

Correlations between molecular and pedigree inbreeding coefficients using the whole population are shown in Table 2. All correlations were greater than 0.84. The correlations between FPED and molecular inbreeding coefficients were between 0.84 (FPED-FVR1) and 0.89 (FPED-FROH). The strongest correlations between molecular inbreeding coefficients were detected between FNEJ and FL&H (1.00). The FYAN showed correlations greater than 0.95 with all other molecular inbreeding coefficients, and the correlation with FPED was 0.88. When only a few generations were used to calculate the correlation coefficients, these correlations were lower between FPED and the molecular inbreeding coefficients (Supplementary Table S2). In the initial generations, FPED showed the highest correlation with FROH (0.25) and the lowest correlation with FVR1 (0.11). However, in the intermediate generations, FPED showed the highest correlations with FNEJ and FL&H (0.39), and the lowest correlation was observed with FVR2 (0.22). Moreover, in the most recent generations, the highest correlations were also observed between FPED-FNEJ and FPED-FL&H (0.22), and the lowest correlation was observed between FPED-FVR2 (0.13). Finally, in general, FYAN showed the highest correlations with other molecular inbreeding coefficients in the initial, intermediate, and final generations.

TABLE 2
www.frontiersin.org

TABLE 2. Correlation coefficients between molecular and pedigree inbreeding coefficients of all genotyped individuals. FPED: inbreeding coefficient computed from pedigree. FNEJ: inbreeding coefficient computed from the diagonal of the Nejati-Javaremi matrix. FL&H i: inbreeding coefficient computed from the diagonal of Li and Horvitz Matrix. FVR1: inbreeding coefficient computed from the diagonal of the Van Raden type 1 matrix. FVR2: inbreeding coefficient computed from the diagonal of the Van Raden type 2 matrix. FYAN: inbreeding coefficient computed from the diagonal of the Yang matrix. FROH: inbreeding coefficient computed using runs of homozygosity.

The evolution of individual inbreeding increases across generations calculated by pedigree and molecular inbreeding coefficients are shown in Figure 3. The ∆FPEDt (values between −0.03 and 0.06, Figure 3A) showed an irregular evolution across generations unlike ∆FPEDi (0.00 and 0.02, Figure 3B), which increased during fifteen generations and became stable thereafter. Individual increases in inbreeding calculated from molecular inbreeding coefficients showed a similar trend, with higher variability among individuals until the 10th generation, after which the trend was stabilised. The ∆FL&H0 (values between −0.13 and 0.13, Figure 3E) performed more like the other molecular increases in inbreeding (∆FNEJ0, ∆FVR1, ∆FVR2, ∆FYAN and ∆FROH0) than ∆FL&H (values between −0.02 and 0.02, Figure 3D) in the first five generations of selection.

FIGURE 3
www.frontiersin.org

FIGURE 3. Evolution across selection generations of molecular and pedigree increases in inbreeding coefficients (∆F) of all genotyped individuals. (A): ∆FPEDt derived from pedigree referenced against the inbreeding coefficient of the individual’s parents. (B): ∆FPEDi derived from pedigree referred to the number of generations of the individual. (C): ∆FNEJ0 derived from inbreeding coefficients (F) of the Nejati-Javaremi matrix (FNEJ) adjusted for the mean FNEJ of the first selection generation. (D): ∆FL&H derived from the diagonal of the Li and Horvitz matrix. (E): ∆FL&H0 derived from F of the Li and Horvitz matrix (FL&H) adjusted for the mean FL&H of the first selection generation. (F): ∆FVR1 derived from the diagonal of the Van Raden type 1 matrix. (G): ∆FVR2 derived from the diagonal of the Van Raden type 2 matrix. (H): ∆FYAN derived from the diagonal of the Yang matrix. (I): ∆FFROH0 derived from F of runs of homozygosity (FROH) adjusted for the mean FROH of the first selection generation.

Table 3 shows the correlations between increases in inbreeding calculated using all generations. Correlations greater than 0.90 were identified between ∆FL&H0—∆FNEJ0 (1.00), ∆FL&H—∆FYAN (0.92), ∆FROH0—∆FNEJ0 (0.91), and ∆FROH0—∆FL&H0 (0.91). Moreover, ∆FYAN was strongly correlated with ∆FVR1 (0.89) and ∆FVR2 (0.86). The lowest correlations were detected between ∆FPEDt and all other inbreeding increases being between 0.38 for ∆FPEDi and 0.18 for ∆FVR1 and ∆FVR2. Interestingly, higher correlations were detected between ∆FL&H and ∆FYAN, ∆FVR1 and ∆FVR2 than between ∆FL&H0 and ∆FYAN, ∆FVR1 and ∆FVR2. When only groups of generations (initial, intermediate and most recent generations) were considered, correlations between individual ∆F values were generally lower than when all generations were used (Supplementary Table S3). The strongest correlations were detected between ∆FL&H0 and ∆FNEJ0(1.00) in all generation groups. In summary, the molecular ∆F with the strongest correlations among the others in the different generation groups was ∆FYAN, followed by ∆FROH0. Nevertheless, when only a few generations were used to calculate ∆F, correlations between molecular and pedigree measures were very low.

TABLE 3
www.frontiersin.org

TABLE 3. Correlation coefficients between molecular and pedigree increases in inbreeding coefficients (∆F) of all genotyped individuals. ∆FPEDi was derived from pedigree referenced against the number of generations of the individual. ∆FPEDt was derived from the pedigree referring to the inbreeding coefficient of the individual’s parents. ∆FNEJ was derived from the diagonal of the Nejati-Javaremi matrix. ∆FNEJ0 was derived from inbreeding coefficients (F) obtained using the Nejati-Javaremi matrix (FNEJ) adjusted for the mean FNEJ of the first selection generation. ∆FL&H was derived from the diagonal of the Li and Horvitz matrix. ∆FL&H0 was derived from F using the Li and Horvitz matrix (FL&H) adjusted for the mean FL&H of the first selection generation. ∆FVR1 was derived from the diagonal of the Van Raden type 1 matrix. ∆FVR2 was derived from the diagonal of the Van Raden type 2 matrix. ∆FYAN was derived from the diagonal of the Yang matrix. ∆FROH was computed using runs of homozygosity. ∆FFROH0 was derived from F obtained using runs of homozygosity (FROH) adjusted for the mean FROH of the first selection generation.

Table 4 shows Ne calculated via different methodologies with molecular and pedigree data across the selection generations. The Ne trend stabilises after the 10th generation of selection except for NePEDt, which presented high variability across generations. In fact, NePEDt showed high values in generations 11 (1,482.27) and 26 (1,138.82) and a negative value in the second generation (−80.03). The NePEDi showed neither values as high as NePEDt nor negative values. The NePEDi showed the highest values in the first generations (94.72 in the first, 143.44 in the second, and 129.46 in the third). The lowest value of NePEDi was identified in the 25th generation (34.37). The mean and standard deviation across generations for NePEDt were 125.94 and 346.28, respectively, and 53.58 and 28.38, respectively, for NePEDi. Regarding molecular measurements, Ne performance differed across methodologies. Extremely negative values were detected with NeL&H in the first five generations of selection, the most negative being in the fifth generation (−4,081.60). These extremely negative values translated into a negative mean value (−228.68) with a high standard deviation (988.62). The NeNEJ0, NeL&H0, NeROH0, and NeYAN also presented negative values in the first and second generations. The NeVR2 had lower values compared with other molecular approaches, with a mean of 23.26 compared with 35.87 for NeNEJ0, 35.88 for NeL&H0, 32.66 for NeVR1, 40.76 for NeYAN and 30.48 for NeROH0. The standard deviation of NeVR1 was the lowest (6.11). The evolution of NeNEJ0 and NeL&H0 was almost identical. The evolution of NeNEJ0, NeL&H0, NeVR1, NeYAN, and NeROH0 as of the 10th generation was similar to that of NePEDi, while NeL&H0 showed greater values than NePEDi, and NeVR2 showed lower values than NePEDi. The evolution of NeEXP was more irregular than that of NeNEJ0, with extremely negative and positive values, such as 1852.44 in the 16th generation and −464.38 in the 26th generation (Supplementary Figure S1). When selected lines were analysed separately, NeEXPH presented a more regular evolution than NeEXPL. However, extremely negative and positive values were also observed in both (Supplementary Figure S1).

TABLE 4
www.frontiersin.org

TABLE 4. Evolution across selection generations, total mean, and standard deviation (SD) of the effective population size (Ne) using pedigree and molecular information. NePEDi derived from the individual increase in inbreeding (∆F) refers to the number of individual generations. NePEDt derived from ∆F based on the increase in pedigree inbreeding between two successive generations. NeNEJ0 derived from the inbreeding coefficient (F) of the Nejati-Javaremi genomic matrix (FNEJ) adjusted for the mean FNEJ of the first generation of selection. NeL&H. derived from ∆F of the Li and Horvitz genomic matrix. NeL&H0 derived from F of the Li and Horvitz genomic matrix (FL&H) adjusted for mean FL&H of the first generation of selection. NeVR1 derived from ∆F of the Van Raden type 1 genomic matrix. NeVR2 derived from ∆F of the Van Raden type 2 genomic matrix. NeYAN derived from ∆F of the Yang genomic matrix. NeROH0 derived from the F of runs of homozygosity (FROH) adjusted for mean FROH of the first generation of selection.

The Nes was 85.33, which was used to infer the number of generations that had passed from the founder population from the hypothetical foundation to the first generation of selection using FNEJ (tNEJ) and FROH(tROH) as described above: tNEJ was 162 generations, and tROH was 133 generations.

When pedigree data was used to calculate the number of generations elapsed if the breeding method of the selection experiment would be applied from the founder population of the recorded pedigree to the first generation of the experiment, tPED was two generations.

4 Discussion

This population, formed by two divergent lines selected for birth weight variability during twenty-six discrete generations, allowed us to test different inbreeding measures using molecular and pedigree data and their application to compute Ne.

The inbreeding coefficient has been defined in terms of correlations between parents uniting gametes (Wright, 1931) and as a probability that two genes randomly sampled in the parents’ gametes are IBD (Malecot, 1948). Because of FPED is completely defined in terms of IBD probability, the performance of the molecular inbreeding coefficients are going to be compared with FPED, although the use of FPED as the golden standard is not straightforward as assumes that the founders of the pedigree are unrelated (Keller et al., 2011) and provide the expected proportion of inbreeding while ignoring the effects of other forces such as Mendelian sampling, linkage, or natural selection against homozygous alleles (Wright, 1931). The FL&H, FVR1, FVR2, and FYAN behave more like correlations and have negative values in some individuals (Supplementary Table S1), as reported by other authors (Alemu et al., 2021; Villanueva et al., 2021). In the ROH definition, a long stretch of homozygosity is expected to be inherited from a common ancestor (Ceballos et al., 2018); FROH values range from 0 to 1 (Alemu et al., 2021) and present the highest correlations with FPED when all individuals from these populations are used. Therefore, among all the different measures of molecular inbreeding analysed in this population, FROH better fits the IBD definition when FPED is used as the reference coefficient. However, ROH analysis is an empirical, rule-based approach, and detection relies on the definition of several parameters that must be adjusted as a function of several constraints, such as marker density or the number of genotyping errors (Ferenčaković et al., 2013; Meyermans et al., 2020). Hence, ROH detection varied across different studies. Moreover, FROH is more representative of IBD than FNEJ, defined as the proportion of homozygous SNPs in the individual and, thus, unable to differentiate between IBD and identity by state (IBS) (Toro et al., 2014); however, FROH theoretically represents the proportion of long IBS regions in the genome. Furthermore, if FROH and FNEJ were used to infer the number of generations that had passed since the founder population, tNEJ would be expected to present the highest value due to FNEJ being fully IBS, while FROH is expected to be more representative of IBD. The high values of tPED and tROH reflected the history of the founder population originating from the genetic contribution of three inbred mouse lines: BALB/c, C57BL and CBA. This mixed starting population was maintained under panmixia for 40 generations before the start of the selection experiment (Formoso-Rafferty et al., 2016b). The lowest correlations between FPED and molecular inbreeding coefficients were detected with FVR1 and FVR2. In addition, other authors (Caballero et al., 2022) reported that the poorest performance with regard to IBD was detected for FVR2.

The number of generations that had passed before starting the experiment inferred from pedigree (tPED) was very low. Pedigrees, however, are usually limited to a few past generations and ignore some of the ancestral relationships (Keller et al., 2011). In this case, FPED was calculated using five generations before initiating the experiment, and the mating method was changed at this start; therefore, tPED was lower than the number of generations that had passed since the founder population in the pedigree.

When all data was used, correlations between different inbreeding measures were high. In other studies in Iberian pigs (Saura et al., 2015), cattle (Alemu et al., 2021; Lozada-Soto et al., 2022) or humans (McQuillan et al., 2012), strong correlations were also detected. However, when a low number of generations was used, the correlations between different pedigrees and genomic inbreeding measures were lower. Therefore, the high correlations observed when using all animals across generations were probably because of the large number of genotypes distributed across generations, resulting in a wide range of inbreeding values in all coefficients. This demonstrates that this experimental population was of particular interest to check the relationship between the different parameters assayed. Moreover, other authors reported that the incompleteness of the pedigree should lead to poor correlations (Schiavo et al., 2022). Furthermore, the correlations between the different increases in inbreeding were lower than between the inbreeding coefficients because, firstly, ∆FPEDi and ∆F molecular are standardised by the depth of pedigree information. Therefore, the evolution of ∆F was expected to be stable if the mating design remains stable across generations, generating less variation than in the case of inbreeding coefficients and lower correlations. In addition, ΔFPEDi was also expected to correlate better with the other coefficients than ∆FPEDt, because the estimating molecular ∆Fs were designed following ∆FPEDi. We could not fit the ∆FPEDt formula to the molecular information because we did not genotype the trios.

If the frequencies of the base population are known, the performances of FL&H, FVR1, FVR2, and FYAN are in terms of IBD (Caballero et al., 2022). However, other authors have reported that FL&H, FVR1, FVR2, and FYAN presented values outside of Malecot and Wright’s inbreeding definitions and that FNEJ better fitted these definitions, as its values ranged between 0 and 1 (Saura et al., 2015; Villanueva et al., 2021). Moreover, Villanueva et al. (2021) observed that only FL&H could be interpreted in terms of loss or gain of genetic variability, while FVR1, FVR2, and FYAN presented some inconsistencies concerning loss or gain of genetic variability. Regarding the application of inbreeding measures, FROH and FPED should be preferred to measure the increase in whole genome homozygosity, and FROH should be preferred over FPED when pedigrees are not deep enough or present many errors (Keller et al., 2011; Alemu et al., 2021).

It has been reported that FYAN and FVR2, giving higher weight to rare alleles (Villanueva et al., 2021), presented the highest correlation with homozygosity at SNPs with moderate to high MAFs (Alemu et al., 2021). However, alleles with low MAFs are more representative of kinship because of the higher probability of being transmitted when individuals belong to the same family. Hence, FPED and FROH. did not reflect the segregation of these low-frequency alleles that could be an indicator of kinship. Therefore, to measure global inbreeding, FPED and FROH were more useful (Alemu et al., 2021; Caballero et al., 2022), and FYAN or FVR2 could be better for measuring population kinship.

The FL&H represents a modification of FNEJ adjusted for the expected homozygosity of the reference population. Hence, these inbreeding coefficients have negative values in some individuals (Supplementary Table S3). For FL&H, this is clearly represented in Figure 1, where expected homozygosity was greater than FNEJ during the first ten generations. This was reflected as negative values of the individual increases in inbreeding and high negative and positive NeL&H values during the first eleven generations. The FL&H is designed using the allele frequencies of a reference population to express homozygosity in terms of IBD, such as FVR1, FVR2, and FYAN. However, ∆FL&H presented more variability during the first five generations of selection (Figure 3) than FVR1, FVR2, and FYAN. This was probably because FL&H had a higher influence on the change in the mating design during the first selection generations and was more dependent on Hardy-Weinberg disequilibrium. This problem was solved when FL&H was adjusted for the mean inbreeding coefficient of the first generation of selection (NeL&H0). Arias et al. (2023) described that the adjustment of the inbreeding coefficients by the mean of the reference population did not allow a complete correction of the estimates, however, the use of the allele frequencies of the RP to estimate the IBD evolution in the population studied hides some assumptions, as the reference population is in Hardy-Weinberg equilibrium or that the RP have no molecular kindship, as it is assumed to be the real founder population. Therefore, adjusting the molecular inbreeding estimates for the mean inbreeding of the RP brings these estimates closer to the IBD definition. In this context, FROH and FNEJ presented the advantage of not being influenced by the allele frequencies of the RP in their estimation and of being easily corrected by the mean FROH and FNEJ of the RP. Furthermore, the first two generations presented negative values for NeNEJ0, NeL&H0, and NeROH0. due to a negative individual increase in inbreeding produced by a negative trend of inbreeding coefficients in the first two generations, which was probably produced by the abrupt change in mating design. Thus, this negative trend could also be because in the five previous generations before the experiment, individuals were randomly mated, and from the first generation of selection, matings were designed to avoid sharing grandparents to avoid inbreeding. In addition, none of the effective population sizes started to present stable values until the 10th generation. Toro et al. (2020) had previously pointed out that when optimal management of genetic diversity is implemented in a population, molecular estimates of Ne could be meaningless because the increase in genetic diversity resulted in negative Ne. In addition, none of the effective population sizes began to present stable values until the 10th generation. Therefore, in this population, ten generations of selection were needed to reflect the change in mating method by the effective population size, as expected, since Ne is a diversity indicator strongly influenced by the change in mating policy, as seen in the results of this study (Table 3).

The NeNEJ0 and NeL&H0 presented almost identical values in this population and similar values to NePEDi. NeNEJ0 presented the advantage of not using the frequencies of a reference population such as NeFROH0 and, therefore, are less dependent on allele frequencies.

Comparing NePEDi and NePEDt, the latter presented a higher standard deviation, as other authors reported that NePEDt presented high fluctuations over time due to the high impact of breeding method changes, errors in pedigree registration, or sampling effects (Gutiérrez et al., 2008; Gutiérrez et al., 2009; Cervantes et al., 2011; Leroy et al., 2013) and could lead to negative values (Groeneveld et al., 2009; Leroy et al., 2009; Nagy et al., 2010). The NePEDi values were high during the 10th–12th first generations. After that, the evolution of NePEDi tended to be stable and reflected the history of mating. Moreover, no negative values were observed, and the standard deviation was lower. Furthermore, NePEDi evolved to be similar to the number of males of females mated in the selection experiment, accounting for the intensity of selection that could be influencing the evolution of homozygosity in this population, when the weight of the random mating previous to the experiment decreased. Another study performed in French Angora rabbits divergently selected for total fleece weight, showed that the low line had high Ne values calculated from pedigree data. However, this line always showed lower levels of inbreeding than the high line over the years of selection, even though the mating policy was the same in both lines (Rafat et al., 2009). This did not occur in our selection lines, both of which had almost the same level of inbreeding since the beginning of the experiment (results not shown).

The minimum viable population size thresholds have been traditionally defined as 50 for the short term and 500 for the long term (Harmon and Braude, 2010). All Ne values of the last generations obtained for this population were lower than the defined thresholds. Low Ne levels relate to a decrease in genetic variability, allele fixation, and a reduction in the selection response (Domínguez-Viveros et al., 2020). Nevertheless, the use of a strict critical level for Ne is not straightforward, as different factors could affect its interpretation, such as the method used, the species, or population structure (Leroy et al., 2013; Mokhtari et al., 2015). Moreover, in other studies, effective population sizes below 50 were reported in some species that did not present viability problems (Leroy et al., 2013). Therefore, this threshold should be revised and adapted to the particularity of each population.

5 Conclusion

The results obtained in these studies allowed to understand better the performance of different inbreeding coefficients and effective population sizes applied in a population with many discrete generations and under strong selection. Moreover, these results could be used as a reference for the study of inbreeding and genetic diversity in other populations. The FROH and FPED presented the strongest correlations and were more representative of inbreeding in terms of IBD. However, other measures such as FYAN were useful in terms of rare alleles, which are more likely inherited from common ancestors. Adjusting the inbreeding estimates to the mean inbreeding of the reference population allowed it to fit more closely to IBD. The NePEDi presented a more reliable performance than NePEDt. In general, Ne obtained by molecular data showed a similar trend when comparing molecular approaches to each other and NePEDi, being more similar to NeL&H0 and NeNEJ0. However, data from several generations was necessary to reach a stable trend for Ne, both with pedigree and molecular data.

Data availability statement

The data presented in the study are deposited in the Figshare repository, accession link: https://figshare.com/s/632fbaefccc501da66e1.

Ethics statement

The animal study was approved by Animal Experimentation Committee (PROEX 224/18). The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

CO-M: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Software, Visualization, Writing–original draft, Writing–review and editing. IC: Conceptualization, Investigation, Project administration, Resources, Supervision, Validation, Writing–review and editing. NF-R: Data curation, Supervision, Writing–review and editing. JPG: Conceptualization, Formal Analysis, Investigation, Methodology, Software, Supervision, Writing–review and editing.

Funding

The authors declare financial support was received for the research, authorship, and/or publication of this article. This study was funded by a grant from the Ministry of Science, Innovation, and Universities; grant number: PGC 2018-096198-A-I00.

Acknowledgments

The genotyping service was carried out at CEGEN-PRB3-ISCIII; it is supported by grant PT17/0019 of the PE I+D+i 2013-2016 and funded by ISCII and ERDF.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2023.1303748/full#supplementary-material

References

Alemu, S. W., Kadri, N. K., Harland, C., Faux, P., CharlierCaballero, A., Druet, T., et al. (2021). An evaluation of inbreeding measures using a whole-genome sequenced cattle pedigree. Hered. (Edinb). 126, 410–423. doi:10.1038/s41437-020-00383-9

CrossRef Full Text | Google Scholar

Arias, K. D., Gutiérrez, J. P., Fernández, I., Álvarez, I., and Goyache, F. (2023). Approaching autozygosity in a small pedigree of Gochu Asturcelta pigs. Genet. Sel. Evol. 25, 74. doi:10.1186/s12711-023-00846-7

CrossRef Full Text | Google Scholar

Caballero, A., Fernández, A., Villanueva, B., and Toro, M. A. (2022). A comparison of marker-based estimators of inbreeding and inbreeding depression. Genet. Sel. Evol. 54 (1), 82. doi:10.1186/s12711-022-00772-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Ceballos, F. C., Joshi, P. K., Clark, D. W., Ramsay, M., and Wilson, J. F. (2018). Runs of homozygosity: windows into population history and trait architecture. Nat. Rev. Genet. 19 (4), 220–234. doi:10.1038/nrg.2017.109

PubMed Abstract | CrossRef Full Text | Google Scholar

Cervantes, I., Goyache, F., Molina, A., Valera, M., and Gutiérrez, J. P. (2008). Application of individual increase in inbreeding to estimate realized effective sizes from real pedigrees. J. Anim. Breed. Genet. 125 (5), 301–310. doi:10.1111/j.1439-0388.2008.00755.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Cervantes, I., Goyache, F., Molina, A., Valera, M., and Gutiérrez, J. P. (2011). Estimation of effective population size from the rate of coancestry in pedigreed populations. J. Anim. Breed. Genet. 128 (1), 56–63. doi:10.1111/j.1439-0388.2010.00881.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Chang, C. C., Chow, C. C., Tellier, L. C., Vattikuti, S., Purcell, S. M., and Lee, J. J. (2015). Second-generation PLINK: rising to the challenge of larger and richer datasets. GigaScience 4, 7. doi:10.1186/s13742-015-0047-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Domínguez-Viveros, J., Rodríguez-Almeida, F. A., Medellín-Cázares, A., and Gutiérrez, J. P. (2020). Análisis del pedigrí en diez poblaciones mexicanas de ovinos. Rev. Mex. Cienc. Pecu. 11 (4), 1071–1086. doi:10.22319/rmcp.v11i4.5457

CrossRef Full Text | Google Scholar

Druet, T., and Gautier, M. (2017). A model-based approach to characterize individual inbreeding at both global and local genomic scales. Mol. Ecol. 26 (20), 5820–5841. doi:10.1111/mec.14324

PubMed Abstract | CrossRef Full Text | Google Scholar

Falconer, D. S., and MacKay, T. F. C. (1996). Introduction to quantitative genetics. Harlow, England: Prentice Hall.

Google Scholar

Ferenčaković, M., Sölkner, J., and Curik, I. (2013). Estimating autozygosity from high-throughput information: effects of SNP density and genotyping errors. Genet. Sel. Evol. 45 (1), 42. doi:10.1186/1297-9686-45-42

PubMed Abstract | CrossRef Full Text | Google Scholar

Formoso-Rafferty, N., Cervantes, I., Ibáñez-Escriche, N., and Gutiérrez, J. P. (2016a). Correlated genetic trends for production and welfare traits in a mouse population divergently selected for birth weight environmental variability. Animal 10, 1770–1777. doi:10.1017/S1751731116000860

PubMed Abstract | CrossRef Full Text | Google Scholar

Formoso-Rafferty, N., Cervantes, I., Ibáñez-Escriche, N., and Gutiérrez, J. P. (2016b). Genetic control of the environmental variance for birth weight in seven generations of a divergent selection experiment in mice. J. Anim. Breed. Genet. 133 (3), 227–237. doi:10.1111/jbg.12174

PubMed Abstract | CrossRef Full Text | Google Scholar

Formoso-Rafferty, N., Cervantes, I., Ibáñez-Escriche, N., and Gutiérrez, J. P. (2017). Modulating birth weight heritability in mice. J. Anim. Sci. 95, 531–537. doi:10.2527/jas2016.1169

PubMed Abstract | CrossRef Full Text | Google Scholar

Formoso-Rafferty, N., Cervantes, I., Sánchez, J. P., Gutiérrez, J. P., and Bodin, L. (2019). Effect of feed restriction on the environmental variability of birth weight in divergently selected lines of mice. Genet. Sel. Evol. 51 (1), 27. doi:10.1186/s12711-019-0471-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Formoso-Rafferty, N., Chavez, K. N., Ojeda, C., Cervantes, I., and Gutiérrez, J. P. (2020). Selection response in a divergent selection experiment for birth weight variability in mice compared with a control line. Anim. (Basel) 10 (6), 920. doi:10.3390/ani10060920

CrossRef Full Text | Google Scholar

Formoso-Rafferty, N., El-Ouazizi El-Kahia, L., Arias-Álvarez, M., Gutiérrez, J. P., and Cervantes, I. (2023). Embryo survival and fertility differ in lines divergently selected for birth weight homogeneity in mice. J. Anim. Breed. Genet. 140, 549–557. doi:10.1111/jbg.12778

PubMed Abstract | CrossRef Full Text | Google Scholar

Formoso-Rafferty, N., Gutiérrez, J. P., García-Álvarez, A., Pérez, T., and Cervantes, I. (2022). Impact of selection for birth weight variability on reproductive longevity: a mice model. J. Anim. Breed. Genet. 139 (4), 370–379. doi:10.1111/jbg.12676

PubMed Abstract | CrossRef Full Text | Google Scholar

Frankham, R., Ballou, J. D., and Briscoe, D. A. (2010). Introduction to conservation genetics. New York: Cambridge University Press.

Google Scholar

González-Recio, O., and López de MaturanaGutiérrez, E. J. P. (2007). Inbreeding depression on female fertility and calving ease in Spanish dairy cattle. J. Dairy. Sci. 90, 5744–5752. doi:10.3168/jds.2007-0203

PubMed Abstract | CrossRef Full Text | Google Scholar

Groeneveld, E., Westhuizen, B. V. D., Maiwashe, A., Voordewind, F., and Ferraz, J. B. S. (2009). POPREP: a generic report for population management. Genet. Mol. Res. 8 (3), 1158–1178. doi:10.4238/vol8-3gmr648

PubMed Abstract | CrossRef Full Text | Google Scholar

Gutiérrez, J. P., Cervantes, I., and Goyache, F. (2009). Improving the estimation of realized effective population sizes in farm animals. J. Anim. Breed. Genet. 126 (4), 327–332. doi:10.1111/j.1439-0388.2009.00810.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Gutiérrez, J. P., Cervantes, I., Molina, A., Valera, M., and Goyache, F. (2008). Individual increase in inbreeding allows estimating effective sizes from pedigrees. Genet. Sel. Evol. 40 (4), 359–378. doi:10.1186/1297-9686-40-4-359

PubMed Abstract | CrossRef Full Text | Google Scholar

Gutierrez, J. P., and Goyache, F. (2005). A note on ENDOG: a computer program for analysing pedigree information. J. Anim. Breed. Genet. 122 (3), 172–176. doi:10.1111/j.1439-0388.2005.00512.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Harmon, L., and Braude, S. (2010). “Conservation of small populations: effective population sizes, inbreeding, and the 50/500 rule,” in An introduction to methods and models in ecology, evolution, and conservation biology. Editors S. Braude, and S. B. Low (Princeton, New Jersey, USA: Princeton University Press), 125–138.

CrossRef Full Text | Google Scholar

Hedrick, P. W. (2012). What is the evidence for heterozygote advantage selection? Trends. Ecol. Evol. 27 (12), 698–704. doi:10.1016/j.tree.2012.08.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Hill, W. G. (1981). Estimation of effective population size from data on linkage disequilibrium. Genet. Res. 38, 209–216. doi:10.1017/s0016672300020553

CrossRef Full Text | Google Scholar

Howard, J. T., Pryce, J. E., Baes, C., and Maltecca, C. (2017). Invited review: inbreeding in the genomics era: inbreeding, inbreeding depression, and management of genomic variability. J. Dairy. Sci. 100 (8), 6009–6024. doi:10.3168/jds.2017-12787

PubMed Abstract | CrossRef Full Text | Google Scholar

Keller, M. C., Visscher, P. M., and Goddard, M. E. (2011). Quantification of inbreeding due to distant ancestors and its detection using dense single nucleotide polymorphism data. Genetics 189 (1), 237–249. doi:10.1534/genetics.111.130922

PubMed Abstract | CrossRef Full Text | Google Scholar

Krimbas, C. B., and Tsakas, S. (1971). The genetics of Dacus oleae V. Changes of esterase polymorphism in a naturas population following insecticide control: selection or drift? Evolution 25, 565–578. doi:10.1111/j.1558-5646.1971.tb01904.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Leroy, G., Mary-Huard, T., Verrier, E., Danvy, S., Charvolin, E., and Danchin-Burge, C. (2013). Methods to estimate effective population size using pedigree data: examples in dog, sheep, cattle and horse. Genet. Sel. Evol. 45 (1), 1. doi:10.1186/1297-9686-45-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Leroy, G., Verrier, E., Meriaux, J. C., and Rognon, X. (2009). Genetic diversity of dog breeds: within-breed diversity comparing genealogical and molecular data. Anim. Genet. 40 (3), 323–332. doi:10.1111/j.1365-2052.2008.01842.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Leutenegger, A. L., Prum, B., Génin, E., Verny, C., Lemainque, A., Clerget-Darpoux, F., et al. (2003). Estimation of the inbreeding coefficient through use of genomic data. Am. J. Hum. Genet. 73 (3), 516–523. doi:10.1086/378207

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, C. C., and Horvitz, D. G. (1953). Some methods of estimating the inbreeding coefficient. Am. J. Hum. Genet. 5 (2), 107–117.

PubMed Abstract | Google Scholar

Lozada-Soto, E. A., Tiezzi, F., Jiang, J., Cole, J. B., VanRaden, P. M., and Maltecca, C. (2022). Genomic characterization of autozygosity and recent inbreeding trends in all major breeds of US dairy cattle. J. Dairy. Sci. 105 (11), 8956–8971. doi:10.3168/jds.2022-22116

PubMed Abstract | CrossRef Full Text | Google Scholar

Malecot, G. (1948). Malécot G. Les mathématiques de l’hérédité. Paris: Masson et Cie.

Google Scholar

McQuillan, R., Eklund, N., Pirastu, N., Kuningas, M., McEvoy, B. P., Esko, T., et al. (2012). Evidence of inbreeding depression on human height. PLOS Genet. 8 (7), e1002655. doi:10.1371/journal.pgen.1002655

PubMed Abstract | CrossRef Full Text | Google Scholar

Meuwissen, T., and Luo, Z. (1992). Computing inbreeding coefficients in large populations. Genet. Sel. Evol. 24 (4), 305–313. doi:10.1186/1297-9686-24-4-305

CrossRef Full Text | Google Scholar

Meyermans, R., Gorssen, W., Buys, N., and Janssens, S. (2020). How to study runs of homozygosity using PLINK? A guide for analyzing medium density SNP data in livestock and pet species. BMC Genomics 21 (1), 94. doi:10.1186/s12864-020-6463-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Mokhtari, M. S., Miraei-Ashtiani, S. R., Jafaroghli, M., and Gutiérrez, J. P. (2015). Studying genetic diversity in moghani sheep using pedigree analysis. J. Agr. Sci. Tech. 17, 1151–1160.

Google Scholar

Nagy, I., Curik, I., Radnai, I., Cervantes, I., Gyovai, P., Baumung, R., et al. (2010). Genetic diversity and population structure of the synthetic Pannon White rabbit revealed by pedigree analyses. J. Anim. Sci. 88 (4), 1267–1275. doi:10.2527/jas.2009-2273

PubMed Abstract | CrossRef Full Text | Google Scholar

Nei, M., and Tajima, F. (1981). Genetic drift and estimation of effective population size. Genetics 98 (3), 625–640. doi:10.1093/genetics/98.3.625

PubMed Abstract | CrossRef Full Text | Google Scholar

Nejati-Javaremi, A., Smith, C., and Gibson, J. P. (1997). Effect of total allelic relationship on accuracy of evaluation and response to selection. J. Anim. Sci. 75 (7), 1738–1745. doi:10.2527/1997.7571738x

PubMed Abstract | CrossRef Full Text | Google Scholar

Peripolli, E., Munari, D. P., Silva, M. V. G. B., Lima, A. L. F., Irgang, R., and Baldi, F. (2017). Runs of homozygosity: current knowledge and applications in livestock. Anim. Genet. 48 (3), 255–271. doi:10.1111/age.12526

PubMed Abstract | CrossRef Full Text | Google Scholar

Pollak, E. (1983). A new method for estimating the effective population size from allele frequency changes. Genetics 104 (3), 531–548. doi:10.1093/genetics/104.3.531

PubMed Abstract | CrossRef Full Text | Google Scholar

Pudovkin, A. I., Zaykin, D. V., and Hedgecock, D. (1996). On the potential for estimating the effective number of breeders from heterozygote-excess in progeny. Genetics 144 (1), 383–387. doi:10.1093/genetics/144.1.383

PubMed Abstract | CrossRef Full Text | Google Scholar

Rafat, S. A., Allain, D., and de Rochambeau and, H. (2009). Genetic description of a divergent selection experiment in Angora rabbits with overlapping generations. J. Anim. Breed. Genet. 126 (3), 189–197. doi:10.1111/j.1439-0388.2008.00769.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Rizopoulos, D. (2007). Ltm: AnRPackage for latent variable modeling and item response theory analyses. J. Stat. Soft 17, 1–25. doi:10.18637/jss.v017.i05

CrossRef Full Text | Google Scholar

Rodríguez-Ramilo, S. T., Elsen, J. M., and Legarra, A. (2019). Inbreeding and effective population size in French dairy sheep: comparison between genomic and pedigree estimates. J. Dairy. Sci. 102 (5), 4227–4237. doi:10.3168/jds.2018-15405

PubMed Abstract | CrossRef Full Text | Google Scholar

Roff, D. A. (1997). Evolutionary quantitative genetics. New York: Chapman & Hall.

Google Scholar

Saura, M., Fernández, A., Varona, L., Fernández, A. I., de Cara, M., Barragán, C., et al. (2015). Detecting inbreeding depression for reproductive traits in Iberian pigs using genome-wide data. Genet. Sel. Evol. 47, 1. doi:10.1186/s12711-014-0081-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Schiavo, G., Bovo, S., Ribani, A., Moscatelli, G., Bonacini, M., Prandi, M., et al. (2022). Comparative analysis of inbreeding parameters and runs of homozygosity islands in 2 Italian autochthonous cattle breeds mainly raised in the Parmigiano-Reggiano cheese production region. J. Dairy. Sci. 105 (3), 2408–2425. doi:10.3168/jds.2021-20915

PubMed Abstract | CrossRef Full Text | Google Scholar

Sjödin, P., Kaj, I., Krone, S., Lascoux, M., and Nordborg, M. (2005). On the meaning and existence of an effective population size. Genetics 169 (2), 1061–1070. doi:10.1534/genetics.104.026799

PubMed Abstract | CrossRef Full Text | Google Scholar

Tallmon, D. A., Koyuk, A., Luikart, G., and Beaumont, M. A. (2008). COMPUTER PROGRAMS: onesamp: a program to estimate effective population size using approximate Bayesian computation. Mol. Ecol. Res. 8 (2), 299–301. doi:10.1111/j.1471-8286.2007.01997.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Toro, M. A., Villanueva, B., and Fernández, J. (2014). Genomics applied to management strategies in conservation programmes. Livest. Sci. 166, 48–53. doi:10.1016/j.livsci.2014.04.020

CrossRef Full Text | Google Scholar

Toro, M. A., Villanueva, B., and Fernández, J. (2020). The concept of effective population size loses its meaning in the context of optimal management of diversity using molecular markers. J. Anim. Breed. Genet. 137 (4), 345–355. doi:10.1111/jbg.12455

PubMed Abstract | CrossRef Full Text | Google Scholar

VanRaden, P. M. (2008). Efficient methods to compute genomic predictions. J. Dairy. Sci. 91 (11), 4414–4423. doi:10.3168/jds.2007-0980

PubMed Abstract | CrossRef Full Text | Google Scholar

VanRaden, P. M., Olson, K. M., Wiggans, G. R., Cole, J. B., and Tooker, M. E. (2011). Genomic inbreeding and relationships among Holsteins, Jerseys, and Brown Swiss. J. Dairy. Sci. 94 (11), 5673–5682. doi:10.3168/jds.2011-4500

PubMed Abstract | CrossRef Full Text | Google Scholar

Villanueva, B., Fernández, A., Saura, M., Caballero, A., Fernández, J., Morales-González, E., et al. (2021). The value of genomic relationship matrices to estimate levels of inbreeding. Genet. Sel. Evol. 53 (1), 42. doi:10.1186/s12711-021-00635-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Santiago, E., and Caballero, A. (2016). Prediction and estimation of effective population size. Heredity 117 (4), 193–206. doi:10.1038/hdy.2016.43

PubMed Abstract | CrossRef Full Text | Google Scholar

Wright, S. (1931). Evolution in mendelian populations. Genetics 16 (2), 97–159. doi:10.1093/genetics/16.2.97

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J., Benyamin, B., McEvoy, B. P., Gordon, S., Henders, A. K., Nyholt, D. R., et al. (2010). Common SNPs explain a large proportion of the heritability for human height. Nat. Genet. 42 (7), 565–569. doi:10.1038/ng.608

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: genomic inbreeding, effective population size, divergent selection, birth weight environmental variability, increase in inbreeding

Citation: Ojeda-Marín C, Cervantes I, Formoso-Rafferty N and Gutiérrez JP (2023) Genomic inbreeding measures applied to a population of mice divergently selected for birth weight environmental variance. Front. Genet. 14:1303748. doi: 10.3389/fgene.2023.1303748

Received: 28 September 2023; Accepted: 27 November 2023;
Published: 14 December 2023.

Edited by:

Luis Varona, University of Zaragoza, Spain

Reviewed by:

Ricardo Pong-Wong, University of Edinburgh, United Kingdom
Thiruvenkadan Aranganoor Kannan, Tamil Nadu Veterinary and Animal Sciences University, India

Copyright © 2023 Ojeda-Marín, Cervantes, Formoso-Rafferty and Gutiérrez. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Candela Ojeda-Marín, Y2FuZGVsYW9AdWNtLmVz

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.