Interpreting SNP heritability in admixed populations

Abstract

AbstractSNP heritabilityis defined as the proportion of phenotypic variance explained by genotyped SNPs and is believed to be a lower bound of heritability (h2), being equal to it if all causal variants are known. Despite the simple intuition behind, its interpretation and equivalence toh2is unclear, particularly in the presence of population structure and assortative mating. It is well known that population structure can lead to inflation inestimates. Here we use analytical theory and simulations to demonstrate thatestimated with genome-wide restricted maximum likelihood (GREML) can be biased in admixed populations, even in the absence of confounding and even if all causal variants are known. This is because admixture generates linkage disequilibrium (LD), which contributes to the genetic variance, and therefore to heritability. GREML implicitly assumes this component is zero, which may not be true, particularly for traits under divergent or stabilizing selection in the source populations, leading under- or over-estimates ofrelative toh2. For the same reason, GREML estimates of local ancestry heritabilitywill also be biased. We describe the bias inandas a function of admixture history and the genetic architecture of the trait and discuss its implications for genome-wide association and polygenic prediction.