The Genetic Architecture of Gliomagenesis-Genetic Risk Variants Linked to Specific Molecular Subtypes

Structural and biochemical characterization of the C-terminal region of the human RTEL1 helicase

RTEL1 is an essential DNA helicase which plays an important role in various aspects of genome stability, from telomere metabolism to DNA replication, repair and recombination. RTEL1 has been implicated in a number of genetic diseases and cancer development, including glioma, breast, lung and gastrointestinal tumors. RTEL1 is a FeS helicase but, in addition to the helicase core, it comprises a long C-terminal region which includes a number of folded domains connected by intrinsically disordered loops and mediates RTEL1 interaction with factors involved in pivotal cellular pathways. However, information on the architecture and the function of this region is still limited. We expressed and purified a variety of fragments encompassing the folded domains and the unstructured regions. We determined the crystal structure of the second repeat, confirming that it has a fold similar to the harmonin homology domains. SAXS data provide low-resolution information on all the fragments and suggest that the presence of the RING domain affects the overall architecture of the C-terminal region, making the structure significantly more compact. NMR data provide experimental information on the interaction between PCNA and the RTEL1 C-terminal region, revealing a putative low-affinity additional site of interaction. A biochemical analysis shows that the C-terminal region, in addition to a preference for telomeric RNA and DNA G-quadruplexes, has a high affinity for R-loops and D-loops, consistent with the role played by the RTEL1 helicase in homologous recombination, telomere maintenance and preventing replication-transcription conflicts. We further dissected the contribution of each domain in binding different substrates.

Multi-ancestry genome-wide association study of 4069 children with glioma identifies 9p21.3 risk locus

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.

The Long Non-Coding RNA ANRIL in Cancers

ANRIL (Antisense Noncoding RNA in the INK4 Locus), a long non-coding RNA encoded in the human chromosome 9p21 region, is a critical factor for regulating gene expression by interacting with multiple proteins and miRNAs. It has been found to play important roles in various cellular processes, including cell cycle control and proliferation. Dysregulation of ANRIL has been associated with several diseases like cancers and cardiovascular diseases, for instance. Understanding the oncogenic role of ANRIL and its potential as a diagnostic and prognostic biomarker in cancer is crucial. This review provides insights into the regulatory mechanisms and oncogenic significance of the 9p21 locus and ANRIL in cancer.

Low-grade glioma risk SNP rs11706832 is associated with type I interferon response pathway genes in cell lines

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.

Prediagnostic biomarkers for early detection of glioma—using case–control studies from cohorts as study approach

Background

Understanding the trajectory and development of disease is important and the knowledge can be used to find novel targets for therapy and new diagnostic tools for early diagnosis.

Methods

Large cohorts from different parts of the world are unique assets for research as they have systematically collected plasma and DNA over long-time periods in healthy individuals, sometimes even with repeated samples. Over time, the population in the cohort are diagnosed with many different diseases, including brain tumors.

Results

Recent studies have detected genetic variants that are associated with increased risk of glioblastoma and lower grade gliomas specifically. The impact for genetic markers to predict disease in a healthy population has been deemed low, and a relevant question is if the genetic variants for glioma are associated with risk of disease or partly consist of genes associated to survival. Both metabolite and protein spectra are currently being explored for early detection of cancer.

Conclusions

We here present a focused review of studies of genetic variants, metabolomics, and proteomics studied in prediagnostic glioma samples and discuss their potential in early diagnostics.

LncRNA CDKN2B-AS1 contributes to glioma development by regulating the miR-199a-5p/DDR1 axis

BACKGROUND

Although cyclin-dependent kinase inhibitor 2B antisense RNA 1 (CDKN2B-AS1) is upregulated in glioma, its function and potential mechanism in glioma remain unclear.

METHODS

CDKN2B-AS1 level in glioma tissues and cell lines LN229, U251, and U87 was measured by qRT-PCR. Loss-of-function assays using short hairpin RNA for CDKN2B-AS1 (sh-CDKN2B-AS1) were performed to evaluate the effect of CDKN2B-AS1 on cell invasion, migration, proliferation, and apoptosis. The relationship among CDKN2B-AS1, miR-199a-5p, and DDR1 was determined by bioinformatics analysis and luciferase reporter assay. Rescue experiments were conducted to explore the function of CDKN2B-AS1 and miR-199a-5p in glioma. An in vivo animal model of lentivirally transduced U87 glioma xenografts in mice was established to confirm the role of CDKN2B-AS1.

RESULTS

CDKN2B-AS1 is significantly upregulated in glioma tissues and cell lines. CDKN2B-AS1 knockdown significantly inhibits cell proliferation, invasion, and migration, while promoting apoptosis of glioma cell lines U251 and U87. Further, a miR-199a-5p inhibitor attenuates the inhibitory effects of sh-CDKN2B-AS1 on these cell phenotypes. CDKN2B-AS1 positively regulates DDR1 expression by directly sponging miR-199a-5p. Moreover, CDKN2B-AS1 knockdown efficiently inhibits U87 tumor xenograft growth in mice.

CONCLUSION

Our study reveals that CDKN2B-AS1 promotes glioma development by regulating the miR-199a-5p/DDR1 axis, suggesting that this lncRNA might be a potential therapeutic target.

A Genome-Wide Profiling of Glioma Patients with an IDH1 Mutation Using the Catalogue of Somatic Mutations in Cancer Database

Simple Summary

Glioma patients that present a somatic mutation in the isocitrate dehydrogenase 1 (IDH1) gene have a significantly better prognosis and overall survival than patients with the wild-type genotype. An IDH1 mutation is hypothesized to occur early during cellular transformation and leads to further genetic instability. A genome-wide profiling of glioma patients in the Catalogue of Somatic Mutations in Cancer (COSMIC) database was performed to classify the genetic differences in IDH1-mutant versus IDH1-wildtype patients. This classification will aid in a better understanding of how this specific mutation influences the genetic make-up of glioma and the resulting prognosis. Key differences in co-mutation and gene expression levels were identified that correlate with an improved prognosis.

Abstract

Gliomas are differentiated into two major disease subtypes, astrocytoma or oligodendroglioma, which are then characterized as either IDH (isocitrate dehydrogenase)-wild type or IDH-mutant due to the dramatic differences in prognosis and overall survival. Here, we investigated the genetic background of IDH1-mutant gliomas using the Catalogue of Somatic Mutations in Cancer (COSMIC) database. In astrocytoma patients, we found that IDH1 is often co-mutated with TP53, ATRX, AMBRA1, PREX1, and NOTCH1, but not CHEK2, EGFR, PTEN, or the zinc finger transcription factor ZNF429. The majority of the mutations observed in these genes were further confirmed to be either drivers or pathogenic by the Cancer-Related Analysis of Variants Toolkit (CRAVAT). Gene expression analysis showed down-regulation of DRG2 and MSN expression, both of which promote cell proliferation and invasion. There was also significant over-expression of genes such as NDRG3 and KCNB1 in IDH1-mutant astrocytoma patients. We conclude that IDH1-mutant glioma is characterized by significant genetic changes that could contribute to a better prognosis in glioma patients.

Meta-Analyses of Splicing and Expression Quantitative Trait Loci Identified Susceptibility Genes of Glioma

Background

The functions of most glioma risk alleles are unknown. Very few studies had evaluated expression quantitative trait loci (eQTL), and insights of susceptibility genes were limited due to scarcity of available brain tissues. Moreover, no prior study had examined the effect of glioma risk alleles on alternative RNA splicing.

Objective

This study explored splicing quantitative trait loci (sQTL) as molecular QTL and improved the power of QTL mapping through meta-analyses of both cis eQTL and sQTL.

Methods

We first evaluated eQTLs and sQTLs of the CommonMind Consortium (CMC) and Genotype-Tissue Expression Project (GTEx) using genotyping, or whole-genome sequencing and RNA-seq data. Alternative splicing events were characterized using an annotation-free method that detected intron excision events. Then, we conducted meta-analyses by pooling the eQTL and sQTL results of CMC and GTEx using the inverse variance-weighted model. Afterward, we integrated QTL meta-analysis results (Q < 0.05) with the Glioma International Case Control Study (GICC) GWAS meta-analysis (case:12,496, control:18,190), using a summary statistics-based mendelian randomization (SMR) method.

Results

Between CMC and GTEx, we combined the QTL data of 354 unique individuals of European ancestry. SMR analyses revealed 15 eQTLs in 11 loci and 32 sQTLs in 9 loci relevant to glioma risk. Two loci only harbored sQTLs (1q44 and 16p13.3). In seven loci, both eQTL and sQTL coexisted (2q33.3, 7p11.2, 11q23.3 15q24.2, 16p12.1, 20q13.33, and 22q13.1), but the target genes were different for five of these seven loci. Three eQTL loci (9p21.3, 20q13.33, and 22q13.1) and 4 sQTL loci (11q23.3, 16p13.3, 16q12.1, and 20q13.33) harbored multiple target genes. Eight target genes of sQTLs (C2orf80, SEC61G, TMEM25, PHLDB1, RP11-161M6.2, HEATR3, RTEL1-TNFRSF6B, and LIME1) had multiple alternatively spliced transcripts.

Conclusion

Our study revealed that the regulation of transcriptome by glioma risk alleles is complex, with the potential for eQTL and sQTL jointly affecting gliomagenesis in risk loci. QTLs of many loci involved multiple target genes, some of which were specific to alternative splicing. Therefore, quantitative trait loci that evaluate only total gene expression will miss many important target genes.

A functional variant on 20q13.33 related to glioma risk alters enhancer activity and modulates expression of multiple genes

Genome‐wide association studies (GWAS) have identified single‐nucleotide polymorphisms (SNPs) associated with glioma risk on 20q13.33, but the biological mechanisms underlying this association are unknown. We tested the hypothesis that a functional SNP on 20q13.33 impacted the activity of an enhancer, leading to an altered expression of nearby genes. To identify candidate functional SNPs, we identified all SNPs in linkage disequilibrium with the risk‐associated SNP rs2297440 that mapped to putative enhancers. Putative enhancers containing candidate functional SNPs were tested for allele‐specific effects in luciferase enhancer activity assays against glioblastoma multiforme (GBM) cell lines. An enhancer containing SNP rs3761124 exhibited allele‐specific effects on activity. Deletion of this enhancer by CRISPR‐Cas9 editing in GBM cell lines correlated with an altered expression of multiple genes, including STMN3, RTEL1, RTEL1‐TNFRSF6B, GMEB2, and SRMS. Expression quantitative trait loci (eQTL) analyses using nondiseased brain samples, isocitrate dehydrogenase 1 (IDH1) wild‐type glioma, and neurodevelopmental tissues showed STMN3 to be a consistent significant eQTL with rs3761124. RTEL1 and GMEB2 were also significant eQTLs in the context of early CNS development and/or in IDH1 wild‐type glioma. We provide evidence that rs3761124 is a functional variant on 20q13.33 related to glioma/GBM risk that modulates the expression of STMN3 and potentially other genes across diverse cellular contexts.