Genome-wide association study identifies susceptibility loci for open angle glaucoma at TMCO1 and CDKN2B-AS1

We report a genome-wide association study for open-angle glaucoma (OAG) blindness using a discovery cohort of 590 individuals with severe visual field loss (cases) and 3,956 controls. We identified associated loci at TMCO1 (rs4656461[G] odds ratio (OR) = 1.68, P = 6.1 × 10(-10)) and CDKN2B-AS1 (rs4977756[A] OR = 1.50, P = 4.7 × 10(-9)). We replicated these associations in an independent cohort of cases with advanced OAG (rs4656461 P = 0.010; rs4977756 P = 0.042) and two additional cohorts of less severe OAG (rs4656461 combined discovery and replication P = 6.00 × 10(-14), OR = 1.51, 95% CI 1.35-1.68; rs4977756 combined P = 1.35 × 10(-14), OR = 1.39, 95% CI 1.28-1.51). We show retinal expression of genes at both loci in human ocular tissues. We also show that CDKN2A and CDKN2B are upregulated in the retina of a rat model of glaucoma.     

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.     

Genome-wide association analysis in primary sclerosing cholangitis identifies two non-HLA susceptibility loci

Primary sclerosing cholangitis (PSC) is a chronic bile duct disease affecting 2.4-7.5% of individuals with inflammatory bowel disease. We performed a genome-wide association analysis of 2,466,182 SNPs in 715 individuals with PSC and 2,962 controls, followed by replication in 1,025 PSC cases and 2,174 controls. We detected non-HLA associations at rs3197999 in MST1 and rs6720394 near BCL2L11 (combined P = 1.1 × 10?1? and P = 4.1 × 10??, respectively).     

Genome-wide association study identifies susceptibility loci for dengue shock syndrome at MICB and PLCE1

Hypovolemic shock (dengue shock syndrome (DSS)) is the most common life-threatening complication of dengue. We conducted a genome-wide association study of 2,008 pediatric cases treated for DSS and 2,018 controls from Vietnam. Replication of the most significantly associated markers was carried out in an independent Vietnamese sample of 1,737 cases and 2,934 controls. SNPs at two loci showed genome-wide significant association with DSS. We identified a susceptibility locus at MICB (major histocompatibility complex (MHC) class I polypeptide-related sequence B), which was within the broad MHC region on chromosome 6 but outside the class I and class II HLA loci (rs3132468, P(meta) = 4.41 × 10(-11), per-allele odds ratio (OR) = 1.34 (95% confidence interval: 1.23-1.46)). We identified associated variants within PLCE1 (phospholipase C, epsilon 1) on chromosome 10 (rs3765524, P(meta) = 3.08 × 10(-10), per-allele OR = 0.80 (95% confidence interval: 0.75-0.86)). We identify two loci associated with susceptibility to DSS in people with dengue, suggesting possible mechanisms for this severe complication of dengue.     

SNPs in KCNQ1 are associated with susceptibility to type 2 diabetes in East Asian and European populations

We conducted a genome-wide association study using 207,097 SNP markers in Japanese individuals with type 2 diabetes and unrelated controls, and identified KCNQ1 (potassium voltage-gated channel, KQT-like subfamily, member 1) to be a strong candidate for conferring susceptibility to type 2 diabetes. We detected consistent association of a SNP in KCNQ1 (rs2283228) with the disease in several independent case-control studies (additive model P = 3.1 x 10(-12); OR = 1.26, 95% CI = 1.18-1.34). Several other SNPs in the same linkage disequilibrium (LD) block were strongly associated with type 2 diabetes (additive model: rs2237895, P = 7.3 x 10(-9); OR = 1.32, 95% CI = 1.20-1.45, rs2237897, P = 6.8 x 10(-13); OR = 1.41, 95% CI = 1.29-1.55). The association of these SNPs with type 2 diabetes was replicated in samples from Singaporean (additive model: rs2237895, P = 8.5 x 10(-3); OR = 1.14, rs2237897, P = 2.4 x 10(-4); OR = 1.22) and Danish populations (additive model: rs2237895, P = 3.7 x 10(-11); OR = 1.24, rs2237897, P = 1.2 x 10(-4); OR = 1.36).     

A genome-wide association study of Hodgkin's lymphoma identifies new susceptibility loci at 2p16.1 (REL), 8q24.21 and 10p14 (GATA3)

To identify susceptibility loci for classical Hodgkin's lymphoma (cHL), we conducted a genome-wide association study of 589 individuals with cHL (cases) and 5,199 controls with validation in four independent samples totaling 2,057 cases and 3,416 controls. We identified three new susceptibility loci at 2p16.1 (rs1432295, REL, odds ratio (OR) = 1.22, combined P = 1.91 × 10(-8)), 8q24.21 (rs2019960, PVT1, OR = 1.33, combined P = 1.26 × 10(-13)) and 10p14 (rs501764, GATA3, OR = 1.25, combined P = 7.05 × 10(-8)). Furthermore, we confirmed the role of the major histocompatibility complex in disease etiology by revealing a strong human leukocyte antigen (HLA) association (rs6903608, OR = 1.70, combined P = 2.84 × 10(-50)). These data provide new insight into the pathogenesis of cHL.     

A genome-wide association study in Han Chinese identifies multiple susceptibility loci for IgA nephropathy

We performed a two-stage genome-wide association study of IgA nephropathy (IgAN) in Han Chinese, with 1,434 affected individuals (cases) and 4,270 controls in the discovery phase and follow-up of the top 61 SNPs in an additional 2,703 cases and 3,464 controls. We identified associations at 17p13 (rs3803800, P = 9.40 × 10(-11), OR = 1.21; rs4227, P = 4.31 × 10(-10), OR = 1.23) and 8p23 (rs2738048, P = 3.18 × 10(-14), OR = 0.79) that implicated the genes encoding tumor necrosis factor (TNFSF13) and α-defensin (DEFA) as susceptibility genes. In addition, we found multiple associations in the major histocompatibility complex (MHC) region (rs660895, P = 4.13 × 10(-20), OR = 1.34; rs1794275, P = 3.43 × 10(-13), OR = 1.30; rs2523946, P = 1.74 × 10(-11), OR = 1.21) and confirmed a previously reported association at 22q12 (rs12537, P = 1.17 × 10(-11), OR = 0.78). We also found that rs660895 was associated with clinical subtypes of IgAN (P = 0.003), proteinuria (P = 0.025) and IgA levels (P = 0.047). Our findings show that IgAN is associated with variants near genes involved in innate immunity and inflammation.     

A genome-wide association study identifies three loci associated with susceptibility to uterine fibroids

Uterine fibroids are a common benign tumor of the female genital tract. We conducted a genome-wide association study in which 457,044 SNPs were analyzed in 1,607 individuals with clinically diagnosed uterine fibroids and 1,428 female controls. SNPs showing suggestive associations (P < 5 × 10(-5)) were further genotyped in 3,466 additional cases and 3,245 female controls. Three loci on chromosomes 10q24.33, 22q13.1 and 11p15.5 revealed genome-wide significant associations with uterine fibroids. The SNPs showing the most significant association in a combination analysis at each of these loci were rs7913069 (P = 8.65 × 10(-14), odds ratio (OR) = 1.47), rs12484776 (P = 2.79 × 10(-12), OR = 1.23) and rs2280543 (P = 3.82 × 10(-12), OR = 1.39), respectively. Subsequent fine mapping of these regions will be necessary to pinpoint the causal variants. Our findings should shed light on the pathogenesis of uterine fibroids.     

Genome-wide association scan identifies a colorectal cancer susceptibility locus on 11q23 and replicates risk loci at 8q24 and 18q21

In a genome-wide association study to identify loci associated with colorectal cancer (CRC) risk, we genotyped 555,510 SNPs in 1,012 early-onset Scottish CRC cases and 1,012 controls (phase 1). In phase 2, we genotyped the 15,008 highest-ranked SNPs in 2,057 Scottish cases and 2,111 controls. We then genotyped the five highest-ranked SNPs from the joint phase 1 and 2 analysis in 14,500 cases and 13,294 controls from seven populations, and identified a previously unreported association, rs3802842 on 11q23 (OR = 1.1; P = 5.8 x 10(-10)), showing population differences in risk. We also replicated and fine-mapped associations at 8q24 (rs7014346; OR = 1.19; P = 8.6 x 10(-26)) and 18q21 (rs4939827; OR = 1.2; P = 7.8 x 10(-28)). Risk was greater for rectal than for colon cancer for rs3802842 (P < 0.008) and rs4939827 (P < 0.009). Carrying all six possible risk alleles yielded OR = 2.6 (95% CI = 1.75-3.89) for CRC. These findings extend our understanding of the role of common genetic variation in CRC etiology.     

Meta-analysis of three genome-wide association studies identifies susceptibility loci for colorectal cancer at 1q41, 3q26.2, 12q13.13 and 20q13.33

Genome-wide association studies (GWAS) have identified ten loci harboring common variants that influence risk of developing colorectal cancer (CRC). To enhance the power to identify additional CRC risk loci, we conducted a meta-analysis of three GWAS from the UK which included a total of 3,334 affected individuals (cases) and 4,628 controls followed by multiple validation analyses including a total of 18,095 cases and 20,197 controls. We identified associations at four new CRC risk loci: 1q41 (rs6691170, odds ratio (OR) = 1.06, P = 9.55 × 10?1? and rs6687758, OR = 1.09, P = 2.27 × 10??, 3q26.2 (rs10936599, OR = 0.93, P = 3.39 × 10??), 12q13.13 (rs11169552, OR = 0.92, P = 1.89 × 10?1? and rs7136702, OR = 1.06, P = 4.02 × 10??) and 20q13.33 (rs4925386, OR = 0.93, P = 1.89 × 10?1?). In addition to identifying new CRC risk loci, this analysis provides evidence that additional CRC-associated variants of similar effect size remain to be discovered.     

Genome-wide association identifies three new susceptibility loci for Paget's disease of bone

Paget's disease of bone (PDB) is a common disorder characterized by focal abnormalities of bone remodeling. We previously identified variants at the CSF1, OPTN and TNFRSF11A loci as risk factors for PDB by genome-wide association study. Here we extended this study, identified three new loci and confirmed their association with PDB in 2,215 affected individuals (cases) and 4,370 controls from seven independent populations. The new associations were with rs5742915 within PML on 15q24 (odds ratio (OR) = 1.34, P = 1.6 × 10(-14)), rs10498635 within RIN3 on 14q32 (OR = 1.44, P = 2.55 × 10(-11)) and rs4294134 within NUP205 on 7q33 (OR = 1.45, P = 8.45 × 10(-10)). Our data also confirmed the association of TM7SF4 (rs2458413, OR = 1.40, P = 7.38 × 10(-17)) with PDB. These seven loci explained ～13% of the familial risk of PDB. These studies provide new insights into the genetic architecture and pathophysiology of PDB.     

A genome-wide association meta-analysis identifies new childhood obesity loci

Multiple genetic variants have been associated with adult obesity and a few with severe obesity in childhood; however, less progress has been made in establishing genetic influences on common early-onset obesity. We performed a North American, Australian and European collaborative meta-analysis of 14 studies consisting of 5,530 cases (≥95th percentile of body mass index (BMI)) and 8,318 controls (<50th percentile of BMI) of European ancestry. Taking forward the eight newly discovered signals yielding association with P < 5 × 10(-6) in nine independent data sets (2,818 cases and 4,083 controls), we observed two loci that yielded genome-wide significant combined P values near OLFM4 at 13q14 (rs9568856; P = 1.82 × 10(-9); odds ratio (OR) = 1.22) and within HOXB5 at 17q21 (rs9299; P = 3.54 × 10(-9); OR = 1.14). Both loci continued to show association when two extreme childhood obesity cohorts were included (2,214 cases and 2,674 controls). These two loci also yielded directionally consistent associations in a previous meta-analysis of adult BMI(1).     

Genome-wide association study identifies two new susceptibility loci for atopic dermatitis in the Chinese Han population

Atopic dermatitis is a chronic, relapsing form of inflammatory skin disorder that is affected by genetic and environmental factors. We performed a genome-wide association study of atopic dermatitis in a Chinese Han population using 1,012 affected individuals (cases) and 1,362 controls followed by a replication study in an additional 3,624 cases and 12,197 controls of Chinese Han ethnicity, as well as 1,806 cases and 3,256 controls from Germany. We identified previously undescribed susceptibility loci at 5q22.1 (TMEM232 and SLC25A46, rs7701890, P(combined) = 3.15 × 10(-9), odds ratio (OR) = 1.24) and 20q13.33 (TNFRSF6B and ZGPAT, rs6010620, P(combined) = 3.0 × 10(-8), OR = 1.17) and replicated another previously reported locus at 1q21.3 (FLG, rs3126085, P(combined) = 5.90 × 10(-12), OR = 0.82) in the Chinese sample. The 20q13.33 locus also showed evidence for association in the German sample (rs6010620, P = 2.87 × 10(-5), OR = 1.25). Our study identifies new genetic susceptibility factors and suggests previously unidentified biological pathways in atopic dermatitis.     

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.     

Discovery of common variants associated with low TSH levels and thyroid cancer risk

To search for sequence variants conferring risk of nonmedullary thyroid cancer, we focused our analysis on 22 SNPs with a P < 5 × 10(-8) in a genome-wide association study on levels of thyroid stimulating hormone (TSH) in 27,758 Icelanders. Of those, rs965513 has previously been shown to associate with thyroid cancer. The remaining 21 SNPs were genotyped in 561 Icelandic individuals with thyroid cancer (cases) and up to 40,013 controls. Variants suggestively associated with thyroid cancer (P < 0.05) were genotyped in an additional 595 non-Icelandic cases and 2,604 controls. After combining the results, three variants were shown to associate with thyroid cancer: rs966423 on 2q35 (OR = 1.34; P(combined) = 1.3 × 10(-9)), rs2439302 on 8p12 (OR = 1.36; P(combined) = 2.0 × 10(-9)) and rs116909374 on 14q13.3 (OR = 2.09; P(combined) = 4.6 × 10(-11)), a region previously reported to contain an uncorrelated variant conferring risk of thyroid cancer. A strong association (P = 9.1 × 10(-91)) was observed between rs2439302 on 8p12 and expression of NRG1, which encodes the signaling protein neuregulin 1, in blood.     

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.     

Genome-wide association study of prostate cancer identifies a second risk locus at 8q24

Recently, common variants on human chromosome 8q24 were found to be associated with prostate cancer risk. While conducting a genome-wide association study in the Cancer Genetic Markers of Susceptibility project with 550,000 SNPs in a nested case-control study (1,172 cases and 1,157 controls of European origin), we identified a new association at 8q24 with an independent effect on prostate cancer susceptibility. The most significant signal is 70 kb centromeric to the previously reported SNP, rs1447295, but shows little evidence of linkage disequilibrium with it. A combined analysis with four additional studies (total: 4,296 cases and 4,299 controls) confirms association with prostate cancer for rs6983267 in the centromeric locus (P = 9.42 x 10(-13); heterozygote odds ratio (OR): 1.26, 95% confidence interval (c.i.): 1.13-1.41; homozygote OR: 1.58, 95% c.i.: 1.40-1.78). Each SNP remained significant in a joint analysis after adjusting for the other (rs1447295 P = 1.41 x 10(-11); rs6983267 P = 6.62 x 10(-10)). These observations, combined with compelling evidence for a recombination hotspot between the two markers, indicate the presence of at least two independent loci within 8q24 that contribute to prostate cancer in men of European ancestry. We estimate that the population attributable risk of the new locus, marked by rs6983267, is higher than the locus marked by rs1447295 (21% versus 9%).     

A genome-wide association study identifies two new susceptibility loci for lung adenocarcinoma in the Japanese population

Lung adenocarcinoma is the most common histological type of lung cancer, and its incidence is increasing worldwide. To identify genetic factors influencing risk of lung adenocarcinoma, we conducted a genome-wide association study and two validation studies in the Japanese population comprising a total of 6,029 individuals with lung adenocarcinoma (cases) and 13,535 controls. We confirmed two previously reported risk loci, 5p15.33 (rs2853677, P(combined) = 2.8 × 10(-40), odds ratio (OR) = 1.41) and 3q28 (rs10937405, P(combined) = 6.9 × 10(-17), OR = 1.25), and identified two new susceptibility loci, 17q24.3 (rs7216064, P(combined) = 7.4 × 10(-11), OR = 1.20) and 6p21.3 (rs3817963, P(combined) = 2.7 × 10(-10), OR = 1.18). These data provide further evidence supporting a role for genetic susceptibility in the development of lung adenocarcinoma.     

Common variants on 8p12 and 1q24.2 confer risk of schizophrenia

Schizophrenia is a severe mental disorder affecting ～1% of the world population, with heritability of up to 80%. To identify new common genetic risk factors, we performed a genome-wide association study (GWAS) in the Han Chinese population. The discovery sample set consisted of 3,750 individuals with schizophrenia and 6,468 healthy controls (1,578 cases and 1,592 controls from northern Han Chinese, 1,238 cases and 2,856 controls from central Han Chinese, and 934 cases and 2,020 controls from the southern Han Chinese). We further analyzed the strongest association signals in an additional independent cohort of 4,383 cases and 4,539 controls from the Han Chinese population. Meta-analysis identified common SNPs that associated with schizophrenia with genome-wide significance on 8p12 (rs16887244, P = 1.27 × 10(-10)) and 1q24.2 (rs10489202, P = 9.50 × 10(-9)). Our findings provide new insights into the pathogenesis of schizophrenia.