92169042419Nat GenetNat. Genet.Nature genetics1061-40361546-171828346443555824610.1038/ng.3823NIHMS891712ArticleGenome-wide association study of glioma subtypes identifies specific differences in genetic susceptibility to glioblastoma and non-glioblastoma tumorsMelinBeatrice S141Barnholtz-SloanJill S241WrenschMargaret R3441JohansenChristoffer541Il’yasovaDora67841KinnersleyBen941OstromQuinn T2LabrecheKarim910ChenYanwen2ArmstrongGeorgina11LiuYanhong11Eckel-PassowJeanette E12DeckerPaul A12LabussièreMarianne10IdbaihAhmed1013Hoang-XuanKhe1013Di StefanoAnna-Luisa1013MokhtariKarima1013DelattreJean-Yves1013BroderickPeter9GalanPilar14GousiasKonstantinos15SchrammJohannes15SchoemakerMinouk J9FlemingSarah J16HermsStefan16HeilmannStefanie17NöthenMarkus M17WichmannHeinz-Erich181920SchreiberStefan21SwerdlowAnthony922LathropMark23SimonMatthias15SansonMarc1013AnderssonUlrika1RajaramanPreetha24ChanockStephen24LinetMartha24WangZhaoming24YeagerMeredith24GliomaScan Consortium25WienckeJohn K34HansenHelen3McCoyLucie3RiceTerri3KoselMatthew L12SicotteHugues12AmosChristopher I26BernsteinJonine L27DavisFaith28LachanceDan29LauChing30MerrellRyan T31ShildkrautJoellen78Ali-OsmanFrancis732SadetzkiSiegal3334ScheurerMichael30SheteSanjay35LaiRose K3642ClausElizabeth B373842OlsonSara H2742JenkinsRobert B3942HoulstonRichard S94042BondyMelissa L1142 Department of Radiation Sciences, Umeå University, Umeå, Sweden Case Comprehensive Cancer Center, School of Medicine, Case Western Reserve University, Cleveland, Ohio, USA Department of Neurological Surgery, School of Medicine, University of California, San Francisco, San Francisco, California, USA Institute of Human Genetics, University of California, San Francisco, San Francisco, California, USA Institute of Cancer Epidemiology, Danish Cancer Society, Copenhagen, Denmark and Rigshospitalet, University of Copenhagen, Copenhagen, Denmark Department of Epidemiology and Biostatistics, School of Public Health, Georgia State University, Atlanta, Georgia, USA Duke Cancer Institute, Duke University Medical Center, Durham, North Carolina, USA Cancer Control and Prevention Program, Department of Community and Family Medicine, Duke University Medical Center, Durham, North Carolina, USA Division of Genetics and Epidemiology, Institute of Cancer Research, London, UK Sorbonne Universités UPMC Univ Paris 06, INSERM CNRS, U1127, UMR 7225, ICM, Paris, France Department of Medicine, Dan L. Duncan Comprehensive Cancer Center, Baylor College of Medicine, Houston, Texas, USA Division of Biomedical Statistics and Informatics, Mayo Clinic College of Medicine, Rochester, Minnesota, USA AP-HP, Groupe Hospitalier Pitié-Salpêtrière, Service de neurologie 2-Mazarin, Paris, France Université Paris 13 Sorbonne Paris Cité, INSERM U557, INRA U1125, CNAM, Paris, France Department of Neurosurgery, University of Bonn Medical Center, Bonn, Germany Centre for Epidemiology and Biostatistics, Faculty of Medicine and Health, University of Leeds, Leeds, UK Institute of Human Genetics, University of Bonn, Bonn, Germany Helmholtz Center Munich, Institute of Epidemiology I, Munich, Germany Institute of Medical Informatics, Biometry and Epidemiology, Ludwig Maximilians University, Munich, Germany Institute of Medical Statistics and Epidemiology, Technical University Munich, Munich, Germany 1st Medical Department, University Clinic Schleswig–Holstein, Campus Kiel, Kiel, Germany Division of Breast Cancer Research, Institute of Cancer Research, London, UK Génome Québec, Department of Human Genetics, McGill University, Montreal, Quebec, Canada Division of Cancer Epidemiology and Genetics, National Cancer Institute, Bethesda, Maryland, USA Department of Biomedical Data Science, Geisel School of Medicine at Dartmouth, Hanover, New Hampshire, USA Department of Epidemiology and Biostatistics, Memorial Sloan Kettering Cancer Center, New York, New York, USA School of Public Health, University of Alberta, Edmonton, Alberta, Canada Department of Neurology, Mayo Clinic Comprehensive Cancer Center, Mayo Clinic, Rochester, Minnesota, USA Department of Pediatrics, Dan L. Duncan Comprehensive Cancer Center, Baylor College of Medicine, Houston, Texas, USA Department of Neurology, NorthShore University HealthSystem, Evanston, Illinois, USA Department of Surgery, Duke University Medical Center, Durham, North Carolina, USA Cancer and Radiation Epidemiology Unit, Gertner Institute, Chaim Sheba Medical Center, Tel Hashomer, Israel Department of Epidemiology and Preventive Medicine, School of Public Health, Sackler Faculty of Medicine, Tel-Aviv University, Tel-Aviv, Israel Department of Biostatistics, University of Texas Maryland Anderson Cancer Center, Houston, Texas, USA Departments of Neurology and Preventive Medicine, Keck School of Medicine, University of Southern California, Los Angeles, California, USA School of Public Health, Yale University, New Haven, Connecticut, USA Department of Neurosurgery, Brigham and Women’s Hospital, Boston, Massachusetts, USA Department of Laboratory Medicine and Pathology, Mayo Clinic Comprehensive Cancer Center, Mayo Clinic, Rochester, Minnesota, USA Division of Molecular Pathology, Institute of Cancer Research, London, UKCorrespondence should be addressed to B.S.M. (beatrice.melin@umu.se), R.S.H. (richard.houlston@icr.ac.uk) or M.L.B. (mbondy@bcm.edu)

A full list of members and affiliations appears in the Acknowledgments

These authors contributed equally to this work

These authors jointly directed this work

12720172732017520171682017495789794Reprints and permissions information is available online at http://www.nature.com/reprints/index.html.

Genome-wide association studies (GWAS) have transformed our understanding of glioma susceptibility, but individual studies have had limited power to identify risk loci. We performed a meta-analysis of existing GWAS and two new GWAS, which totaled 12,496 cases and 18,190 controls. We identified five new loci for glioblastoma (GBM) at 1p31.3 (rs12752552; P = 2.04 × 10−9, odds ratio (OR) = 1.22), 11q14.1 (rs11233250; P = 9.95 × 10−10, OR = 1.24), 16p13.3 (rs2562152; P = 1.93 × 10−8, OR = 1.21), 16q12.1 (rs10852606; P = 1.29 × 10−11, OR = 1.18) and 22q13.1 (rs2235573; P = 1.76 × 10−10, OR = 1.15), as well as eight loci for non-GBM tumors at 1q32.1 (rs4252707; P = 3.34 × 10−9, OR = 1.19), 1q44 (rs12076373; P = 2.63 × 10−10, OR = 1.23), 2q33.3 (rs7572263; P = 2.18 × 10−10, OR = 1.20), 3p14.1 (rs11706832; P = 7.66 × 10−9, OR = 1.15), 10q24.33 (rs11598018; P = 3.39 × 10−8, OR = 1.14), 11q21 (rs7107785; P = 3.87 × 10−10, OR = 1.16), 14q12 (rs10131032; P = 5.07 × 10−11, OR = 1.33) and 16p13.3 (rs3751667; P = 2.61 × 10−9, OR = 1.18). These data substantiate that genetic susceptibility to GBM and non-GBM tumors are highly distinct, which likely reflects different etiology.

Glioma accounts for around 27% of all primary brain tumors and is responsible for approximately 13,000 cancer-related deaths in the United States each year1, 2. Gliomas can be broadly classified into GBM and lower-grade non-GBM tumors3. Gliomas typically have a poor prognosis irrespective of medical care, with the most common form, GBM, having a five-year survival rate of only 5% (ref. 4).

So far, no environmental exposures have been robustly linked to the risk of developing glioma, except for moderate to high doses of ionizing radiation, which accounts for a small proportion of cases5. Evidence for an inherited predisposition to glioma is provided by a number of rare inherited cancer syndromes, such as Turcot’s and Li–Fraumeni syndromes, as well as neurofibromatosis. Even collec-tively, however, these account for little of the twofold familial risk of glioma6. Our understanding of the heritability of glioma has been transformed by recent GWAS, which have identified single-nucleotide polymorphisms (SNPs) at 13 loci influencing risk7–14.

Previous individual studies have had limited statistical power for the additional discovery of new glioma risk loci15. Therefore, to gain more comprehensive insight into glioma etiology, we performed a meta-analysis of previously published GWAS and two new GWAS, which allowed us to identify 13 new risk loci for glioma.

We analyzed GWAS SNP data that passed quality control for 12,496 cases (6,191 classified as GBM and 5,819 classified as non-GBM tumors) and 18,190 controls from eight studies with individuals of European ancestry, a new GWAS of 4,572 cases and 3,286 controls performed by the Glioma International Case Control Consortium (GICC) (Supplementary Table 1), a new GWAS of 1,591 cases and 804 controls from the University of California, San Francisco (UCSF)-Mayo, and six previously reported GWAS9, 10, 13 totaling 6,405 cases and 14,100 controls (Supplementary Table 2). To increase genomic resolution, we imputed >10 million SNPs. Quantile–quantile (Q-Q) plots for SNPs with a minor allele frequency (MAF) >1% after imputation did not show evidence of substantive overdispersion (λ = 1.02–1.10, λ90 = 1.02–1.05; Supplementary Fig. 1). We derived joint ORs and 95% confidence intervals (CIs) under a fixed-effects model for each SNP with MAF >1% and associated per-allele principal component (PCA) corrected P-values for all glioma, GBM and non-GBM cases versus those for the controls (Fig. 1).

In the combined meta-analysis, among previously published glioma risk SNPs, those for all glioma at 17p13.1 (TP53), for GBM at 5p15.33 (TERT), 7p11.2 (EGFR), 9p21.3 (CDKN2B–AS1) and 20q13.33 (RTEL1), and for non-GBM tumors at 8q24.21 (CCDC26), 11q23.2, 11q23.3 (PHLDB1) and 15q24.2 (ETFA) showed even greater evidence for association (Supplementary Fig. 2 and Supplementary Table 3). SNPs at 10q25.2 and 12q12.1 for non-GBM tumors retained genome-wide significance (i.e., P < 5.0 × 10−8). Associations at the previously reported 3q26.2 (near TERC)11 and 12q23.33 (POLR3B)10 loci for GBM did not retain statistical significance (P values for the most associated SNPs are 2.68 × 10−5 and 1.60 × 10−5, respectively; Supplementary Table 3).

In addition to previously reported loci, we identified genome-wide significant associations marking new risk loci (Table 1, Supplementary Fig. 3 and Supplementary Data 1) for GBM at 1p31.3 (rs12752552; P = 2.04 × 10−9), 11q14.1 (rs11233250; P = 9.95 × 10−10), 16p13.3 (rs2562152; P = 1.93 × 10−8), 16q12.1 (rs10852606; P = 1.29 × 10−11) and 22q13.1 (rs2235573; P = 1.76 × 10−10) and for non-GBM tumors at 1q32.1 (rs4252707; P = 3.34 × 10−9), 1q44 (rs12076373; P = 2.63 × 10−10), 2q33.3 (rs7572263; P = 2.18 × 10−10), 3p14.1 (rs11706832; P = 7.66 × 10−9), 10q24.33 (rs11598018; P = 3.39 × 10−8), 11q21 (rs7107785; P = 3.87 × 10−10), 14q12 (rs10131032; P = 5.07 × 10−11) and 16p13.3 (rs3751667; P = 2.61 × 10−9). Conditional analysis confirmed the existence of two independent association signals at 7p11.2 (EGFR) as previously reported7 but did not provide evidence for additional signals at any of the other established identified risk loci or at the 13 newly identified loci. Case-only analyses con-firmed the specificity of 11q14.1, 16p13.3 and 22q13.1 associations for GBM and of 1q44, 2q33.3, 3p14.1, 11q21 and 14q12 associations for non-GBM tumors (Fig. 2 and Supplementary Table 4). Collectively, our findings provide strong evidence for specific associations for the different glioma subtypes, consistent with their previously described distinctive molecular profiles, presumably resulting from different etiological pathways.

Across the new and known risk loci, we found a significant enrichment of overlap with enhancers in H9-Derived neuronal progenitor cells (P = 8.2 × 10−5; Supplementary Data 2). These observations support the assertion that the loci identified in the GWAS influence glioma risk through effects on neural cis regulatory networks and that they are strongly involved in transcriptional initiation and enhancement. To gain further insight into the biological basis for associations at the 13 new risk loci, we performed an expression quantitative trait loci (eQTL) analysis using RNA-seq data on ten regions of normal human brain from up to 103 individuals from the Genotype–Tissue Expression (GTEx) project16 and blood eQTL data on 5,311 individuals from Westra et al.17. We used summary-level mendelian randomization (SMR)18 analysis to test for a concordance between signals from GWAS and cis eQTL for genes within 1 Mb of the sentinel and correlated SNPs (r2 > 0.8) at each locus (Supplementary Data 3) and derived bXY statistics, which estimate the effect of gene expression on glioma risk. Additionally, for each of the risk SNPs at the 13 new loci (as well as the correlated variants), we examined published data19, 20 and made use of the online resources HaploRegv4, RegulomeDB and SeattleSeq for evidence of functional effects (Supplementary Table 5).

At 16q12.1, the GBM association signal was significantly associated with HEATR3 expression in nine of ten regions of the brain (PSMR = 3.38 × 10−6 to 6.55 × 10−10; bXY = 0.14–0.24; Supplementary Fig. 4 and Supplementary Data 3). The risk allele ‘C’ of rs10852606 that was associated with reduced HEATR3 expression was consistent with differential expression of HEATR3 being the functional basis of the 16q12.1 association. The observation that variation at 16q12.1 is associated with risk of testicular21 (rs8046148; pairwise r2 and D′ with rs10852606 of 0.67 and 1.0, respectively) and esophageal22 (rs4785204; pairwise r2 and D′ with rs10852606 of 0.16 and 1.0, respectively) cancer suggests that the locus has pleiotropic effects on tumor risk, which are compatible with generic effects as shown by the observation of a HEATR3 eQTL signal in blood (PSMR = 5.84 × 10−11; bXY = 0.30).

Similarly, significant associations between gene expression and glioma risk were observed at the GBM loci 1p31.3 (JAK1, brain cortex and cerebellar hemisphere), 16p13.3 (POLR3K, whole blood) and 22q13.1 (CTA-228A9.3, brain cerebellum; PICK1, brain hippocampus) (Supplementary Fig. 4 and Supplementary Data 3). The non-GBM association at 1q32.1 marked by rs4252707 (Supplementary Fig. 3) maps to intron 8 of the gene encoding MDM4, a p53-binding protein. The SNP rs4252707 is in strong linkage disequilibrium (LD) with rs12031912 and rs12028476 (r2 = 0.92), both of which map to the MDM4 promoter. Although no significant eQTL was shown in any brain tissue, an association with MDM4 was seen in blood (PSMR = 4.74 × 10−6; bXY = 0.31; Supplementary Fig. 4 and Supplementary Data 3). Overexpression of MDM4 is a feature in glioma tumors containing wild-type TP53 and no amplification of the MDM2 gene, consistent with MDM4 amplification being a mechanism by which the p53-dependent growth control is inactivated23.

The 1q44 association with non-GBM that is marked by rs12076373 maps to intron 8 of AKT3, whose encoded product is one of the major downstream effectors of phosphatidylinositol 3-kinase (PI3K) and is highly expressed during active neurogenesis, with haploinsufficiency causing postnatal microcephaly and agenesis of the corpus callosum24. Notably, AKT3 is hyper-expressed in glioma, thus having a role in tumor viability by activating DNA repair25. Although rs12076373 does not map to a regulatory element, the correlated SNPs rs12124113 (r2 = 0.94) and rs59953491 (r2 = 0.90) locate within an enhancer element in brain cells and tissues, including H9-derived neuronal progenitor cultured cells, cortex-derived primary cultured neurospheres and NH-A astrocytes.

The 3p14.1 association with non-GBM that is marked by rs11706832 localizes to intron 2 of LRIG1. Although we did not identify an eQTL in this gene, LRIG1 is highly expressed in the brain and is a pan-negative regulator of the epidermal growth factor receptor (EGFR) signaling pathway, which inhibits hypoxia-induced vasculogenic mimicry via EGFR–PI3K–AKT pathway suppression and epithelial-to-mesenchymal transition26. Reduced LRIG1 expression is linked to tumor aggressiveness, temozolomide resistance and radio-resistance27, 28. We have previously shown an association for glioma at EGFR (7p11.2)7, which is well established to be pivotal in both the initiation of primary GBM and the progression of lower-grade glioma to grade IV. Although speculative, our new findings now suggest a more extensive pathway involving variation at LRIG1 and AKT3.

Of particular interest is rs7572263, which maps to 2q33.3, localizes within intron 3 of C2orf80 and is 50 kb telomeric to IDH1. Mutation of IDH1 is a driver for gliomagenesis29, 30 and is responsible for the CpG island methylator (G-CIMP) phenotype31, 32. Mutations in IDH1 predominate in non-GBM glioma33, 34; therefore, the association at 2q33.3 is plausible as a basis for susceptibility to non-GBM glioma. In the absence of convincing eQTL or other functional support, this does not preclude C2orf80 or another gene mapping to the region of LD as being the functional basis for the 2q33.3 association.

The maintenance of telomeres is central to cell immortalization, and it has a central role in gliomagenesis35. We have previously shown that the risk of GBM is strongly linked to genetic variation in the telomere-related genes TERT (5p15.33) and RTEL1 (20q13.33), and possibly also TERC (3q26.2)8, 9, 11. The 10q24.33 association with non-GBM that is marked by rs11598018 lies intronic to OBFC1, which functions in a telomere-associated complex that protects telomeres independently of POT1 (ref. 36). The CST complex, whose components are encoded by OBFC1, CTC1, and TEN1, competes with shelterin for telomeric-DNA-inhibiting telomerase-based telomere extension37. The significant association between the risk of non-GBM tumors and OBFC1 variation is particularly of note in light of our recent exome-sequencing report demonstrating that rare germline loss-of-function mutations in genes that encode components of the shelterin complex are a cause of familial oligodendroglioma38. The glioma risk alleles at TERT, TERC and OBFC1 are associated with increased leukocyte telomere length, thereby supporting a relationship between genotype and biology (Supplementary Table 6)35, 39, 40. However, the RTEL1 locus is not consistent with such a postulate, and recent data that have not shown a relationship between mutations in the TERT promoter and telomere length in glioma41 raise the possibility of a role for extratelomeric effects.

The deregulation of pathways involved in telomere length and EGFR signaling are thus consistent with glioma risk being governed by pathways that are important in the longevity of glial cells, and they substantiate early observations that genetic susceptibility to GBM and non-GBM tumors is highly distinct, presumably reflecting different etiologies between GBM and non-GBM tumors (Fig. 2).

The other associations we identified mark genes with varying degrees of plausibility for having a role in glioma oncogenesis. The GBM association at 16p13.33 marked by rs2562152 localizes 3 kb telomeric to MPG, which encodes a N-methylpurine DNA glycosylase whose expression is linked to temozolomide resistance in glioma42. Although attractive as a candidate, the only genes for which there was found to be a significant association between expression and glioma risk were POLR3K and C16ORF33 in blood (Supplementary Fig. 4 and Supplementary Data 3). At 1p31.3, only JAK1 provided convincing evidence for a significant eQTL with glioma risk SNPs in brain tissue. The strongest association was shown in the cortex (PSMR = 1.61 × 10−6; bXY = 0.22; Supplementary Fig. 4 and Supplementary Data 3), with the risk allele ‘T’ of rs12752552 showing increased JAK1 expression. The cis-eQTL signal for JAK1 in the cortex maps from 65.3 Mb to 65.35 Mb and shows a consistent direction of effect with the glioma-associated SNPs. JAK1–STAT6 signaling is increasingly being recognized to be relevant in glioma progression43. Hence, although JAK1 remains an attractive candidate mechanistic basis for the glioma association at 1p31.3, we cannot exclude the possibility that the cluster of SNPs between 65.3 Mb and 65.35 Mb contains the true causal variant. In the absence of functional data, potential target genes for associations at 11q14.1 (GBM), 16p13.3 (non-GBM), 11q21 (non-GBM) and 14q12 (non-GBM) remain to be elucidated.

In conclusion, we have performed the largest glioma GWAS to date and have identified 13 new glioma risk loci, thereby providing further evidence for a polygenic basis of genetic susceptibility to glioma. Histological classification of glioma is, in part, being superseded by molecular profiling34, 44; hence, it is important to understand the biology behind these risk variants in the context of molecularly defined glioma subtypes. Currently identified risk SNPs for glioma account for, at best, ~27% and ~37% of the familial risk of GBM and non-GBM tumors, respectively (Supplementary Table 7). Therefore, further GWAS-based analyses in concert with functional analyses should lead to additional insights into the biology and etiological basis of the different glioma histologies. Notably, such information can inform gene discovery initiatives and thus have a measurable effect on the successful development of new therapeutic agents.

ONLINE METHODSEthics

Collection of patient samples and associated clinico-pathological information was undertaken with written informed consent and relevant ethical review board approval at the respective study centers in accordance with the tenets of the Declaration of Helsinki. Specifically informed consent and ethical board approval was obtained from the South-East Multicentre Research Ethics Committee (MREC) (UK), the Scottish MREC (UK), the APHP ethical committee-CPP (Comité de Protection des Personnes) (France), the Ethics Commission of the Medical Faculty of the University of Bonn (Germany), the University of Texas MD Anderson Cancer Institutional Review Board (USA), the Mayo Clinic Office for Human Research Protection (USA), the UCSF Committee on Human Research (USA), the University Hospitals of Cleveland Institutional Review Board (USA) and the Cleveland Clinic Institutional Review Board (board for the Case Comprehensive Cancer Center) (USA). The diagnosis of glioma (ICDO-3 codes 9380-9480 or equivalent) was established through histology in all cases in accordance with World Health Organization guidelines. Every effort was made to classify tumors as GBM or non-GBM.

GWAS data setsGICC, UK, French, German, MDA, SFAGS and GliomaScan

Studies participating in GICC are described in Amirian et al.46 and in Supplementary Table 1. Briefly, they comprise 5,189 glioma cases and 3,827 controls that were ascertained through centers in the USA, Denmark, Sweden and the UK. Cases had newly diagnosed glioma, and controls had no personal history of central nervous system tumor at the time of ascertainment. Detailed information regarding recruitment protocol is given in Amirian et al.46. Cases and controls were genotyped using the Illumina Oncoarray according to the manufacturer’s recommendations (Illumina Inc.). Individuals with a call rate <99%, as well as all individuals evaluated to be of non-European ancestry (<80% estimated European ancestry using the FastPop47 procedure developed by the GAMEON consortium with HapMap version 2 CEU, JPT/CHB and YRI populations as a reference; Supplementary Fig. 5), were excluded. For pairs of apparent first-degree relatives, we removed the control from a case–control pair; otherwise, we excluded the individual with the lower call rate. SNPs with a call rate <95% were excluded as were those with a MAF <0.01 or those displaying significant deviation from the Hardy–Weinberg equilibrium (HWE) (i.e., P < 10−5). After performing these quality-control measures, there were 4,572 cases and 3,286 controls remaining for downstream analyses.

The UK, French, German, MDA, SFAGS and GliomaScan GWAS of non-overlapping case–control series of Northern European ancestry have been the subject of previous studies. Briefly, the UK GWAS7, 8, 10 was based on 636 cases (401 males; mean age 46 years) who were ascertained through the INTERPHONE study48. Individuals from the 1958 Birth Cohort (n = 2,930) served as a source of controls. The French GWAS7, 10 comprised 1,495 patients with glioma who were ascertained through the Service de Neurologie Mazarin, Groupe Hospitalier Pitié-Salpêtrière Paris. The controls (n = 1,213) were ascertained from the SU.VI.MAX (Supplementation en Vitamines et MinerauxAntioXydants) study of 12,735 healthy subjects (women aged 35–60 years; men aged 45–60 years)49. The German GWAS10 comprised 880 patients who had undergone surgery for a glioma at the Department of Neurosurgery, University of Bonn Medical Center, between 1996 and 2008. Control subjects were taken from three population studies: KORA (Co-operative Health Research in the Region of Augsburg; n = 488)50; POPGEN (Population Genetic Cohort; n = 678)51 and the Heinz Nixdorf Recall study (n = 380)52. Standard quality-control measures were applied to the UK, French and German GWAS and have previously been reported. The MDA GWAS8 was based on 1,281 cases (786 males; mean age 47 years) who were ascertained through the MD Anderson Cancer Center, Texas, between 1990 and 2008. Individuals from the Cancer Genetic Markers of Susceptibility (CGEMS, n = 2,245) studies served as controls53, 54. Quality-control measures were applied as per the primary GWAS. The UCSF adult glioma case–control study (SFAGS–GWAS) included participants of the San Francisco Bay Area Adult Glioma Study (AGS). Details of subject recruitment for AGS have been reported previously9, 12, 34, 55, 56. Briefly, cases were adults (>18 years of age) with newly diagnosed, histologically confirmed glioma. Population-based cases who were diagnosed between 1991 and 2009 (series 1–4) and who were residing in the six San Francisco Bay area counties were ascertained using the Cancer Prevention Institute of California’s early-case ascertainment system. Clinic-based cases who were diagnosed between 2002 and 2012 (series 3–5) were recruited from the UCSF Neuro-oncology Clinic, regardless of the place of residence. From 1991 to 2010, population-based controls from the same residential area as the population-based cases were identified using random digit-dialing and were frequency matched to population-based cases for age, gender and ethnicity. Between 2010 and 2012, all controls were selected from the UCSF general medicine phlebotomy clinic. Clinic-based controls were matched to clinic-based glioma cases for age, gender and ethnicity. Consenting participants provided blood, buccal and/or saliva specimens, and information, during in-person or telephone interviews. A total of 677 cases and 3,940 controls (including 3,347 Illumina iControlDB iControls) were used in the current analysis. For the GliomaScan GWAS13, in addition to the published analysis, we excluded samples from the ATBC (Finnish study) and controls from NSHDS due to exhibiting outlying population ancestry after manual inspection of PCA plots. In total 1,653 cases and 2,725 controls were used in the current study.

GWAS data from the seven studies were imputed to >10 million SNPs with IMPUTE2 (v2.3)57 software using a merged reference panel consisting of data from the 1000 Genomes Project (phase 1 integrated release 3, March 2012)58 and UK10K (ALSPAC, EGAS00001000090 and EGAD00001000195, and TwinsUK EGAS00001000108 and EGAS00001000194 studies). Genotypes were aligned to the positive strand in both imputation and genotyping. Imputation was conducted separately for each study, and in each the data were pruned to a common set of SNPs between cases and controls before imputation. We set thresholds for imputation quality to retain potential risk variants with MAF > 0.01. Poorly imputed SNPs, defined by an information measure <0.40 with IMPUTE2, were excluded, as were SNPs exhibiting a significant deviation from Hardy–Weinberg equilibrium (P < 1 × 10−8) in controls. Test of association between imputed SNPs and glioma was performed using SNPTEST (v2.5)59 under an additive frequentist model. The adequacy of the case–control matching and the possibility of differential genotyping of cases and controls were formally evaluated using Q-Q plots of test statistics (Supplementary Fig. 1). Where appropriate, principal components, generated using common SNPs, were included in the analysis to limit the effects of cryptic population stratification that otherwise might cause inflation of test statistics. Principal components, based on genotyped SNPs, were generated for the GICC, GliomaScan, MDA-GWAS and SFAGS studies using PLINK60. Eigenvectors for the German GWAS were inferred using smartpca (part of EIGENSOFTv2.4)61 by merging cases and controls with Phase II HapMap samples10. PCA plots for all studies are provided in Supplementary Figure 4.

UCSF-Mayo GWAS

The UCSF-Mayo study comprised Mayo cases (n = 945) and UCSF cases (n = 574) and Mayo Clinic Biobank control (n = 806) data. The Mayo Clinic case–control study has been described previously9, 34, 62. Briefly, adult cases (>18 years of age) were identified at diagnosis (diagnosed at Mayo Clinic) or at pathologic confirmation (diagnosed elsewhere and treated at Mayo Clinic), and the patients had a surgical resection or biopsy between 1973 and 2014. Consenting participants provided blood, buccal and/or saliva specimens, and information, during in-person or telephone interviews. This analysis used 574 non-overlapping cases from the UCSF Adult Glioma Study described above. Mayo Clinic and UCSF cases were genotyped using the Illumina Oncoarray. The Mayo Clinic Biobank controls comprised volunteers who donated biological specimens and provided risk factor data, access to clinical data obtained from the medical record and consent to participate in any study approved by the Access Committee. Recruitment for the Mayo Clinic Biobank took place from April 2009 through December 2015. Although participants could be unselected volunteers, the vast majority of participants were contacted as part of a pre-scheduled medical examination in the Department of Medicine, Divisions of Community Internal Medicine, Family Medicine and General Internal Medicine at Mayo Clinic sites in Rochester (Minnesota), Jacksonville (Florida), and the Mayo Clinic Health System sites in La Crosse and Onalaska (Wisconsin). All individuals were aged 18 years and older at the time of consent. Illumina Omni Express genotyping arrays were run on the 806 Mayo Clinic Biobank participants.

Quality-control analyses were performed on each cohort separately (Mayo cases, UCSF cases and Mayo Clinic Biobank controls). SNPs with call rates <95% were removed, followed by removal of subjects with call rates <95%. Concordance of replicate samples was assessed, and the sample with the higher call rate was retained. Subject’s sex was verified using the sex check option in PLINK. Relationship checking was performed by estimating the proportion of alleles shared identical by descent (IBD) for all pairs of subjects in PLINK60. STRUCTURE63 was used to assess population admixture with 1000 Genomes as a reference. Subjects indicated to be non-Caucasian were excluded. Prior to imputation, SNPs were tested for HWE, and SNPs with HWE P < 10−6 were removed. Mayo Clinic, UCSF and Mayo Clinic Biobank SNP data were each phased and imputed using the Michigan Imputation Server with the Haplotype Reference Consortium (release 1; http://www.haplotype-reference-consortium.org) as reference. Genotypes were forward-strand-aligned to the 1000 Genomes reference, and for ambiguous SNPs the Browning strand checking utility was used (http://faculty.washington.edu/sguy/beagle/strand_switching/strand_switching.html). PCA was used to correct for population stratification using SNPs common to cases and controls. The first three principal components were significantly (P < 0.05) associated with case–control status. An additive logistic regression model was used to assess the association between each SNP and disease status, with genotype being coded as 0, 1 or 2 copies of the minor allele, adjusted for age, sex and the first three principal components.

Meta-analysis and additional statistical analyses

Meta-analyses were performed using the fixed-effects inverse-variance method based on the β-esti-mates and standard errors from each study using META (v1.6)64. Cochran’s Q-statistic was used to test for heterogeneity, and the I2 statistic was used to quantify the proportion of the total variation due to heterogeneity65, taking I2 values >75 to indicate significant heterogeneity. Using the meta-analysis summary statistics and LD correlations from a reference panel of the 1000 Genomes Project combined with UK10K, we used GCTA66, 67 to perform conditional association analysis. Association statistics were calculated for all SNPs, conditioning on the top SNP in each locus showing genome-wide significance. This was carried out in a step-wise fashion. We performed a case-only analysis to test for differences in SNP-risk-allele frequency between GBM and non-GBM tumors.

ENCODE and chromatin state dynamics

Risk SNPs and their proxies (i.e., r2 > 0.8 in the 1000 Genomes EUR reference panel) were annotated for putative functional effect using HaploReg (v4)68, RegulomeDB69 and SeattleSeq Annotation70. These servers make use of data from ENCODE, genomic evolutionary rate profiling (GERP) conservation metrics, combined annotation-dependent depletion (CADD) scores and PolyPhen scores. We searched for overlap of associated SNPs with enhancers defined by the FANTOM5 enhancer atlas19, annotating by overlap with ubiquitous, permissive and robust enhancers, as well as enhancer–promoter correlations and enhancers specifically expressed in astrocytes, neuronal stem cells and brain tissue. Similarly, we searched for overlap with ‘super-enhancer’ regions, as defined by Hnisz et al.20, restricting analysis to data from U87 GBM cells, astrocyte cells and brain tissue. We additionally made use of 15-state chromHMM data from H1- and H9-derived neuronal progenitor cells available from the Epigenome Roadmap Project71. Enhancer enrichment analysis was carried out using HaploReg (v4.0)68. Briefly, from a query list of variants, the overlap with enhancers in each of 107 cell types, as predicted from the Roadmap Epigenomics Project71 chromatin-state segmentations, was calculated. A binomial test for enrichment was performed against a background set of all (i) 1000 Genomes variants with MAF > 0.05 and (ii) all unique GWAS loci in the European population. We applied a cutoff of P < 3.94 × 10−4 corresponding to a Bonferroni correction for 127 cell lines and tissues.

Expression quantitative trait loci (eQTL) analysis

To examine the relationship between SNP genotype and gene expression, we carried out summary-data-based mendelian randomization (SMR) analysis as per Zhu et al.18 (at http://cnsgenomics.com/software/smr/index.html). We used publicly available brain tissue data from the GTEx16 (http://www.gtexportal.org) v6p release. Briefly, GWAS summary statistics files were generated from the meta-analysis. Reference files were generated from merging 1000 Genomes phase 3 and UK10K (ALSPAC and TwinsUK) vcfs. Summary eQTL files for GTEx samples were generated from downloaded v6p “all_snpgene_pairs” files. Besd files were generated from these summary eQTL files using the –make-besd command. Additionally, we analyzed downloaded whole-blood eQTL data from Westra et al.17. Results from the SMR test for each of the 13 new glioma loci are reported in Supplementary Data 3. As previously advocated18, only probes with at least one eQTL P value <5.0 × 10−8 were considered for SMR analysis. We set a threshold for the SMR test of PSMR < 1.06 × 10−4 corresponding to a Bonferroni correction for 473 tests (473 probes with a top eQTL P < 5.0 × 10−8 across the 13 loci, 10 brain regions and Westra data set). For all genes passing this threshold, we generated plots of the eQTL and GWAS associations at the locus, as well as plots of GWAS and eQTL effect sizes (i.e., corresponding to input for the HEIDI heterogeneity test). HEIDI test P values <0.05 were taken to indicate significant heterogeneity. Respective SMR plots for significant eQTLs are shown in Supplementary Figure 4.

Additional statistical and bioinformatics analysis

Estimates of individual variance in risk associated with glioma risk SNPs was carried out using the method described in Pharoah et al.72, assuming the familial risk of high-grade and low-grade glioma to be 1.76 and 1.54, respectively, from analysis of the Swedish series in Scheurer et al.73. Briefly, for a single allele (i) of frequency p, relative risk R and ln risk r, the variance (Vi) of the risk distribution due to that allele is given by: Vi=(1−p)2E2+2p(1−p)(r−E)2+p2(2r−E)2

Where E is the expected value of r given by: E=2p(1−p)r+2p2r

For multiple risk alleles, the distribution of risk in the population tends toward the normal with variance: V=∑Vi

The total genetic variance (V) for all susceptibility alleles has been estimated to be √1.77. Thus, the fraction of the genetic risk explained by a single allele is given by: Vi/V

LD metrics were calculated in vcftools (v0.1.12b)74 using UK10K data and plotted using visPIG75. LD blocks were defined on the basis of HapMap recombination rate (cM/Mb), as defined using the Oxford recombination hotspots and on the basis of distribution of confidence intervals defined by Gabriel et al.76.

Data availability

Genotype data from the GICC GWAS are available from the database of Genotypes and Phenotypes (dbGaP) under accession phs001319.v1.p1. Additionally, genotypes from the GliomaScan GWAS can be accessed through dbGaP accession phs000652.v1.p1. Data from the other studies are available upon request.

Supplementary Material

Note: Any Supplementary Information and Source Data files are available in the online version of the paper.

AUTHOR CONTRIBUTIONS

M.L.B., B.S.M., R.S.H. and J.S.B.-S. managed the project; R.S.H., M.L.B., B.S.M., J.S.B.-S., R.B.J., Q.T.O., B.K. and M.R.W. drafted the manuscript; Q.T.O., K.L., B.K., J.E.E.-P. and P.A.D. performed statistical analyses; Y.C., K.L., Y.L. and B.K. performed bioinformatics analyses; B.S.M., J.S.B.-S., M.R.W., J.K.W., C.J., D.I., R.K.L., G.A., P.A.D., U.A., T.R., H.H., L.M., M.L.K., H.S., J.L.B., F.D., D.L, C.I.A, C.L., R.T.M., J. Shildkraut, F.A.-O., S. Sadetski, M. Scheurer, S. Shete, E.B.C., S.H.O., R.B.J., R.S.H. and M.L.B. developed the GICC protocol and performed sample acquisition; and P.R., S.C., M. Linet, Z.W. and M.Y. provided the National Cancer Institute (NCI) data. In the UK, P.B., A.S., M.J.S., S.J.F. and R.S.H. developed patient recruitment, performed sample acquisition and performed sample collection of cases; P.B. oversaw DNA isolation and storage, and performed case and control ascertainment, and supervision of DNA extractions. In Germany, M. Simon, M.M.N., H.-E.W., S. Schreiber and J. Schramm developed patient recruitment, and oversaw performed blood sample collection; M. Simon oversaw DNA isolation and storage and performed case and control ascertainment, and supervision of DNA extractions; and S. Herms, S. Heilmann and K.G. performed experimental work. In France, M. Sanson and J.-Y.D. developed patient recruitment; M. Labussière, A.-L.D.S., P.G., K.M., A.I. and K.H.-X. performed patient ascertainment. M. Lathrop performed laboratory management and oversaw genotyping of the French samples. All authors contributed to the final manuscript.

COMPETING FINANCIAL INTERESTS

The authors declare no competing financial interests.

We are grateful to all of the patients and individuals for their participation, and we would also like to thank the clinicians and other hospital staff members, cancer registries and the study staff members in the respective centers who contributed to the blood sample and data collection.

The GICC was supported by grants from the US National Institutes of Health (NIH) (R01CA139020 (M.L.B. and B.S.M.), R01CA52689 (M.R.W.), R01CA52689 (M.L.B.) and P30CA125123 (M. Scheurer). Additional support was provided by the McNair Medical Institute (M. Scheurer) and the Population Sciences Biorepository at Baylor College of Medicine (M. Scheurer).

In Sweden, work was additionally supported by Acta Oncologica through the Royal Swedish Academy of Science (B.S.M.’s salary) and by the Swedish Research Council (B.S.M.) and the Swedish Cancer Foundation (B.S.M.). We are grateful to the National Clinical Brain Tumor Group and to all of the clinicians and research nurses throughout Sweden who identified all of the cases.

In the UK, funding was provided by Cancer Research UK (C1298/A8362 supported by the Bobby Moore Fund (R.S.H., B.K. and P.B.), the Wellcome Trust (R.S.H., B.K. and P.B.) and the DJ Fielding Medical Research Trust (R.S.H., B.K. and P.B.). The National Brain Tumor Study is supported by the National Cancer Research Network, and we acknowledge all clinicians and healthcare professionals who contributed to this initiative. The UK INTERPHONE study was supported by the European Union Fifth Framework Program ‘Quality of Life and Management of Living Resources’ (QLK4-CT-1999-01563) (A.S., M.J.S. and S.J.F.) and the International Union against Cancer (UICC) (A.S., M.J.S. and S.J.F.). The UICC received funds from the Mobile Manufacturers’ Forum and the GSM Association. Provision of funds via the UICC was governed by agreements that guaranteed INTERPHONE’s scientific independence (http://www.iarc.fr/ENG/Units/RCAd.html), and the views expressed in the paper are not necessarily those of the funders. The UK centers were also supported by the Mobile Telecommunications and Health Research (MTHR) Programme, and the Northern UK Centre (A.S., M.J.S. and S.J.F.) was supported by the Health and Safety Executive, Department of Health and Safety Executive and the UK Network Operators.

In France, funding was provided by the Ligue Nationale Contre le Cancer (J.-Y.D.), the Fondation ARC (M. Sanson), the Institut National du Cancer (INCa; PL046; (M. Sanson)), the French Ministry of Higher Education and Research and the program “Investissements d’avenir” ANR-10-IAIHU-06 (M. Sanson, J.-Y.D., M. Labussière, A.-L.D.S., P.G., K.M., A.I., K.H.-X. and K.L.). This study was additionally supported by a grant from Génome Québec, le Ministère de l’Enseignement supérieur, de la Recherche, de la Science et de la Technologie (MESRST) Québec (M. Lathrop) and McGill University (M. Lathrop).

In Germany, funding was provided to M. Simon and J. Schramm by the Deutsche Forschungsgemeinschaft (Si552, Schr285), the Deutsche Krebshilfe (70-2385-Wi2, 70-3163-Wi3, 10-6262) and BONFOR. Funding for the WTCCC was provided by the Wellcome Trust (076113 and 085475; M. Simon and J. Schramm). The KORA Ausburg studies are supported by grants from the German Federal Ministry of Education and Research (BMBF) and were mainly financed by the Helmholtz Zentrum München, German Research Center for Environmental Health, Neuherberg. This work was financed by the German National Genome Research Network (NGFN) (S. Schreiber and H.E.-W.) and supported within the Munich Center of Health Sciences (MC Health) as part of LMUinnovativ (S. Schreiber and H.E.-W.). Generation of the German control data was partially supported by a grant of the German Federal Ministry of Education and Research (BMBF) through the Integrated Network IntegraMent (Integrated Understanding of Causes and Mechanisms in Mental Disorders), under the auspices of the e:Med research and funding concept (01ZX1314A) (M.M.N., S. Herms and S. Heilmann). M.M.N. is a member of the DFG-funded Excellence Cluster ImmunoSensation and received support from the Alfried Krupp von Bohlen und Halbach-Stiftung.

For the UK GWAS, we acknowledge the funders, organizations and individuals who contributed to the blood sample and data collection as listed in Hepworth et al.45. MD Anderson acknowledges the work of P. Adatto, F. Morice, H. Zhang, V. Levin, A.W.K. Yung, M. Gilbert, R. Sawaya, V. Puduvalli, C. Conrad, F. Lang and J. Weinberg from the Brain and Spine Center for the MDA GWAS. For the French study, we are indebted to A. Rahimian (Onconeurotek), A.M. Lekieffre and M. Brandel for help in collecting data and to Y. Marie for database support. For the German study, we are indebted to B. Harzheim (Bonn), S. Ott and A. Müller-Erkwoh (Bonn) for help with the acquisition of clinical data and to R. Mahlberg (Bonn), who provided technical support. The UK study made use of control genotyping data generated by the Wellcome Trust Case–Control Consortium. A full list of the investigators who contributed to the generation of the data is available from http://www.wtccc.org.uk. The MDA GWAS made use of control genotypes from the CGEMS prostate and breast cancer studies. A full list of the investigators who contributed to the generation of the data is available from http://cgems.cancer.gov/. French controls were taken from the SU.VI.MAX study. The German GWAS made use of genotyping data from three population control sources: KORA-gen39, the Heinz-Nixdorf RECALL study and POPGEN. The HNR cohort was established with the support of the Heinz-Nixdorf Foundation. F.D. received support from the BONFOR Programme of the University of Bonn, Germany.

The UCSF Adult Glioma Study was supported by the NIH (grant numbers R01CA52689 (M.R.W. and J.K.W.), P50CA097257 (M.R.W. and J.K.W.), R01CA126831 (J.K.W.) and R01CA139020 (M.R.W.)), the Loglio Collective (M.R.W. and J.K.W.), the National Brain Tumor Foundation (M.R.W.), the Stanley D. Lewis and Virginia S. Lewis Endowed Chair in Brain Tumor Research (M.R.W.), the Robert Magnin Newman Endowed Chair in Neuro-oncology (J.K.W.) and by donations from the families and friends of J. Berardi, H. Glaser, E. Olsen, R.E. Cooper and W. Martinusen. This project also was supported by the National Center for Research Resources and the National Center for Advancing Translational Sciences, NIH, through UCSF–CTSI grant UL1 RR024131 (UCSF CTSI). The contents of this work are solely the responsibility of the authors and do not necessarily represent the official views of the NIH. The collection of cancer incidence data used in this study was supported by the California Department of Public Health as part of the statewide cancer reporting program mandated by California Health and Safety Code section 103885, the National Cancer Institute’s Surveillance, Epidemiology and End Results Program under contract HHSN261201000140C (awarded to the Cancer Prevention Institute of California), contract HHSN261201000035C (awarded to the University of Southern California) and contract HHSN261201000034C (awarded to the Public Health Institute), and the Centers for Disease Control and Prevention’s National Program of Cancer Registries under agreement # U58DP003862-01 (awarded to the California Department of Public Health). The ideas and opinions expressed herein are those of the author(s), and endorsement by the State of California Department of Public Health, the National Cancer Institute and the Centers for Disease Control and Prevention, or their contractors and subcontractors, is not intended nor should be inferred. Other significant contributors for the UCSF Adult Glioma Study include M. Berger, P. Bracci, S. Chang, J. Clarke, A. Molinaro, A. Perry, M. Pezmecki, M. Prados, I. Smirnov, T. Tihan, K. Walsh, J. Wiemels and S. Zheng.

At Mayo, the authors wish to acknowledge the study participants and the clinicians and research staff at the participating medical centers, the Mayo Clinic Biobank and Biospecimens Accessioning and Processing Shared Resource (in particular its manager, M. Cicek). Work at the Mayo Clinic beyond the GICC was also supported by the NIH (grants P50CA108961 (B. O’Neill) and P30CA15083 (R. Diasio)), the National Institute of Neurological Disorders and Stroke (grant RC1NS068222Z (R.B.J.)), the Bernie and Edith Waterman Foundation (R.B.J.) and the Ting Tsung and Wei Fong Chao Family Foundation (R.B.J.).

The GliomaScan Consortium comprised (apart from authors listed in the author list): L.E.B. Freeman (Division of Cancer Epidemiology and Genetics, National Cancer Institute, Rockville, Maryland, USA), S. Koutros (Division of Cancer Epidemiology and Genetics, National Cancer Institute, Rockville, Maryland, USA), D. Albanes (Division of Cancer Epidemiology and Genetics, National Cancer Institute, Rockville, Maryland, USA), K. Visvanathan (Johns Hopkins Bloomberg School of Public Health, Baltimore, Maryland, USA and Johns Hopkins Sidney Kimmel Comprehensive Cancer Center, Baltimore, Maryland, USA), V.L. Stevens (Epidemiology Research Program, American Cancer Society, Atlanta, Georgia, USA), R. Henriksson (Department of Radiation Sciences, Oncology, Umea University, Umea, Sweden), D.S. Michaud (Department of Public Health and Community Medicine, Tufts University Medical School, Boston, Massachusetts, USA), M. Feychting (Unit of Epidemiology, Institute of Environmental Medicine, Karolinska Institutet, Stockholm, Sweden), A. Ahlbom (Unit of Epidemiology, Institute of Environmental Medicine, Karolinska Institutet, Stockholm, Sweden), G.G. Giles (Cancer Epidemiology Centre, Cancer Council of Victoria, Melbourne, Victoria, Australia and Centre for Molecular, Environmental, Genetic and Analytic Epidemiology, University of Melbourne, Melbourne, Victoria, Australia), R. Milne (Cancer Epidemiology Centre, Cancer Council of Victoria, Melbourne, Victoria, Australia and Centre for Molecular, Environmental, Genetic and Analytic Epidemiology, University of Melbourne, Melbourne, Victoria, Australia), R. McKean-Cowdin (Department of Preventive Medicine, Keck School of Medicine, University of Southern California, Los Angeles, California, USA), L. Le Marchand (Cancer Research Center, University of Hawaii, Honolulu, Hawaii, USA), M. Stampfer (Channing Laboratory, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, USA and Departments of Epidemiology and Nutrition, Harvard School of Public Health, Boston, Massachusetts, USA), A.M. Ruder (National Institute for Occupational Safety and Health, Centers for Disease Control and Prevention, Cincinnati, Ohio, USA), T. Carreon (National Institute for Occupational Safety and Health, Centers for Disease Control and Prevention, Cincinnati, Ohio, USA), G. Hallmans (Department of Public Health and Clinical Medicine/Nutritional Research, Umea University, Umea, Sweden), A. Zeleniuch-Jacquotte (Division of Epidemiology, Department of Environmental Medicine, New York University School of Medicine, New York, New York, USA), J.M. Gaziano (Massachusetts Veteran’s Epidemiology, Research and Information Center, Geriatric Research Education and Clinical Center, VA Boston Healthcare System, Boston, Massachusetts, USA), H.D. Sesso (Division of Preventive Medicine, Brigham and Women’s Hospital and Harvard Medical School, Boston, Massachusetts, USA), M.P. Purdue (Division of Cancer Epidemiology and Genetics, National Cancer Institute, Rockville, Maryland, USA), E. White (Fred Hutchinson Cancer Research Center, Seattle, Washington, USA and Department of Epidemiology, University of Washington, Seattle, Washington, USA) and J. Buring (Division of Preventive Medicine, Brigham and Women’s Hospital and Harvard Medical School, Boston, Massachusetts, USA).

UK10K data generation and access was organized by the UK10K consortium and funded by the Wellcome Trust.

BondyMLBrain tumor epidemiology: consensus from the Brain Tumor Epidemiology ConsortiumCancer11319531968200818798534OstromQTCBTRUS statistical report: primary brain and central nervous system tumors diagnosed in the United States in 2008–2012Neuro-oncol17Suppl. 4iv1iv62201526511214LouisDNThe 2007 WHO classification of tumors of the central nervous systemActa Neuropathol11497109200717618441OstromQTCBTRUS statistical report: primary brain and central nervous system tumors diagnosed in the United States in 2007–2011Neuro-oncol16iv1iv63201425304271OstromQTThe epidemiology of glioma in adults: a ‘state of the science’ reviewNeuro-oncol16896913201424842956HemminkiKTretliSSundquistJJohannesenTBGranstromCFamilial risks in nervous-system tumors: a histology-specific analysis from Sweden and NorwayLancet Oncol10481488200919356978SansonMChromosome 7p11.2 (EGFR) variation influences glioma riskHum. Mol. Genet2028972904201121531791SheteSGenome-wide association study identifies five susceptibility loci for gliomaNat. Genet41899904200919578367WrenschMVariants in the CDKN2B and RTEL1 regions are associated with high-grade glioma susceptibilityNat. Genet41905908200919578366KinnersleyBGenome-wide association study identifies multiple susceptibility loci for gliomaNat. Commun68559201526424050WalshKMVariants near TERT and TERC influencing telomere length are associated with high-grade glioma riskNat. Genet46731735201424908248JenkinsRBA low-frequency variant at 8q24.21 is strongly associated with risk of oligodendroglial tumors and astrocytomas with IDH1 or IDH2 mutationNat. Genet4411221125201222922872RajaramanPGenome-wide association study of glioma and meta-analysisHum. Genet13118771888201222886559StaceySNA germline variant in the TP53 polyadenylation signal confers cancer susceptibilityNat. Genet4310981103201121946351KinnersleyBQuantifying the heritability of glioma using genome-wide complex trait analysisSci. Rep517267201526625949LonsdaleJThe Genotype–Tissue Expression (GTEx) projectNat. Genet45580585201323715323WestraHJSystematic identification of trans eQTLs as putative drivers of known disease associationsNat. Genet4512381243201324013639ZhuZIntegration of summary data from GWAS and eQTL studies predicts complex-trait gene targetsNat. Genet48481487201627019110AnderssonRAn atlas of active enhancers across human cell types and tissuesNature507455461201424670763HniszDSuper-enhancers in the control of cell identity and diseaseCell155934947201324119843RuarkEIdentification of nine new susceptibility loci for testicular cancer, including variants near DAZL and PRDM14Nat. Genet45686689201323666240WuCGenome-wide association analyses of esophageal squamous cell carcinoma in Chinese identify multiple susceptibility loci and gene–environment interactionsNat. Genet4410901097201222960999RiemenschneiderMJAmplification and overexpression of the MDM4 (MDMX) gene from 1q32 in a subset of malignant gliomas without TP53 mutation or MDM2 amplificationCancer Res5960916096199910626796BolandEMapping of deletion and translocation breakpoints in 1q44 implicates the serine–threonine kinase AKT3 in postnatal microcephaly and agenesis of the corpus callosumAm. J. Hum. Genet81292303200717668379TurnerKMGenomically amplified Akt3 activates DNA repair pathway and promotes glioma progressionProc. Natl. Acad. Sci. USA11234213426201525737557GurGLRIG1 restricts growth factor signaling by enhancing receptor ubiquitylation and degradationEMBO J2332703281200415282549YangJALRIG1 enhances the radio-sensitivity of radio-resistant human glioblastoma U251 cells via attenuation of the EGFR–AKT signaling pathwayInt. J. Clin. Exp. Pathol835803590201526097540WeiJmiR-20a mediates temozolomide resistance in glioblastoma cells via negatively regulating LRIG1 expressionBiomed. Pharmacother71112118201525960225WatanabeTNobusawaSKleihuesPOhgakiHIDH1 mutations are early events in the development of astrocytomas and oligodendrogliomasAm. J. Pathol17411491153200919246647YanHIDH1 and IDH2 mutations in gliomasN. Engl. J. Med360765773200919228619NoushmehrHIdentification of a CpG island methylator phenotype that defines a distinct subgroup of gliomaCancer Cell17510522201020399149ChristensenBCDNA methylation, isocitrate dehydrogenase mutation and survival in gliomaJ. Natl. Cancer Inst103143153201121163902SansonMIsocitrate dehydrogenase 1 codon 132 mutation is an important prognostic biomarker in gliomasJ. Clin. Oncol2741504154200919636000Eckel-PassowJEGlioma groups based on 1p/19q, IDH and TERT promoter mutations in tumorsN. Engl. J. Med37224992508201526061753WalshKMTelomere maintenance and the etiology of adult gliomaNeuro-oncol1714451452201526014050MiyakeYRPA-like mammalian Ctc1–Stn1–Ten1 complex binds to single-stranded DNA and protects telomeres independently of the Pot1 pathwayMol. Cell36193206200919854130ChenLYRedonSLingnerJThe human CST complex is a terminator of telomerase activityNature488540544201222763445BainbridgeMNGermline mutations in shelterin complex genes are associated with familial gliomaJ. Natl. Cancer Inst107384201425482530ZhangCGenetic determinants of telomere length and risk of common cancers: a mendelian randomization studyHum. Mol. Genet2453565366201526138067WalshKMLonger genotypically estimated leukocyte telomere length is associated with increased adult glioma riskOncotarget64246842477201526646793CeccarelliMMolecular profiling reveals biologically discrete subsets and pathways of progression in diffuse gliomaCell164550563201626824661XipellEEndoplasmic-reticulum-stress-inducing drugs sensitize glioma cells to temozolomide through downregulation of MGMT, MPG and Rad51Neuro-oncol1811091119201626951384NicolasCSThe role of JAK–STAT signaling within the CNSJAK-STAT2e22925201324058789Cancer Genome Atlas Research NetworkComprehensive, integrative genomic analysis of diffuse lower-grade gliomasN. Engl. J. Med37224812498201526061751HepworthSJMobile phone use and risk of glioma in adults: case–control studyBMJ332883887200616428250AmirianESThe Glioma International Case–Control Study: a report from the Genetic Epidemiology of Glioma International ConsortiumAm. J. Epidemiol1838591201626656478LiYFastPop: a rapid principal component–derived method to infer intercontinental ancestry using genetic dataBMC Bioinformatics17122201626961892CardisEThe INTERPHONE study: design, epidemiological methods and description of the study populationEur. J. Epidemiol22647664200717636416HercbergSThe SU.VI.MAX study: a randomized, placebo-controlled trial of the health effects of antioxidant vitamins and mineralsArch. Intern. Med16423352342200415557412WichmannHEGiegerCIlligTMONICA–KORA Study GroupKORA-gen—resource for population genetics, controls and a broad spectrum of disease phenotypesGesundheitswesen67Suppl. 1S26S30200516032514KrawczakMPopGen: population-based recruitment of patients and controls for the analysis of complex genotype–phenotype relationshipsCommunity Genet95561200616490960SchmermundAAssessment of clinically silent atherosclerotic disease, and established and novel risk factors for predicting myocardial infarction and cardiac death in healthy middle-aged subjects: rationale and design of the Heinz Nixdorf RECALL study. Risk factors, evaluation of coronary calcium and lifestyleAm. Heart J144212218200212177636HunterDJA genome-wide association study identifies alleles in FGFR2 associated with risk of sporadic postmenopausal breast cancerNat. Genet39870874200717529973YeagerMGenome-wide association study of prostate cancer identifies a second risk locus at 8q24Nat. Genet39645649200717401363WiemelsJLHistory of allergies among adults with glioma and controlsInt. J. Cancer98609615200211920623FeliniMJReproductive factors and hormone use risk of adult gliomasCancer Causes Control208796200918766447HowieBNDonnellyPMarchiniJA flexible and accurate genotype imputation method for the next generation of genome-wide association studiesPLoS Genet5e10005292009195433731000 Genomes Project ConsortiumA map of human genome variation from population-scale sequencingNature46710611073201020981092MarchiniJHowieBMyersSMcVeanGDonnellyPA new multipoint method for genome-wide association studies by imputation of genotypesNat. Genet39906913200717572673PurcellSPLINK: a tool set for whole-genome association and population-based linkage analysesAm. J. Hum. Genet81559575200717701901PattersonNPriceALReichDPopulation structure and eigenanalysisPLoS Genet2e190200617194218JenkinsRBDistinct germline polymorphisms underlie glioma morphologic heterogeneityCancer Genet2041318201121356187PritchardJKStephensMDonnellyPInference of population structure using multilocus genotype dataGenetics155945959200010835412LiuJZMeta-analysis and imputation refines the association of 15q25 with smoking quantityNat. Genet42436440201020418889HigginsJPThompsonSGDeeksJJAltmanDGMeasuring inconsistency in meta-analysesBr. Med. J327557560200312958120YangJLeeSHGoddardMEVisscherPMGCTA: a tool for genome-wide complex-trait analysisAm. J. Hum. Genet887682201121167468YangJConditional and joint multiple-SNP analysis of GWAS summary statistics identifies additional variants influencing complex traitsNat. Genet44369375201222426310WardLDKellisMHaploReg: a resource for exploring chromatin states, conservation and regulatory motif alterations within sets of genetically linked variantsNucleic Acids Res40D930D934201222064851BoyleAPAnnotation of functional variation in personal genomes using RegulomeDBGenome Res2217901797201222955989NgSBTargeted capture and massively parallel sequencing of 12 human exomesNature461272276200919684571Roadmap Epigenomics ConsortiumIntegrative analysis of 111 reference human epigenomesNature518317330201525693563PharoahPDAntoniouACEastonDFPonderBAPolygenes, risk prediction and targeted prevention of breast cancerN. Engl. J. Med35827962803200818579814ScheurerMEFamilial aggregation of glioma: a pooled analysisAm. J. Epidemiol17210991107201020858744DanecekPThe variant call format and VCFtoolsBioinformatics2721562158201121653522ScalesMJagerRMiglioriniGHoulstonRSHenrionMYvisPIG—a web tool for producing multiregion, multitrack, multiscale plots of genetic dataPLoS One9e107497201425208325GabrielSBThe structure of haplotype blocks in the human genomeScience29622252229200212029063

Genome-wide discovery-phase meta-analysis P-values (−log10P) plotted against their chromosomal positions. (a) All glioma. (b) GBM. (c) Non-GBM tumors. The red horizontal line corresponds to a significance threshold of P = 5.0 × 10−8. New and known loci are labeled in red and blue, respectively.

Relative impact of SNP associations at known and newly identified risk loci for GBM and non-GBM tumors. Odds ratios (ORs) derived with respect to the risk allele. Asterisks denote SNPs showing a significant difference between GBM and non-GBM tumors from the case-only analysis as detailed in Supplementary Table 4.

Association statistics for the top SNP at each of the newly-reported glioma risk loci

All glioma
GBM glioma
Non-GBM glioma
LocusSubtypeSNPPositionAllelesRAFINFOPOR (95% CI)POR (95% CI)POR (95% CI)
1p31.3GBMrs1275255265229299T/C0.8700.9924.07 × 10−91.18 (1.11–1.24)2.04 × 10−91.22 (1.15–1.31)4.78 × 10−31.11 (1.03–1.18)
1q32.1Non-GBMrs4252707204508147G/A0.2200.9922.97 × 10−71.12 (1.07–1.17)0.0151.07 (1.01–1.13)3.34 × 10−91.19 (1.12–1.26)
1q44Non-GBMrs12076373243851947G/C0.8370.9964.97 × 10−41.09 (1.04–1.15)0.8460.99 (0.94–1.06)2.63 × 10−101.23 (1.16–1.32)
2q33.3Non-GBMrs7572263209051586A/G0.7560.9972.58 × 10−61.11 (1.06–1.15)0.0191.06 (1.01–1.12)2.18 × 10−101.20 (1.13–1.26)
3p14.1Non-GBMrs1170683266502981A/C0.4560.9971.06 × 10−51.08 (1.05–1.12)0.1581.03 (0.99–1.08)7.66 × 10−91.15 (1.09–1.20)
10q24.33Non-GBMrs11598018105661315C/A0.4620.9603.07 × 10−71.10 (1.06–1.14)0.01031.06 (1.01–1.11)3.39 × 10−81.14 (1.09–1.20)
11q14.1GBMrs1123325082397014C/T0.8680.9905.40 × 10−61.14 (1.08–1.21)9.95 × 10−101.24 (1.16–1.33)0.5920.98 (0.91–1.05)
11q21Non-GBMrs710778595747337T/C0.4790.9972.96 × 10−41.07 (1.03–1.11)0.8441.00 (0.95–1.04)3.87 × 10−101.16 (1.11–1.21)
14q12Non-GBMrs1013103233250081G/A0.9160.9912.33 × 10−61.17 (1.09–1.24)0.2471.05 (0.97–1.13)5.07 × 10−111.33 (1.22–1.44)
16p13.3GBMrs2562152123896A/T0.8500.9371.18 × 10−31.09 (1.04–1.15)1.93 × 10−81.21 (1.13–1.29)0.9481.00 (0.93–1.07)
16p13.3Non-GBMrs37516671004554C/T0.2080.9858.75 × 10−101.14 (1.09–1.19)5.95 × 10−61.13 (1.07–1.19)2.61 × 10−91.18 (1.12–1.25)
16q12.1GBMrs1085260650128872T/C0.7130.9903.66 × 10−111.14 (1.10–1.19)1.29 × 10−111.18 (1.13–1.24)2.42 × 10−31.08 (1.03–1.14)
22q13.1GBMrs223557338477930G/A0.5070.9958.64 × 10−71.09 (1.06–1.13)1.76 × 10−101.15 (1.10–1.20)0.3251.02 (0.97–1.07)

Associations at P < 5 × 10−8 are highlighted in bold. Odds ratios (ORs) were derived with respect to the risk allele underlined and highlighted in bold. Risk allele frequency (RAF) is according to European samples from the 1000 Genomes project. The INFO column indicates the average imputation info score across all studies, with a score of 1 indicating that the SNP is directly genotyped in all studies. CI, confidence interval.