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

Abstract

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00772-0.