On individual-specific prediction of hidden inbreeding depression load

Inbreeding depression and runs of homozygosity islands in Asturiana de los Valles cattle breed after 30 years of selection

Inbreeding depression results in a decrease in the average phenotypic values of affected traits. It has been traditionally estimated from pedigree-based inbreeding coefficients. However, with the development of single-nucleotide polymorphism arrays, novel methods were developed for calculating the inbreeding coefficient, and consequently, inbreeding depression. The aim of the study was to analyse inbreeding depression in 6 growth and 2 reproductive traits in the Asturiana de los Valles cattle breed using both genealogical and molecular information. The pedigree group comprised 225,848 records and an average equivalent number of complete generations of 2.3. The molecular data comprised genotypes of 2693 animals using the Affymetrix medium-density chip. Using the pedigree information, three different inbreeding coefficients were estimated for the genotyped animals: the full pedigree coefficient (FPED), and the recent and ancient inbreeding coefficients based on the information of the last three generations (FPED<3G) and until the last three generations (FPED>3G), respectively. Using the molecular data, seven inbreeding coefficients were calculated. Four of them were estimated based on runs of homozygosity (ROH), considering (1) the total length (FROH), (2) segments shorter than 4 megabases (FROH<4), (3) between 4 and 17 megabases (FROH4-17), and (4) longer than 17 Mb (FROH>17). Additionally, the three inbreeding coefficients implemented in the Plink software (FHAT1-3) were estimated. Inbreeding depression was estimated using linear mixed-effects model with inbreeding coefficients used as covariates. All analysed traits (birth weight, preweaning average daily gain, weaning weight adjusted at 180 days, carcass weight, calving ease, age at first calving, calving interval) showed a statistically significant non-zero effect of inbreeding depression estimated from the pedigree group, except for the Postweaning Average Daily Gain trait. When inbreeding coefficients were based on the genomic group, statistically significant inbreeding depression was observed for two traits, Preweaning Average Daily Gain and Weaning Weight based on FROH, FROH>17, and FHAT3 inbreeding coefficients. Nevertheless, similar to inbreeding depression estimated based on pedigree information, estimates of inbreeding depression based on genomic information had no relevant economic impact. Despite this, from a long-term perspective, genotyped data could be included to maximize genetic progress in genetic programs following an optimal genetic contribution strategy and to consider individual inbreeding load instead global inbreeding. ROH islands were identified on chromosomes 2, 3, 8, 10, and 16. Such regions contain several candidate genes for growth development, intramuscular fat, body weight and lipid metabolism that are related to production traits selected in Asturiana de los Valles breed.

Genetic inbreeding depression load for fertility traits in Pura Raza Española mares

Fertility is a key factor in the economic success of horse farms. However, it has received little attention due to the difficulty of measuring fertility objectively. Since its studbook creation (1912), the Pura Raza Española (PRE) breed has been a closed population and become high in-bred resulting in inbreeding depression (poor phenotypic values). Nevertheless, heterogeneous effects of inbreeding depression have been detected among founders and nonfounders. The aims of this study were (1) to analyze the genetic parameters for reproductive traits in mares of the PRE horse breed and (2) to estimate, for the first time, the inbreeding depression load associated with common ancestors of the breed. A total of 22,799 mares were analyzed. Heritability estimates ranged from 0.05 (interval between first and second foaling) to 0.16 (age at first foaling), whereas inbreeding depression load ratios ranged from 0.06 (parturition efficiency at 6th foaling) to 0.17 (age at first foaling), for a partial inbreeding coefficient of 10%. Although heritability is related to the variability expressed in the population, inbreeding depression load ratios measure the potential variability, whether expressed in the population or not. Most correlations between additive and inbreeding depression load genetic values were significant (P < 0.001) and of low to moderate magnitude. Our results confirm that individual inbreeding depression loads allow us to select horses that have a genetic value resistant to the deleterious effects of inbreeding.

Bayesian inference of the inbreeding load variance for fertility traits in Brown Swiss cattle

Our study investigated the inbreeding load for fertility traits in the Italian Brown Swiss dairy cattle breed. Fertility traits included continuous traits (i.e., interval from calving to first service, days open, and calving interval) and categorical traits (i.e., calving rate at first insemination and nonreturn date at d 56). We included only records of the first 3 parities of cows that calved between 2010 and 2018. We traced up the pedigree of the cows with records as far as possible, ending up with a total of 73,246 animals. The final data set consisted of 59,864 records from 34,921 cows. We analyzed all models using a Bayesian approach that included a covariate with total inbreeding in addition to systematic, permanent environment, additive genetic, and inbreeding load effects. We then evaluated the trends in heritabilities and ratios of the inbreeding load using a continuum of partial inbreeding coefficients from 0.001 to 0.100 as reference. Posterior estimates of heritabilities tended to decrease across the continuum, whereas ratios of the inbreeding load tended to increase, more noticeably in categorical traits (calving rate at first insemination and nonreturn date at d 56). From the results obtained, we confirmed the presence of heterogeneity in inbreeding depression. We then predicted the inbreeding load effects, which had a low reliability of prediction, explained by having only 513 ancestors generating inbreeding. However, reliability of prediction was high enough for some of the individuals, obtaining a favorable prediction of inbreeding load for a relevant percentage, which improved the phenotypic performance of their inbred descendants. These results make it feasible to implement breeding and management strategies that select ancestors with a favorable inbreeding load prediction. In addition, it opens the possibility to define a global index for the expected consequences of the inbreeding generated by each individual.

Genomic and pedigree estimation of inbreeding depression for semen traits in the Basco-Béarnaise dairy sheep breed

Inbreeding depression is associated with a decrease in performance and fitness of the animals. The goal of this study was to evaluate pedigree-based and genomic methods to estimate the level of inbreeding and inbreeding depression for 3 semen traits (volume, concentration, and motility score) in the Basco-Béarnaise sheep breed. Data comprised 16,196 (or 15,071) phenotypic records from 620 rams (of which 533 rams had genotypes of 36,464 SNPs). The pedigree included 8,266 animals, composed of the 620 rams and their ancestors. The number of equivalent complete generations for the 620 rams was 7.04. Inbreeding coefficients were estimated using genomic and pedigree-based information. Genomic inbreeding coefficients were estimated from individual SNP and using segments of homozygous SNP (runs of homozygosity, ROH). Short ROH are of old origin, whereas long ROH are due to recent inbreeding. Considering that the equivalent number of generations in Basco-Béarnaise was 6, inbreeding coefficients for ROH with a length >4 Mb refer to all (recent + old) inbreeding, those with a length >17 Mb correspond to recent inbreeding, and the difference between them indicates old inbreeding. Pedigree-based inbreeding coefficients were also estimated classically, or accounting for nonzero relationships for unknown parents, or including metafounder relationships (estimated using markers) to account for missing pedigree information. Finally, inbreeding coefficients combining genotyped and nongenotyped animal information were computed from matrix H of the single-step approach, also including metafounders. Inbreeding depression was estimated differently depending on the approach used to compute inbreeding coefficients. These 8 estimators of inbreeding coefficients were included as covariates in different animal models. No inbreeding depression was detected for sperm volume or sperm concentration. Inbreeding depression was significant for the motility of spermatozoa. The effect of old and recent inbreeding on motility was null and negative, respectively, demonstrating the existence of purging by selection of deleterious recessive alleles affecting motility. A 10% increase in inbreeding would result in a reduction in mean motility ranging between 0.09 and 0.22 points in the score (from 0 to 5). Motility is unfavorably affected by increasing recent inbreeding but the impact is very small. Runs of homozygosity and metafounders allow us to accurately estimate inbreeding depression and detect recent inbreeding.

Genetic inbreeding depression load for morphological traits and defects in the Pura Raza Española horse

Background

Inbreeding is caused by mating between related individuals and is associated with reduced fitness and performance (inbreeding depression). Several studies have detected heterogeneity in inbreeding depression among founder individuals. Recently, a procedure was developed to predict hidden inbreeding depression load that is associated with founders using the Mendelian sampling of non-founders. The objectives of this study were to: (1) analyse the population structure and general inbreeding, and (2) test this recent approach for predicting hidden inbreeding depression load for four morphological traits and two morphology defects in the Pura Raza Española (PRE) horse breed.

Results

The regression coefficients that were calculated between trait performances and inbreeding coefficients demonstrated the existence of inbreeding depression. In total, 58,772,533 partial inbreeding coefficients (F_ij) were estimated for the whole PRE population (328,706 horses). We selected the descendants of horses with a F_ij ≥ 6.25% that contributed to at least four offspring and for which morphological traits were measured for the subsequent analysis of inbreeding depression load (639 horses). A pedigree was generated with the last five generations (5026 animals) used as the reference population (average inbreeding coefficient of 8.39% and average relatedness coefficient of 10.76%). Heritability estimates ranged from 0.08 (cresty neck) to 0.80 (height at withers), whereas inbreeding depression load ratios ranged from 0.01 (knock knee) to 0.40 (length of shoulder), for an inbreeding coefficient of 10%. Most of the correlations between additive and inbreeding depression load genetic values and correlations between inbreeding depression load genetic values for the different traits were positive or near 0.

Conclusions

Although the average inbreeding depression loads presented negative values, a certain percentage of the animals showed neutral or even positive values. Thus, high levels of inbreeding do not always lead to a decrease in mean phenotypic value or an increase in morphological defects. Hence, individual inbreeding depression loads could be used as a tool to select the most appropriate breeding animals. The possibility of selecting horses that have a high genetic value and are more resistant to the deleterious effects of inbreeding should help improve selection outcomes.

Genetic Structure Analysis of the Pura Raza Español Horse Population through Partial Inbreeding Coefficient Estimation

The aim of this work was to analyze genetic parameters such as the inbreeding coefficient (F), relatedness coefficient (AR) and partial inbreeding coefficient (Fij) of the whole PRE population, and the ancestors which account for 50% of the total genetic variability of the current population, from genealogical information. The average F of the whole PRE population (328,706 animals) has decreased from 8.45% to 7.51% in the least 20 years. The Fij was estimated for the whole PRE population, resulting in a database of 58,772,533 records containing one record for each Fij that each animal receives from a certain common ancestor (CA). A total of 10,244 CAs contributed to the Fij with an average of 5370 descendants, with each descendant having an average of 170 CAs. Over the generations, the number of CAs has increased, while the proportion of Fij by each one has decreased. In addition, the contributions of the more influential ancestors have changed. The increased census, the limited use of artificial insemination and our increased knowledge about inbreeding depression and the animals' breeding values allow breeders to select horses more for their functionality and conformation than for their pedigree reputation, which is the cause of all these changes.

A multivariate analysis with direct additive and inbreeding depression load effects

Background

Inbreeding is caused by mating between related individuals and its most common consequence is inbreeding depression. Several studies have detected heterogeneity in inbreeding depression among founder individuals, and recently a procedure for predicting hidden inbreeding depression loads associated with founders and the Mendelian sampling of non-founders has been developed. The objectives of our study were to expand this model to predict the inbreeding loads for all individuals in the pedigree and to estimate the covariance between the inbreeding loads and the additive genetic effects for the trait of interest. We tested the proposed approach with simulated data and with two datasets of records on weaning weight from the Spanish Pirenaica and Rubia Gallega beef cattle breeds.

Results

The posterior estimates of the variance components with the simulated datasets did not differ significantly from the simulation parameters. In addition, the correlation between the predicted and simulated inbreeding loads were always positive and ranged from 0.27 to 0.82. The beef cattle datasets comprised 35,126 and 75,194 records on weights between 170 and 250 days of age, and pedigrees of 308,836 and 384,434 individual-sire-dam entries for the Pirenaica and Rubia Gallega breeds, respectively. The posterior mean estimates of the variance of inbreeding depression loads were 29,967.8 and 28,222.4 for the Pirenaica and Rubia Gallega breeds, respectively. They were larger than those of the additive variance (695.0 and 439.8 for Pirenaica and Rubia Gallega, respectively), because they should be understood as the variance of the inbreeding depression achieved by a fully inbred (100%) descendant. Therefore, the inbreeding loads have to be rescaled for smaller inbreeding coefficients. In addition, a strong negative correlation (− 0.43 ± 0.10) between additive effects and inbreeding loads was detected in the Pirenaica, but not in the Rubia Gallega breed.

Conclusions

The results of the simulation study confirmed the ability of the proposed procedure to predict inbreeding depression loads for all individuals in the populations. Furthermore, the results obtained from the two real datasets confirmed the variability in the inbreeding depression loads in both breeds and suggested a negative correlation of the inbreeding loads with the additive genetic effects in the Pirenaica breed.

Inbreeding depression load for litter size in Entrepelado and Retinto Iberian pig varieties1

Individual-specific hidden inbreeding depression load (IDL) can be accounted for in livestock populations by appropriate best linear unbiased prediction approaches. This genetic effect has a recessive pattern and reveals when inherited in terms of identity-by-descent. Nevertheless, IDL inherits as a pure additive genetic background and can be selected using standard breeding values. The main target of this research was to evaluate IDL for litter size in 2 Iberian pig varieties (Entrepelado and Retinto) from a commercial breeding-stock. Analyses were performed on the total number of piglets born (both alive and dead) and used data from 3,200 (8.02 ± 0.04 piglets/litter) Entrepelado and 4,744 Retinto litters (8.40 ± 0.03 piglets/litter). Almost 50% of Entrepelado sows were inbred (1.7% to 25.0%), whereas this percentage reduced to 37.4% in the Retinto variety (0.2% to 25.0%). The analytical model was solved by Bayesian inference and accounted for 2 systematic effects (sow age and breed/variety of the artificial insemination boar), 2 permanent environmental effects (herd-year-season and sow), and 2 genetic effects (IDL and infinitesimal additive). In terms of posterior means (PM), additive genetic and IDL variances were similar in the Entrepelado variety (PM, 0.68 vs. 0.76 piglets2, respectively) and their 95% credibility intervals (95CI) overlapped, although without including zero (0.38 to 0.94 vs. 0.15 to 1.31 piglets2, respectively). The same pattern revealed in the Retinto variety, with IDL variance (PM, 0.41 piglets2; 95CI, 0.07 to 0.88 piglets2) slightly larger than the additive genetic variance (PM, 0.37 piglets2; 95CI, 0.16 to 0.59 piglets2). The relevance of IDL was also checked by a Bayes factor and the deviance information criterion, the model including this effect being clearly favored in both cases. Although the analysis assumed null genetic covariance between IDL and infinitesimal additive effects, a moderate negative correlation (-0.31) was suggested when plotting the PM of breeding values in the Entrepelado variety; a negative genetic trend for IDL was also revealed in this Iberian pig variety (-0.25 piglets for 100% inbred offspring of individuals born in 2014), whereas no trend was detected in Retinto breeding-stock. Those were the first estimates of IDL in a commercial livestock population, they giving evidence of a relevant genetic background with potential consequences on the reproductive performance of Iberian sows.