Umeå University's logo

Genome-wide association studies (GWAS) have contributed to our understanding of glioma susceptibility. To date, 25 risk loci for development of any of the glioma subtypes are known. However, GWAS studies reveal little about the molecular processes that lead to increased risk, especially for non-coding single nucleotide polymorphisms (SNP). A particular SNP in intron 2 of LRIG1, rs11706832, has been shown to increase the susceptibility for IDH1 mutated low-grade gliomas (LGG). Leucine-rich repeats and immunoglobulin-like domains protein 1 (LRIG1) is important in cancer development as it negatively regulates the epidermal growth factor receptor (EGFR); however, the mechanism responsible for this particular risk SNP and its potential effect on LRIG1 are not known. Using CRISPR-CAS9, we edited rs11706832 in HEK293T cells. Four HEK293T clones with the risk allele were compared to four clones with the non-risk allele for LRIG1 and SLC25A26 gene expression using RT-qPCR, for global gene expression using RNA-seq, and for metabolites using gas chromatography-mass spectrometry (GC–MS). The experiment did not reveal any significant effect of the SNP on the expression levels or splicing patterns of LRIG1 or SLC25A26. The global gene expression analysis revealed that the risk allele C was associated with upregulation of several mitochondrial genes. Gene enrichment analysis of 74 differentially expressed genes in the genome revealed a significant enrichment of type I interferon response genes, where many genes were downregulated for the risk allele C. Gene expression data of IDH1 mutated LGGs from the cancer genome atlas (TCGA) revealed a similar under expression of type I interferon genes associated with the risk allele. This study found the expression levels and splicing patterns of LRIG1 and SLC25A26 were not affected by the SNP in HEK293T cells. However, the risk allele was associated with a downregulation of genes involved in the innate immune response both in the HEK293T cells and in the LGG data from TCGA.

BACKGROUND: Although recent sequencing studies have revealed that 10% of childhood gliomas are caused by rare germline mutations, the role of common variants is undetermined and no genome-wide significant risk loci for pediatric central nervous system tumors have been identified to date.

METHODS: Meta-analysis of 3 population-based genome-wide association studies comprising 4069 children with glioma and 8778 controls of multiple genetic ancestries. Replication was performed in a separate case-control cohort. Quantitative trait loci analyses and a transcriptome-wide association study were conducted to assess possible links with brain tissue expression across 18 628 genes.

RESULTS: Common variants in CDKN2B-AS1 at 9p21.3 were significantly associated with astrocytoma, the most common subtype of glioma in children (rs573687, P-value of 6.974e-10, OR 1.273, 95% CI 1.179-1.374). The association was driven by low-grade astrocytoma (P-value of 3.815e-9) and exhibited unidirectional effects across all 6 genetic ancestries. For glioma overall, the association approached genome-wide significance (rs3731239, P-value of 5.411e-8), while no significant association was observed for high-grade tumors. Predicted decreased brain tissue expression of CDKN2B was significantly associated with astrocytoma (P-value of 8.090e-8).

CONCLUSIONS: In this population-based genome-wide association study meta-analysis, we identify and replicate 9p21.3 (CDKN2B-AS1) as a risk locus for childhood astrocytoma, thereby establishing the first genome-wide significant evidence of common variant predisposition in pediatric neuro-oncology. We furthermore provide a functional basis for the association by showing a possible link to decreased brain tissue CDKN2B expression and substantiate that genetic susceptibility differs between low- and high-grade astrocytoma.