Genetic inbreeding load and its individual prediction for milk yield in French dairy sheep

The magnitude of inbreeding depression depends on the recessive burden of the individual, which can be traced back to the hidden (recessive) inbreeding load among ancestors. However, these ancestors carry different alleles at potentially deleterious loci and therefore there is individual variability of this inbreeding load. Estimation of the additive genetic value for inbreeding load is possible using a decomposition of inbreeding in partial inbreeding components due to ancestors. Both the magnitude of variation in partial inbreeding components and the additive genetic variance of inbreeding loads are largely unknown. Our study had three objectives. First, based on substitution effect under non-random matings, we showed analytically that inbreeding load of an ancestor can be expressed as an additive genetic effect. Second, we analysed the structure of individual inbreeding by examining the contributions of specific ancestors/founders using the concept of partial inbreeding coefficients in three French dairy sheep populations (Basco-Béarnaise, Manech Tête Noire and Manech Tête Rousse). Third, we included these coefficients in a mixed model as random regression covariates, to predict genetic variance and breeding values of the inbreeding load for milk yield in the same breeds.

Pedigrees included 190,276, 166,028 and 633,655 animals of Basco-Béarnaise, Manech Tête Noire and Manech Tête Rousse, respectively, born between 1985 and 2021. A fraction of 99.1% of the partial inbreeding coefficients were lower than 0.01 in all breeds, meaning that in practice inbreeding occurs in pedigree loops that span several generations backwards. Less than 5% ancestors generate inbreeding, because mating is essentially between unrelated individuals. Inbreeding load estimations involved 658,731, 541,180 and 2,168,454 records of yearly milk yield from 178,123, 151,863 and 596,586 females in Basco-Béarnaise, Manech Tête Noire and Manech Tête Rousse, respectively. Adding the inbreeding load effect to the model improved the fitting (values of the statistic Likelihood Ratio Test between 132 and 383) for milk yield in the three breeds. The inbreeding load variances were equal to 11,804 and 9435 L squared of milk yield for a fully inbred (100%) descendant in Manech Tête Noire and Manech Tête Rousse. In Basco-Béarnaise, the estimate of the inbreeding load variance (11,804) was not significantly different from zero. The correlations between (direct effect) additive genetic and inbreeding load effects were - 0.09, - 0.08 and - 0.12 in Basco-Béarnaise, Manech Tête Noire and Manech Tête Rousse.

The decomposition of inbreeding in partial coefficients in these populations shows that inbreeding is mostly due to several small contributions of ancestors (lower than 0.001) going back several generations (5 to 7 generations), which is according to the policy of avoiding close matings. There is variation of inbreeding load among animals, although its magnitude does not seem enough to warrant selection based on this criterion.