PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Nat Genet. Author manuscript; available in PMC 2010 August 3.
Published in final edited form as:
Published online 2008 January 13. doi:  10.1038/ng.74
PMCID: PMC2914680
NIHMSID: NIHMS92527

Common variants in the GDF5-BFZB region are associated with variation in human height

Abstract

Identifying genetic variants that influence human height will further our understanding of skeletal growth and development. A number of rare genetic variants have been convincingly and reproducibly associated with height in Mendelian syndromes, and common variants in HMGA2 were recently found to be associated with variation in height in the general population1. Here, we report genome-wide association analyses of 6,669 individuals from Finland and Sardinia, using genotyped and imputed markers, and follow-up in an additional 28,801 individuals. We show that common variants in the osteoarthritis-associated2 GDF5-BFZB locus are responsible for variation in height (estimated additive effect of 0.44 cm, overall p<10−15). Our results suggest a link between the genetic basis of height and osteoarthritis, potentially mediated through alterations in bone growth and development.

Height is a complex trait influenced by genes and a variety of environmental factors, including diet and the prenatal environment3. Heritability estimates suggest that ≥ 80% of variation in height may be genetically determined3,4. Rare mutations with large effects on height in Mendelian syndromes have been identified in FBN1, FGFR3, GH1, EVC1, and GPC3 (MIM 154700, 134934, 262400, 604831, and 312870). Despite the high heritability, numerous candidate gene and linkage studies to identify loci influencing height in individuals of “normal stature” have been inconclusive5. Overall, variation in human height is likely to be polygenic and heterogeneous. Recently, the first GWAS of height1 identified common variants in HMGA2 associated with normal variation in height, both in adults (p=4 × 10−16) and children (p=1 × 10−6); the variants account for a small fraction (~0.3%) of the overall variation in height.

To identify additional genetic variants associated with height, we analyzed genome-wide SNP data on 2,371 Finns from the FUSION study6 and 4,298 Sardinians from the SardiNIA study7 (Table 1). The two samples were originally genotyped with distinct sets of markers. We used genotype imputation methods6,8 to facilitate comparison of the two studies and evaluate association between height and ~2.28 million common genetic variants. After verifying the overall accuracy of imputed genotypes in a subset of markers, we conducted within-study analyses using a rapid variance components-based association test9 and then carried out a meta-analysis of the two studies (Supplementary Figure 1).

Table 1
Subject characteristics. Values are reported as means (standard deviation)

Our results provided confirmatory evidence of association with rs1042725 and rs7968682, the recently reported common variants in HMGA2 (p=.031 and .0093, respectively, both in the same direction as the original report: Supplementary Table 1)1. The five loci showing the most significant evidence of association in our study are shown in Supplementary Table 2. To our knowledge, common variants in these loci have not been associated with height previously.

The genes near our strongest signal (p<2 × 10−7), located on chromosome 20, have a plausible biological role in human height. Rare variants in growth differentiation factor 5 (GDF5) have been associated with disorders of skeletal development (see below), and BFZB (also known as C20orf44 and UQCC) encodes a ZIC-binding protein repressed by basic fibroblast growth factor10.

We pursued the chromosome 20 signal because it was the single best result in our initial scan, the surrounding region accounts for 40 of the 50 lowest p-values in our meta-analysis, and it overlaps an osteoarthritis susceptibility locus2. We focused our follow-up efforts for this locus on SNP rs6060369 (p = 9.7 × 10−7 in the initial meta-analysis, Table 2) because it exhibited the strongest evidence for association among all Affymetrix Mapping 500K SNPs, which had been genotyped in the majority of our GWAS samples; we favored genotyped over imputed SNPs for initial follow-up. The SNP most significantly associated with height in the meta-analysis of imputed HapMap SNPs was rs725908, which is in strong LD with rs6060369 (r2 = 0.96, Supplementary Table 2). SNP rs6060369 was initially imputed in the FUSION GWAS, but direct genotyping sustained strong evidence of association (meta-analysis p=1.5 × 10−6). In the GWAS samples, each copy of the C allele at rs6060369 was associated with an increase in height of 0.3 to 0.7 cm (accounting for 0.3–0.6% of the variance in height, after adjustment for age and sex).

Table 2
Association between rs6060369 and height.

Motivated by previous reports of sex differences at osteoarthritis-associated loci11, we carried out an analysis stratified by sex. The stratified analyses show strong evidence of association in both males and females with no evidence of heterogeneity (Supplementary Table 3). We did not detect significant association of the SNP with other anthropometric measures (p>.05 for weight, body mass index, waist circumference, hip circumference, and waist-to-hip ratio).

We next tested the association of rs6060369 with height in nine follow-up samples, comprising 23,684 individuals of European ancestry and 3,860 African American individuals (Table 2). Six of the samples provided significant evidence of association (p<.05) and the other three showed a trend (p<.20) in the same direction (Table 2). The p-value for association in all 27,544 follow-up samples was 1.05 × 10−11, and in all 34,213 GWAS and follow-up samples combined was 2.22 × 10−16 (Figure 1). In the follow-up samples, each copy of the C allele at rs6060369 was associated with an increase in height of 0.2 – 0.6 cm (Table 2), and overall we estimate the SNP explains 0.3–0.5% of the variance in height, both in males and females (Supplementary Table 3). Further independent evidence for association between rs6060369 and adult height comes from the 1958 British Birth Cohort for which rs6060369 is associated with height measured at 44–45 years of age (p = .0046, explaining 0.5% of the variance) (http://www.b58cgene.sgul.ac.uk/, accessed October 2007).

Figure 1
Evidence of (A) association with height and (B) linkage disequilibrium around GDF5 and BFZB/C20orf44/UQCC

Our association signal lies in a 136 kb stretch of linkage disequilibrium (LD) that contains two genes, GDF5 and BFZB (Figure 1, Supplementary Table 4). BFZB is expressed in differentiating chondrocytes12, starting at early stages of mesenchymal cell proliferation13. Studies in mouse embryonic stem cells have shown that BFZB is down-regulated on addition of FGF2 (bFGF)10, which acts in concert with bone morphogenic proteins and several Hox gene products to initiate and promote morphogenesis and growth of the skeleton. Thus, BFZB appears to be involved in a network of FGF2-regulated growth control. GDF5 is a member of the TGF-beta superfamily involved in bone growth and differentiation, both in adult and embryonic tissues14,15. GDF5 expression is typically restricted to the primordial cartilage of long bones, with little expression in the vertebrae and ribs14. Mutations in GDF5 are associated with several disorders of skeletal development (MIM 201250, 200700, 112600, 113100, 228900, 185800 and 186500). Recombinant human GDF5 has been shown to restore vertebral disc height in a rabbit disc degeneration model, perhaps through enhanced extracellular matrix production15. Other nearby genes do not appear to be involved in chondrocyte differentiation, bone growth, or development, but the locus that includes GDF5 and BFZB is also interesting because it was highlighted as showing very strong evidence for selection in a genome-wide search for regions that underwent recent selection16. The target of selection is presently unknown.

A SNP located in the 5’ untranslated region of GDF5, rs143383, is strongly associated with osteoarthritis2,17 and is estimated to be in very strong LD with rs6060369 in the HapMap, FUSION, and SardiNIA samples (r2=.83–.90). The SNP appears to influence GDF5 expression2,17 and, we reasoned, might be a causal variant. Thus, we analyzed this SNP in our screening samples and a subset of our follow-up samples. The rs143383 A allele previously associated with increased risk of osteoarthritis was associated with decreased height in our studies (p=5.01 × 10−12 versus p=4.08 × 10−9 for rs6060369 in the same subset of samples; Figure 1, Table 3). Analysis stratified by gender shows strong association in both males and females (Supplementary Table 5).

Table 3
Association between rs143383 and height.

Interestingly, the ARIC African American samples that showed only a trend toward association with rs6060369 (p=.17) showed significant evidence of association with rs143383 (p=.011), illustrating the utility of studying different ancestral groups in the fine-mapping of complex disease genes18,19. In the ARIC African American samples, even when rs6060369 was included in a regression model, rs143383 remained marginally associated with height (p=.034, estimated additive effect = 0.579 cm), while the association of rs6060369 disappeared (p=.92, estimated additive effect = −0.019 cm) when conditioning on rs143383. These results suggest that GDF5 is more likely to influence height, although other nonsynonymous SNPs present in GDF5 and BFZB may affect function instead or as well.

Miyamoto and colleagues demonstrated that the A allele of rs143383 is associated with decreased GDF5 transcriptional activity in chondrogenic cells2. Decreased expression of GDF5 would logically lead to decreased limb bone growth, consistent with decreased height, as we observed. Decreased transcription of GDF5 may influence the amount of cartilage of the vertebrae, limb proportions, or joint angles, resulting both in a modest decrease in stature and susceptibility to osteoarthritis.

To evaluate the impact of osteoarthritis as a confounding factor, we repeated the association analysis restricted to younger individuals (age < 40). In SardiNIA, we analyzed 1,964 individuals and confirmed the association (p=.0018 for rs6060369, and p=.015 for rs143383) with effect size similar to estimates as for the combined sample (0.70 cm per copy of the C allele for rs6060369). In the Old Order Amish, the younger subgroup of 891 individuals showed a trend toward association of allele C with increased height (0.60 cm per copy of the C allele at rs6060369), but no significant association (p = .86), likely due to low statistical power.

To compare the evidence of association with length of long bones compared to vertebrae and skull, we tested rs6060369 and rs143383 for evidence of association with sitting height, which was measured in the ARIC and BWHHS studies. In ARIC European Americans and the BWHHS British sample, similar evidence of association was observed for both standing and sitting height. In ARIC African Americans, only rs143383 was significantly associated (p<.05) with height, and it was associated only with standing height, not with sitting height (Table 3), perhaps suggesting a stronger effect on long bones than on vertebrae.

Multiple regression analysis of our data suggests that a single common variant in the region may underlie the evidence of association. Specifically, multiple regression analysis in GWAS samples showed that after including rs6060369, rs143383, or rs725908 as a covariate, other association signals in the region become non-significant. One of these common variants or another nearby unmeasured variant in LD may influence height through expression of GDF5 2,17. However, either or both GDF5 and BFZB could be affected. Thoroughly evaluating the contribution of this locus to human height will require re-sequencing to catalog all genetic variants in the region and genotyping to evaluate their effects.

Combined, the variants identified here and previously reported in HMGA2 account for <1% of the variance in height, so that most of the 80% of variation in height estimated to be genetic remains unexplained. Our GWAS provides suggestive evidence for several other loci influencing height. For example, after excluding SNPs within 250 kb from GDF5, we observed a slight excess of SNPs with p-value <10−5 (38 observed versus 23 expected, Supplementary Figure 2). Still, it appears likely that many of the common variants influencing height will have only small effects. Follow-up of additional SNPs in larger meta-analyses will be necessary to define these variants20, which may also be relevant not only to normal variation in height but also to musculoskeletal or other diseases.

Methods

Study Subjects

Informed consent was obtained from all study participants and, in addition, ethics approval was obtained from the participating institutions.

FUSION GWAS

The Finland-United States Investigation of NIDDM Genetics (FUSION) study GWAS included 1,161 Finnish type 2 diabetes (T2D) cases, 1,174 normal glucose tolerant (NGT) controls, and 122 offspring of case/control pairs6. Cases and controls were approximately frequency matched, taking into account age, sex, and birth province within Finland6. For the height analysis, our sample consisted of 1,084 T2D individuals and 1,287 NGT individuals with height measurements from clinical exams. Samples were genotyped with the Illumina Infinium II HumanHap300 BeadChip6 and with an Illumina GoldenGate Custom Panel designed to improve genomic coverage around T2D candidate genes. The two imputed SNPs selected for additional follow-up were subsequently genotyped using the TaqMan allelic discrimination assay (Applied Biosystems, Foster City, CA).

SardiNIA GWAS

The SardiNIA GWAS examined a total of 4,305 related individuals participating in a longitudinal study of aging-related quantitative traits in the Ogliastra region of Sardinia, Italy. These individuals are distributed across 1,133 inter-related sibships, each with an average of 3.9 phenotyped individuals. For this study, we analyzed phenotypes for 4,298 individuals, excluding 4 cases of short stature due to Morquio Syndrome (MIM 253000) and 3 individuals for whom height measurements were not available. Among the individuals examined, 1,412 were genotyped with the Affymetrix Mapping 500K Array Set. Taking advantage of the relatedness among individuals in the SardiNIA sample, we conducted a second round of computational analysis to impute genotypes for analysis in an additional 2,893 individuals who were genotyped only with the Affymetrix Mapping 10K Array. In this second round, we identified large stretches of chromosome shared within each family and probabilistically “filled-in” genotypes within each stretch whenever one or more of its carriers was genotyped with the 500K Array Set9,21. These 2,893 individuals were mostly offspring and siblings of the 1,412 individuals genotyped at high density. For computational efficiency, the second round of imputation was applied to sub-pedigrees, each including no more than 20–25 individuals.

Follow-up Samples and Genotyping

Based on analysis of the combined SardiNIA and FUSION results, SNPs rs6060369 and rs143383 were each examined in up to 28,801 additional individuals. These included individuals of European ancestry from the FUSION study stage 2 samples (N = 2,466), the Diabetes Genetics Initiative22 (N = 2,985), the Old Order Amish23,24 (N = 2,711), the Atherosclerosis Risk in Communities (ARIC) Study25 (N = 11,370), the Caerphilly Study26,27 (1,389 men), and the British Women’s Heart and Health Study28 (3,685 women). 4,195 African American individuals from the ARIC study were also examined. Additional details of follow-up samples and genotyping methods are included in Supplementary Methods. Within each follow-up sample, SNP genotype success rates were >90% and genotype counts were consistent with Hardy-Weinberg equilibrium (p>.05).

Genotype Imputation

Our initial screen was based on the meta-analysis of two genome-wide association studies. Because the studies used two different genotyping platforms, we imputed genotypes for all polymorphic HapMap SNPs in each study, using a Hidden Markov Model programmed in MACH6,8. Details are provided in the Supplementary Methods.

GWAS Analysis

Within the FUSION and SardiNIA study samples, we carried out association analyses to relate observed and estimated genotypes to height. At each SNP, height was related to allele counts for a reference allele in a regression model that also included sex, age, and age2 as covariates; FUSION covariates also included birth province and study6. For SNPs genotyped in the laboratory, allele counts were discrete (0, 1, or 2), whereas for imputed SNPs, allele counts were fractional (between 0.0 and 2.0, depending on the expected number of copies of the allele for each individual). For FUSION, T2D and control individuals were analyzed separately and results combined using the meta-analytic techniques described below. To allow for relatedness, regression coefficients were estimated in the context of a variance components model that can handle imputed genotypes and accounts for background polygenic effects9. When evaluating significance, we applied quantile normalization to trait values (SardiNIA) or to residuals after adjustment for covariates (FUSION), by ranking all height values and then converting them to z-scores according to quantiles of the standard normal distribution. In parallel to the analysis of quantile normalized data, we also analyzed untransformed height (in centimeters) to estimate effect sizes.

Meta-Analysis

To summarize results for the three initial scans (1,084 T2D cases and 1,287 controls from FUSION, and 4,298 individuals from Sardinia), we carried out a meta-analysis. We used meta-analysis rather than an analysis of pooled data to avoid an increase in false positive rates due to population stratification. The Sardinian and Finnish populations are genetically and geographically distinct, with an average Fst of .025 among the 45,284 autosomal SNPs genotyped in both samples, and clear differences in height. The genomic control parameter for our meta-analysis, which estimates inflation in test statistics due to the combined effects of population stratification, cryptic relatedness, and genotyping error29, was 1.02, suggesting both that population stratification and unmodeled relatedness had a negligible impact on our association results and that our meta-analysis of imputed genotypes resulted in appropriate control of false-positive rates.

For each marker, we selected an arbitrary reference allele and calculated a z-statistic characterizing the evidence for association in each study (summarizing both the p-value, in its magnitude, and the direction of effect, in its sign). We then calculated an overall z-statistic as a weighted average of the three individual statistics and calculated the corresponding p-value20. Weights were proportional to the square-root of the number of individuals examined in each sample and were selected such that the squared weights summed to 1.0. An analogous strategy was used to summarize results of follow-up genotyping.

Supplementary Material

Acknowledgements

The authors thank the many individuals who generously participated in this study. We thank Yun Li for computational analysis, Kyle Gaulton, Amy Swift, Michael Erdos, and Narisu Narisu for FUSION genotyping and technical expertise, Center for Inherited Disease Research investigators for FUSION GWAS genotyping, the staff and participants of the ARIC study for their important contributions, and the Amish Research Clinic Staff for their energetic efforts in study subject recruitment and characterization. The SardiNIA team thanks Monsignore Piseddu, Bishop of Ogliastra, the mayors and citizens of the 4 Sardinian towns (Lanusei, Ilbono, Arzana, and Elini), and the head of the Public Health Unit ASL4 for their volunteerism and cooperation; the team also thanks Andrea Maschio, Fabio Busonero, Antonella Mulas, Natascia Sestu and Maria Grazia Piras for genotyping and technical expertise, and all the physicians, nurses and the recruitment personnel of the ProgeNIA group in Lanusei. This research was supported (in part) by the intramural Research Program of the NIH, National Institute on Aging [contracts NO1-AG-1-2109 to the SardiNIA (“Progenia”) team and 263-MA-410953 to the University of Michigan (G.R.A)]. Other support was from NIH grants DK072193 (K.L.M.), HL084729 (G.R.A.), HG002651 (G.R.A.), DK062370 (M.B.), DK54361 (A.R.S.), HL72515 (A.R.S.), AG18728 (A.R.S.), AR046838 (A.R.S.); National Human Genome Research Institute intramural project number 1 Z01 HG000024 (F.S.C.); the American Diabetes Association, including a post-doctoral fellowship award (C.J.W.); and March of Dimes research grant 6-FY04-61 (J.N.H). K.L.M. and G.R.A are Pew Scholars in the Biomedical Sciences. The Atherosclerosis Risk in Communities Study is carried out as a collaborative study supported by National Heart, Lung, and Blood Institute (NHLBI) contracts N01-HC-55015, N01-HC-55016, N01-HC-55018, N01-HC-55019, N01-HC-55020, N01-HC-55021, and N01-HC-55022. FUSION GWAS genotyping was performed with support from CIDR NIH Contract Number N01-HG-65403 and the JHU GRCF SNP Center. The whole-genome genotyping and analysis in the DGI genome scan was supported by Novartis Institutes for BioMedical Research (to D. Altshuler), with additional support from The Richard and Susan Smith Family Foundation/American Diabetes Association Pinnacle Program Project Award (to D. Altshuler, J.N.H. and M.J. Daly). Funding and support were also provided by the University of Maryland General Clinical Research Center (M01 RR 16500), the NIDDK Clinical Nutrition Research Unit of Maryland (NIH P30 DK072488), the Department of Veterans Affairs, and Veterans Affairs Medical Center Baltimore Geriatric Research, Education and Clinical Center (GRECC). The Caerphilly study was funded by the Medical Research Council (UK). The Caerphilly study was undertaken by the former MRC Epidemiology Unit (South Wales) and was funded by the Medical Research Council of the United Kingdom. The data archive is maintained by the Department of Social Medicine, University of Bristol. We acknowledge use of genotype data from the British 1958 Birth Cohort DNA collection, funded by the Medical Research Council grant G0000934 and the Wellcome Trust grant 068545/Z/02.

Footnotes

Accession numbers. GenBank mRNA sequence for growth differentiation factor 5 (GDF5), NM_000557 and for basic FGF-repressed Zic binding protein (BFZB, UQCC, C20orf44) is NM_018244.

References

1. Weedon MN, et al. A common variant of HMGA2 is associated with adult and childhood height in the general population. Nat Genet. 2007;39:1245–1250. [PMC free article] [PubMed]
2. Miyamoto Y, et al. A functional polymorphism in the 5' UTR of GDF5 is associated with susceptibility to osteoarthritis. Nat Genet. 2007;39:529–533. [PubMed]
3. Silventoinen K, et al. Heritability of adult body height: a comparative study of twin cohorts in eight countries. Twin Res. 2003;6:399–408. [PubMed]
4. Pilia G, et al. Heritability of cardiovascular and personality traits in 6,148 Sardinians. PLoS Genet. 2006;2:e132. [PMC free article] [PubMed]
5. Perola M, et al. Combined genome scans for body stature in 6,602 European twins: evidence for common Caucasian loci. PLoS Genet. 2007;3:e97. [PubMed]
6. Scott LJ, et al. A genome-wide association study of type 2 diabetes in Finns detects multiple susceptibility variants. Science. 2007;316:1341–1345. [PMC free article] [PubMed]
7. Scuteri A, et al. Genome-wide association scan shows genetic variants in the FTO gene are associated with obesity-related traits. PLoS Genet. 2007;3:e115. [PubMed]
8. Li Y, Willer CJ, Ding J, Scheet P, Abecasis GR. Markov model for rapid haplotyping and genotype imputation in genome wide studies. Nature Genetics. 2007 Submitted. http://www.sph.umich.edu/csg/abecasis/MACH/
9. Chen WM, Abecasis GR. Family based association tests for genome wide association scans. Am J Hum Genet. 2007;81:913–926. [PubMed]
10. Vetter K, Wurst W. Expression of a novel mouse gene 'mbFZb' in distinct regions of the developing nervous system and the adult brain. Mech Dev. 2001;100:123–125. [PubMed]
11. Valdes AM, et al. Sex and ethnic differences in the association of ASPN, CALM1, COL2A1, COMP, and FRZB with genetic susceptibility to osteoarthritis of the knee. Arthritis Rheum. 2007;56:137–146. [PubMed]
12. Imabayashi H, et al. Redifferentiation of dedifferentiated chondrocytes and chondrogenesis of human bone marrow stromal cells via chondrosphere formation with expression profiling by large-scale cDNA analysis. Exp Cell Res. 2003;288:35–50. [PubMed]
13. Goldring MB, Tsuchimochi K, Ijiri K. The control of chondrogenesis. J Cell Biochem. 2006;97:33–44. [PubMed]
14. Chang SC, et al. Cartilage-derived morphogenetic proteins. New members of the transforming growth factor-beta superfamily predominantly expressed in long bones during human embryonic development. J Biol Chem. 1994;269:28227–28234. [PubMed]
15. Chujo T, et al. Effects of growth differentiation factor-5 on the intervertebral disc--in vitro bovine study and in vivo rabbit disc degeneration model study. Spine. 2006;31:2909–2917. [PubMed]
16. Voight BF, Kudaravalli S, Wen X, Pritchard JK. A map of recent positive selection in the human genome. PLoS Biol. 2006;4:e72. [PMC free article] [PubMed]
17. Southam L, et al. An SNP in the 5'-UTR of GDF5 is associated with osteoarthritis susceptibility in Europeans and with in vivo differences in allelic expression in articular cartilage. Hum Mol Genet. 2007;16:2226–2232. [PubMed]
18. McKenzie CA, et al. Trans-ethnic fine mapping of a quantitative trait locus for circulating angiotensin I-converting enzyme (ACE) Hum Mol Genet. 2001;10:1077–1084. [PubMed]
19. Frere C, et al. Fine mapping of quantitative trait nucleotides underlying thrombin-activatable fibrinolysis inhibitor antigen levels by a transethnic study. Blood. 2006;108:1562–1568. [PubMed]
20. Skol AD, Scott LJ, Abecasis GR, Boehnke M. Joint analysis is more efficient than replication-based analysis for two-stage genome-wide association studies. Nat Genet. 2006;38:209–213. [PubMed]
21. Burdick JT, Chen WM, Abecasis GR, Cheung VG. In silico method for inferring genotypes in pedigrees. Nat Genet. 2006;38:1002–1004. [PMC free article] [PubMed]
22. Diabetes Genetics Initiative. Genome-wide association analysis identifies loci for type 2 diabetes and triglyceride levels. Science. 2007;316:1331–1336. [PubMed]
23. Hsueh WC, et al. Diabetes in the Old Order Amish: characterization and heritability analysis of the Amish Family Diabetes Study. Diabetes Care. 2000;23:595–601. [PubMed]
24. Streeten EA, et al. Reduced incidence of hip fracture in the Old Order Amish. J Bone Miner Res. 2004;19:308–313. [PubMed]
25. ARIC Investigators. The Atherosclerosis Risk in Communities (ARIC) Study: design and objectives. Am J Epidemiol. 1989;129:687–702. [PubMed]
26. The Caerphilly and Speedwell Collaborative Group. Caerphilly and Speedwell collaborative heart disease studies. J Epidemiol Community Health. 1984;38:259–262. [PMC free article] [PubMed]
27. Fehily AM, Butland BK, Yarnell JW. Body fatness and frame size: the Caerphilly study. Eur J Clin Nutr. 1990;44:107–111. [PubMed]
28. Lawlor DA, Bedford C, Taylor M, Ebrahim S. Geographical variation in cardiovascular disease, risk factors, and their control in older women: British Women's Heart and Health Study. J Epidemiol Community Health. 2003;57:134–140. [PMC free article] [PubMed]
29. Devlin B, Roeder K. Genomic control for association studies. Biometrics. 1999;55:997–1004. [PubMed]
30. The International HapMap Consortium. A second generation human haplotype map of over 3.1 million SNPs. Nature. 2007;449:851–861. [PMC free article] [PubMed]