Background: Inbreeding decreases animal performance (inbreeding depression), but not all inbreeding is expected to be equally harmful. Recent inbreeding is expected to be more harmful than ancient inbreeding, because selection decreases the frequency of deleterious alleles over time. Selection efficiency is increased by inbreeding, a process called purging. Our objective was to investigate effects of recent and ancient inbreeding on yield, fertility and udder health traits in Dutch Holstein–Friesian cows. Methods: In total, 38,792 first-parity cows were included. Pedigree inbreeding ( $$F_{PED}$$ F PED ) was computed and 75 k genotype data were used to compute genomic inbreeding, among others based on regions of homozygosity (ROH) in the genome ( $$F_{ROH}$$ F ROH ). Results: Inbreeding depression was observed, e.g. a 1% increase in $$F_{ROH}$$ F ROH was associated with a 36.3 kg (SE = 2.4) decrease in 305-day milk yield, a 0.48 day (SE = 0.15) increase in calving interval and a 0.86 unit (SE = 0.28) increase in somatic cell score for day 150 through to 400. These effects equalled − 0.45, 0.12 and 0.05% of the trait means, respectively. When $$F_{PED}$$ F PED was split into generation-based components, inbreeding on recent generations was more harmful than inbreeding on more distant generations for yield traits. When $$F_{PED}$$ F PED was split into new and ancestral components, based on whether alleles were identical-by-descent for the first time or not, new inbreeding was more harmful than ancestral inbreeding, especially for yield traits. For example, a 1% increase in new inbreeding was associated with a 2.42 kg (SE = 0.41) decrease in 305-day fat yield, compared to a 0.03 kg (SE = 0.71) increase for ancestral inbreeding. There were no clear differences between effects of long ROH (recent inbreeding) and short ROH (ancient inbreeding). Conclusions: Inbreeding depression was observed for yield, fertility and udder health traits. For yield traits and based on pedigree, inbreeding on recent generations was more harmful than inbreeding on distant generations and there was evidence of purging. Across all traits, long and short ROH contributed to inbreeding depression. In future work, inbreeding depression and purging should be assessed in more detail at the genomic level, using higher density information and genomic time series.

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.

In recent decades, Holstein–Friesian (HF) selection schemes have undergone profound changes, including the introduction of optimal contribution selection (OCS; around 2000), a major shift in breeding goal composition (around 2000) and the implementation of genomic selection (GS; around 2010). These changes are expected to have influenced genetic diversity trends. Our aim was to evaluate genome-wide and region-specific diversity in HF artificial insemination (AI) bulls in the Dutch-Flemish breeding program from 1986 to 2015. Pedigree and genotype data (~ 75.5 k) of 6280 AI-bulls were used to estimate rates of genome-wide inbreeding and kinship and corresponding effective population sizes. Region-specific inbreeding trends were evaluated using regions of homozygosity (ROH). Changes in observed allele frequencies were compared to those expected under pure drift to identify putative regions under selection. We also investigated the direction of changes in allele frequency over time. Effective population size estimates for the 1986–2015 period ranged from 69 to 102. Two major breakpoints were observed in genome-wide inbreeding and kinship trends. Around 2000, inbreeding and kinship levels temporarily dropped. From 2010 onwards, they steeply increased, with pedigree-based, ROH-based and marker-based inbreeding rates as high as 1.8, 2.1 and 2.8% per generation, respectively. Accumulation of inbreeding varied substantially across the genome. A considerable fraction of markers showed changes in allele frequency that were greater than expected under pure drift. Putative selected regions harboured many quantitative trait loci (QTL) associated to a wide range of traits. In consecutive 5-year periods, allele frequencies changed more often in the same direction than in opposite directions, except when comparing the 1996–2000 and 2001–2005 periods. Genome-wide and region-specific diversity trends reflect major changes in the Dutch-Flemish HF breeding program. Introduction of OCS and the shift in breeding goal were followed by a drop in inbreeding and kinship and a shift in the direction of changes in allele frequency. After introduction of GS, rates of inbreeding and kinship increased substantially while allele frequencies continued to change in the same direction as before GS. These results provide insight in the effect of breeding practices on genomic diversity and emphasize the need for efficient management of genetic diversity in GS schemes.

To test for the presence of purging in populations, the classical pedigree-based inbreeding coefficient (F) can be decomposed into Kalinowski’s ancestral (F_ANC) and new (F_NEW) inbreeding coefficients. The F_ANC and F_NEW can be calculated by a stochastic approach known as gene dropping. However, the only publicly available algorithm for the calculation of F_ANC and F_NEW, implemented in GRain v 2.1 (and also incorporated in the PEDIG software package), has produced biased estimates. The F_ANC was systematically underestimated and consequently, F_NEW was overestimated. To illustrate this bias, we calculated F_ANC and F_NEW by hand for simple example pedigrees. We revised the GRain program so that it now provides unbiased estimates. Correlations between the biased and unbiased estimates of F_ANC and F_NEW, obtained for example data sets of Hungarian Pannon White rabbits (22,781 individuals) and Dutch Holstein Friesian cattle (37,061 individuals), were high, i.e., >0.96. Although the magnitude of bias appeared to be small, results from studies based on biased estimates should be interpreted with caution. The revised GRain program (v 2.2) is now available online and can be used to calculate unbiased estimates of F_ANC and F_NEW.

Background: Inbreeding depression refers to the decrease in mean performance due to inbreeding. Inbreeding depression is caused by an increase in homozygosity and reduced expression of (on average) favourable dominance effects. Dominance effects and allele frequencies differ across loci, and consequently inbreeding depression is expected to differ along the genome. In this study, we investigated differences in inbreeding depression across the genome of Dutch Holstein Friesian cattle, by estimating dominance effects and effects of regions of homozygosity (ROH). Methods: Genotype (75 k) and phenotype data of 38,792 cows were used. For nine yield, fertility and udder health traits, GREML models were run to estimate genome-wide inbreeding depression and estimate additive, dominance and ROH variance components. For this purpose, we introduced a ROH-based relationship matrix. Additive, dominance and ROH effects per SNP were obtained through back-solving. In addition, a single SNP GWAS was performed to identify significant additive, dominance or ROH associations. Results: Genome-wide inbreeding depression was observed for all yield, fertility and udder health traits. For example, a 1% increase in genome-wide homozygosity was associated with a decrease in 305-d milk yield of approximately 99 kg. For yield traits only, including dominance and ROH effects in the GREML model resulted in a better fit ( P < 0.05) than a model with only additive effects. After correcting for the effect of genome-wide homozygosity, dominance and ROH variance explained less than 1% of the phenotypic variance for all traits. Furthermore, dominance and ROH effects were distributed evenly along the genome. The most notable region with a favourable dominance effect for yield traits was on chromosome 5, but overall few regions with large favourable dominance effects and significant dominance associations were detected. No significant ROH-associations were found. Conclusions: Inbreeding depression was distributed quite equally along the genome and was well captured by genome-wide homozygosity. These findings suggest that, based on 75 k SNP data, there is little benefit of accounting for region-specific inbreeding depression in selection schemes.

National gene bank collections for Holstein Friesian (HF) dairy cattle were set up in the 1990s. In this study, we assessed the value of bulls from the Dutch HF germplasm collection, also known as cryobank bulls, to increase genetic variability and improve genetic merit in the current bull population (bulls born in 2010-2015). Genetic variability was defined as 1 minus the mean genomic similarity (SIMSNP) or as 1 minus the mean pedigree-based kinship (fPED). Genetic merit was defined as the mean estimated breeding value for the total merit index or for 1 of 3 subindices (yield, fertility, and udder health). Using optimal contribution selection, we minimized relatedness (maximized variability) or maximized genetic merit at restricted levels of relatedness. We compared breeding schemes with only bulls from 2010 to 2015 with schemes in which cryobank bulls were also included. When we minimized relatedness, inclusion of genotyped cryobank bulls decreased mean SIMSNP by 0.7% and inclusion of both genotyped and nongenotyped cryobank bulls decreased mean fPED by 2.6% (in absolute terms). When we maximized merit at restricted levels of relatedness, inclusion of cryobank bulls provided additional merit at any level of mean SIMSNP or mean fPED except for the total merit index at high levels of mean SIMSNP. Additional merit from cryobank bulls depended on (1) the relative emphasis on genetic variability and (2) the selection criterion. Additional merit was higher when more emphasis was put on genetic variability. For fertility, for example, it was 1.74 SD at a mean SIMSNP restriction of 64.5% and 0.37 SD at a mean SIMSNP restriction of 67.5%. Additional merit was low to nonexistent for the total merit index and higher for the subindices, especially for fertility. At a mean SIMSNP of 64.5%, for example, it was 0.60 SD for the total merit index and 1.74 SD for fertility. In conclusion, Dutch HF cryobank bulls can be used to increase genetic variability and improve genetic merit in the current population, although their value is very limited when selecting for the current total merit index. Anticipating changes in the breeding goal in the future, the germplasm collection is a valuable resource for commercial breeding populations.