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

Reproductive traits in Pura Raza Española mares manifest inbreeding depression from low levels of homozygosity

Inbreeding depression is a genetic phenomenon associated with the loss of fitness and mean phenotypic performance due to mating between relatives. Historically, inbreeding coefficients have been estimated from pedigree information. However, the onset of genomic selection programs provides large datasets of individuals genotyped using SNP arrays, enabling more precise assessment of an individual's genomic-level inbreeding using genomic data. One of the traits most sensitive to issues stemming from increased inbreeding is reproduction. This is particularly important in equine, in which fertility is only moderate compared to other livestock species. To explore this further, we evaluated the effect of inbreeding on five reproductive traits (age at first foaling (AFF), average interval between foalings (AIF), total number of foalings (NF), productive life (PL) and reproductive efficiency (RE)) in Pura Raza Español mares using genomic data. Residual predicted phenotypes were obtained by purging these traits through the REML (wgResidual) and ssGREML (gResidual) approaches in reproductive data of 29,847 PRE mares using the BLUPF90+ program. Next, we used pedigree-based (Fped) and ROH-based genomic (FROH) inbreeding coefficients derived from 1018 animals genotyped with 61,271 SNPs to estimate the inbreeding depression (linear regression). Our results indicated significant levels of inbreeding depression for all reproductive traits, with the exception of the AIF trait when Fped was used. However, all traits were negatively affected by the increase in genomic inbreeding, and FROH was found to capture more inbreeding depression than Fped. Likewise, REML models (ssGREML) using genomic data for estimated predicted residual phenotypes resulted in higher variance explained by the model compared with the models not using genomics (REML). Finally, a segmented regression analysis was conducted to evaluate the effect of inbreeding depression, revealing that the levels of genealogical and genomic homozygosity do not manifest uniformly in reproductive traits. In contrast, the levels of inbreeding depression ranged from low to high as homozygosity increased. This analysis also showed that reproductive traits are very sensitive to inbreeding depression, even with relatively low levels of homozygosity.

Impact of inbreeding on production, fertility, and health traits in German Holstein dairy cattle utilizing various inbreeding estimators

In dairy cattle production, it is important to understand how inbreeding affects production, fertility, and health traits. However, there is still limited use of genomic information to estimate inbreeding, despite advancements in genotyping technologies. To address this gap, we investigated the effect of inbreeding on German Holstein dairy cattle using both pedigree-based and genomic-based inbreeding estimators. We employed one method based on pedigree information (Fped) together with 6 genomic-based methods, including 3 genome-wide complex trait analysis software estimators (Fhat1, Fhat2, Fhat3), VanRaden's first method (FVR1, with observed allele frequencies, and FVR0.5, when allele frequencies are set to 0.5), and one based on runs of homozygosity (Froh). Data from 24,489 cows with both phenotypes and genotypes were used, with a pedigree including 232,780 animals born between 1970 and 2018. We analyzed the effects of inbreeding depression on production, fertility, and health traits separately, using single-trait linear animal models as well as threshold models to account for the binary nature of the health traits. For the health traits, we transformed solutions from the liability scale to a probability scale for easier interpretation. Our results showed that the mean inbreeding coefficients from all estimators ranged from -0.003 to 0.243, with negative values observed for most genomic-based methods. We found out that a 1% increase in inbreeding caused a depression ranging from 25.94 kg (Fhat1) to 40.62 kg (Fhat3), 1.18 kg (Fhat2) to 1.70 kg (Fhat3), 0.90 kg (Fhat2) to 1.45 kg (Froh and Fhat3), 0.19 (Fped) to 0.34 d (Fhat3) for 305-d milk yield, fat, protein, and calving interval, respectively. The health traits showed very slight gradual changes when inbreeding was increased steadily from 0% to 50%, with digital dermatitis showing a rather contrasting trend to that of mastitis, which increased the more an animal was inbred. Overall, our study highlights the importance of considering both pedigree-based and genomic-based inbreeding estimators when assessing the impact on inbreeding, emphasizing that not all inbreeding is harmful.

Characterizing and Eliminating the Inbreeding Load

The authors evaluated the relevant literature related to purging, which is the interaction between selection and inbreeding in which the population may eliminate its inbreeding load at least partially. According to the relevant literature, the inbreeding load and the process of purging were evaluated via pedigree methods based on ancestral inbreeding, the inbreeding-purging model, and expressed opportunity of purging, along with genomic methods. Most ancestral inbreeding-related studies were performed in zoos, where only a small proportion of the studied populations show signs of purging. The inbreeding-purging model was developed with Drosophila, and it was used to evaluate different zoo ungulates and Pannon white rabbits. Purging was detected in both studies. The expressed opportunity of purging was applied in Jersey cattle and Pannon white rabbits. In the Jersey cattle, it had an effect of 12.6% for fitness, while in the Pannon white rabbits, the inbreeding load was between 40% and 80% of its original value. The genomic studies also signalled purging, but they also made it clear that, contrary to the detected purging, the evaluated populations still suffered from inbreeding depression. Therefore, especially for domesticated animals, it can be concluded that deliberate inbreeding with the purpose of generating purging is not advocated.

Effect of Individual Rate of Inbreeding, Recent and Ancestral Inbreeding on Wool Quality, Muscling Conformation and Exterior in German Sheep Breeds

This study provides comprehensive results on the current status of inbreeding depression for traits upon which sheep are selected for the herdbook in Germany. A total of 30 sheep breeds from the OviCap national database met the inclusion criteria for the present analysis regarding the depth and completeness of pedigrees and the number of animals with phenotypic data. We analyzed heritabilities and inbreeding depression for the three breeding objective traits of wool quality, muscling conformation and exterior. Heritabilities were across all breeds of moderate size, with estimates of 0.18 for wool quality and muscling conformation and of 0.14 for exterior. The models employed to estimate linear regression slopes for individual and ancestral inbreeding rates also account for non-genetic effects and the additive genetic effect of the animal. Inbreeding depression was obvious for all three traits when we averaged the estimates across all 30 sheep breeds. Inbreeding depression was significant for wool quality for only a few breeds, whereas for muscling conformation, 14/30 breeds achieved significant estimates. A 1% increase in inbreeding decreased the mean of all three traits across all sheep breeds by 0.33% of their standard deviation. Positive effects due to ancestral inbreeding were only significant in very few breeds in the different traits. Across all 30 sheep breeds, there were indications that purging effects (a reduction in negative effects of inbreeding depression due to selection for heterozygotes) may play a role for the exterior. The results of this study should help for reviewing breeding programs, particularly for sheep breeds with critical effective population sizes and increasing rates of inbreeding, with regard to the selection policy and selection intensity applied.

Genome-wide assessment and mapping of inbreeding depression identifies candidate genes associated with semen traits in Holstein bulls

BACKGROUND

The reduction in phenotypic performance of a population due to mating between close relatives is called inbreeding depression. The genetic background of inbreeding depression for semen traits is poorly understood. Thus, the objectives were to estimate the effect of inbreeding and to identify genomic regions underlying inbreeding depression of semen traits including ejaculate volume (EV), sperm concentration (SC), and sperm motility (SM). The dataset comprised ~ 330 K semen records from ~ 1.5 K Holstein bulls genotyped with 50 K single nucleotide polymorphism (SNP) BeadChip. Genomic inbreeding coefficients were estimated using runs of homozygosity (i.e., F_ROH > 1 Mb) and excess of SNP homozygosity (F_SNP). The effect of inbreeding was estimated by regressing phenotypes of semen traits on inbreeding coefficients. Associated variants with inbreeding depression were also detected by regressing phenotypes on ROH state of the variants.

RESULTS

Significant inbreeding depression was observed for SC and SM (p < 0.01). A 1% increase in F_ROH reduced SM and SC by 0.28% and 0.42% of the population mean, respectively. By splitting F_ROH into different lengths, we found significant reduction in SC and SM due to longer ROH, which is indicative of more recent inbreeding. A genome-wide association study revealed two signals positioned on BTA 8 associated with inbreeding depression of SC (p < 0.00001; FDR < 0.02). Three candidate genes of GALNTL6, HMGB2, and ADAM29, located in these regions, have established and conserved connections with reproduction and/or male fertility. Moreover, six genomic regions on BTA 3, 9, 21 and 28 were associated with SM (p < 0.0001; FDR < 0.08). These genomic regions contained genes including PRMT6, SCAPER, EDC3, and LIN28B with established connections to spermatogenesis or fertility.

CONCLUSIONS

Inbreeding depression adversely affects SC and SM, with evidence that longer ROH, or more recent inbreeding, being especially detrimental. There are genomic regions associated with semen traits that seems to be especially sensitive to homozygosity, and evidence to support some from other studies. Breeding companies may wish to consider avoiding homozygosity in these regions for potential artificial insemination sires.

Genetic Diversity and Trends of Ancestral and New Inbreeding in German Sheep Breeds by Pedigree Data

In Germany, many autochthonous sheep breeds have developed, adapted to mountain, heath, moorland, or other marginal sites, but breeds imported from other countries have also contributed to the domestic breeds, particularly improving wool and meat quality. Selective breeding and the intense use of rams may risk losing genetic diversity and increasing rates of inbreeding. On the other hand, breeds with a low number of founder animals and only regional popularity may not leave their endangered status, as the number of breeders interested in the breed is limited. The objective of the present study was to determine demographic measures of genetic diversity and recent as well as ancestral trends of inbreeding in all autochthonous German sheep breeds and sheep of all breeding directions, including wool, meat, and milk. We used pedigree data from 1,435,562 sheep of 35 different breeds and a reference population of 981,093 sheep, born from 2010 to 2020. The mean number of equivalent generations, founders, effective founders, effective ancestors, and effective founder genomes were 5.77, 1669, 123.2, 63.5, and 33.0, respectively. Genetic drift accounted for 69% of the loss of genetic diversity, while loss due to unequal founder contributions was 31%. The mean inbreeding coefficient, individual rate of inbreeding (∆F_i), and realized effective population size across breeds were 0.031, 0.0074, and 91.4, respectively, with a significantly decreasing trend in ∆F_i in 11/35 breeds. New inbreeding, according to Kalinowski, contributed to 71.8% of individual inbreeding, but ancestral inbreeding coefficients showed an increasing trend in all breeds. In conclusion, in our study, all but one of the mountain-stone sheep breeds and the country sheep breed Wald were the most vulnerable populations, with N_e < 50. The next most endangered breeds are exotic, country, and heath breeds, with average N_e of 66, 83, and 89, respectively. The wool, meat, and milk breeds showed the highest genetic diversity, with average N_e of 158, 120, and 111, respectively. The results of our study should help strengthen conservation program efforts for the most endangered sheep breeds and maintain a high genetic diversity in all sheep breeds.

Analysis of runs of homozygosity of cattle living in different climate zones

Aim of this study was to analyze the distribution and characteristics of runs of homozygosity in Bos taurus taurus and Bos taurus indicus breeds, as well as their crosses, farmed all around the world. With this aim in view, we used single-nucleotide polymorphisms (SNP) genotypes for 3,263 cattle belonging to 204 different breeds. After quality control, 23,311 SNPs were retained for the analysis. Animals were divided into seven different groups: 1) continental taurus, 2) temperate taurus, 3) temperate indicus, 4) temperate composite, 5) tropical taurus, 6) tropical indicus, and 7) tropical composite. The climatic zones were created according to the latitude of the breeds' country of origin: i) continental, latitude ≥ 45°; ii) temperate, 45°< Latitude >23.26°; iii) tropics, latitude ≤ 23.26°. Runs of homozygosity were computed as 15 SNPs spanning in at least 2 Mb; number of ROH per animal (nROH), average ROH length (meanMb), and ROH-based inbreeding coefficients (FROH) were also computed. Temperate indicus showed the largest nROH, whereas Temperate taurus the lowest value. Moreover, the largest meanMb was observed for Temperate taurus, whereas the lowest value for Tropics indicus. Temperate indicus breeds showed the largest FROH values. Genes mapped in the identified ROH were reported to be associated with the environmental adaptation, disease resistance, coat color determinism, and production traits. Results of the present study confirmed that runs of homozygosity could be used to identify genomic signatures due to both artificial and natural selection.

A comparison of marker-based estimators of inbreeding and inbreeding depression

BACKGROUND

The availability of genome-wide marker data allows estimation of inbreeding coefficients (F, the probability of identity-by-descent, IBD) and, in turn, estimation of the rate of inbreeding depression (ΔID). We investigated, by computer simulations, the accuracy of the most popular estimators of inbreeding based on molecular markers when computing F and ΔID in populations under random mating, equalization of parental contributions, and artificially selected populations. We assessed estimators described by Li and Horvitz (F_LH1 and F_LH2), VanRaden (F_VR1 and F_VR2), Yang and colleagues (F_YA1 and F_YA2), marker homozygosity (F_HOM), runs of homozygosity (F_ROH) and estimates based on pedigree (F_PED) in comparison with estimates obtained from IBD measures (F_IBD).

RESULTS

If the allele frequencies of a base population taken as a reference for the computation of inbreeding are known, all estimators based on marker allele frequencies are highly correlated with F_IBD and provide accurate estimates of the mean ΔID. If base population allele frequencies are unknown and current frequencies are used in the estimations, the largest correlation with F_IBD is generally obtained by F_LH1 and the best estimator of ΔID is F_YA2. The estimators F_VR2 and F_LH2 have the poorest performance in most scenarios. The assumption that base population allele frequencies are equal to 0.5 results in very biased estimates of the average inbreeding coefficient but they are highly correlated with F_IBD and give relatively good estimates of ΔID. Estimates obtained directly from marker homozygosity (F_HOM) substantially overestimated ΔID. Estimates based on runs of homozygosity (F_ROH) provide accurate estimates of inbreeding and ΔID. Finally, estimates based on pedigree (F_PED) show a lower correlation with F_IBD than molecular estimators but provide rather accurate estimates of ΔID. An analysis of data from a pig population supports the main findings of the simulations.

CONCLUSIONS

When base population allele frequencies are known, all marker-allele frequency-based estimators of inbreeding coefficients generally show a high correlation with F_IBD and provide good estimates of ΔID. When base population allele frequencies are unknown, F_LH1 is the marker frequency-based estimator that is most correlated with F_IBD, and F_YA2 provides the most accurate estimates of ΔID. Estimates from F_ROH are also very precise in most scenarios. The estimators F_VR2 and F_LH2 have the poorest performances.

Fine-tuning genomic and pedigree inbreeding rates in equine population with a deep and reliable stud book: the case of the Pura Raza Española horse

BACKGROUND

Estimating inbreeding, which is omnipresent and inevitable in livestock populations, is a primary goal for management and animal breeding especially for those interested in mitigating the negative consequences of inbreeding. Inbreeding coefficients have been historically estimated by using pedigree information; however, over the last decade, genome-base inbreeding coefficients have come to the forefront in this field. The Pura Raza Española (PRE) horse is an autochthonous Spanish horse breed which has been recognised since 1912. The total PRE population (344,718 horses) was used to estimate Classical (F), Ballou's ancestral, Kalinowski's ancestral, Kalinowski's new and the ancestral history coefficient values. In addition, genotypic data from a selected population of 805 PRE individuals was used to determine the individual inbreeding coefficient using SNP-by-SNP-based techniques (methods of moments -F_HOM-, the diagonal elements of the genomic -F_G-, and hybrid matrixes -F_H-) and ROH measures (F_RZ). The analyse of both pedigree and genomic based inbreeding coefficients in a large and robust population such as the PRE horse, with proven parenteral information for the last 40 years and a high degree of completeness (over 90% for the last 70 years) will allow us to understand PRE genetic variability better and the correlations between the estimations will give the data greater reliability.

RESULTS

The mean values of the pedigree-based inbreeding coefficients ranged from 0.01 (F for the last 3 generations -F3-) to 0.44 (ancestral history coefficient) and the mean values of genomic-based inbreeding coefficients varied from 0.05 (F_RZ for three generations, F_H and F_HOM) to 0.11 (F_RZ for nine generations). Significant correlations were also found between pedigree and genomic inbreeding values, which ranged between 0.58 (F3 with F_HOM) and 0.79 (F with F_RZ). In addition, the correlations between F_RZ estimated for the last 20 generations and the pedigree-based inbreeding highlight the fact that fewer generations of genomic data are required when comparing total inbreeding values, and the opposite when ancient values are calculated.

CONCLUSIONS

Ultimately, our results show that it is still useful to work with a deep and reliable pedigree in pedigree-based genetic studies with very large effective population sizes. Obtaining a satisfactory parameter will always be desirable, but the approximation obtained with a robust pedigree will allow us to work more efficiently and economically than with massive genotyping.

A comprehensive genome-wide scan detects genomic regions related to local adaptation and climate resilience in Mediterranean domestic sheep

Background

The management of farm animal genetic resources and the adaptation of animals to climate change will probably have major effects on the long-term sustainability of the livestock sector. Genomic data harbour useful relevant information that needs to be harnessed for effectively managing genetic resources. In this paper, we report the genome characterization of the highly productive Mediterranean Chios dairy sheep and focus on genetic diversity measures related with local adaptation and selection and the genetic architecture of animal resilience to weather fluctuations as a novel adaptative trait linked to climate change.

Results

We detected runs of homozygosity (ROH) and heterozygosity (ROHet) that revealed multiple highly homozygous and heterozygous hotspots across the Chios sheep genome. A particularly highly homozygous region was identified on chromosome 13 as a candidate of directional genetic selection associated with milk traits, which includes annotated genes that were previously shown to be linked to local adaptation to harsh environmental conditions. Favourable heterozygosity related with a potentially protective role against livestock diseases and enhanced overall fitness was revealed in heterozygous-rich regions on sheep chromosomes 3, 10, 13 and 19. Furthermore, genomic analyses were conducted on sheep resilience phenotypes that display changes in milk production in response to weather variation. Sheep resilience to heat stress was a significantly heritable trait (h² = 0.26) and genetically antagonistic to milk production. Genome-wide association and regional heritability mapping analyses revealed novel genomic markers and regions on chromosome 5 that were significantly associated with sheep resilience to climate change. Subsequently, an annotation analysis detected a set of genes on chromosome 5 that were associated with olfactory receptor complexes that could participate in heat stress mitigation through changes in respiration rate and respiratory evaporation. Other genes were grouped in previously reported biological processes relevant to livestock heat dissipation, including stress and immune response.

Conclusions

Our results may contribute to the optimal management of sheep genetic resources and inform modern selective breeding programmes that aim at mitigating future environmental challenges towards sustainable farming, while better balancing animal adaptation and productivity. Our results are directly relevant to the studied breed and the respective environmental conditions; however, the methodology may be extended to other livestock species of interest.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00682-7.

Estimating inbreeding depression for growth and reproductive traits using pedigree and genomic methods in Argentinean Brangus cattle

Inbreeding depression reduces the mean phenotypic value of important traits in livestock populations. The goal of this work was to estimate the level of inbreeding and inbreeding depression for growth and reproductive traits in Argentinean Brangus cattle, in order to obtain a diagnosis and monitor breed management. Data comprised 359,257 (from which 1,990 were genotyped for 40,678 single nucleotide polymorphisms [SNPs]) animals with phenotypic records for at least one of three growth traits: birth weight (BW), weaning weight (WW), and finishing weight (FW). For scrotal circumference (SC), 52,399 phenotypic records (of which 256 had genotype) were available. There were 530,938 animals in pedigree. Three methods to estimate inbreeding coefficients were used. Pedigree-based inbreeding coefficients were estimated accounting for missing parents. Inbreeding coefficients combining genotyped and nongenotyped animal information were also computed from matrix H of the single-step approach. Genomic inbreeding coefficients were estimated using homozygous segments obtained from a Hidden Markov model (HMM) approach. Inbreeding depression was estimated from the regression of the phenotype on inbreeding coefficients in a multiple-trait mixed model framework, either for the whole dataset or for the dataset of genotyped animals. All traits were unfavorably affected by inbreeding depression. A 10% increase in pedigree-based or combined inbreeding would result in a reduction of 0.34 to 0.39 kg in BW, 2.77 to 3.28 kg in WW, and 0.23 cm in SC. For FW, a 10% increase in pedigree-based, genomic, or combined inbreeding would result in a decrease of 8.05 to 11.57 kg. Genomic inbreeding based on the HMM was able to capture inbreeding depression, even in such a compressed genotyped dataset.

How Depressing Is Inbreeding? A Meta-Analysis of 30 Years of Research on the Effects of Inbreeding in Livestock

Inbreeding depression has been widely documented for livestock and other animal and plant populations. Inbreeding is generally expected to have a stronger unfavorable effect on fitness traits than on other traits. Traditionally, the degree of inbreeding depression in livestock has been estimated as the slope of the linear regression of phenotypic values on pedigree-based inbreeding coefficients. With the increasing availability of SNP-data, pedigree inbreeding can now be replaced by SNP-based measures. We performed a meta-analysis of 154 studies, published from 1990 to 2020 on seven livestock species, and compared the degree of inbreeding depression (1) across different trait groups, and (2) across different pedigree-based and SNP-based measures of inbreeding. Across all studies and traits, a 1% increase in pedigree inbreeding was associated with a median decrease in phenotypic value of 0.13% of a trait’s mean, or 0.59% of a trait’s standard deviation. Inbreeding had an unfavorable effect on all sorts of traits and there was no evidence for a stronger effect on primary fitness traits (e.g., reproduction/survival traits) than on other traits (e.g., production traits or morphological traits). p-values of inbreeding depression estimates were smaller for SNP-based inbreeding measures than for pedigree inbreeding, suggesting more power for SNP-based measures. There were no consistent differences in p-values for percentage of homozygous SNPs, inbreeding based on runs of homozygosity (ROH) or inbreeding based on a genomic relationship matrix. The number of studies that directly compares these different measures, however, is limited and comparisons are furthermore complicated by differences in scale and arbitrary definitions of particularly ROH-based inbreeding. To facilitate comparisons across studies in future, we provide the dataset with inbreeding depression estimates of 154 studies and stress the importance of always reporting detailed information (on traits, inbreeding coefficients, and models used) along with inbreeding depression estimates.

Association of inbreeding and regional equine leucocyte antigen homozygosity with the prevalence of insect bite hypersensitivity in Old Kladruber horse

Inbreeding depression is the reduction of performance caused by mating of close relatives. In livestock populations, inbreeding depression has been traditionally estimated by regression of phenotypes on pedigree inbreeding coefficients. This estimation can be improved by utilising genomic inbreeding coefficients. Here we estimate inbreeding depression for insect bite hypersensitivity (IBH) prevalence, the most common allergic horse disease worldwide, in Old Kladruber horse. In a deep pedigree with 3214 horses (187 genotyped), we used a generalised linear mixed model with IBH phenotype from 558 horses examined between 1996 and 2009 (1368 records). In addition to the classical pedigree information, we used the single-step approach that enabled joint use of pedigree and genomic information to estimate inbreeding depression overall genome and equine leucocyte antigen (ELA) class II region. Significant inbreeding depression was observed in all models fitting overall inbreeding coefficients (odds ratio between 1.018 and 1.074, P < 0.05) with the exception of Kalinowski's new inbreeding (P = 0.0516). The increase of ELA class II inbreeding was significantly associated with increased prevalence of IBH (odds ratio 1.018; P = 0.027). However, when fitted jointly with the overall inbreeding coefficient, the effect of ELA class II inbreeding was not significant (odds ratio 1.016; P = 0.062). Overall, the higher ELA class II and/or overall inbreeding (pedigree or genomic) was associated with increased prevalence of IBH in Old Kladruber horses. The single-step approach provides an efficient use of all the available pedigree, genomic, and phenotype information for estimation of overall and regional inbreeding effects.