PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Hum Genet. Author manuscript; available in PMC 2013 October 1.
Published in final edited form as:
PMCID: PMC3432754
NIHMSID: NIHMS395516

Heritability in the genome-wide association era

Noah Zaitlen1,2,3,4 and Peter Kraft1,2,3,4

Abstract

Heritability, the fraction of phenotypic variation explained by genetic variation, has been estimated for many phenotypes in a range of populations, organisms, and time points. The recent development of efficient genotyping and sequencing technology has led researchers to attempt to identify the genetic variants responsible for the genetic component of phenotype directly via GWAS. The gap between the phenotypic variance explained by GWAS results and those estimated by from classical heritability methods has been termed the “missing heritability problem”. In this work, we examine modern methods for estimating heritability, which use the genotype and sequence data directly. We discuss them in the context of classical heritability methods, the missing heritability problem, and describe their implications for understanding the genetic architecture of complex phentoypes.

Keywords: heritability, genome-wide association study, prediction, polygenic models, linear mixed models

Introduction

Since their debut in 2005 genome-wide associations studies (GWAS) have identified thousands of single nucleotide polymorphisms (SNPs) associated with hundreds of different phenotypes (Hindorff et al. 2009). Despite this success, the total faction of the phenotypic variation explained for most phenotypes remains small relative to the published heritability estimates, which are estimated using the trait covariance among relatives (Eichler et al. 2010; Maher 2008; Manolio et al. 2009). This “missing heritability problem” raises questions about the methods used to estimate heritability as well as the genetic architecture of complex phenotypes.

Many explanations for the sources of missing heritability have been proposed including structural variations, gene-environment interactions, epistatic interactions, parent of origin effects, and errors in narrow-sense heritability estimates (Eichler et al. 2010; Manolio et al. 2009; Zuk et al. 2012). Of particular interest is the distribution of causal variants along the genome, their number, and their frequency spectrum. GWAS are particularly suited to capture common variants and so violation of the common disease common variant model may lead to missing heritability. In Fisher’s infinitesimal model there are expected to be a large number of rare variants associated with disease. The rare-allele model proposes that rare variants of large effect account for a significant fraction of phenotypic variation, and it has been proposed that these can give rise to synthetic association in common variants (Dickson et al. 2010; Gibson 2011).

Determining which combination of these hypotheses is correct and where the majority of phenotypic variation lays has significant implications for the future success of association studies as well as the clinical utility of genetic risk prediction. It is possible to decouple some of these proposed genetic architectures without directly identifying the causal variants themselves. For example, Wray et al. show the potential for comparing heritability and sibling relative risk estimates to determine the validity of a rare-variant model (Gibson 2011; Wray and Goddard 2010; Wray et al. 2011).

Recently, Yang et al. (2010) proposed using linear mixed models (LMM) to estimate a lower bound on the total narrow-sense heritability estimation from GWAS data as well determining how much of the phenotypic variation is due to SNPs in LD with those on genotyping platforms. The results of this approach have broad implications for the genetic architecture of phenotypes as well as the future success of GWAS.

In this work we examine the problem heritability estimation in the GWAS era and how it relates to the missing heritability problem. We briefly review the classical methods of heritability estimation and contrast them with relatively recent use of genotype data to estimate the component of heritability explained by common SNPs via the LMM approach. We discuss the relative merits of the different methods in terms of potential confounding factors as well as what they tell us about the distribution of causal variants and the potential returns of future GWAS. Finally, we discuss the prospects for using LMM to predict human traits, including disease risk.

Background

Heritability is a measure of the contribution of genetics to phenotype. Wright and Fisher formalized the concept by writing phenotypic variance as the sum of genetic variance and environmental variance, σP2=σG2+σε2. Broad sense heritability H2 is the ratio of total genetic variance to phenotypic variance H2=σG2σP2. This measure includes the effects of gene-gene interactions (epistatic effects) σI2, dominance effects σD2, and additive effects σg2 such that σG2=σg2+σD2+σI2. Narrow-sense heritability h2 measures just the additive contribution of genetic variation to phenotype h2=σg2σP2 (Falconer 1989; Lynch and Walsh 1998).

In this work we discuss estimates of narrow-sense heritability h2 unless stated otherwise. This is done because we focus on GWAS and the missing heritability problem. Most traditional estimates of heritability using the correlations among related individuals are presumed to estimate h2, although these estimates can be biased. For example, the classical estimate involving the regression of offspring trait values on the mean parental values does not include the dominance component of variance, but the epistatic component does contribute to the estimate. The epistatic component is typically (and perhaps incorrectly) assumed to be 0 for identifiability purposes (Falconer 1989; Zuk et al. 2012) GWAS estimates of individual-marker effect sizes are generally measured marginally, ignoring dominance and interaction effects, so the “bottom up” heritability estimates from GWAS (defined below) are narrow-sense estimates.

The additive model

In a GWAS, we are given a set of Ns SNPs S = {s1, s2, …, sNs} genotyped on Ni individuals with phenotypesY = y1, y2, …, yNi. Each genotype has value (0,1,2) and the genotypes of the jth individual are Gj = g1j, g2j, …, gNsj with minor allele frequencies p1, p2, …, pNs. Let C be the set of Nc causal SNPs, which along with environmental factors determine the phenotype of each individual. The ability of GWAS to identify the genetic contribution to trait variance will depend on the proportion of SNPs in C that are in S or in linkage disequilibrium with one or more SNPs in S.

In an additive model, the phenotype of each individual is defined by a sum of linear effects

yj=m+iCzijαi+εj
[Eqn. 1]

where Zij=gij-2pi2pi(1-pi) are the normalized genotypes, αi is the effect size of SNP si, εj is the environmental contribution, and Y is normalized to have variance 1. The environmental contribution is assumed to be normally distributed εj~N(0,σε2), and εj and εk are independently distributed for j≠k.

Marginal GWAS “Bottom Up” heritability estimation

The genetic variance in an additive model is computed by the sum of the squared effect sizes of the normalized genotypes σg2=iαi2 and the heritability is the ratio of the genetic variance to the total phenotypic variance h2=σg2σg2+σε2=σg2, where σε2 is environmental contribution to phenotype and σg2+σε2=σY2=1.

Given a GWAS one can compute an estimate of the genetic variance σ^g2 by using the effect size estimates from the markers with a pre-specified genome-wide significance level. This can be used to compute an estimate of the heritability hGWAS2=σ^g2σY2, which is defined as “bottom up” heritability estimation by Zuk et al. (2012).

Unfortunately, the full set of casual variants and their effect sizes are not known, so hGWAS2 will typically underestimate the total heritability. (The winner’s curse (Ioannidis 2007, 2008; Kraft 2008) and the inclusion of false-positive markers in the bottom-up estimate of genetic variance could in principle lead to an overestimate of heritability.) The difference between h2 and hGWAS2 is known as the “missing heritability”. It is the additive genetic variance not yet captured with GWAS or other methods of identifying associated variants.

Classical “Top Down” heritability estimation

The classical methods of heritability estimation are based on an intuitive concept. Phenotypes that are highly correlated amongst relatives in patterns consistent with Mendelian inheritance are more heritable than those that are weakly correlated amongst relatives. The formalization of this idea by Fisher (1918) and Wright (1921) is the foundation of heritability estimation.

Consider the correlation between the phenotype of two individuals in the additive model above:

cor(yj,yk)=cov(yj,yk)=cov(iCzijαi,iCzikαi)=ag2NciCcov(zij,zik)var(zi)=σg2KCausal[j,k].

KCausal is the genetic covariance matrix (Kang et al. 2010; Price et al. 2006; Yang et al. 2010) defined at the causal SNPs. The entry for element j,k in the matrix is:

KCausal,jk=1NCiC(gij-2p^i)(gik-2p^i)2p^i(1-p^i).

Until recently, the genotypes of individuals were unavailable and even now the set of causal variants is unknown, so alternative means of estimating KCausal are required. The classical and still widely used approach is to collect sets of related individuals from known pedigrees. The estimate of KCausal,jk is twice the kinship coefficient or 2 Φjk. Here Φjk is the probability that an allele drawn at random from j is identical by descent to a randomly drawn allele from k, and can be calculated from the known pedigree structure (Lange 2002). Many of the familiar values for Φjk such as Φjk= ¼ for full siblings assume founders share no alleles identical by descent, which may not be true in the presence of inbreeding or population substructure (Lange 2002; Powell et al. 2010). We call the matrix estimated from these pedigree-based estimates KPed, and it serves as an estimate of KCausal. Given this matrix, the problem of heritability estimation is reduced to estimating σ^g2 from the observed covariance of the phenotypes of the related individuals.

It is worth stressing that the entries in KPed are the sums of the expected cross products E[Zij Zik], while the actual covariance KCausal depends on the observed cross products zij zik. The actual covariance will vary around its expected value for most relative pairs. Visscher et al. (2006) proposed using an estimate of KCausal based on observed genotype data as a more accurate method for estimating heritiability using related individuals. For a sample of unrelated siblings pairs, the value of the entries in KPed are ½ for siblings and 0 otherwise. It follows (making the questionable assumption that the dominance, epistatic, and shared environmental components of variance are 0) that σ^g2 is twice the average correlation of the phenotype across the sib-pairs (Falconer 1989). If the average correlation amongst the normalized height of siblings in a population is 0.4 then the heritability estimate for height is 0.8.

Pedigree-based linear-mixed-model-estimates of heritability

When multiple classes of relationship are measured, as is the case in extended pedigrees, one can take advantage of all the relationships simultaneously via a linear mixed model (Lange 2002; Shaw 1987), where the 1×Nsubjects phenotype vector Y is distributed as a multivariate normal random variable with mean M and variance-covariance matrix Σ. The mean vector M captures the fixed effects of observed covariates (e.g. sex, age, or principal components of genetic variation). The variance-covariance matrix is:

=var(iCziαi)+var(ε)=KCausalσg2+Iσε2.

To estimate heritability via a linear mixed model the restricted maximum likelihood (REML) estimate of σ^g2 is computed and the heritability estimate is h^2=σ^g2σY2. REML is used to estimate the components of variance instead of maximum likelihood to avoid a bias introduced by the fixed effects (Shaw 1987). Since KCausal is not known KPed serves as an estimate. There are several algorithms for REML estimation (Kang et al. 2010; Lange 2002; Shaw 1987), but most are computationally expensive due to the cost of matrix inversion. Lippert et al. (2011) recently developed a fast method when the number of individuals exceeds the number of markers. When only one type of relationship is available (e.g. only sibs) then the REML estimator will give the same estimate as the covariance based approach described above.

There are many extensions to this linear mixed model approach that allow estimation different components of heritability. These include dominance effects (Lynch and Walsh 1998), gene-gene interaction (Yang et al. 2011a), the shared genetic basis of multiple phenotypes (cross heritability) (Boehnke et al. 1986; Deary et al. 2012; Lange and Boehnke 1983; Macgregor et al. 2006; Price et al. 2011), heritability from different genomic regions (Yang et al. 2011b), and the effects of shared environment (Lynch and Walsh 1998).

Top down heritability estimates are susceptible to a range of confounding factors, which can bias estimates. These include gene-environment correlations, selection, non-random mating, and inbreeding (Lynch and Walsh 1998; Visscher et al. 2008). Recently, Zuk et al. (2012) showed that certain types of epistatic interactions can inflate estimates of narrow-sense heritability.

Heritability in the GWAS Era

The availability of genotype data over large collections of individuals has opened up new approaches to estimating heritability. These methods apply the same linear mixed model method described above, but replace the KPed estimate of KCausal with estimates based on genotype data. We examine only the simple additive estimate of heritability, but each of the extensions listed above may be utilized for each estimator of KCausal.

Heritability using realized IBD

When genetic data are collected over the set of individuals in the study it is possible to estimate the total fraction of the genome shared identical by descent (IBD). Siblings for example do not share exactly 50% of their genome with each other (Visscher et al. 2006). Using the genetic data to estimate the fraction of genome shared IBD gives another means of estimating KCausal, which we call [K with circumflex]IBD. Provided that the IBD estimates are accurate, this matrix will be a better estimate of KCausal than KPed and therefore require fewer individuals to achieve a robust estimate of the heritability.

To illustrate this approach, Visscher et al. (2007) used the software package Merlin (Abecasis et al. 2002) to estimate IBD for a collection of twins to generate [K with circumflex]IBD from 791 autosomal markers and estimate several components of the heritability of height.

Heritability of common variants using observed genetic covariance

Recently, linear mixed models (LMMs) have been applied to GWAS data in an attempt to partition the “missing” heritability into variants tagged by GWAS SNPs (mostly common) and those that are not (mostly rare) (Yang et al. 2010). This use of the LMM links modern statistical approaches for high-dimensional data analysis (penalized regression) with classical models in statistical genetics (de los Campos et al. 2010).

This LMM approach uses the same REML based estimate of σ^g2 given above, but the matrix used is an empirical estimate of the genetic covariance (KGCV) instead of [K with circumflex]IBD or KPed (Yang et al. 2010; Yang et al. 2011a). This is similar to the “pseudo-heritability” estimate proposed by Kang et al. (2010). The relationship matrix KGCV is computed in a nearly identical way to KCausal, but because the set of causal variants C is unknown, the full set of genotyped SNPs in the GWAS is used directly as a proxy for KCausal.

This approach—which we refer to as the Yang-Visscher or LMM-KGCV approach—relies on the equivalence between the LMM,

yj=α+gj+εj,

with cov(gj, gk)= KGCV,jk σg2 and cov (εj, εk), and the random effects model,

yj=α+iSβizij+εj,

with the βi i.i.d. N(0, σg2/Ns). This equivalency provides the motivation for the claim that the LMM using KGCV estimates the proportion of additive genetic variance tagged by the GWAS markers. For unrelated individuals, causal variants that are not correlated with any of the zij for iεS (e.g. many rare variants) do not contribute to the estimate of σg2.

Note that hidden relatedness between individuals would bias σg2 since untagged causal variants would still tend to have the same correlation structure (related to KIBD) as the causal variants that are tagged thereby inflating the estimate of the portion of variability explained by the measured SNPs. Moreover, epistatic effects may also confound estimates of the additive genetic component σg2 if the sample contains closely related individuals. For distantly related individuals, such as those in most GWAS, the effects of dominance and epistasis on heritability estimates are much more attenuated.

The LMM just described is a special case of a general class of regression models defined by any similarity matrix K calculated from the GWAS data, with cov(gj, gk)= Kjk σg2 (de los Campos et al. 2010). Other choices of K may yield improved trait prediction, as they implicitly include non-additive effects (dominance, epistasis). However, precisely because the genetic variance includes non-additive components, the heritability estimates from these more general models can be difficult to interpret.

One of the advantages of this LMM approach using KGCV is that individuals may be selected randomly with respect to their environmental exposures preventing confounding from shared environments that can affect pedigree-based estimates. In addition, they can inform researchers about the potential success of future GWAS conducted on the phenotype of interest. The KGCV based estimates provide an upper bound on the total fraction of phenotypic variance explained by future GWAS on the same set of SNPs. They also provide a lower bound on the total narrow-sense heritability of the phenotype.

The application of KGCV to heritability estimation was proposed by Hayes et al. (2009) in the context of related individuals. In this case, KKGV will serve as an estimate of KIBD and therefore give an estimate of the total narrow-sense heritability. The crucial difference in the Yang-Visscher approach is the assumption that when the individuals are distantly related, the KGGV matrix provides no information apart from that contained in the genotyped SNPs. Thus the heritability estimate in this situation will be the narrow-sense heritability due exclusively to the SNPs in LD with those on the genotyping platform.

Violations of Model Assumptions

Each of the heritability estimation methods described above make different assumptions about the model generating phenotype. The estimates of heritability may be biased when these assumptions are broken.

While pedigree-based estimates of heritability have been examined for decades, the Yang-Visscher approach is a very recent development and there are many open questions about the factors that can affect these estimates of heritability. Here we give several examples of such factors and perform some simple experiments to examine their effects. These are in no way meant to be exhaustive or conclusive, but rather to inform the reader of potential issues.

Violations of Additivity

Zuk et al. (2012) show that when certain types of epistatic (gene-gene) interactions exist the estimates of heritability found from pedigree estimates, such as MZ versus DZ twins, will be upwardly biased. In this situation, bottom up estimates will never reach the top down estimate of heritability. They propose that this is a possible element of the “missing heritability problem”, and that the true narrow-sense heritability maybe substantially lower than current estimates for certain phenotypes (Zuk et al. 2012).

To examine this problem in the context of Yang-Visscher heritability estimates we simulated data sets using the epistatic “limiting pathway” models of Zuk et al. (Zuk et al. 2012), LP(1), LP(3), and LP(4). We simulated case-control genotypes and phenotypes of 2000 randomly ascertained unrelated individuals with 200 causal variants in each pathway, an effect size of 0.1, a minor allele frequency of 0.5, and prevalence of 50%. We computed a bottom up adjusted h2 estimate via linear regression as well as Yang-Visscher estimate of heritability, using all causal variants to estimate KGCV. The results are shown in Table 1 and demonstrate that the Yang-Visscher approach is not susceptible to confounding from epistatic interaction under the LP model of interaction. If closely related individuals were used then the Yang-Visscher estimate would be upwardly biased from the epistatic component of variance.

Table 1
Yang-Visscher and bottom-up estimates of heritability (and their standard error over 1000 replications) under three limiting pathway models of phenotype. For K>1 the pedigree-based top down estimates of heritability will be inflated. An LP(4) ...

Thus, the LMM estimates of heritability from unrelated individuals provide a benchmark to assess how much of the total narrow-sense heritability currently-known GWAS-identified trait markers explain—a benchmark that is not influenced by “phantom heritability” due to epistatic interactions. The ratio of the bottom-up additive genetic variance estimated using GWAS-identified markers to the LMM estimate of the additive genetic variance estimates the proportion of GWAS-identifiable markers that have been identified to date.

Violations of exchangeability

The Yang-Vissher approach assumes a polygenic model of disease in which many markers of small effect contribute to variance in genetic risk. Specifically, it assumes marker effect sizes are all drawn from the same normal distribution, β~N(0,σg2/Ns). There are however many diseases where there are outlier markers with strikingly different effects. For example, GWAS have identified dozens of markers associated with Type 1 Diabetes and rheumatoid arthritis, most of which have very small effects relative to the long-established risk variants in the MHC; for both of these diseases, the variants in the MHC have per-allele relative risks roughly three times larger than the relative risks for the GWAS-identified risk variants (Barrett et al. 2009; Stahl et al. 2010).

To examine the effect of such extreme variants this we simulated 1000 GWAS of 1500 individuals with a single causal variant. The genotypes at 1,000 marker loci (including the causal locus) were generated by random binomials with minor allele frequencies drawn uniformly between 0.05 and 0.5. The true heritability of the phenotype was 0.5 and the average estimate over the 1000 GWAS was 0.50, suggesting that violations of the infinitesimal model do not strongly effect estimates of heritability.

Addition of non-causal variants

For many phenotypes KGCV will contain a large number of variants unlinked to any causal variants. To examine the effect of these variants on the estimates of heritability we repeated the experiment above with 10 causal variants and 102, 103, and 5×103 additional independent (i.e. non-causal) variants. The true heritability of the phenotype was 0.5 and the mean heritability across the 1000 simulated GWAS was 0.50 in all studies. However, the standard deviations were 0.018, 0.025, and 0.067, showing that the effect of additional variants is to increase to standard error of the heritability estimates. The results did not change qualitatively for other values of h2. Other factors that affect the standard error of heritability estimates include the study sample size as well as the true heritability. Alternative disease models, such as mixtures of infinitesimals described by Park et al. (2011) have not yet been investigated in this context; the possibility that they lead to biased heritability estimates remains open.

Sample size considerations

To investigate the precision of LMM estimates of h2 using KGCV in real-world situations, we used GWAS data on 10,503 individuals from two European-ancestry cohorts, the Nurses’ Health Study and Health Professionals Follow-up Study. We simulated continuous phenotypes as a function of 500 SNPs, according to Eqn. 1, constraining the SNP effects so that the resulting phenotype had the desired heritability (h2=0.50, 0.25, 0.10). We also simulated a binary phenotype using the liability threshold model, with liability given by Eqn. 1, and prevalence 10%. We estimated h2 using the LMM approach as implemented in GCTA (Yang et al. 2011a), applied to a set of 151,019 markers (including the 500 causal variants) chosen to have low linkage disequilibrium (r2<0.2), varying the sample size from 1,000 to 10,503.

Results from single replicates are shown in Table 2. Precision increases roughly linearly with increasing log sample size. For sample sizes under 2,000, the 95% confidence intervals are wide (>0.40), and, for modest heritabilities (under 25%, consistent with the observed heritabilities for many complex traits), they include 0. This suggests that accurate estimation of narrow-sense heritabilities will require large sample sizes, on the order of 5,000 to 10,000 or more, at least as big as those needed to identify individual markers with modest effects. Published studies using the LMM-KGCV approach to estimate the narrow-sense heritability due to GWAS markers for continuous traits like height and body mass index used between 4,000 and 11,500 subjects (Yang et al. 2010; Yang et al. 2011b). Care must be taken when combining studies to reach such large sample sizes as this may introduce population substructure and corresponding environmental variation of non-genetic risk factors, potentially biasing estimates of heritability

Table 2
LMM estimates of narrow-sense heritability using KGCV and their standard errors for phenotypes simulated conditional on empirical GWAS data (described in text, under “Sample size considerations”). For binary phenotypes, heritabilities ...

Addition of markers in LD with the causal variants

The additive model assumes that all of the tested variants are independent. In reality, there is extensive LD between causal and non-causal variants in the genome. To examine the potential for LD to affect heritability estimates we repeated the experiment above with 4 causal variants, and 1 additional causal variant repeated 100 times simulating extensive LD for a particular SNP, and 104 non-causal variants. The true heritability was 0.5 and the average estimated heritability was 0.40 showing that LD patterns can significantly affect heritability estimates. We note that this is an extreme example meant to demonstrate the potential for bias. Yang et al simulated phenotypes over real GWAS data (i.e. with real LD patterns) and found estimates within two standard errors of the true heritability (Yang et al. 2010).

Distant/Cryptic Relatedness in the study

Provided that the individuals in a GWAS are unrelated, the matrix KGCV contains no information about SNPs out of LD with the genotyped SNPs. If the study contains related, individuals, however, the LMM estimate of heritability will contain some additional genetic variance due to variants not tagged by the GWAS SNPs. This is because KGCV is an unbiased estimate of KIBD, as we illustrate below. Since there are no truly “unrelated” individuals, any GWAS will contain a range of distantly related individuals. Yang et al suggest removing individuals with KGCV > 0.025 in the case of quantitative phenotypes and 0.05 for dichotomous phenotypes (Yang et al. 2011a).

We simulated 1000 pairs of individuals that shared 0.5, 0.1, 0.05, 0.025 of their genome IBD and compute KGCV for each pair. We repeated this experiment using 104, 105, and 106 SNPs. The results are presented in Table 3. In each case the mean estimate of IBD is close to the true IBD showing the KGCV is a good estimate of KIBD. The standard error is independent of the true IBD and decreases as a function of the number of independent SNPs.

Table 3
The genetic covariance between pairs of individuals with a range of IBDs, estimate from Ns SNPs. GCV is an unbiased estimate of IBD and the variance of the esimate in ()’s is function of the number of available SNPs.

For distantly related individuals, the signal from IBD will typically be small relative to the signal from the causal variants. Here a concern is confounding due to cryptic relatedness, where more closely related individuals tend to have similar trait values for non-genetic reasons (Kang et al. 2010). The influence of low levels of IBD in the Yang-Visscher approach remains an open question. It is possible to test explicitly for inflation due to relatedness, by simulating phenotypes over odd chromosomes and estimating heritability over even chromosomes (Visscher et al. 2010).

Population Substructure

Individuals from different populations have different minor allele frequencies as well different environmental exposures. In a case control study this can lead to significant confounding if there is a difference in the phenotypic mean between the populations, and is usually corrected with a principal component adjustment. Browning and Browning (2011) show that under certain extreme population differences this can lead to biases in heritability estimates. Yang et al show that using PC adjustment will mitigate this inflation. They also propose to estimate the effects of population stratification and cryptic relatedness by performing heritability estimation over each chromosome (Yang et al. 2011b). This procedure has not yet been examined in detail in the published literature.

Another type of population stratification arises when there is a difference in the phenotypic variance (but not necessarily mean phenotype) between the populations. In this case PCA will not adequately adjust for population substructure leading to inflation in standard GWAS (McPeek and Abney, 2008). Furthermore, the interpretation of heritability may ambiguous in this scenario, since each of the sub populations will likely have different heritability estimates.

Heritability is defined with respect to a population at a particular time. The heritability of lung cancer will be dramatically different between a population where some people smoke and a population of only non-smokers. Thus bottom-up GWAS heritability estimates and those from published heritability studies can only be compared if they come from the same population and are conducted at similar times.

Imputation and Rare Variants

Currently the Yang-Visscher approach has been performed using observed SNPs, genotyped using the same platform (Yang et al. 2010; Yang et al. 2011b). Given the success of imputation in the GWAS community one open questions is the possibility of leveraging external reference panels such as the HapMap to determine if additional signal lies within the additional SNPs genotypes in the panel. High-throughput sequencing data are becoming available and with them a large number of rare variants. The proper way to include dense maps of common markers and rare variants in heritability estimation—notably in light of the discussion of the impact of linkage disequilibrium patterns, above—is an area of current research.

Ascertainment and Case Control Phenotypes

For binary traits, the percent of trait variance captured by KGCV when analyzing a discontinuous 1-0 case-control phenotype in the LMM framework is not directly comparable to commonly-quoted heritabilites from some family studies (e.g. MZ-DZ twin comparisons), which are measures of the percentage of the underlying liability captured by inherited factors (Dempster and Lerner 1949; Visscher et al. 2008). Nor is there a simple link between individual-locus odds ratios and bottom-up estimates of the genetic variance and liability-scale heritability. (There is a simple relationship between the familial recurrence risk and the additive genetic component of variance on the log relative risk scale (Pharoah et al. 2002). Moreover, individual-marker allele frequencies and relative risk estimates from GWAS can be directly related to heritability on the log relative risk scale (Pharoah et al. 2008). See Wray and Goddard (2010) for a thorough discussion of the relationship between individual-marker relative risks and heritability measured on different scales.)

For the linear mixed model, the phenotypic variance captured by KGCV depends on disease prevalence and sampling scheme. By construction, the heritability of liability is independent of prevalence. When estimating the heritability of case-control phenotypes, the ascertainment strategy and prevalence of disease will affect the final heritability estimate. To address this issue, it is possible to transform the disease scale heritability estimate to a liability scale heritability estimate, which accounts for both ascertainment and prevalence (Dempster and Lerner 1949; Lee et al. 2011):

hliability2=hObs2F(1-F)ϕ(Φ-1(F))2F(1-F)P(1-P).

F is the prevalence ϕ is the normal pdf and Φ is the normal cdf and P is the proportion of cases in the sample. The justification for this elegant adjustment depends on a rather simple model for ascertainment, namely, that selection for inclusion is independent of all other covariates conditional on disease status. This will not be the case in many practical situations (e.g. matched case-control studies), where ascertainment depends on other factors that are usually associated with disease risk and may also be associated with genotype. The impact of violations of this assumption is unclear.

Phenotypic Prediction

The LMM using KGCV also offers a means of phenotypic prediction using the best linear unbiased predictors or BLUPs (Lynch and Walsh 1998). The expected trait value for a new individual (who did not contribute to the data set used to fit the LMM) is given by:

y^=α+g^=α+iSβ^izi.

This is similar to the “polygenic” models proposed by Purcell et al. (2009) and Evans et al. (Evans et al. 2009; Purcell et al. 2009), in that the predictor uses information contained in SNPs that do not reach the genome-wide significance threshold. But where the “polygenic” perform feature selection, only building predictors using markers with single-SNP (marginal) p-values below some threshold (often much larger than the stringent GWAS threshold), the LMM approach builds predictors using all available SNPs simultaneously. The LMM predictor is closely related to ridge regression, a penalized regression procedure that often outperforms variable selection procedures in terms of minimizing prediction error in new data sets (Harrell 2001; Hastie et al. 2001).

The accuracy of the LMM predictor is a function of narrow-sense heritability, the number of markers included in the LMM, the true genetic architecture, and the sample size in the data set used to fit the LMM. The sample size determines the accuracy with which βi can be estimated. The squared correlation between the LMM predictor and trait values in new observations is typically far smaller than the heritability estimate from the LMM (the theoretical maximum of the squared correlation); this is because of the variability in the estimated βis (Daetwyler et al. 2008; Visscher et al. 2010).

Summary

The Yang-Visscher approach to heritability estimation provides a means of estimating the contribution of SNPs in LD with those on genotyping platforms to the total phenotypic variation. In the context of GWAS these estimates answer questions about the genetic architecture of complex phenotypes. The growing number of GWAS identified loci, as well as their small effect sizes, has led to speculation about genetic models of disease.

There has been significant recent debate about the success or failure of GWAS (Eichler et al. 2010; Gibson 2011; Visscher et al. 2012). This has in turn reinvigorated the debate about the distribution of causal variants. Goldstein demonstrated the possibility for rare variants to induce synthentic associations (Dickson et al. 2010), and there have been several recent works discussing the common disease common variant, strong and weak rare variants, the infinitesimal, and other disease models (Gibson 2011).

There has also been speculation about the location of the “missing heritability” with discussions of parent of origin effects, epistatic interactions, gene-environment interactions, structural variation, and other cache’s of genetic variation not well captured by current GWAS or their analysis methods(Eichler et al. 2010; Visscher et al. 2012; Zuk et al. 2012).

The work of Yang and Visscher discussed here as well as other GWAS-based approaches (Lango Allen et al. 2010; So et al. 2011a; So et al. 2011b; Yang et al. 2011c) provide insights relevant to these questions. They estimate heritability restricted to a certain class of SNPs (i.e. those in LD with genotyped SNPs), are not confounded by many of the factors biasing traditional methods of heritability estimation, and are fundamentally different than bottom up methods. In principle these procedures could also be used to build phenotype prediction algorithms incorporating markers beyond the small number identified at genome-wide significance levels. However, very large sample sizes will be needed to obtain accurate estimates and precise prediction algorithms.

Acknowledgments

The authors thank Alkes Price, Eli Stahl and Dan Stram for helpful discussions, and Poorva Mudgal for programming support. NZ was supported by NIH fellowship 5T32ES007142-27, PK by NIH grant R21 DK084529.

References

  • Abecasis GR, Cherny SS, Cookson WO, Cardon LR. Merlin--rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet. 2002;30:97–101. [PubMed]
  • Barrett JC, Clayton DG, Concannon P, Akolkar B, Cooper JD, Erlich HA, Julier C, Morahan G, Nerup J, Nierras C, Plagnol V, Pociot F, Schuilenburg H, Smyth DJ, Stevens H, Todd JA, Walker NM, Rich SS. Genome-wide association study and meta-analysis find that over 40 loci affect risk of type 1 diabetes. Nat Genet. 2009;41:703–7. [PMC free article] [PubMed]
  • Boehnke M, Moll PP, Lange K, Weidman WH, Kottke BA. Univariate and bivariate analyses of cholesterol and triglyceride levels in pedigrees. Am J Med Genet. 1986;23:775–92. [PubMed]
  • Browning SR, Browning BL. Population structure can inflate SNP-based heritability estimates. Am J Hum Genet. 2011;89:191–3. author reply 193–5. [PubMed]
  • Daetwyler HD, Villanueva B, Woolliams JA. Accuracy of predicting the genetic risk of disease using a genome-wide approach. PLoS One. 2008;3:e3395. [PMC free article] [PubMed]
  • de los Campos G, Gianola D, Allison DB. Predicting genetic predisposition in humans: the promise of whole-genome markers. Nat Rev Genet. 2010;11:880–6. [PubMed]
  • Deary IJ, Yang J, Davies G, Harris SE, Tenesa A, Liewald D, Luciano M, Lopez LM, Gow AJ, Corley J, Redmond P, Fox HC, Rowe SJ, Haggarty P, McNeill G, Goddard ME, Porteous DJ, Whalley LJ, Starr JM, Visscher PM. Genetic contributions to stability and change in intelligence from childhood to old age. Nature 2012 [PubMed]
  • Dempster E, Lerner I. Heritability of threshold characters. Genetics. 1949;35:212–236. [PubMed]
  • Dickson SP, Wang K, Krantz I, Hakonarson H, Goldstein DB. Rare variants create synthetic genome-wide associations. PLoS Biol. 2010;8:e1000294. [PMC free article] [PubMed]
  • Eichler EE, Flint J, Gibson G, Kong A, Leal SM, Moore JH, Nadeau JH. Missing heritability and strategies for finding the underlying causes of complex disease. Nat Rev Genet. 2010;11:446–50. [PMC free article] [PubMed]
  • Evans DM, Visscher PM, Wray NR. Harnessing the information contained within genome-wide association studies to improve individual prediction of complex disease risk. Hum Mol Genet. 2009;18:3525–31. [PubMed]
  • Falconer DS. Introduction to quantitative genetics. 3. Longman Wiley; Burnt Mill, Harlow, Essex, England New York: 1989.
  • Fisher R. The correlation among relatives on the supposition of Mendelian inheritance. Trans Roy Soc Edinburgh. 1918;52:399–433.
  • Gibson G. Rare and common variants: twenty arguments. Nat Rev Genet. 2011;13:135–45. [PubMed]
  • Harrell F. Regression modeling strategies. Springer; New York: 2001.
  • Hastie T, Tibshirani R, Friedman J. The Elements of Statistical Learning. Springer; New York: 2001.
  • Hayes BJ, Visscher PM, Goddard ME. Increased accuracy of artificial selection by using the realized relationship matrix. Genet Res (Camb) 2009;91:47–60. [PubMed]
  • Hindorff LA, Sethupathy P, Junkins HA, Ramos EM, Mehta JP, Collins FS, Manolio TA. Potential etiologic and functional implications of genome-wide association loci for human diseases and traits. Proc Natl Acad Sci U S A. 2009;106:9362–7. [PubMed]
  • Ioannidis JP. Non-replication and inconsistency in the genome-wide association setting. Hum Hered. 2007;64:203–13. [PubMed]
  • Ioannidis JP. Why most discovered true associations are inflated. Epidemiology. 2008;19:640–8. [PubMed]
  • Kang HM, Sul JH, Service SK, Zaitlen NA, Kong SY, Freimer NB, Sabatti C, Eskin E. Variance component model to account for sample structure in genome-wide association studies. Nat Genet. 2010;42:348–54. [PMC free article] [PubMed]
  • Kraft P. Curses--winner’s and otherwise--in genetic epidemiology. Epidemiology. 2008;19:649–51. discussion 657–8. [PubMed]
  • Lange K. Mathematical and statistical methods for genetic analysis. Springer; New York: 2002.
  • Lange K, Boehnke M. Extensions to pedigree analysis. IV. Covariance components models for multivariate traits. Am J Med Genet. 1983;14:513–24. [PubMed]
  • Lango Allen H, Estrada K, Lettre G, Berndt SI, Weedon MN, Rivadeneira F, Willer CJ, Jackson AU, Vedantam S, Raychaudhuri S, Ferreira T, Wood AR, Weyant RJ, Segre AV, Speliotes EK, Wheeler E, Soranzo N, Park JH, Yang J, Gudbjartsson D, Heard-Costa NL, Randall JC, Qi L, Vernon Smith A, Magi R, Pastinen T, Liang L, Heid IM, Luan J, Thorleifsson G, Winkler TW, Goddard ME, Sin Lo K, Palmer C, Workalemahu T, Aulchenko YS, Johansson A, Zillikens MC, Feitosa MF, Esko T, Johnson T, Ketkar S, Kraft P, Mangino M, Prokopenko I, Absher D, Albrecht E, Ernst F, Glazer NL, Hayward C, Hottenga JJ, Jacobs KB, Knowles JW, Kutalik Z, Monda KL, Polasek O, Preuss M, Rayner NW, Robertson NR, Steinthorsdottir V, Tyrer JP, Voight BF, Wiklund F, Xu J, Hua Zhao J, Nyholt DR, Pellikka N, Perola M, Perry JR, Surakka I, Tammesoo ML, Altmaier EL, Amin N, Aspelund T, Bhangale T, Boucher G, Chasman DI, Chen C, Coin L, Cooper MN, Dixon AL, Gibson Q, Grundberg E, Hao K, Juhani Junttila M, Kaplan LM, Kettunen J, Konig IR, Kwan T, Lawrence RW, Levinson DF, Lorentzon M, McKnight B, Morris AP, Muller M, Suh Ngwa J, Purcell S, Rafelt S, Salem RM, Salvi E, et al. Hundreds of variants clustered in genomic loci and biological pathways affect human height. Nature. 2010;467:832–8. [PMC free article] [PubMed]
  • Lee SH, Wray NR, Goddard ME, Visscher PM. Estimating missing heritability for disease from genome-wide association studies. Am J Hum Genet. 2011;88:294–305. [PubMed]
  • Lippert C, Listgarten J, Liu Y, Kadie CM, Davidson RI, Heckerman D. FaST linear mixed models for genome-wide association studies. Nat Methods. 2011;8:833–5. [PubMed]
  • Lynch M, Walsh B. Genetics and analysis of quantitative traits. Sinauer; Sunderland, Mass: 1998.
  • Macgregor S, Cornes BK, Martin NG, Visscher PM. Bias, precision and heritability of self-reported and clinically measured height in Australian twins. Hum Genet. 2006;120:571–80. [PubMed]
  • Maher B. Personal genomes: The case of the missing heritability. Nature. 2008;456:18–21. [PubMed]
  • Manolio TA, Collins FS, Cox NJ, Goldstein DB, Hindorff LA, Hunter DJ, McCarthy MI, Ramos EM, Cardon LR, Chakravarti A, Cho JH, Guttmacher AE, Kong A, Kruglyak L, Mardis E, Rotimi CN, Slatkin M, Valle D, Whittemore AS, Boehnke M, Clark AG, Eichler EE, Gibson G, Haines JL, Mackay TF, McCarroll SA, Visscher PM. Finding the missing heritability of complex diseases. Nature. 2009;461:747–53. [PMC free article] [PubMed]
  • Park JH, Gail MH, Weinberg CR, Carroll RJ, Chung CC, Wang Z, Chanock SJ, Fraumeni JF, Jr, Chatterjee N. Distribution of allele frequencies and effect sizes and their interrelationships for common genetic susceptibility variants. Proc Natl Acad Sci U S A. 2011;108:18026–31. [PubMed]
  • Pharoah PD, Antoniou A, Bobrow M, Zimmern RL, Easton DF, Ponder BA. Polygenic susceptibility to breast cancer and implications for prevention. Nat Genet. 2002;31:33–6. [PubMed]
  • Pharoah PD, Antoniou AC, Easton DF, Ponder BA. Polygenes, risk prediction, and targeted prevention of breast cancer. N Engl J Med. 2008;358:2796–803. [PubMed]
  • Powell JE, Visscher PM, Goddard ME. Reconciling the analysis of IBD and IBS in complex trait studies. Nat Rev Genet 2010 [PubMed]
  • Price AL, Helgason A, Thorleifsson G, McCarroll SA, Kong A, Stefansson K. Single-tissue and cross-tissue heritability of gene expression via identity-by-descent in related or unrelated individuals. PLoS Genet. 2011;7:e1001317. [PMC free article] [PubMed]
  • Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D. Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006;38:904–9. [PubMed]
  • Purcell SM, Wray NR, Stone JL, Visscher PM, O’Donovan MC, Sullivan PF, Sklar P. Common polygenic variation contributes to risk of schizophrenia and bipolar disorder. Nature. 2009;460:748–52. [PubMed]
  • Shaw R. Maximum-likelihood approaches applied to quantitative genetics of natural populations. Evolution. 1987;41:812–826.
  • So HC, Gui AH, Cherny SS, Sham PC. Evaluating the heritability explained by known susceptibility variants: a survey of ten complex diseases. Genet Epidemiol. 2011a;35:310–7. [PubMed]
  • So HC, Li M, Sham PC. Uncovering the total heritability explained by all true susceptibility variants in a genome-wide association study. Genet Epidemiol. 2011b;35:447–56. [PubMed]
  • Stahl EA, Raychaudhuri S, Remmers EF, Xie G, Eyre S, Thomson BP, Li Y, Kurreeman FA, Zhernakova A, Hinks A, Guiducci C, Chen R, Alfredsson L, Amos CI, Ardlie KG, Barton A, Bowes J, Brouwer E, Burtt NP, Catanese JJ, Coblyn J, Coenen MJ, Costenbader KH, Criswell LA, Crusius JB, Cui J, de Bakker PI, De Jager PL, Ding B, Emery P, Flynn E, Harrison P, Hocking LJ, Huizinga TW, Kastner DL, Ke X, Lee AT, Liu X, Martin P, Morgan AW, Padyukov L, Posthumus MD, Radstake TR, Reid DM, Seielstad M, Seldin MF, Shadick NA, Steer S, Tak PP, Thomson W, van der Helm-van Mil AH, van der Horst-Bruinsma IE, van der Schoot CE, van Riel PL, Weinblatt ME, Wilson AG, Wolbink GJ, Wordsworth BP, Wijmenga C, Karlson EW, Toes RE, de Vries N, Begovich AB, Worthington J, Siminovitch KA, Gregersen PK, Klareskog L, Plenge RM. Genome-wide association study meta-analysis identifies seven new rheumatoid arthritis risk loci. Nat Genet. 2010;42:508–14. [PubMed]
  • Visscher PM, Brown MA, McCarthy MI, Yang J. Five years of GWAS discovery. Am J Hum Genet. 2012;90:7–24. [PubMed]
  • Visscher PM, Hill WG, Wray NR. Heritability in the genomics era--concepts and misconceptions. Nat Rev Genet. 2008;9:255–66. [PubMed]
  • Visscher PM, Macgregor S, Benyamin B, Zhu G, Gordon S, Medland S, Hill WG, Hottenga JJ, Willemsen G, Boomsma DI, Liu YZ, Deng HW, Montgomery GW, Martin NG. Genome partitioning of genetic variation for height from 11,214 sibling pairs. Am J Hum Genet. 2007;81:1104–10. [PubMed]
  • Visscher PM, Medland SE, Ferreira MA, Morley KI, Zhu G, Cornes BK, Montgomery GW, Martin NG. Assumption-free estimation of heritability from genome-wide identity-by-descent sharing between full siblings. PLoS Genet. 2006;2:e41. [PubMed]
  • Visscher PM, Yang J, Goddard ME. A commentary on ’common SNPs explain a large proportion of the heritability for human height’ by Yang et al. (2010) Twin Res Hum Genet. 2010;13:517–24. [PubMed]
  • Wray NR, Goddard ME. Multi-locus models of genetic risk of disease. Genome Med. 2010;2:10. [PMC free article] [PubMed]
  • Wray NR, Purcell SM, Visscher PM. Synthetic associations created by rare variants do not explain most GWAS results. PLoS Biol. 2011;9:e1000579. [PMC free article] [PubMed]
  • Wright S. Systems of Mating. Genetics. 1921;6:111–78. [PubMed]
  • Yang J, Benyamin B, McEvoy BP, Gordon S, Henders AK, Nyholt DR, Madden PA, Heath AC, Martin NG, Montgomery GW, Goddard ME, Visscher PM. Common SNPs explain a large proportion of the heritability for human height. Nat Genet. 2010;42:565–9. [PMC free article] [PubMed]
  • Yang J, Lee SH, Goddard ME, Visscher PM. GCTA: A Tool for Genome-wide Complex Trait Analysis. Am J Hum Genet. 2011a;88:76–82. [PubMed]
  • Yang J, Manolio TA, Pasquale LR, Boerwinkle E, Caporaso N, Cunningham JM, de Andrade M, Feenstra B, Feingold E, Hayes MG, Hill WG, Landi MT, Alonso A, Lettre G, Lin P, Ling H, Lowe W, Mathias RA, Melbye M, Pugh E, Cornelis MC, Weir BS, Goddard ME, Visscher PM. Genome partitioning of genetic variation for complex traits using common SNPs. Nat Genet. 2011b;43:519–25. [PubMed]
  • Yang J, Weedon MN, Purcell S, Lettre G, Estrada K, Willer CJ, Smith AV, Ingelsson E, O’Connell JR, Mangino M, Magi R, Madden PA, Heath AC, Nyholt DR, Martin NG, Montgomery GW, Frayling TM, Hirschhorn JN, McCarthy MI, Goddard ME, Visscher PM. Genomic inflation factors under polygenic inheritance. Eur J Hum Genet. 2011c;19:807–12. [PMC free article] [PubMed]
  • Zuk O, Hechter E, Sunyaev SR, Lander ES. The mystery of missing heritability: Genetic interactions create phantom heritability. Proc Natl Acad Sci U S A 2012 [PubMed]