Genome-wide association study reveals three susceptibility loci for common migraine in the general population

Migraine is a common, heterogeneous and heritable neurological disorder. Its pathophysiology is incompletely understood, and its genetic influences at the population level are unknown. In a population-based genome-wide analysis including 5,122 migraineurs and 18,108 non-migraineurs, rs2651899 (1p36.32, PRDM16), rs10166942 (2q37.1, TRPM8) and rs11172113 (12q13.3, LRP1) were among the top seven associations (P < 5 × 10(-6)) with migraine. These SNPs were significant in a meta-analysis among three replication cohorts and met genome-wide significance in a meta-analysis combining the discovery and replication cohorts (rs2651899, odds ratio (OR) = 1.11, P = 3.8 × 10(-9); rs10166942, OR = 0.85, P = 5.5 × 10(-12); and rs11172113, OR = 0.90, P = 4.3 × 10(-9)). The associations at rs2651899 and rs10166942 were specific for migraine compared with non-migraine headache. None of the three SNP associations was preferential for migraine with aura or without aura, nor were any associations specific for migraine features. TRPM8 has been the focus of neuropathic pain models, whereas LRP1 modulates neuronal glutamate signaling, plausibly linking both genes to migraine pathophysiology.     

A genome-wide association study identifies new susceptibility loci for non-cardia gastric cancer at 3q13.31 and 5p13.1

Gastric cancer, including the cardia and non-cardia types, is the second leading cause of cancer-related deaths worldwide. To identify genetic risk variants for non-cardia gastric cancer, we performed a genome-wide association study in 3,279 individuals (1,006 with non-cardia gastric cancer and 2,273 controls) of Chinese descent. We replicated significant associations in an additional 6,897 subjects (3,288 with non-cardia gastric cancer and 3,609 controls). We identified two new susceptibility loci for non-cardia gastric cancer at 5p13.1 (rs13361707 in the region including PTGER4 and PRKAA1; odds ratio (OR) = 1.41; P = 7.6 × 10(-29)) and 3q13.31 (rs9841504 in ZBTB20; OR = 0.76; P = 1.7 × 10(-9)). Imputation analyses also confirmed previously reported associations of rs2294008 and rs2976392 on 8q24, rs4072037 on 1q22 and rs13042395 on 20p13 with non-cardia gastric cancer susceptibility in the Han Chinese population.     

A genome-wide association study identifies two new lung cancer susceptibility loci at 13q12.12 and 22q12.2 in Han Chinese

Lung cancer is the leading cause of cancer-related deaths worldwide. To identify genetic factors that modify the risk of lung cancer in individuals of Chinese ancestry, we performed a genome-wide association scan in 5,408 subjects (2,331 individuals with lung cancer (cases) and 3,077 controls) followed by a two-stage validation among 12,722 subjects (6,313 cases and 6,409 controls). The combined analyses identified six well-replicated SNPs with independent effects and significant lung cancer associations (P < 5.0 × 10(-8)) located in TP63 (rs4488809 at 3q28, P = 7.2 × 10(-26)), TERT-CLPTM1L (rs465498 and rs2736100 at 5p15.33, P = 1.2 × 10(-20) and P = 1.0 × 10(-27), respectively), MIPEP-TNFRSF19 (rs753955 at 13q12.12, P = 1.5 × 10(-12)) and MTMR3-HORMAD2-LIF (rs17728461 and rs36600 at 22q12.2, P = 1.1 × 10(-11) and P = 6.2 × 10(-13), respectively). Two of these loci (13q12.12 and 22q12.2) were newly identified in the Chinese population. These results suggest that genetic variants in 3q28, 5p15.33, 13q12.12 and 22q12.2 may contribute to the susceptibility of lung cancer in Han Chinese.     

Association analyses identify multiple new lung cancer susceptibility loci and their interactions with smoking in the Chinese population

To find additional susceptibility loci for lung cancer, we tested promising associations from our previous genome-wide association study (GWAS) of lung cancer in the Chinese population in an extended validation sample size of 7,436 individuals with lung cancer (cases) and 7,483 controls. We found genome-wide significant (P < 5.0 × 10(-8)) evidence for three additional lung cancer susceptibility loci at 10p14 (rs1663689, close to GATA3, P = 2.84 × 10(-10)), 5q32 (rs2895680 in PPP2R2B-STK32A-DPYSL3, P = 6.60 × 10(-9)) and 20q13.2 (rs4809957 in CYP24A1, P = 1.20 × 10(-8)). We also found consistent associations for rs247008 at 5q31.1 (IL3-CSF2-P4HA2, P = 7.68 × 10(-8)) and rs9439519 at 1p36.32 (AJAP1-NPHP4, P = 3.65 × 10(-6)). Four of these loci showed evidence for interactions with smoking dose (P = 1.72 × 10(-10), P = 5.07 × 10(-3), P = 6.77 × 10(-3) and P = 4.49 × 10(-2) for rs2895680, rs4809957, rs247008 and rs9439519, respectively). These results advance our understanding of lung cancer susceptibility and highlight potential pathways that integrate genetic variants and smoking in the development of lung cancer.     

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.     

Genome-wide association study identifies a susceptibility locus for schizophrenia in Han Chinese at 11p11.2

To identify susceptibility loci for schizophrenia, we performed a two-stage genome-wide association study (GWAS) of schizophrenia in the Han Chinese population (GWAS: 746 individuals with schizophrenia and 1,599 healthy controls; validation: 4,027 individuals with schizophrenia and 5,603 healthy controls). We identified two susceptibility loci for schizophrenia at 6p21-p22.1 (rs1233710 in an intron of ZKSCAN4, P(combined) = 4.76 × 10(-11), odds ratio (OR) = 0.79; rs1635 in an exon of NKAPL, P(combined) = 6.91 × 10(-12), OR = 0.78; rs2142731 in an intron of PGBD1, P(combined) = 5.14 × 10(-10), OR = 0.79) and 11p11.2 (rs11038167 near the 5' UTR of TSPAN18, P(combined) = 1.09 × 10(-11), OR = 1.29; rs11038172, P(combined) = 7.21 × 10(-10), OR = 1.25; rs835784, P(combined) = 2.73 × 10(-11), OR = 1.27). These results add to previous evidence of susceptibility loci for schizophrenia at 6p21-p22.1 in the Han Chinese population. We found that NKAPL and ZKSCAN4 were expressed in postnatal day 0 (P0) mouse brain. These findings may lead to new insights into the pathogenesis of schizophrenia.     

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.     

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 Han Chinese identifies new susceptibility loci for ankylosing spondylitis

To identify susceptibility loci for ankylosing spondylitis, we performed a two-stage genome-wide association study in Han Chinese. In the discovery stage, we analyzed 1,356,350 autosomal SNPs in 1,837 individuals with ankylosing spondylitis and 4,231 controls; in the validation stage, we analyzed 30 suggestive SNPs in an additional 2,100 affected individuals and 3,496 controls. We identified two new susceptibility loci between EDIL3 and HAPLN1 at 5q14.3 (rs4552569; P = 8.77 × 10(-10)) and within ANO6 at 12q12 (rs17095830; P = 1.63 × 10(-8)). We also confirmed previously reported associations in Europeans within the major histocompatibility complex (MHC) region (top SNP, rs13202464; P < 5 × 10(-324)) and at 2p15 (rs10865331; P = 1.98 × 10(-8)). We show that rs13202464 within the MHC region mainly represents the risk effect of HLA-B*27 variants (including HLA-B*2704, HLA-B*2705 and HLA-B*2715) in Chinese. The two newly discovered loci implicate genes related to bone formation and cartilage development, suggesting their potential involvement in the etiology of ankylosing spondylitis.     

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.     

Common variation at 3p22.1 and 7p15.3 influences multiple myeloma risk

To identify risk variants for multiple myeloma, we conducted a genome-wide association study of 1,675 individuals with multiple myeloma and 5,903 control subjects. We identified risk loci for multiple myeloma at 3p22.1 (rs1052501 in ULK4; odds ratio (OR) = 1.32; P = 7.47 × 10(-9)) and 7p15.3 (rs4487645, OR = 1.38; P = 3.33 × 10(-15)). In addition, we observed a promising association at 2p23.3 (rs6746082, OR = 1.29; P = 1.22 × 10(-7)). Our study identifies new genomic regions associated with multiple myeloma risk that may lead to new etiological insights.     

Genome-wide association study identifies susceptibility loci for polycystic ovary syndrome on chromosome 2p16.3, 2p21 and 9q33.3

Polycystic ovary syndrome (PCOS) is a common metabolic disorder in women. To identify causative genes, we conducted a genome-wide association study (GWAS) of PCOS in Han Chinese. The discovery set included 744 PCOS cases and 895 controls; subsequent replications involved two independent cohorts (2,840 PCOS cases and 5,012 controls from northern Han Chinese; 498 cases and 780 controls from southern and central Han Chinese). We identified strong evidence of associations between PCOS and three loci: 2p16.3 (rs13405728; combined P-value by meta-analysis P(meta) = 7.55 × 10?21, odds ratio (OR) 0.71); 2p21 (rs13429458, P(meta) = 1.73 × 10?23, OR 0.67); and 9q33.3 (rs2479106, P(meta) = 8.12 × 10?1?, OR 1.34). These findings provide new insight into the pathogenesis of PCOS. Follow-up studies of the candidate genes in these regions are recommended.     

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 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.     

A genome-wide association study identifies two susceptibility loci for duodenal ulcer in the Japanese population

Through a genome-wide association analysis with a total of 7,035 individuals with duodenal ulcer and 25,323 controls from Japan, we identified two susceptibility loci at the PSCA gene (encoding prostate stem cell antigen) at 8q24 and at the ABO blood group locus at 9q34. The C allele of rs2294008 at PSCA was associated with increased risk of duodenal ulcer (odds ratio (OR) = 1.84; P = 3.92 × 10(-33)) in a recessive model but was associated with decreased risk of gastric cancer (OR = 0.79; P = 6.79 × 10(-12)), as reported previously. The T allele of rs2294008 encodes a translation initiation codon upstream of the reported site and changes protein localization from the cytoplasm to the cell surface. rs505922 at ABO was also associated with duodenal ulcer in a recessive model (OR = 1.32; P = 1.15 × 10(-10)). Our findings demonstrate a role for genetic variants in the pathogenesis of duodenal ulcer.     

Common variation near CDKN1A, POLD3 and SHROOM2 influences colorectal cancer risk

We performed a meta-analysis of five genome-wide association studies to identify common variants influencing colorectal cancer (CRC) risk comprising 8,682 cases and 9,649 controls. Replication analysis was performed in case-control sets totaling 21,096 cases and 19,555 controls. We identified three new CRC risk loci at 6p21 (rs1321311, near CDKN1A; P = 1.14 × 10(-10)), 11q13.4 (rs3824999, intronic to POLD3; P = 3.65 × 10(-10)) and Xp22.2 (rs5934683, near SHROOM2; P = 7.30 × 10(-10)) This brings the number of independent loci associated with CRC risk to 20 and provides further insight into the genetic architecture of inherited susceptibility to CRC.     

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.     

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.     

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.