The genomic landscape of 2,023 colorectal cancers.

Colorectal carcinoma (CRC) is a common cause of mortality1, but a comprehensive description of its genomic landscape is lacking2-9. Here we perform whole-genome sequencing of 2,023 CRC samples from participants in the UK 100,000 Genomes Project, thereby providing a highly detailed somatic mutational landscape of this cancer. Integrated analyses identify more than 250 putative CRC driver genes, many not previously implicated in CRC or other cancers, including several recurrent changes outside the coding genome. We extend the molecular pathways involved in CRC development, define four new common subgroups of microsatellite-stable CRC based on genomic features and show that these groups have independent prognostic associations. We also characterize several rare molecular CRC subgroups, some with potential clinical relevance, including cancers with both microsatellite and chromosomal instability. We demonstrate a spectrum of mutational profiles across the colorectum, which reflect aetiological differences. These include the role of Escherichia colipks+ colibactin in rectal cancers10 and the importance of the SBS93 signature11-13, which suggests that diet or smoking is a risk factor. Immune-escape driver mutations14 are near-ubiquitous in hypermutant tumours and occur in about half of microsatellite-stable CRCs, often in the form of HLA copy number changes. Many driver mutations are actionable, including those associated with rare subgroups (for example, BRCA1 and IDH1), highlighting the role of whole-genome sequencing in optimizing patient care.

Germline MBD4 deficiency causes a multi-tumor predisposition syndrome.

We report an autosomal recessive, multi-organ tumor predisposition syndrome, caused by bi-allelic loss-of-function germline variants in the base excision repair (BER) gene MBD4. We identified five individuals with bi-allelic MBD4 variants within four families and these individuals had a personal and/or family history of adenomatous colorectal polyposis, acute myeloid leukemia, and uveal melanoma. MBD4 encodes a glycosylase involved in repair of G:T mismatches resulting from deamination of 5'-methylcytosine. The colorectal adenomas from MBD4-deficient individuals showed a mutator phenotype attributable to mutational signature SBS1, consistent with the function of MBD4. MBD4-deficient polyps harbored somatic mutations in similar driver genes to sporadic colorectal tumors, although AMER1 mutations were more common and KRAS mutations less frequent. Our findings expand the role of BER deficiencies in tumor predisposition. Inclusion of MBD4 in genetic testing for polyposis and multi-tumor phenotypes is warranted to improve disease management.

Genetic predisposition to mosaic Y chromosome loss in blood.

Mosaic loss of chromosome Y (LOY) in circulating white blood cells is the most common form of clonal mosaicism1-5, yet our knowledge of the causes and consequences of this is limited. Here, using a computational approach, we estimate that 20% of the male population represented in the UK Biobank study (n = 205,011) has detectable LOY. We identify 156 autosomal genetic determinants of LOY, which we replicate in 757,114 men of European and Japanese ancestry. These loci highlight genes that are involved in cell-cycle regulation and cancer susceptibility, as well as somatic drivers of tumour growth and targets of cancer therapy. We demonstrate that genetic susceptibility to LOY is associated with non-haematological effects on health in both men and women, which supports the hypothesis that clonal haematopoiesis is a biomarker of genomic instability in other tissues. Single-cell RNA sequencing identifies dysregulated expression of autosomal genes in leukocytes with LOY and provides insights into why clonal expansion of these cells may occur. Collectively, these data highlight the value of studying clonal mosaicism to uncover fundamental mechanisms that underlie cancer and other ageing-related diseases.

Searching for causal relationships of glioma: a phenome-wide Mendelian randomisation study.

BACKGROUND: The aetiology of glioma is poorly understood. Summary data from genome-wide association studies (GWAS) can be used in a Mendelian randomisation (MR) phenome-wide association study (PheWAS) to search for glioma risk factors. METHODS: We performed an MR-PheWAS analysing 316 phenotypes, proxied by 8387 genetic variants, and summary genetic data from a GWAS of 12,488 glioma cases and 18,169 controls. Causal effects were estimated under a random-effects inverse-variance-weighted (IVW-RE) model, with robust adjusted profile score (MR-RAPS), weighted median and mode-based estimates computed to assess the robustness of findings. Odds ratios per one standard deviation increase in each phenotype were calculated for all glioma, glioblastoma (GBM) and non-GBM tumours. RESULTS: No significant associations (P < 1.58 × 10-4) were observed between phenotypes and glioma under the IVW-RE model. Suggestive associations (1.58 × 10-4 < P < 0.05) were observed between leukocyte telomere length (LTL) with all glioma (ORSD = 3.91, P = 9.24 × 10-3) and GBM (ORSD = 4.86, P = 3.23 × 10-2), but the association was primarily driven by the TERT variant rs2736100. Serum low-density lipoprotein cholesterol and plasma HbA1C showed suggestive associations with glioma (ORSD = 1.11, P = 1.39 × 10-2 and ORSD = 1.28, P = 1.73 × 10-2, respectively), both associations being reliant on single genetic variants. CONCLUSIONS: Our study provides further insight into the aetiological basis of glioma for which published data have been mixed.

Transcriptome-wide association study of multiple myeloma identifies candidate susceptibility genes.

BACKGROUND: While genome-wide association studies (GWAS) of multiple myeloma (MM) have identified variants at 23 regions influencing risk, the genes underlying these associations are largely unknown. To identify candidate causal genes at these regions and search for novel risk regions, we performed a multi-tissue transcriptome-wide association study (TWAS). RESULTS: GWAS data on 7319 MM cases and 234,385 controls was integrated with Genotype-Tissue Expression Project (GTEx) data assayed in 48 tissues (sample sizes, N = 80-491), including lymphocyte cell lines and whole blood, to predict gene expression. We identified 108 genes at 13 independent regions associated with MM risk, all of which were in 1 Mb of known MM GWAS risk variants. Of these, 94 genes, located in eight regions, had not previously been considered as a candidate gene for that locus. CONCLUSIONS: Our findings highlight the value of leveraging expression data from multiple tissues to identify candidate genes responsible for GWAS associations which provide insight into MM tumorigenesis. Among the genes identified, a number have plausible roles in MM biology, notably APOBEC3C, APOBEC3H, APOBEC3D, APOBEC3F, APOBEC3G, or have been previously implicated in other malignancies. The genes identified in this TWAS can be explored for follow-up and validation to further understand their role in MM biology.

Age-specific genome-wide association study in glioblastoma identifies increased proportion of 'lower grade glioma'-like features associated with younger age.

Glioblastoma (GBM) is the most common malignant brain tumor in the United States. Incidence of GBM increases with age, and younger age-at-diagnosis is significantly associated with improved prognosis. While the relationship between candidate GBM risk SNPs and age-at-diagnosis has been explored, genome-wide association studies (GWAS) have not previously been stratified by age. Potential age-specific genetic effects were assessed in autosomal SNPs for GBM patients using data from four previous GWAS. Using age distribution tertiles (18-53, 54-64, 65+) datasets were analyzed using age-stratified logistic regression to generate p values, odds ratios (OR), and 95% confidence intervals (95%CI), and then combined using meta-analysis. There were 4,512 total GBM cases, and 10,582 controls used for analysis. Significant associations were detected at two previously identified SNPs in 7p11.2 (rs723527 [p54-63 = 1.50x10-9 , OR54-63 = 1.28, 95%CI54-63 = 1.18-1.39; p64+ = 2.14x10-11 , OR64+ = 1.32, 95%CI64+ = 1.21-1.43] and rs11979158 [p54-63 = 6.13x10-8 , OR54-63 = 1.35, 95%CI54-63 = 1.21-1.50; p64+ = 2.18x10-10 , OR64+ = 1.42, 95%CI64+ = 1.27-1.58]) but only in persons >54. There was also a significant association at the previously identified lower grade glioma (LGG) risk locus at 8q24.21 (rs55705857) in persons ages 18-53 (p18-53 = 9.30 × 10-11 , OR18-53 = 1.76, 95%CI18-53 = 1.49-2.10). Within The Cancer Genome Atlas (TCGA) there was higher prevalence of 'LGG'-like tumor characteristics in GBM samples in those 18-53, with IDH1/2 mutation frequency of 15%, as compared to 2.1% [54-63] and 0.8% [64+] (p = 0.0005). Age-specific differences in cancer susceptibility can provide important clues to etiology. The association of a SNP known to confer risk for IDH1/2 mutant glioma and higher prevalence of IDH1/2 mutation within younger individuals 18-53 suggests that more younger individuals may present initially with 'secondary glioblastoma.'

Influence of obesity-related risk factors in the aetiology of glioma.

BACKGROUND: Obesity and related factors have been implicated as possible aetiological factors for the development of glioma in epidemiological observation studies. We used genetic markers in a Mendelian randomisation framework to examine whether obesity-related traits influence glioma risk. This methodology reduces bias from confounding and is not affected by reverse causation. METHODS: Genetic instruments were identified for 10 key obesity-related risk factors, and their association with glioma risk was evaluated using data from a genome-wide association study of 12,488 glioma patients and 18,169 controls. The estimated odds ratio of glioma associated with each of the genetically defined obesity-related traits was used to infer evidence for a causal relationship. RESULTS: No convincing association with glioma risk was seen for genetic instruments for body mass index, waist-to-hip ratio, lipids, type-2 diabetes, hyperglycaemia or insulin resistance. Similarly, we found no evidence to support a relationship between obesity-related traits with subtypes of glioma-glioblastoma (GBM) or non-GBM tumours. CONCLUSIONS: This study provides no evidence to implicate obesity-related factors as causes of glioma.

Validation of Recently Proposed Colorectal Cancer Susceptibility Gene Variants in an Analysis of Families and Patients-a Systematic Review.

High-throughput sequencing analysis has accelerated searches for genes associated with risk for colorectal cancer (CRC); germline mutations in NTHL1, RPS20, FANCM, FAN1, TP53, BUB1, BUB3, LRP6, and PTPN12 have been recently proposed to increase CRC risk. We attempted to validate the association between variants in these genes and development of CRC in a systematic review of 11 publications, using sequence data from 863 familial CRC cases and 1604 individuals without CRC (controls). All cases were diagnosed at an age of 55 years or younger and did not carry mutations in an established CRC predisposition gene. We found sufficient evidence for NTHL1 to be considered a CRC predisposition gene-members of 3 unrelated Dutch families were homozygous for inactivating p.Gln90Ter mutations; a Canadian woman with polyposis, CRC, and multiple tumors was reported to be heterozygous for the inactivating NTHL1 p.Gln90Ter/c.709+1G>A mutations; and a man with polyposis was reported to carry p.Gln90Ter/p.Gln287Ter; whereas no inactivating homozygous or compound heterozygous mutations were detected in controls. Variants that disrupted RPS20 were detected in a Finnish family with early-onset CRC (p.Val50SerfsTer23), a 39-year old individual with metachronous CRC (p.Leu61GlufsTer11 mutation), and a 41-year-old individual with CRC (missense p.Val54Leu), but not in controls. We therefore found published evidence to support the association between variants in NTHL1 and RPS20 with CRC, but not of other recently reported CRC susceptibility variants. We urge the research community to adopt rigorous statistical and biological approaches coupled with independent replication before making claims of pathogenicity.

Impact of atopy on risk of glioma: a Mendelian randomisation study.

BACKGROUND: An inverse relationship between allergies with glioma risk has been reported in several but not all epidemiological observational studies. We performed an analysis of genetic variants associated with atopy to assess the relationship with glioma risk using Mendelian randomisation (MR), an approach unaffected by biases from temporal variability and reverse causation that might have affected earlier investigations. METHODS: Two-sample MR was undertaken using genome-wide association study data. We used single nucleotide polymorphisms (SNPs) associated with atopic dermatitis, asthma and hay fever, IgE levels, and self-reported allergy as instrumental variables. We calculated MR estimates for the odds ratio (OR) for each risk factor with glioma using SNP-glioma estimates from 12,488 cases and 18,169 controls, using inverse-variance weighting (IVW), maximum likelihood estimation (MLE), weighted median estimate (WME) and mode-based estimate (MBE) methods. Violation of MR assumptions due to directional pleiotropy were sought using MR-Egger regression and HEIDI-outlier analysis. RESULTS: Under IVW, MLE, WME and MBE methods, associations between glioma risk with asthma and hay fever, self-reported allergy and IgE levels were non-significant. An inverse relationship between atopic dermatitis and glioma risk was found by IVW (OR 0.96, 95% confidence interval (CI) 0.93-1.00, P = 0.041) and MLE (OR 0.96, 95% CI 0.94-0.99, P = 0.003), but not by WME (OR 0.96, 95% CI 0.91-1.01, P = 0.114) or MBE (OR 0.97, 95% CI 0.92-1.02, P = 0.194). CONCLUSIONS: Our investigation does not provide strong evidence for relationship between atopy and the risk of developing glioma, but findings do not preclude a small effect in relation to atopic dermatitis. Our analysis also serves to illustrate the value of using several MR methods to derive robust conclusions.

Germ line mutations in shelterin complex genes are associated with familial chronic lymphocytic leukemia.

Chronic lymphocytic leukemia (CLL) can be familial; however, thus far no rare germ line disruptive alleles for CLL have been identified. We performed whole-exome sequencing of 66 CLL families, identifying 4 families where loss-of-function mutations in protection of telomeres 1 (POT1) co-segregated with CLL. The p.Tyr36Cys mutation is predicted to disrupt the interaction between POT1 and the telomeric overhang. The c.1164-1G>A splice-site, p.Gln358SerfsTer13 frameshift, and p.Gln376Arg missense mutations are likely to impact the interaction between POT1 and adrenocortical dysplasia homolog (ACD), which is a part of the telomere-capping shelterin complex. We also identified mutations in ACD (c.752-2A>C) and another shelterin component, telomeric repeat binding factor 2, interacting protein (p.Ala104Pro and p.Arg133Gln), in 3 CLL families. In a complementary analysis of 1083 cases and 5854 controls, the POT1 p.Gln376Arg variant, which has a global minor allele frequency of 0.0005, conferred a 3.61-fold increased risk of CLL (P = .009). This study further highlights telomere dysregulation as a key process in CLL development.

Diffuse gliomas classified by 1p/19q co-deletion, TERT promoter and IDH mutation status are associated with specific genetic risk loci.

Recent genome-wide association studies of glioma have led to the discovery of single nucleotide polymorphisms (SNPs) at 25 loci influencing risk. Gliomas are heterogeneous, hence to investigate the relationship between risk SNPs and glioma subtype we analysed 1659 tumours profiled for IDH mutation, TERT promoter mutation and 1p/19q co-deletion. These data allowed definition of five molecular subgroups of glioma: triple-positive (IDH mutated, 1p/19q co-deletion, TERT promoter mutated); TERT-IDH (IDH mutated, TERT promoter mutated, 1p/19q-wild-type); IDH-only (IDH mutated, 1p/19q wild-type, TERT promoter wild-type); triple-negative (IDH wild-type, 1p/19q wild-type, TERT promoter wild-type) and TERT-only (TERT promoter mutated, IDH wild-type, 1p/19q wild-type). Most glioma risk loci showed subtype specificity: (1) the 8q24.21 SNP for triple-positive glioma; (2) 5p15.33, 9p21.3, 17p13.1 and 20q13.33 SNPs for TERT-only glioma; (3) 1q44, 2q33.3, 3p14.1, 11q21, 11q23.3, 14q12, and 15q24.2 SNPs for IDH mutated glioma. To link risk SNPs to target candidate genes we analysed Hi-C and gene expression data, highlighting the potential role of IDH1 at 2q33.3, MYC at 8q24.21 and STMN3 at 20q13.33. Our observations provide further insight into the nature of susceptibility to glioma.

Identification of susceptibility loci for colorectal cancer in a genome-wide meta-analysis.

To identify common variants influencing colorectal cancer (CRC) risk, we performed a meta-analysis of five genome-wide association studies, comprising 5626 cases and 7817 controls of European descent. We conducted replication of top ranked single nucleotide polymorphisms (SNPs) in additional series totalling 14 037 cases and 15 937 controls, identifying a new CRC risk locus at 10q24.2 [rs1035209; odds ratio (OR) = 1.13, P = 4.54 × 10(-11)]. We also performed meta-analysis of our studies, with previously published data, of several recently purported CRC risk loci. We failed to find convincing evidence for a previously reported genome-wide association at rs11903757 (2q32.3). Of the three additional loci for which evidence of an association in Europeans has been previously described we failed to show an association between rs59336 (12q24.21) and CRC risk. However, for the other two SNPs, our analyses demonstrated new, formally significant associations with CRC. These are rs3217810 intronic in CCND2 (12p13.32; OR = 1.19, P = 2.16 × 10(-10)) and rs10911251 near LAMC1 (1q25.3; OR = 1.09, P = 1.75 × 10(-8)). Additionally, we found some evidence to support a relationship between, rs647161, rs2423297 and rs10774214 and CRC risk originally identified in East Asians in our European datasets. Our findings provide further insights into the genetic and biological basis of inherited genetic susceptibility to CRC.

Deciphering the 8q24.21 association for glioma.

We have previously identified tagSNPs at 8q24.21 influencing glioma risk. We have sought to fine-map the location of the functional basis of this association using data from four genome-wide association studies, comprising a total of 4147 glioma cases and 7435 controls. To improve marker density across the 700 kb region, we imputed genotypes using 1000 Genomes Project data and high-coverage sequencing data generated on 253 individuals. Analysis revealed an imputed low-frequency SNP rs55705857 (P = 2.24 × 10(-38)) which was sufficient to fully capture the 8q24.21 association. Analysis by glioma subtype showed the association with rs55705857 confined to non-glioblastoma multiforme (non-GBM) tumours (P = 1.07 × 10(-67)). Validation of the non-GBM association was shown in three additional datasets (625 non-GBM cases, 2412 controls; P = 1.41 × 10(-28)). In the pooled analysis, the odds ratio for low-grade glioma associated with rs55705857 was 4.3 (P = 2.31 × 10(-94)). rs55705857 maps to a highly evolutionarily conserved sequence within the long non-coding RNA CCDC26 raising the possibility of direct functionality. These data provide additional insights into the aetiological basis of glioma development.

Radiogenomic biomarkers for immunotherapy in glioblastoma: A systematic review of magnetic resonance imaging studies.

Background: Immunotherapy is an effective "precision medicine" treatment for several cancers. Imaging signatures of the underlying genome (radiogenomics) in glioblastoma patients may serve as preoperative biomarkers of the tumor-host immune apparatus. Validated biomarkers would have the potential to stratify patients during immunotherapy clinical trials, and if trials are beneficial, facilitate personalized neo-adjuvant treatment. The increased use of whole genome sequencing data, and the advances in bioinformatics and machine learning make such developments plausible. We performed a systematic review to determine the extent of development and validation of immune-related radiogenomic biomarkers for glioblastoma. Methods: A systematic review was performed following PRISMA guidelines using the PubMed, Medline, and Embase databases. Qualitative analysis was performed by incorporating the QUADAS 2 tool and CLAIM checklist. PROSPERO registered: CRD42022340968. Extracted data were insufficiently homogenous to perform a meta-analysis. Results: Nine studies, all retrospective, were included. Biomarkers extracted from magnetic resonance imaging volumes of interest included apparent diffusion coefficient values, relative cerebral blood volume values, and image-derived features. These biomarkers correlated with genomic markers from tumor cells or immune cells or with patient survival. The majority of studies had a high risk of bias and applicability concerns regarding the index test performed. Conclusions: Radiogenomic immune biomarkers have the potential to provide early treatment options to patients with glioblastoma. Targeted immunotherapy, stratified by these biomarkers, has the potential to allow individualized neo-adjuvant precision treatment options in clinical trials. However, there are no prospective studies validating these biomarkers, and interpretation is limited due to study bias with little evidence of generalizability.

Genetic landscape of interval and screen detected breast cancer.

Interval breast cancers (IBCs) are cancers diagnosed between screening episodes. Understanding the biological differences between IBCs and screen-detected breast-cancers (SDBCs) has the potential to improve mammographic screening and patient management. We analysed and compared the genomic landscape of 288 IBCs and 473 SDBCs by whole genome sequencing of paired tumour-normal patient samples collected as part of the UK 100,000 Genomes Project. Compared to SDBCs, IBCs were more likely to be lobular, higher grade, and triple negative. A more aggressive clinical phenotype was reflected in IBCs displaying features of genomic instability including a higher mutation rate and number of chromosomal structural abnormalities, defective homologous recombination and TP53 mutations. We did not however, find evidence to indicate that IBCs are associated with a significantly different immune response. While IBCs do not represent a unique molecular class of invasive breast cancer they exhibit a more aggressive phenotype, which is likely to be a consequence of the timing of tumour initiation. This information is relevant both with respect to treatment as well as informing the screening interval for mammography.

Analysis of 10,478 cancer genomes identifies candidate driver genes and opportunities for precision oncology.

Tumor genomic profiling is increasingly seen as a prerequisite to guide the treatment of patients with cancer. To explore the value of whole-genome sequencing (WGS) in broadening the scope of cancers potentially amenable to a precision therapy, we analysed whole-genome sequencing data on 10,478 patients spanning 35 cancer types recruited to the UK 100,000 Genomes Project. We identified 330 candidate driver genes, including 74 that are new to any cancer. We estimate that approximately 55% of patients studied harbor at least one clinically relevant mutation, predicting either sensitivity or resistance to certain treatments or clinical trial eligibility. By performing computational chemogenomic analysis of cancer mutations we identify additional targets for compounds that represent attractive candidates for future clinical trials. This study represents one of the most comprehensive efforts thus far to identify cancer driver genes in the real world setting and assess their impact on informing precision oncology.

Whole genome sequencing refines stratification and therapy of patients with clear cell renal cell carcinoma.

Clear cell renal cell carcinoma (ccRCC) is the most common form of kidney cancer, but a comprehensive description of its genomic landscape is lacking. We report the whole genome sequencing of 778 ccRCC patients enrolled in the 100,000 Genomes Project, providing for a detailed description of the somatic mutational landscape of ccRCC. We identify candidate driver genes, which as well as emphasising the major role of epigenetic regulation in ccRCC highlight additional biological pathways extending opportunities for therapeutic interventions. Genomic characterisation identified patients with divergent clinical outcome; higher number of structural copy number alterations associated with poorer prognosis, whereas VHL mutations were independently associated with a better prognosis. The observations that higher T-cell infiltration is associated with better overall survival and that genetically predicted immune evasion is not common supports the rationale for immunotherapy. These findings should inform personalised surveillance and treatment strategies for ccRCC patients.

Therapeutic Targeting of Glioblastoma and the Interactions with Its Microenvironment.

Glioblastoma (GBM) is the most common primary malignant brain tumour, and it confers a dismal prognosis despite intensive multimodal treatments. Whilst historically, research has focussed on the evolution of GBM tumour cells themselves, there is growing recognition of the importance of studying the tumour microenvironment (TME). Improved characterisation of the interaction between GBM cells and the TME has led to a better understanding of therapeutic resistance and the identification of potential targets to block these escape mechanisms. This review describes the network of cells within the TME and proposes treatment strategies for simultaneously targeting GBM cells, the surrounding immune cells, and the crosstalk between them.

Whole genome sequencing refines stratification and therapy of patients with clear cell renal cell carcinoma.

Clear cell renal cell carcinoma (ccRCC) is the most common form of kidney cancer, but a comprehensive description of its genomic landscape is lacking. We report the whole genome sequencing of 778 ccRCC patients enrolled in the 100,000 Genomes Project, providing the most detailed somatic mutational landscape to date. We identify new driver genes, which as well as emphasising the major role of epigenetic regulation in ccRCC highlight additional biological pathways extending opportunities for drug repurposing. Genomic characterisation identified patients with divergent clinical outcome; higher number of structural copy number alterations associated with poorer prognosis, whereas VHL mutations were independently associated with a better prognosis. The twin observations that higher T-cell infiltration is associated with better outcome and that genetically predicted immune evasion is not common supports the rationale for immunotherapy. These findings should inform personalised surveillance and treatment strategies for ccRCC patients.

Relationship between genetically determined telomere length and glioma risk.

BACKGROUND: Telomere maintenance is increasingly recognized as being fundamental to glioma oncogenesis with longer leukocyte telomere length (LTL) reported to increase risk of glioma. To gain further insight into the relationship between telomere genetics and risk of glioma, we conducted several complementary analyses, using genome-wide association studies data on LTL (78 592 individuals) and glioma (12 488 cases and 18 169 controls). METHODS: We performed both classical and summary Mendelian randomization (SMR), coupled with heterogeneity in dependent instruments tests, at genome-wide significant LTL loci to examine if an association was mediated by the same causal variant in glioma. To prioritize genes underscoring glioma-LTL associations, we analyzed gene expression and DNA methylation data. RESULTS: Genetically increased LTL was significantly associated with increased glioma risk, random-effects inverse variance weighted ORs per 1 SD unit increase in the putative risk factor (odds ratio [OR]SD) 4.79 (95% confidence interval: 2.11-10.85; P = 1.76 × 10-4). SMR confirmed the previously reported LTL associations at 3q26.2 (TERC; PSMR = 1.33 × 10-5), 5p15.33 (TERT; PSMR = 9.80 × 10-27), 10q24.33 (STN1 alias OBFC1; PSMR = 4.31 × 10-5), and 20q13.3 (STMN3/RTEL1; PSMR = 2.47 × 10-4) glioma risk loci. Our analysis implicates variation at 1q42.12 (PSMR = 1.55 × 10-2), 6p21.3 (PSMR = 9.76 × 10-3), 6p22.2 (PSMR = 5.45 × 10-3), 7q31.33 (PSMR = 6.52 × 10-3), and 11q22.3 (PSMR = 8.89 × 10-4) as risk factors for glioma risk. While complicated by patterns of linkage disequilibrium, genetic variation involving PARP1, PRRC2A, CARMIL1, POT1, and ATM-NPAT1 was implicated in the etiology of glioma. CONCLUSIONS: These observations extend the role of telomere-related genes in the development of glioma.