Common variants at 11p13 are associated with susceptibility to tuberculosis

After imputation of data from the 1000 Genomes Project into a genome-wide dataset of Ghanaian individuals with tuberculosis and controls, we identified a resistance locus on chromosome 11p13 downstream of the WT1 gene (encoding Wilms tumor 1). The strongest signal was obtained at the rs2057178 SNP (P = 2.63 × 10(-9)). Replication in Gambian, Indonesian and Russian tuberculosis case-control study cohorts increased the significance level for the association with this SNP to P = 2.57 × 10(-11).     

Common variants near TARDBP and EGR2 are associated with susceptibility to Ewing sarcoma

Ewing sarcoma, a pediatric tumor characterized by EWSR1-ETS fusions, is predominantly observed in populations of European ancestry. We performed a genome-wide association study (GWAS) of 401 French individuals with Ewing sarcoma, 684 unaffected French individuals and 3,668 unaffected individuals of European descent and living in the United States. We identified candidate risk loci at 1p36.22, 10q21 and 15q15. We replicated these loci in two independent sets of cases and controls. Joint analysis identified associations with rs9430161 (P = 1.4 × 10(-20); odds ratio (OR) = 2.2) located 25 kb upstream of TARDBP, rs224278 (P = 4.0 × 10(-17); OR = 1.7) located 5 kb upstream of EGR2 and, to a lesser extent, rs4924410 at 15q15 (P = 6.6 × 10(-9); OR = 1.5). The major risk haplotypes were less prevalent in Africans, suggesting that these loci could contribute to geographical differences in Ewing sarcoma incidence. TARDBP shares structural similarities with EWSR1 and FUS, which encode RNA binding proteins, and EGR2 is a target gene of EWSR1-ETS. Variants at these loci were associated with expression levels of TARDBP, ADO (encoding cysteamine dioxygenase) and EGR2.     

Common variants in DVWA on chromosome 3p24.3 are associated with susceptibility to knee osteoarthritis

Susceptibility to osteoarthritis, the most common human arthritis, is known to be influenced by genetic factors. Through a genome-wide association study using approximately 100,000 SNPs, we have identified a previously unknown gene on chromosome 3p24.3, DVWA, which is associated with susceptibility to knee osteoarthritis. Expressed specifically in cartilage, DVWA encodes a 276-amino-acid protein with two regions corresponding to the von Willebrand factor type A domain (VWA domain). Several DVWA SNPs are significantly associated with knee osteoarthritis in two independent Japanese case-control cohorts. This association was replicated in a Japanese population cohort and a Han Chinese case-control cohort (combined P = 7.3 x 10(-11)). DVWA protein binds to beta-tubulin, and the binding is influenced by two highly associated missense SNPs (rs11718863 and rs7639618) located in the VWA domain. The Tyr169-Cys260 isoform of DVWA, which is overrepresented in knee osteoarthritis, showed weaker interaction. Our findings reveal a new paradigm for study of osteoarthritis etiology and pathogenesis.     

Common variants at 11q12, 10q26 and 3p11.2 are associated with prostate cancer susceptibility in Japanese

We have previously reported multiple loci associated with prostate cancer susceptibility in a Japanese population using a genome-wide association study (GWAS). To identify additional prostate cancer susceptibility loci, we genotyped nine SNPs that were nominally associated with prostate cancer (P < 1 × 10(-4)) in our previous GWAS in three independent studies of prostate cancer in Japanese men (2,557 individuals with prostate cancer (cases) and 3,003 controls). In a meta-analysis of our previous GWAS and the replication studies, which included a total of 7,141 prostate cancer cases and 11,804 controls from a single ancestry group, three new loci reached genome-wide significance on chromosomes 11q12 (rs1938781; P = 1.10 × 10(-10); FAM111A-FAM111B), 10q26 (rs2252004; P = 1.98 × 10(-8)) and 3p11.2 (rs2055109; P = 3.94 × 10(-8)). We also found suggestive evidence of association at a previously reported prostate cancer susceptibility locus at 2p11 (rs2028898; P = 1.08 × 10(-7)). The identification of three new susceptibility loci should provide additional insight into the pathogenesis of prostate cancer and emphasizes the importance of conducting GWAS in diverse populations.     

Common variants in P2RY11 are associated with narcolepsy

Growing evidence supports the hypothesis that narcolepsy with cataplexy is an autoimmune disease. We here report genome-wide association analyses for narcolepsy with replication and fine mapping across three ethnic groups (3,406 individuals of European ancestry, 2,414 Asians and 302 African Americans). We identify a SNP in the 3' untranslated region of P2RY11, the purinergic receptor subtype P2Y?? gene, which is associated with narcolepsy (rs2305795, combined P = 6.1 × 10?1?, odds ratio = 1.28, 95% CI 1.19-1.39, n = 5689). The disease-associated allele is correlated with reduced expression of P2RY11 in CD8(+) T lymphocytes (339% reduced, P = 0.003) and natural killer (NK) cells (P = 0.031), but not in other peripheral blood mononuclear cell types. The low expression variant is also associated with reduced P2RY11-mediated resistance to ATP-induced cell death in T lymphocytes (P = 0.0007) and natural killer cells (P = 0.001). These results identify P2RY11 as an important regulator of immune-cell survival, with possible implications in narcolepsy and other autoimmune diseases.     

Common variants near ATM are associated with glycemic response to metformin in type 2 diabetes

Metformin is the most commonly used pharmacological therapy for type 2 diabetes. We report a genome-wide association study for glycemic response to metformin in 1,024 Scottish individuals with type 2 diabetes with replication in two cohorts including 1,783 Scottish individuals and 1,113 individuals from the UK Prospective Diabetes Study. In a combined meta-analysis, we identified a SNP, rs11212617, associated with treatment success (n = 3,920, P = 2.9 × 10(-9), odds ratio = 1.35, 95% CI 1.22-1.49) at a locus containing ATM, the ataxia telangiectasia mutated gene. In a rat hepatoma cell line, inhibition of ATM with KU-55933 attenuated the phosphorylation and activation of AMP-activated protein kinase in response to metformin. We conclude that ATM, a gene known to be involved in DNA repair and cell cycle control, plays a role in the effect of metformin upstream of AMP-activated protein kinase, and variation in this gene alters glycemic response to metformin.     

Common variants of FUT2 are associated with plasma vitamin B12 levels

We identified a strong association (P = 5.36 x 10(-17)) between rs492602 in FUT2 and plasma vitamin B(12) levels in a genome-wide scan (n = 1,658) and an independent replication sample (n = 1,059) from the Nurses' Health Study. Women homozygous for the rs492602[G] allele had higher B(12) levels. This allele is in strong linkage disequilibrium with the FUT2 nonsecretor variant encoding W143X, suggesting a plausible mechanism for altered B(12) absorption and plasma levels.     

Common variants at 12q14 and 12q24 are associated with hippocampal volume

Aging is associated with reductions in hippocampal volume that are accelerated by Alzheimer's disease and vascular risk factors. Our genome-wide association study (GWAS) of dementia-free persons (n = 9,232) identified 46 SNPs at four loci with P values of <4.0 × 10(-7). In two additional samples (n = 2,318), associations were replicated at 12q14 within MSRB3-WIF1 (discovery and replication; rs17178006; P = 5.3 × 10(-11)) and at 12q24 near HRK-FBXW8 (rs7294919; P = 2.9 × 10(-11)). Remaining associations included one SNP at 2q24 within DPP4 (rs6741949; P = 2.9 × 10(-7)) and nine SNPs at 9p33 within ASTN2 (rs7852872; P = 1.0 × 10(-7)); along with the chromosome 12 associations, these loci were also associated with hippocampal volume (P < 0.05) in a third younger, more heterogeneous sample (n = 7,794). The SNP in ASTN2 also showed suggestive association with decline in cognition in a largely independent sample (n = 1,563). These associations implicate genes related to apoptosis (HRK), development (WIF1), oxidative stress (MSR3B), ubiquitination (FBXW8) and neuronal migration (ASTN2), as well as enzymes targeted by new diabetes medications (DPP4), indicating new genetic influences on hippocampal size and possibly the risk of cognitive decline and dementia.     

Common variants at 6q22 and 17q21 are associated with intracranial volume

During aging, intracranial volume remains unchanged and represents maximally attained brain size, while various interacting biological phenomena lead to brain volume loss. Consequently, intracranial volume and brain volume in late life reflect different genetic influences. Our genome-wide association study (GWAS) in 8,175 community-dwelling elderly persons did not reveal any associations at genome-wide significance (P < 5 × 10(-8)) for brain volume. In contrast, intracranial volume was significantly associated with two loci: rs4273712 (P = 3.4 × 10(-11)), a known height-associated locus on chromosome 6q22, and rs9915547 (P = 1.5 × 10(-12)), localized to the inversion on chromosome 17q21. We replicated the associations of these loci with intracranial volume in a separate sample of 1,752 elderly persons (P = 1.1 × 10(-3) for 6q22 and 1.2 × 10(-3) for 17q21). Furthermore, we also found suggestive associations of the 17q21 locus with head circumference in 10,768 children (mean age of 14.5 months). Our data identify two loci associated with head size, with the inversion at 17q21 also likely to be involved in attaining maximal brain size.     

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 variants in DGKK are strongly associated with risk of hypospadias

Hypospadias is a common congenital malformation of the male external genitalia. We performed a genome-wide association study using pooled DNA from 436 individuals with hypospadias (cases) and 494 controls of European descent and selected the highest ranked SNPs for individual genotyping in the discovery sample, an additional Dutch sample of 133 cases and their parents, and a Swedish series of 266 cases and 402 controls. Individual genotyping of two SNPs (rs1934179 and rs7063116) in DGKK, encoding diacylglycerol kinase κ, produced compelling evidence for association with hypospadias in the discovery sample (allele-specific odds ratio (OR) = 2.5, P = 2.5 × 10?11 and OR = 2.3, P = 2.9 × 10??, respectively) and in the Dutch (OR = 3.9, P = 2.4 × 10?? and OR = 3.8, P = 3.4 × 10??) and Swedish (OR = 2.5, P = 2.6 × 10?? and OR = 2.2, P = 2.7 × 10??) replication samples. Expression studies showed expression of DGKK in preputial tissue of cases and controls, which was lower in carriers of the risk allele of rs1934179 (P = 0.047). We propose DGKK as a major risk gene for hypospadias.     

Common variants near MBNL1 and NKX2-5 are associated with infantile hypertrophic pyloric stenosis

Infantile hypertrophic pyloric stenosis (IHPS) is a severe condition characterized by hypertrophy of the pyloric sphincter muscle. We conducted a genome-wide association study (GWAS) on 1,001 surgery-confirmed cases and 2,401 controls from Denmark. The six most strongly associated loci were tested in a replication set of 796 cases and 876 controls. Three SNPs reached genome-wide significance. One of these SNPs, rs11712066 (odds ratio (OR) = 1.61; P = 1.5 × 10(-17)) at 3p25.1, is located 150 kb upstream of MBNL1, which encodes a factor that regulates splicing transitions occurring shortly after birth. The second SNP, rs573872 (OR = 1.41; P = 4.3 × 10(-12)), maps to an intergenic region at 3p25.2 approximately 1.3 Mb downstream of MBNL1. The third SNP, rs29784 (OR = 1.42; P = 1.5 × 10(-15)) at 5q35.2, is 64 kb downstream of NKX2-5, which is involved in development of cardiac muscle tissue and embryonic gut development.     

Common variants at CDKAL1 and KLF9 are associated with body mass index in east Asian populations

Obesity is a disorder with a complex genetic etiology, and its epidemic is a worldwide problem. Although multiple genetic loci associated with body mass index, the most common measure of obesity, have been identified in European populations, few studies have focused on Asian populations. Here we report a genome-wide association study and replication studies with 62,245 east Asian subjects, which identified two new body mass index-associated loci in the CDKAL1 locus at 6p22 (rs2206734, P = 1.4 × 10(-11)) and the KLF9 locus at 9q21 (rs11142387, P = 1.3 × 10(-9)), as well as several previously reported loci (the SEC16B, BDNF, FTO, MC4R and GIPR loci, P < 5.0 × 10(-8)). We subsequently performed gene-gene interaction analyses and identified an interaction (P = 2.0 × 10(-8)) between a SNP in the KLF9 locus (rs11142387) and one in the MSTN (also known as GDF8) locus at 2q32 (rs13034723). These findings should provide useful insights into the etiology of obesity.     

Common and rare variants in multifactorial susceptibility to common diseases

Here, we give a historical overview of the search for genetic variants that influence the susceptibility of an individual to a chronic disease, from RA Fisher's seminal work to the current excitement of whole-genome association studies (WGAS). We then discuss the concepts behind the identification of common variants as disease causal factors and contrast them to the basic ideas that underlie the rare variant hypothesis. The identification of rare variants involves the careful selection of candidate genes to examine, the availability of highly efficient resequencing techniques and the appropriate assessment of the functional consequences of the implicated variant. We believe that this strategy can be successfully applied at present in order to unravel the contribution of rare variants to the multifactorial inheritance of common diseases, which could lead to the implementation of much needed preventative screening schemes.     

Common variants in ZNF365 are associated with both mammographic density and breast cancer risk

High-percent mammographic density adjusted for age and body mass index is one of the strongest risk factors for breast cancer. We conducted a meta analysis of five genome-wide association studies of percent mammographic density and report an association with rs10995190 in ZNF365 (combined P = 9.6 × 10(-10)). Common variants in ZNF365 have also recently been associated with susceptibility to breast cancer.     

Common sequence variants on 2p15 and Xp11.22 confer susceptibility to prostate cancer

We conducted a genome-wide SNP association study on prostate cancer on over 23,000 Icelanders, followed by a replication study including over 15,500 individuals from Europe and the United States. Two newly identified variants were shown to be associated with prostate cancer: rs5945572 on Xp11.22 and rs721048 on 2p15 (odds ratios (OR) = 1.23 and 1.15; P = 3.9 x 10(-13) and 7.7 x 10(-9), respectively). The 2p15 variant shows a significantly stronger association with more aggressive, rather than less aggressive, forms of the disease.     

Common variants at CD40 and other loci confer risk of rheumatoid arthritis

To identify rheumatoid arthritis risk loci in European populations, we conducted a meta-analysis of two published genome-wide association (GWA) studies totaling 3,393 cases and 12,462 controls. We genotyped 31 top-ranked SNPs not previously associated with rheumatoid arthritis in an independent replication of 3,929 autoantibody-positive rheumatoid arthritis cases and 5,807 matched controls from eight separate collections. We identified a common variant at the CD40 gene locus (rs4810485, P = 0.0032 replication, P = 8.2 x 10(-9) overall, OR = 0.87). Along with other associations near TRAF1 (refs. 2,3) and TNFAIP3 (refs. 4,5), this implies a central role for the CD40 signaling pathway in rheumatoid arthritis pathogenesis. We also identified association at the CCL21 gene locus (rs2812378, P = 0.00097 replication, P = 2.8 x 10(-7) overall), a gene involved in lymphocyte trafficking. Finally, we identified evidence of association at four additional gene loci: MMEL1-TNFRSF14 (rs3890745, P = 0.0035 replication, P = 1.1 x 10(-7) overall), CDK6 (rs42041, P = 0.010 replication, P = 4.0 x 10(-6) overall), PRKCQ (rs4750316, P = 0.0078 replication, P = 4.4 x 10(-6) overall), and KIF5A-PIP4K2C (rs1678542, P = 0.0026 replication, P = 8.8 x 10(-8) overall).     

Variants in KCNQ1 are associated with susceptibility to type 2 diabetes mellitus

We carried out a multistage genome-wide association study of type 2 diabetes mellitus in Japanese individuals, with a total of 1,612 cases and 1,424 controls and 100,000 SNPs. The most significant association was obtained with SNPs in KCNQ1, and dense mapping within the gene revealed that rs2237892 in intron 15 showed the lowest Pvalue (6.7 x 10(-13), odds ratio (OR) = 1.49). The association of KCNQ1 with type 2 diabetes was replicated in populations of Korean, Chinese and European ancestry as well as in two independent Japanese populations, and meta-analysis with a total of 19,930 individuals (9,569 cases and 10,361 controls) yielded a P value of 1.7 x 10(-42) (OR = 1.40; 95% CI = 1.34-1.47) for rs2237892. Among control subjects, the risk allele of this polymorphism was associated with impairment of insulin secretion according to the homeostasis model assessment of beta-cell function or the corrected insulin response. Our data thus implicate KCNQ1 as a diabetes susceptibility gene in groups of different ancestries.