Adult diffuse glioma GWAS by molecular subtype identifies variants in D2HGDH and FAM20C

Functional germline variants in DNA damage repair pathways are associated with altered survival in adults with glioma treated with temozolomide

Background

Temozolomide (TMZ) treatment has demonstrated, but variable, impact on glioma prognosis. This study examines associations of survival with DNA repair gene germline polymorphisms among glioma patients who did and did not have TMZ treatment. Identifying genetic markers which sensitize tumor cells to TMZ could personalize therapy and improve outcomes.

Methods

We evaluated TMZ-related survival associations of pathogenic germline SNPs and genetically predicted transcript levels within 34 DNA repair genes among 1504 glioma patients from the UCSF Adult Glioma Study and Mayo Clinic whose diagnoses spanned pre- and post-TMZ eras within the major known glioma prognostic molecular subtypes.

Results

Among those who received TMZ, 5 SNPs were associated with overall survival, but not in those who did not receive TMZ. Only rs2308321-G, in MGMT, was associated with decreased survival (HR=1.21, p=0.019) for all glioma subtypes. Rs73191162-T (near UNG), rs13076508-C (near PARP3), rs7840433-A (near NEIL2), and rs3130618-A (near MSH5) were only associated with survival and TMZ treatment for certain subtypes, suggesting subtype-specific germline chemo-sensitization.Genetically predicted elevated compared to normal brain expression of PNKP was associated with dramatically worse survival for TMZ-treated patients with IDH-mutant and 1p/19q non-codeleted gliomas (p=0.015). Similarly, NEIL2 and TDG expressions were associated with altered TMZ-related survival only among certain subtypes.

Conclusions

Functional germline alterations within DNA repair genes were associated with TMZ sensitivity, measured by overall survival, among adults with glioma, these variants should be evaluated in prospective analyses and functional studies.

Genome-wide polygenic risk scores predict risk of glioma and molecular subtypes

BACKGROUND

Polygenic risk scores (PRS) aggregate the contribution of many risk variants to provide a personalized genetic susceptibility profile. Since sample sizes of glioma genome-wide association studies (GWAS) remain modest, there is a need to efficiently capture genetic risk using available data.

METHODS

We applied a method based on continuous shrinkage priors (PRS-CS) to model the joint effects of over 1 million common variants on disease risk and compared this to an approach (PRS-CT) that only selects a limited set of independent variants that reach genome-wide significance (P < 5 × 10-8). PRS models were trained using GWAS stratified by histological (10 346 cases and 14 687 controls) and molecular subtype (2632 cases and 2445 controls), and validated in 2 independent cohorts.

RESULTS

PRS-CS was generally more predictive than PRS-CT with a median increase in explained variance (R2) of 24% (interquartile range = 11-30%) across glioma subtypes. Improvements were pronounced for glioblastoma (GBM), with PRS-CS yielding larger odds ratios (OR) per standard deviation (SD) (OR = 1.93, P = 2.0 × 10-54 vs. OR = 1.83, P = 9.4 × 10-50) and higher explained variance (R2 = 2.82% vs. R2 = 2.56%). Individuals in the 80th percentile of the PRS-CS distribution had a significantly higher risk of GBM (0.107%) at age 60 compared to those with average PRS (0.046%, P = 2.4 × 10-12). Lifetime absolute risk reached 1.18% for glioma and 0.76% for IDH wildtype tumors for individuals in the 95th PRS percentile. PRS-CS augmented the classification of IDH mutation status in cases when added to demographic factors (AUC = 0.839 vs. AUC = 0.895, PΔAUC = 6.8 × 10-9).

CONCLUSIONS

Genome-wide PRS has the potential to enhance the detection of high-risk individuals and help distinguish between prognostic glioma subtypes.

Serum exosomal miRNA promote glioma progression by targeting SOS1 via abscopal effect of radiation

INTRODUCTION

Local exposure to ionizing radiation (IR) can induce changes in biological processes in distant tissues and organs. Exosomes are nanoscale vesicles that transport biomolecules, mediate communication between cells and tissues, and can affect the abscopal effects of radiotherapy.

METHODS

Mice were treated with 8.0 Gy doses of chest and abdomen IR, after which serum samples were taken 24 h after exposure. Their serum exosomes were then isolated via ultracentrifugation and the small RNA portions were extracted for sequencing and bioinformatic analysis. Exosomes were injected intravenously into the mice to assess their ability to cross the blood-brain barrier (BBB). Glioma cells and glioma stem cells (GSCs) were examined for malignant biological behaviors, stemness, and tumorigenic capacity after co-culturing with different groups of exosomes.

RESULTS

We found that serum exosomes crossed the BBB in mice after local IR exposure-which induced decreases in the expression of BBB tight-junction proteins and increased brain endothelial cell apoptosis. Exosomes from the exposed groups promoted malignant biological behaviors, stemness, and tumorigenic capacity in glioma cells and GSCs by upregulating the expression of SOS1. Phospho-MEK1/2 and Phospho-ERK1/2, of the MAPK signaling pathway, were found to be up-regulated in cells that were co-cultured with the exposing groups of the exosomes. Further analyses demonstrated that differentially expressed levels of miR-93-5p in mouse serum exosomes regulated the cellular expression of SOS1.

CONCLUSION

Following local IR exposure, serum exosomes cross the BBB to promote the progression of distant gliomas. Exosomal microRNAs play an important role in this process.

AllergoOncology: Biomarkers and refined classification for research in the allergy and glioma nexus-A joint EAACI-EANO position paper

Epidemiological studies have explored the relationship between allergic diseases and cancer risk or prognosis in AllergoOncology. Some studies suggest an inverse association, but uncertainties remain, including in IgE-mediated diseases and glioma. Allergic disease stems from a Th2-biased immune response to allergens in predisposed atopic individuals. Allergic disorders vary in phenotype, genotype and endotype, affecting their pathophysiology. Beyond clinical manifestation and commonly used clinical markers, there is ongoing research to identify novel biomarkers for allergy diagnosis, monitoring, severity assessment and treatment. Gliomas, the most common and diverse brain tumours, have in parallel undergone changes in classification over time, with specific molecular biomarkers defining glioma subtypes. Gliomas exhibit a complex tumour-immune interphase and distinct immune microenvironment features. Immunotherapy and targeted therapy hold promise for primary brain tumour treatment, but require more specific and effective approaches. Animal studies indicate allergic airway inflammation may delay glioma progression. This collaborative European Academy of Allergy and Clinical Immunology (EAACI) and European Association of Neuro-Oncology (EANO) Position Paper summarizes recent advances and emerging biomarkers for refined allergy and adult-type diffuse glioma classification to inform future epidemiological and clinical studies. Future research is needed to enhance our understanding of immune-glioma interactions to ultimately improve patient prognosis and survival.

Genome-wide Polygenic Risk Scores Predict Risk of Glioma and Molecular Subtypes

Background

Polygenic risk scores (PRS) aggregate the contribution of many risk variants to provide a personalized genetic susceptibility profile. Since sample sizes of glioma genome-wide association studies (GWAS) remain modest, there is a need to find efficient ways of capturing genetic risk factors using available germline data.

Methods

We developed a novel PRS (PRS-CS) that uses continuous shrinkage priors to model the joint effects of over 1 million polymorphisms on disease risk and compared it to an approach (PRS-CT) that selects a limited set of independent variants that reach genome-wide significance (P<5×10-8). PRS models were trained using GWAS results stratified by histological (10,346 cases, 14,687 controls) and molecular subtype (2,632 cases, 2,445 controls), and validated in two independent cohorts.

Results

PRS-CS was consistently more predictive than PRS-CT across glioma subtypes with an average increase in explained variance (R2) of 21%. Improvements were particularly pronounced for glioblastoma tumors, with PRS-CS yielding larger effect sizes (odds ratio (OR)=1.93, P=2.0×10-54 vs. OR=1.83, P=9.4×10-50) and higher explained variance (R2=2.82% vs. R2=2.56%). Individuals in the 95th percentile of the PRS-CS distribution had a 3-fold higher lifetime absolute risk of IDH mutant (0.63%) and IDH wildtype (0.76%) glioma relative to individuals with average PRS. PRS-CS also showed high classification accuracy for IDH mutation status among cases (AUC=0.895).

Conclusions

Our novel genome-wide PRS may improve the identification of high-risk individuals and help distinguish between prognostic glioma subtypes, increasing the potential clinical utility of germline genetics in glioma patient management.

Construction and validation of a risk model based on the key SNARE proteins to predict the prognosis and immune microenvironment of gliomas

Background

Synaptic transmission between neurons and glioma cells can promote glioma progression. The soluble N-ethylmaleimide-sensitive fusion factor attachment protein receptors (SNARE) play a key role in synaptic functions. We aimed to construct and validate a novel model based on the SNARE proteins to predict the prognosis and immune microenvironment of glioma.

Methods

Differential expression analysis and COX regression analysis were used to identify key SRGs in glioma datasets, and we constructed a prognostic risk model based on the key SRGs. The prognostic value and predictive performance of the model were assessed in The Cancer Genome Atlas (TCGA) and Chinese glioma Genome Atlas (CGGA) datasets. Functional enrichment analysis and immune-related evaluation were employed to reveal the association of risk scores with tumor progression and microenvironment. A prognostic nomogram containing the risk score was established and assessed by calibration curves and time-dependent receiver operating characteristic curves. We verified the changes of the key SRGs in glioma specimens and cells by real-time quantitative PCR and Western blot analyses.

Results

Vesicle-associated membrane protein 2 (VAMP2) and vesicle-associated membrane protein 5 (VAMP5) were identified as two SRGs affecting the prognoses of glioma patients. High-risk patients characterized by higher VAMP5 and lower VAMP2 expression had a worse prognosis. Higher risk scores were associated with older age, higher tumor grades, IDH wild-type, and 1p19q non-codeletion. The SRGs risk model showed an excellent predictive performance in predicting the prognosis in TCGA and CGGA datasets. Differentially expressed genes between low- and high-risk groups were mainly enriched in the pathways related to immune infiltration, tumor metastasis, and neuronal activity. Immune score, stromal score, estimate score, tumor mutational burden, and expression of checkpoint genes were positively correlated with risk scores. The nomogram containing the risk score showed good performance in predicting the prognosis of glioma. Low VAMP2 and high VAMP5 were found in different grades of glioma specimens and cell lines.

Conclusion

We constructed and validated a novel risk model based on the expression of VAMP2 and VAMP5 by bioinformatics analysis and experimental confirmation. This model might be helpful for clinically predicting the prognosis and response to immunotherapy of glioma patients.

Genetic predisposition to altered blood cell homeostasis is associated with glioma risk and survival

Glioma is a highly fatal brain tumor comprised of molecular subtypes with distinct clinical trajectories. Observational studies have suggested that variability in immune response may play a role in glioma etiology. However, their findings have been inconsistent and susceptible to reverse causation due to treatment effects and the immunosuppressive nature of glioma. We applied genetic variants associated (p<5×10-8) with blood cell traits to a meta-analysis of 3418 glioma cases and 8156 controls. Genetically predicted increase in the platelet to lymphocyte ratio (PLR) was associated with an increased risk of glioma (odds ratio (OR)=1.25, p=0.005), especially in IDH-mutant (IDHmut OR=1.38, p=0.007) and IDHmut 1p/19q non-codeleted (IDHmut-noncodel OR=1.53, p=0.004) tumors. However, reduced glioma risk was observed for higher counts of lymphocytes (IDHmut-noncodel OR=0.70, p=0.004) and neutrophils (IDHmut OR=0.69, p=0.019; IDHmut-noncodel OR=0.60, p=0.009), which may reflect genetic predisposition to enhanced immune-surveillance. In contrast to susceptibility, there was no association with survival in IDHmut-noncodel; however, in IDHmut 1p/19q co-deleted tumors, we observed higher mortality with increasing genetically predicted counts of lymphocytes (hazard ratio (HR)=1.65, 95% CI: 1.24-2.20), neutrophils (HR=1.49, 1.13-1.97), and eosinophils (HR=1.59, 1.18-2.14). Polygenic scores for blood cell traits were also associated with tumor immune microenvironment features, with heterogeneity by IDH status observed for 17 signatures related to interferon signaling, PD-1 expression, and T-cell/Cytotoxic responses. In summary, we identified novel, immune-mediated susceptibility mechanisms for glioma with potential disease management implications.

Antibodies to varicella-zoster virus and three other herpesviruses and survival in adults with glioma

BACKGROUND

Lifetime exposure to the varicella-zoster virus (VZV) has been consistently inversely associated with glioma risk, however, the relationship of VZV with survival in adults with glioma has not been investigated. In this study, we analyzed the survival of adults with glioma in relation to their antibody measurements to 4 common herpes viral infections, including VZV, measured post-diagnosis.

METHODS

We analyzed IgG antibody measurements to VZV, cytomegalovirus (CMV), herpes simplex virus 1/2 (HSV), and Epstein-Barr virus (EBV) collected from 1378 adults with glioma diagnosed between 1991 and 2010. Blood was obtained a median of 3 months after surgery. Associations of patient IgG levels with overall survival were estimated using Cox models adjusted for age, sex, self-reported race, surgery type, dexamethasone usage at blood draw, and tumor grade. Models were stratified by recruitment series and meta-analyzed to account for time-dependent treatment effects.

RESULTS

VZV antibody seropositivity was associated with improved survival outcomes in adults with glioma (Hazard ratio, HR = 0.70, 95% Confidence Interval 0.54-0.90, P = .006). Amongst cases who were seropositive for VZV antibodies, survival was significantly improved for those above the 25th percentile of continuous reactivity measurements versus those below (HR = 0.76, 0.66-0.88, P = .0003). Antibody seropositivity to EBV was separately associated with improved survival (HR = 0.71, 0.53-0.96, P = .028). Antibody positivity to 2 other common viruses (CMV, HSV) was not associated with altered survival.

CONCLUSIONS

Low levels of VZV or EBV antibodies are associated with poorer survival outcomes for adults with glioma. Differential immune response rather than viral exposure may explain these findings.

Mutational landscape of primary spinal cord astrocytoma

Primary spinal cord astrocytoma (SCA) is a rare disease. Knowledge about the molecular profiles of SCAs mostly comes from intracranial glioma; the pattern of genetic alterations of SCAs is not well understood. Herein, we describe genome-sequencing analyses of primary SCAs, aiming to characterize the mutational landscape of primary SCAs. We utilized whole exome sequencing (WES) to analyze somatic nucleotide variants (SNVs) and copy number variants (CNVs) among 51 primary SCAs. Driver genes were searched using four algorithms. GISTIC2 was used to detect significant CNVs. Additionally, recurrently mutated pathways were also summarized. A total of 12 driver genes were identified. Of those, H3F3A (47.1%), TP53 (29.4%), NF1 (19.6%), ATRX (17.6%), and PPM1D (17.6%) were the most frequently mutated genes. Furthermore, three novel driver genes seldom reported in glioma were identified: HNRNPC, SYNE1, and RBM10. Several germline mutations, including three variants (SLC16A8 rs2235573, LMF1 rs3751667, FAM20C rs774848096) that were associated with risk of brain glioma, were frequently observed in SCAs. Moreover, 12q14.1 (13.7%) encompassing the oncogene CDK4 was recurrently amplified and negatively affected patient prognosis. Besides frequently mutated RTK/RAS pathway and PI3K pathway, the cell cycle pathway controlling the phosphorylation of retinoblastoma protein (RB) was mutated in 39.2% of patients. Overall, a considerable degree of the somatic mutation landscape is shared between SCAs and brainstem glioma. Our work provides a key insight into the molecular profiling of primary SCAs, which might represent candidate drug targets and complement the molecular atlas of glioma. © 2023 The Pathological Society of Great Britain and Ireland.

Epidemiological characteristics and genetic alterations in adult diffuse glioma in East Asian populations

Understanding the racial specificities of diseases-such as adult diffuse glioma, the most common primary malignant tumor of the central nervous system-is a critical step toward precision medicine. Here, we comprehensively review studies of gliomas in East Asian populations and other ancestry groups to clarify the racial differences in terms of epidemiology and genomic characteristics. Overall, we observed a lower glioma incidence in East Asians than in Whites; notably, patients with glioblastoma had significantly younger ages of onset and longer overall survival than the Whites. Multiple genome-wide association studies of various cohorts have revealed single nucleotide polymorphisms associated with overall and subtype-specific glioma susceptibility. Notably, only 3 risk loci-5p15.33, 11q23.3, and 20q13.33-were shared between patients with East Asian and White ancestry, whereas other loci predominated only in particular populations. For instance, risk loci 12p11.23, 15q15-21.1, and 19p13.12 were reported in East Asians, whereas risk loci 8q24.21, 1p31.3, and 1q32.1 were reported in studies in White patients. Although the somatic mutational profiles of gliomas between East Asians and non-East Asians were broadly consistent, a lower incidence of EGFR amplification in glioblastoma and a higher incidence of 1p19q-IDH-TERT triple-negative low-grade glioma were observed in East Asian cohorts. By summarizing large-scale disease surveillance, germline, and somatic genomic studies, this review reveals the unique characteristics of adult diffuse glioma among East Asians, to guide clinical management and policy design focused on patients with East Asian ancestry.

A noncoding single-nucleotide polymorphism at 8q24 drives IDH1-mutant glioma formation

Establishing causal links between inherited polymorphisms and cancer risk is challenging. Here, we focus on the single-nucleotide polymorphism rs55705857, which confers a sixfold greater risk of isocitrate dehydrogenase (IDH)-mutant low-grade glioma (LGG). We reveal that rs55705857 itself is the causal variant and is associated with molecular pathways that drive LGG. Mechanistically, we show that rs55705857 resides within a brain-specific enhancer, where the risk allele disrupts OCT2/4 binding, allowing increased interaction with the Myc promoter and increased Myc expression. Mutating the orthologous mouse rs55705857 locus accelerated tumor development in an Idh1R132H-driven LGG mouse model from 472 to 172 days and increased penetrance from 30% to 75%. Our work reveals mechanisms of the heritable predisposition to lethal glioma in ~40% of LGG patients.

Inherited genetics of adult diffuse glioma and polygenic risk scores-a review

Knowledge about inherited and acquired genetics of adult diffuse glioma has expanded significantly over the past decade. Genomewide association studies (GWAS) stratified by histologic subtype identified six germline variants that were associated specifically with glioblastoma (GBM) and 12 that were associated with lower grade glioma. A GWAS performed using the 2016 WHO criteria, stratifying patients by IDH mutation and 1p/19q codeletion (as well as TERT promoter mutation), discovered that many of the known variants are associated with specific WHO glioma subtypes. In addition, the GWAS stratified by molecular group identified two additional novel regions: variants in D2HGDH that were associated with tumors that had an IDH mutation and a variant near FAM20C that was associated with tumors that had both IDH mutation and 1p/19q codeletion. The results of these germline associations have been used to calculate polygenic risk scores, from which to estimate relative and absolute risk of overall glioma and risk of specific glioma subtypes. We will review the concept of polygenic risk models and their potential clinical utility, as well as discuss the published adult diffuse glioma polygenic risk models. To date, these prior genetic studies have been done on European populations. Using the published glioma polygenic risk model, we show that the genetic associations published to date do not generalize across genetic ancestries, demonstrating that genetic studies need to be done on more diverse populations.

The immunogenetics of viral antigen response is associated with subtype-specific glioma risk and survival

Glioma is a highly fatal cancer with prognostically significant molecular subtypes and few known risk factors. Multiple studies have implicated infections in glioma susceptibility, but evidence remains inconsistent. Genetic variants in the human leukocyte antigen (HLA) region modulate host response to infection and have been linked to glioma risk. In this study, we leveraged genetic predictors of antibody response to 12 viral antigens to investigate the relationship with glioma risk and survival. Genetic reactivity scores (GRSs) for each antigen were derived from genome-wide-significant (p < 5 × 10-8) variants associated with immunoglobulin G antibody response in the UK Biobank cohort. We conducted parallel analyses of glioma risk and survival for each GRS and HLA alleles imputed at two-field resolution by using data from 3,418 glioma-affected individuals subtyped by somatic mutations and 8,156 controls. Genetic reactivity scores to Epstein-Barr virus (EBV) ZEBRA and EBNA antigens and Merkel cell polyomavirus (MCV) VP1 antigen were associated with glioma risk and survival (Bonferroni-corrected p < 0.01). GRSZEBRA and GRSMCV were associated in opposite directions with risk of IDH wild-type gliomas (ORZEBRA = 0.91, p = 0.0099/ORMCV = 1.11, p = 0.0054). GRSEBNA was associated with both increased risk for IDH mutated gliomas (OR = 1.09, p = 0.040) and improved survival (HR = 0.86, p = 0.010). HLA-DQA1∗03:01 was significantly associated with decreased risk of glioma overall (OR = 0.85, p = 3.96 × 10-4) after multiple testing adjustment. This systematic investigation of the role of genetic determinants of viral antigen reactivity in glioma risk and survival provides insight into complex immunogenomic mechanisms of glioma pathogenesis. These results may inform applications of antiviral-based therapies in glioma treatment.

Fam20C Overexpression Predicts Poor Outcomes and is a Diagnostic Biomarker in Lower-Grade Glioma

Glioma is a relatively low aggressive brain tumor. Although the median survival time of patients for lower-grade glioma (LGG) was longer than that of patients for glioblastoma, the overall survival was still short. Therefore, it is urgent to find out more effective molecular prognostic markers. The role of the Fam20 kinase family in different tumors was an emerging research field. However, the biological function of Fam20C and its prognostic value in brain tumors have rarely been reported. This study aimed to evaluate the value of Fam20C as a potential prognostic marker for LGG. A total of 761 LGG samples (our cohort, TCGA and CGGA) were included to investigate the expression and role of Fam20C in LGG. We found that Fam20C was drastically overexpressed in LGG and was positively associated with its clinical progression. Kaplan-Meier analysis and a Cox regression model were employed to evaluate its prognostic value, and Fam20C was found as an independent risk factor in LGG patients. Gene set enrichment analysis also revealed the potential signaling pathways associated with Fam20C gene expression in LGG; these pathways were mainly enriched in extracellular matrix receptor interactions, cell adhesion, cell apoptosis, NOTCH signaling, cell cycle, etc. In summary, our findings provide insights for understanding the potential role of Fam20C and its application as a new prognostic biomarker for LGG.

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⁻⁴). SMR confirmed the previously reported LTL associations at 3q26.2 (TERC; P_SMR = 1.33 × 10⁻⁵), 5p15.33 (TERT; P_SMR = 9.80 × 10⁻²⁷), 10q24.33 (STN1 alias OBFC1; P_SMR = 4.31 × 10⁻⁵), and 20q13.3 (STMN3/RTEL1; P_SMR = 2.47 × 10⁻⁴) glioma risk loci. Our analysis implicates variation at 1q42.12 (P_SMR = 1.55 × 10⁻²), 6p21.3 (P_SMR = 9.76 × 10⁻³), 6p22.2 (P_SMR = 5.45 × 10⁻³), 7q31.33 (P_SMR = 6.52 × 10⁻³), and 11q22.3 (P_SMR = 8.89 × 10⁻⁴) 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.

Partitioned glioma heritability shows subtype-specific enrichment in immune cells

BACKGROUND

Epidemiological studies of adult glioma have identified genetic syndromes and 25 heritable risk loci that modify individual risk for glioma, as well increased risk in association with exposure to ionizing radiation and decreased risk in association with allergies. In this analysis, we assess whether there is a shared genome-wide genetic architecture between glioma and atopic/autoimmune diseases.

METHODS

Using summary statistics from a glioma genome-wide association studies (GWAS) meta-analysis, we identified significant enrichment for risk variants associated with gene expression changes in immune cell populations. We also estimated genetic correlations between glioma and autoimmune, atopic, and hematologic traits using linkage disequilibrium score regression (LDSC), which leverages genome-wide single-nucleotide polymorphism (SNP) associations and patterns of linkage disequilibrium.

RESULTS

Nominally significant negative correlations were observed for glioblastoma (GB) and primary biliary cirrhosis (rg = -0.26, P = .0228), and for non-GB gliomas and celiac disease (rg = -0.32, P = .0109). Our analyses implicate dendritic cells (GB pHM = 0.0306 and non-GB pHM = 0.0186) in mediating both GB and non-GB genetic predisposition, with GB-specific associations identified in natural killer (NK) cells (pHM = 0.0201) and stem cells (pHM = 0.0265).

CONCLUSIONS

This analysis identifies putative new associations between glioma and autoimmune conditions with genomic architecture that is inversely correlated with that of glioma and that T cells, NK cells, and myeloid cells are involved in mediating glioma predisposition. This provides further evidence that increased activation of the acquired immune system may modify individual susceptibility to glioma.