Genome-wide association analysis identifies 20 loci that influence adult height

Adult height is a model polygenic trait, but there has been limited success in identifying the genes underlying its normal variation. To identify genetic variants influencing adult human height, we used genome-wide association data from 13,665 individuals and genotyped 39 variants in an additional 16,482 samples. We identified 20 variants associated with adult height (P < 5 x 10(-7), with 10 reaching P < 1 x 10(-10)). Combined, the 20 SNPs explain approximately 3% of height variation, with a approximately 5 cm difference between the 6.2% of people with 17 or fewer 'tall' alleles compared to the 5.5% with 27 or more 'tall' alleles. The loci we identified implicate genes in Hedgehog signaling (IHH, HHIP, PTCH1), extracellular matrix (EFEMP1, ADAMTSL3, ACAN) and cancer (CDK6, HMGA2, DLEU7) pathways, and provide new insights into human growth and developmental processes. Finally, our results provide insights into the genetic architecture of a classic quantitative trait.     

A common variant of HMGA2 is associated with adult and childhood height in the general population

Human height is a classic, highly heritable quantitative trait. To begin to identify genetic variants influencing height, we examined genome-wide association data from 4,921 individuals. Common variants in the HMGA2 oncogene, exemplified by rs1042725, were associated with height (P = 4 x 10(-8)). HMGA2 is also a strong biological candidate for height, as rare, severe mutations in this gene alter body size in mice and humans, so we tested rs1042725 in additional samples. We confirmed the association in 19,064 adults from four further studies (P = 3 x 10(-11), overall P = 4 x 10(-16), including the genome-wide association data). We also observed the association in children (P = 1 x 10(-6), N = 6,827) and a tall/short case-control study (P = 4 x 10(-6), N = 3,207). We estimate that rs1042725 explains approximately 0.3% of population variation in height (approximately 0.4 cm increased adult height per C allele). There are few examples of common genetic variants reproducibly associated with human quantitativetraits; these results represent, to our knowledge, the first consistently replicated association with adult and childhood height.     

Evidence of widespread selection on standing variation in Europe at height-associated SNPs

Strong signatures of positive selection at newly arising genetic variants are well documented in humans, but this form of selection may not be widespread in recent human evolution. Because many human traits are highly polygenic and partly determined by common, ancient genetic variation, an alternative model for rapid genetic adaptation has been proposed: weak selection acting on many pre-existing (standing) genetic variants, or polygenic adaptation. By studying height, a classic polygenic trait, we demonstrate the first human signature of widespread selection on standing variation. We show that frequencies of alleles associated with increased height, both at known loci and genome wide, are systematically elevated in Northern Europeans compared with Southern Europeans (P < 4.3 × 10(-4)). This pattern mirrors intra-European height differences and is not confounded by ancestry or other ascertainment biases. The systematic frequency differences are consistent with the presence of widespread weak selection (selection coefficients ～10(-3)-10(-5) per allele) rather than genetic drift alone (P < 10(-15)).     

Variation in factor B (BF) and complement component 2 (C2) genes is associated with age-related macular degeneration

Age-related macular degeneration (AMD) is the most common form of irreversible blindness in developed countries. Variants in the factor H gene (CFH, also known as HF1), which encodes a major inhibitor of the alternative complement pathway, are associated with the risk for developing AMD. Here we test the hypothesis that variation in genes encoding other regulatory proteins of the same pathway is associated with AMD. We screened factor B (BF) and complement component 2 (C2) genes, located in the major histocompatibility complex class III region, for genetic variation in two independent cohorts comprising approximately 900 individuals with AMD and approximately 400 matched controls. Haplotype analyses identify a statistically significant common risk haplotype (H1) and two protective haplotypes. The L9H variant of BF and the E318D variant of C2 (H10), as well as a variant in intron 10 of C2 and the R32Q variant of BF (H7), confer a significantly reduced risk of AMD (odds ratio = 0.45 and 0.36, respectively). Combined analysis of the C2 and BF haplotypes and CFH variants shows that variation in the two loci can predict the clinical outcome in 74% of the affected individuals and 56% of the controls. These data expand and refine our understanding of the genetic risk for AMD.     

Common variants at 12q15 and 12q24 are associated with infant head circumference

To identify genetic variants associated with head circumference in infancy, we performed a meta-analysis of seven genome-wide association studies (GWAS) (N = 10,768 individuals of European ancestry enrolled in pregnancy and/or birth cohorts) and followed up three lead signals in six replication studies (combined N = 19,089). rs7980687 on chromosome 12q24 (P = 8.1 × 10(-9)) and rs1042725 on chromosome 12q15 (P = 2.8 × 10(-10)) were robustly associated with head circumference in infancy. Although these loci have previously been associated with adult height, their effects on infant head circumference were largely independent of height (P = 3.8 × 10(-7) for rs7980687 and P = 1.3 × 10(-7) for rs1042725 after adjustment for infant height). A third signal, rs11655470 on chromosome 17q21, showed suggestive evidence of association with head circumference (P = 3.9 × 10(-6)). SNPs correlated to the 17q21 signal have shown genome-wide association with adult intracranial volume, Parkinson's disease and other neurodegenerative diseases, indicating that a common genetic variant in this region might link early brain growth with neurological disease in later life.     

Common deletions and SNPs are in linkage disequilibrium in the human genome

Humans show great variation in phenotypic traits such as height, eye color and susceptibility to disease. Genomic DNA sequence differences among individuals are responsible for the inherited components of these complex traits. Reports suggest that intermediate and large-scale DNA copy number and structural variations are prevalent enough to be an important source of genetic variation between individuals. Because association studies to identify genomic loci associated with particular phenotypic traits have focused primarily on genotyping SNPs, it is important to determine whether common structural polymorphisms are in linkage disequilibrium with common SNPs, and thus can be assessed indirectly in SNP-based studies. Here we examine 100 deletion polymorphisms ranging from 70 bp to 7 kb. We show that common deletions and SNPs ascertained with similar criteria have essentially the same distribution of linkage disequilibrium with surrounding SNPs, indicating that these polymorphisms may share evolutionary history and that most deletion polymorphisms are effectively assayed by proxy in SNP-based association studies.     

Bayesian inference analyses of the polygenic architecture of rheumatoid arthritis

The genetic architectures of common, complex diseases are largely uncharacterized. We modeled the genetic architecture underlying genome-wide association study (GWAS) data for rheumatoid arthritis and developed a new method using polygenic risk-score analyses to infer the total liability-scale variance explained by associated GWAS SNPs. Using this method, we estimated that, together, thousands of SNPs from rheumatoid arthritis GWAS explain an additional 20% of disease risk (excluding known associated loci). We further tested this method on datasets for three additional diseases and obtained comparable estimates for celiac disease (43% excluding the major histocompatibility complex), myocardial infarction and coronary artery disease (48%) and type 2 diabetes (49%). Our results are consistent with simulated genetic models in which hundreds of associated loci harbor common causal variants and a smaller number of loci harbor multiple rare causal variants. These analyses suggest that GWAS will continue to be highly productive for the discovery of additional susceptibility loci for common diseases.     

Population-based resequencing of ANGPTL4 uncovers variations that reduce triglycerides and increase HDL

Resequencing genes provides the opportunity to assess the full spectrum of variants that influence complex traits. Here we report the first application of resequencing to a large population (n = 3,551) to examine the role of the adipokine ANGPTL4 in lipid metabolism. Nonsynonymous variants in ANGPTL4 were more prevalent in individuals with triglyceride levels in the lowest quartile than in individuals with levels in the highest quartile (P = 0.016). One variant (E40K), present in approximately 3% of European Americans, was associated with significantly lower plasma levels of triglyceride and higher levels of high-density lipoprotein cholesterol in European Americans from the Atherosclerosis Risk in Communities Study and in Danes from the Copenhagen City Heart Study. The ratio of nonsynonymous to synonymous variants was higher in European Americans than in African Americans (4:1 versus 1.3:1), suggesting population-specific relaxation of purifying selection. Thus, resequencing of ANGPTL4 in a multiethnic population allowed analysis of the phenotypic effects of both rare and common variants while taking advantage of genetic variation arising from ethnic differences in population history.     

Estimating the proportion of variation in susceptibility to schizophrenia captured by common SNPs

Schizophrenia is a complex disorder caused by both genetic and environmental factors. Using 9,087 affected individuals, 12,171 controls and 915,354 imputed SNPs from the Schizophrenia Psychiatric Genome-Wide Association Study (GWAS) Consortium (PGC-SCZ), we estimate that 23% (s.e. = 1%) of variation in liability to schizophrenia is captured by SNPs. We show that a substantial proportion of this variation must be the result of common causal variants, that the variance explained by each chromosome is linearly related to its length (r = 0.89, P = 2.6 × 10(-8)), that the genetic basis of schizophrenia is the same in males and females, and that a disproportionate proportion of variation is attributable to a set of 2,725 genes expressed in the central nervous system (CNS; P = 7.6 × 10(-8)). These results are consistent with a polygenic genetic architecture and imply more individual SNP associations will be detected for this disease as sample size increases.     

A genome-wide association study in Chinese men identifies three risk loci for non-obstructive azoospermia

Non-obstructive azoospermia (NOA) is one of the most severe forms of male infertility. Its pathophysiology is largely unknown, and few genetic influences have been defined. To identify common variants contributing to NOA in Han Chinese men, we performed a three-stage genome-wide association study of 2,927 individuals with NOA and 5,734 controls. The combined analyses identified significant (P < 5.0 × 10(-8)) associations between NOA risk and common variants near PRMT6 (rs12097821 at 1p13.3: odds ratio (OR) = 1.25, P = 5.7 × 10(-10)), PEX10 (rs2477686 at 1p36.32: OR = 1.39, P = 5.7 × 10(-12)) and SOX5 (rs10842262 at 12p12.1: OR = 1.23, P = 2.3 × 10(-9)). These findings implicate genetic variants at 1p13.3, 1p36.32 and 12p12.1 in the etiology of NOA in Han Chinese men.     

A common polymorphism of the growth hormone receptor is associated with increased responsiveness to growth hormone

Growth hormone is used to increase height in short children who are not deficient in growth hormone, but its efficacy varies largely across individuals. The genetic factors responsible for this variation are entirely unknown. In two cohorts of short children treated with growth hormone, we found that an isoform of the growth hormone receptor gene that lacks exon 3 (d3-GHR) was associated with 1.7 to 2 times more growth acceleration induced by growth hormone than the full-length isoform (P < 0.0001). In transfection experiments, the transduction of growth hormone signaling through d3-GHR homo- or heterodimers was approximately 30% higher than through full-length GHR homodimers (P < 0.0001). One-half of Europeans are hetero- or homozygous with respect to the allele encoding the d3-GHR isoform, which is dominant over the full-length isoform. These observations suggest that the polymorphism in exon 3 of GHR is important in growth hormone pharmacogenetics.     

Comparing strategies to fine-map the association of common SNPs at chromosome 9p21 with type 2 diabetes and myocardial infarction

Noncoding variants at human chromosome 9p21 near CDKN2A and CDKN2B are associated with type 2 diabetes, myocardial infarction, aneurysm, vertical cup disc ratio and at least five cancers. Here we compare approaches to more comprehensively assess genetic variation in the region. We carried out targeted sequencing at high coverage in 47 individuals and compared the results to pilot data from the 1000 Genomes Project. We imputed variants into type 2 diabetes and myocardial infarction cohorts directly from targeted sequencing, from a genotyped reference panel derived from sequencing and from 1000 Genomes Project low-coverage data. Polymorphisms with frequency >5% were captured well by all strategies. Imputation of intermediate-frequency polymorphisms required a higher density of tag SNPs in disease samples than is available on first-generation genome-wide association study (GWAS) arrays. Our association analyses identified more comprehensive sets of variants showing equivalent statistical association with type 2 diabetes or myocardial infarction, but did not identify stronger associations than the original GWAS signals.     

Many sequence variants affecting diversity of adult human height

Adult human height is one of the classical complex human traits. We searched for sequence variants that affect height by scanning the genomes of 25,174 Icelanders, 2,876 Dutch, 1,770 European Americans and 1,148 African Americans. We then combined these results with previously published results from the Diabetes Genetics Initiative on 3,024 Scandinavians and tested a selected subset of SNPs in 5,517 Danes. We identified 27 regions of the genome with one or more sequence variants showing significant association with height. The estimated effects per allele of these variants ranged between 0.3 and 0.6 cm and, taken together, they explain around 3.7% of the population variation in height. The genes neighboring the identified loci cluster in biological processes related to skeletal development and mitosis. Association to three previously reported loci are replicated in our analyses, and the strongest association was with SNPs in the ZBTB38 gene.     

Identification of ten loci associated with height highlights new biological pathways in human growth

Height is a classic polygenic trait, reflecting the combined influence of multiple as-yet-undiscovered genetic factors. We carried out a meta-analysis of genome-wide association study data of height from 15,821 individuals at 2.2 million SNPs, and followed up the strongest findings in >10,000 subjects. Ten newly identified and two previously reported loci were strongly associated with variation in height (P values from 4 x 10(-7) to 8 x 10(-22)). Together, these 12 loci account for approximately 2% of the population variation in height. Individuals with < or =8 height-increasing alleles and > or =16 height-increasing alleles differ in height by approximately 3.5 cm. The newly identified loci, along with several additional loci with strongly suggestive associations, encompass both strong biological candidates and unexpected genes, and highlight several pathways (let-7 targets, chromatin remodeling proteins and Hedgehog signaling) as important regulators of human stature. These results expand the picture of the biological regulation of human height and of the genetic architecture of this classical complex trait.     

Identification of common variants influencing risk of the tauopathy progressive supranuclear palsy

Progressive supranuclear palsy (PSP) is a movement disorder with prominent tau neuropathology. Brain diseases with abnormal tau deposits are called tauopathies, the most common of which is Alzheimer's disease. Environmental causes of tauopathies include repetitive head trauma associated with some sports. To identify common genetic variation contributing to risk for tauopathies, we carried out a genome-wide association study of 1,114 individuals with PSP (cases) and 3,247 controls (stage 1) followed by a second stage in which we genotyped 1,051 cases and 3,560 controls for the stage 1 SNPs that yielded P ≤ 10(-3). We found significant previously unidentified signals (P < 5 × 10(-8)) associated with PSP risk at STX6, EIF2AK3 and MOBP. We confirmed two independent variants in MAPT affecting risk for PSP, one of which influences MAPT brain expression. The genes implicated encode proteins for vesicle-membrane fusion at the Golgi-endosomal interface, for the endoplasmic reticulum unfolded protein response and for a myelin structural component.     

Variation in the DEPDC5 locus is associated with progression to hepatocellular carcinoma in chronic hepatitis C virus carriers

Chronic viral hepatitis is the most important risk factor for progression to hepatocellular carcinoma (HCC). To identify genetic risk factors for progression to HCC in individuals with chronic hepatitis C virus (HCV), we analyzed 467,538 SNPs in 212 Japanese individuals with chronic HCV with HCC and 765 individuals with chronic HCV without HCC. We identified one intronic SNP in the DEPDC5 locus on chromosome 22 associated with HCC risk and confirmed the association using an independent case-control population (710 cases and 1,625 controls). The association was highly significant when we analyzed the stages separately as well as together (rs1012068, P(combined) = 1.27 × 10(-13), odds ratio = 1.75). The significance level of the association further increased after adjustment for gender, age and platelet count (P = 1.35 × 10(-14), odds ratio = 1.96). Our findings suggest that common variants within the DEPDC5 locus affect susceptibility to HCC in Japanese individuals with chronic HCV infection.     

Population genetic structure of variable drug response.

Geographic patterns of genetic variation, including variation at drug metabolizing enzyme (DME) loci and drug targets, indicate that geographic structuring of inter-individual variation in drug response may occur frequently. This raises two questions: how to represent human population genetic structure in the evaluation of drug safety and efficacy, and how to relate this structure to drug response. We address these by (i) inferring the genetic structure present in a heterogeneous sample and (ii) comparing the distribution of DME variants across the inferred genetic clusters of individuals. We find that commonly used ethnic labels are both insufficient and inaccurate representations of the inferred genetic clusters, and that drug-metabolizing profiles, defined by the distribution of DME variants, differ significantly among the clusters. We note, however, that the complexity of human demographic history means that there is no obvious natural clustering scheme, nor an obvious appropriate degree of resolution. Our comparison of drug-metabolizing profiles across the inferred clusters establishes a framework for assessing the appropriate level of resolution in relating genetic structure to drug response.     

Genome-wide association study identifies a susceptibility locus for thoracic aortic aneurysms and aortic dissections spanning FBN1 at 15q21.1

Although thoracic aortic aneurysms and dissections (TAAD) can be inherited as a single-gene disorder, the genetic predisposition in the majority of affected people is poorly understood. In a multistage genome-wide association study (GWAS), we compared 765 individuals who had sporadic TAAD (STAAD) with 874 controls and identified common SNPs at a 15q21.1 locus that were associated with STAAD, with odds ratios of 1.6-1.8 that achieved genome-wide significance. We followed up 107 SNPs associated with STAAD with P < 1 × 10(-5) in the region, in two separate STAAD cohorts. The associated SNPs fall into a large region of linkage disequilibrium encompassing FBN1, which encodes fibrillin-1. FBN1 mutations cause Marfan syndrome, whose major cardiovascular complication is TAAD. This study shows that common genetic variants at 15q21.1 that probably act via FBN1 are associated with STAAD, suggesting a common pathogenesis of aortic disease in Marfan syndrome and STAAD.     

CFH haplotypes without the Y402H coding variant show strong association with susceptibility to age-related macular degeneration

In developed countries, age-related macular degeneration is a common cause of blindness in the elderly. A common polymorphism, encoding the sequence variation Y402H in complement factor H (CFH), has been strongly associated with disease susceptibility. Here, we examined 84 polymorphisms in and around CFH in 726 affected individuals (including 544 unrelated individuals) and 268 unrelated controls. In this sample, 20 of these polymorphisms showed stronger association with disease susceptibility than the Y402H variant. Further, no single polymorphism could account for the contribution of the CFH locus to disease susceptibility. Instead, multiple polymorphisms defined a set of four common haplotypes (of which two were associated with disease susceptibility and two seemed to be protective) and multiple rare haplotypes (associated with increased susceptibility in aggregate). Our results suggest that there are multiple disease susceptibility alleles in the region and that noncoding CFH variants play a role in disease susceptibility.