Interpreting SNP heritability in admixed populations.

SNP heritability (h2snp) is defined as the proportion of phenotypic variance explained by genotyped SNPs and is believed to be a lower bound of heritability (h), being equal to it if all causal variants are genotyped. Despite the simple intuition behind h2snp, its interpretation and equivalence to h2 is unclear, particularly in the presence of admixture and assortative mating. Here we use analytical theory and simulations to describe the behavior of h2 and three widely used random-effect estimators of h2snp -- Genome-wide restricted maximum likelihood, Haseman-Elston regression, and LD score regression -- in admixed populations. We show that h2snp estimates can be biased in admixed populations, even if all causal variants are genotyped and in the absence of confounding due to shared environment. This is largely because admixture generates directional LD, which contributes to the genetic variance, and therefore to heritability. Random-effect estimators of h2snp, because they assume that SNP effects are independent, do not capture the contribution, which can be positive or negative depending on the genetic architecture, leading to under- or over-estimates of h2snp relative to h2. For the same reason, estimates of local ancestry heritability (ĥ2γ) are also biased in the presence of directional LD. We describe this bias in ĥ2snp and ĥ2γ as a function of admixture history and the genetic architecture of the trait, clarifying their interpretation and implication for genome-wide association studies and polygenic prediction in admixed populations.