Volume 5 Supplement 9
Estimating heritability using family and unrelated individuals data
© Shetty et al; licensee BioMed Central Ltd. 2011
Published: 29 November 2011
For the family data from Genetic Analysis Workshop 17, we obtained heritability estimates of quantitative traits Q1 and Q4 using the ASSOC program in the S.A.G.E. software package. ASSOC is a family-based method that estimates heritability through the estimation of variance components. The covariate-adjusted mean heritability was 0.650 for Q1 and 0.745 for Q4. For the unrelated individuals data, we estimated the heritability of Q1 as the proportion of total variance that can be accounted for by all single-nucleotide polymorphisms under an additive model. We examined a novel ordinary least-squares method, a naïve restricted maximum-likelihood method, and a calibrated restricted maximum-likelihood method. We applied the different methods to all 200 replicates for Q1. We observed that the ordinary least-squares method yielded many estimates outside the interval [0, 1]. The restricted maximum-likelihood estimates were more stable than the ordinary least-squares estimates. The naïve restricted maximum-likelihood method yielded an average estimate of 0.462 ± 0.1, and the calibrated restricted maximum-likelihood method yielded an average of 0.535 ± 0.121. Our results demonstrate discrepancies in heritability estimates using the family data and the unrelated individuals data.
The heritability of a trait is usually calculated using family data. The identified genetic variants found through genome-wide association studies account for only a small portion of heritability for most complex traits  compared with the heritability estimated from family data. This discrepancy in the estimates, the missing heritability, is of great interest because the sources of this difference are still unknown . Recently, Yang et al. , using a novel statistical method, suggested that the missing heritability can be recovered using the genome-wide associations of unrelated samples . Because the Genetic Analysis Workshop 17 (GAW17) data set included family data and unrelated individuals data for the same traits , we estimated the “heritability” of Q1 with the unrelated individuals data and estimated the “heritability” of Q1 and Q4 with the family data.
For the family data, the heritability is the narrow sense heritability, estimated with the polygenetic effect model; we conducted a George-Elston transformation  to estimate the heritability. For the unrelated data, the heritability is the proportion of the total variance in a phenotype that can be described by all single-nucleotide polymorphisms (SNPs) under an additive model; we estimated it using the ordinary least-squares (OLS) method suggested by Yang et al. , a naïve restricted maximum-likelihood (REML) method, and a calibrated REML method. In all our analyses, the heritability estimates were obtained after adjustments for age, sex, and smoking status.
PEDINFO and ASSOC
Heritability estimates for Q1 and Q4 using the family data
OLS and REML estimates
where consists of trait values of n unrelated individuals, , where x i = (x i 1, …, x i 3) consists of the sex, age, and smoking status of the ith individual, respectively, consists of the effect sizes of the covariates, summarizes genotype data of m unknown causal variants such that z i = (z i 1, …, z im ), and , or if the genotype of the ith individual at the jth causal variant is aa, aA, or AA, respectively, f j is the frequency of allele A and Here the prime indicates the transpose of a vector or matrix.
where is the residual variance, I n is the identity matrix of order n,
Because this estimate is based on OLS, it does not need iteration. By replacing G with A and A* in the model given by , we can constructed the naïve and calibrated REML estimates by maximizing the likelihood of
Heritability estimates using the family data
In the family data, 697 individuals (202 founders and 495 nonfounders) form eight pedigrees. The pedigrees all have four generations of family members and a mean size of 87.13 individuals (range, 73–128). The pedigrees include 194 sibships with a mean size of 2.55 (range, 1–9). In the four randomly selected replicates, the heritability estimates for Q1 ranged from 0.608 to 0.698 with an average of 0.650; the heritability estimates for Q4 ranged from 0.687 to 0.773 with an average of 0.745 (Table 1).
Heritability estimates using the unrelated individuals data
We were unable to obtain REML estimates for Q4 because the convergence rate of the REML was extremely slow. We found that the convergence of the REML failed because no SNP contributed any phenotypic variation in the simulated model .
Discussion and conclusions
In our analyses, we estimated heritability using both the family data for Q1 and Q4 and the unrelated individuals data for Q1. The heritability estimates for Q1 and Q4 using the family data appeared stable and reasonable. In the simulation, Q1 has a heritability of 0.575, where 0.135 is due to the 39 causal SNPs and 0.440 is due to a polygenic component, and Q4 has a heritability of 0.70 resulting from a polygenic effect. The mean heritability estimates for Q1 and Q4 with the family data were 0.650 and 0.745, respectively.
The heritability estimates using the unrelated individuals data seem less reasonable. The OLS method did not work well for the GAW17 unrelated individuals data because the method was designed for genome-wide common SNPs. In the GAW17 unrelated individuals data, most of the SNPs are rare variants and a few of them are causal variants. The genetic relationships estimated using many rare variants may be unreliable, and this results in the instability of the OLS estimates. The REML approaches appear to be more stable than the OLS method for Q1. We observed that the heritability estimates using the unrelated individuals data were less than those using the family data on average. For example, the mean of the heritability estimates for Q1 for the unrelated individuals data was 0.462 (by naïve REML), which was 0.188 less than the mean for the family data. One possible reason is that the polygenic component (0.440) in Q1 is not due to any SNPs in the GAW17 sequence data set. We should not be able to uncover the polygenic effect using unrelated samples. However, the mean naïve REML estimate (0.462) is much larger than the heritability because of the causal SNPs (0.135). The reason is that we used all 24,487 SNPs to estimate the relationships among individuals. There might be other sources contributing to the heritability estimates.
Finally, we failed to estimate the heritability for Q4 using the unrelated samples because of the convergence problem, which was the result of no genotyped exonic SNPs in the data contributing to the phenotypic variation.
The Genetic Analysis Workshop is supported by National Institutes of Health (NIH) grant R01 GM031575 from the National Institute of General Medical Sciences. This work was supported by National Cancer Institute grant P30 CAD43703 and NIH grants HL074166, HL086718, R01 HG003054 and R01 HG005854. Some of the results of this paper were obtained by using the program package S.A.G.E., which is supported by U.S. Public Health Service Resource Grant RR03655 from the National Center for Research Resources. We thank the other members of Xiaofeng Zhu’s laboratory for their critiques and comments.
This article has been published as part of BMC Proceedings Volume 5 Supplement 9, 2011: Genetic Analysis Workshop 17. The full contents of the supplement are available online at http://www.biomedcentral.com/1753-6561/5?issue=S9.
- Manolio TA, Collins FS, Cox NJ, Goldstein DB, Hindorff LA, Hunter DJ, McCarthy MI, Ramos EM, Cardon LR, Chakravarti A, et al: Finding the missing heritability of complex diseases. Nature. 2009, 461: 747-753. 10.1038/nature08494.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang J, Benyamin B, McEvoy BP, Gordon S, Henders AK, Nyholt DR, Madden PA, Heath AC, Martin NG, Montgomery GW, et al: Common SNPs explain a large proportion of the heritability for human height. Nat Genet. 2010, 42: 565-569. 10.1038/ng.608.PubMed CentralView ArticlePubMedGoogle Scholar
- Almasy LA, Dyer TD, Peralta JM, Kent JW, Charlesworth JC, Curran JE, Blangero J: Genetic Analysis Workshop 17 mini-exome simulation. BMC Proc. 2011, 5 (suppl 9): S2-10.1186/1753-6561-5-S9-S2.PubMed CentralView ArticlePubMedGoogle Scholar
- George V, Elston RC: Generalized modulus power transformations. Commun Stat Theory Meth. 1988, 17: 2933-2952. 10.1080/03610928808829781.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.