Redefining and interpreting genomic relationships of metafounders

Abstract

Abstract Metafounders are a useful concept to characterize relationships within and across populations, and to help genetic evaluations because they help modelling the means and variances of unknown base population animals. Current definitions of metafounder relationships are sensitive to the choice of reference alleles and have not been compared to their counterparts in population genetics—namely, heterozygosities, F ST coefficients, and genetic distances. We redefine the relationships across populations with an arbitrary base of a maximum heterozygosity population in Hardy–Weinberg equilibrium. Then, the relationship between or within populations is a cross-product of the form $${\Gamma }_{\left(b,{b}^{\prime}\right)}=\left(\frac{2}{n}\right)\left(2{\mathbf{p}}_{b}-\mathbf{1}\right)\left(2{\mathbf{p}}_{{b}^{\prime}}-\mathbf{1}\right)^{\prime}$$ Γ b , b ′ = 2 n 2 p b - 1 2 p b ′ - 1 ′ with $$\mathbf{p}$$ p being vectors of allele frequencies at $$n$$ n markers in populations $$b$$ b and $$b^{\prime}$$ b ′ . This is simply the genomic relationship of two pseudo-individuals whose genotypes are equal to twice the allele frequencies. We also show that this coding is invariant to the choice of reference alleles. In addition, standard population genetics metrics (inbreeding coefficients of various forms; F ST differentiation coefficients; segregation variance; and Nei’s genetic distance) can be obtained from elements of matrix $${\varvec{\Gamma}}$$ Γ .