RNA-seq of 272 gliomas revealed a novel, recurrent PTPRZ1-MET fusion transcript in secondary glioblastomas

Case of multifocal glioblastoma with four fusion transcripts of ALK, FGFR2, NTRK2, and NTRK3 genes stresses the need for tumor tissue multisampling for transcriptomic analysis

Glioblastoma multiforme (GBM) is the most malignant brain tumor with patient mortality rate close to 100%, 5-yr survival rate of ∼5%, and a median survival of 14 mo. GBMs have notorious histomorphologic and molecular heterogeneities thus giving hope for development of future personalized therapies. We describe here a case of a 48-yr-old male patient with three-nodular GBM. To address the question of intratumoral molecular heterogeneity, a comparative analysis of gene expression was performed by using multiple samples collected from different tumor sites with the aid of intraoperative magnetic resonance imaging (MRI). Sixteen GBM biosamples from parietal, temporal, and temporo-polar localizations were collected from primary, recurrent, and second recurrent tumors and were obtained and investigated by RNA sequencing. Our investigations revealed that biosamples derived from different tumor sites differ in their gene expression profiles with classical or mesenchymal signatures associated with clinically distinct molecular subtypes of GBM found within the same tumor. The results also showed significant differences in the expression of genes specific for targeted therapeutics. Our investigations have enabled the identification of four novel fusion transcripts—KIF5C-NTRK3, AC016907.2-ALK, CNTNAP3-NTRK2, and ZNF135-FGFR2—each present in only one sample. We found no differences between untreated and recurrent stages in the expression levels of genes involved in fusion transcripts, suggesting the lack of association between fusion transcript and treatment response. In contrast, longitudinal changes in the expression of VEGF and MGMT genes were concordant with the tumor response to bevacizumab and temozolomide. Our study underscores the importance of integrating a multisampling approach and RNA sequencing and demonstrates the predictive merit of an integrated approach for differentiating genomic aberrations associated with untreated or post-treatment recurrent GBMs.

Development and Validation of a Radiosensitivity Prediction Model for Lower Grade Glioma Based on Spike-and-Slab Lasso

Background and Purpose

Lower grade glioma (LGG) is one of the leading causes of death world worldwide. We attempted to develop and validate a radiosensitivity model for predicting the survival of lower grade glioma by using spike-and-slab lasso Cox model.

Methods

In this research, differentially expressed genes based on tumor microenvironment was obtained to further analysis. Log-rank test was used to identify genes in patients who received radiotherapy and patients who did not receive radiotherapy, respectively. Then, spike-and-slab lasso was performed to select genes in patients who received radiotherapy. Finally, three genes (INA, LEPREL1 and PTCRA) were included in the model. A radiosensitivity-related risk score model was established based on overall rate of TCGA dataset in patients who received radiotherapy. The model was validated in TCGA dataset that PFS as endpoint and two CGGA datasets that OS as endpoint. A novel nomogram integrated risk score with age and tumor grade was developed to predict the OS of LGG patients.

Results

We developed and verified a radiosensitivity-related risk score model. The radiosensitivity-related risk score is served as an independent prognostic indicator. This radiosensitivity-related risk score model has prognostic prediction ability. Moreover, the nomogram integrated risk score with age and tumor grade was established to perform better for predicting 1, 3, 5-year survival rate.

Conclusions

This model can be used by clinicians and researchers to predict patient’s survival rates and achieve personalized treatment of LGG.

Enhanced B7-H4 expression in gliomas with low PD-L1 expression identifies super-cold tumors

Background

Characterizing expression profiles of different immune checkpoint molecules are promising for personalized checkpoint inhibitory immunotherapy. Gliomas have been shown as potential targets for immune checkpoint inhibitors recently. Our study was performed to determine coexpression levels of two major B7 immune regulatory molecules programmed death ligand 1 (PD-L1) and B7-H4, both of which have been demonstrated to inhibit antitumor host immunity in gliomas.

Methods

We assessed tumor tissues from stage II–IV primary gliomas (n=505) by immunohistochemistry (IHC) for protein levels of both PD-L1 and B7-H4. Gene coexpression analysis assessing clusters based on extent of PD-L1/B7-H4 classifier genes expression were investigated in two transcriptome datasets (The Cancer Genome Atlas and Chinese Glioma Genome Atlas). In addition, levels of immune cell infiltrates were estimated with IHC and RNA-seq data for assessing the tumor immune microenvironment of PD-L1/B7-H4 subgroups.

Results

High expression of PD-L1 and B7-H4 in gliomas was 23% and 20%, respectively, whereas coexpression of two proteins at high levels was limited to 2% of the cases. Comparable results were seen in RNA-seq datasets where PD-L1 mRNA expression levels negatively correlated with that of B7-H4. Gene coexpression modules clustered within each grade of gliomas demonstrated lack of double-high modules (cluster with high expression of both PD-L1 and B7-H4 classifier genes). B7-H4 mRNA expression levels showed negative correlation with extent of immune cell infiltration and High-B7-H4 module gliomas (high B7-H4 but low PD-L1 classifier genes expression) had less tumor-infiltrating lymphocytes (TILs) and tumor-associated macrophages (TAMs). IHC assessment also showed few TILs and TAMs in High-B7-H4 subgroup gliomas.

Conclusions

The majority of gliomas express PD-L1 or B7-H4, however, coexpression of both at high levels is minimal. The high-B7-H4 patients could be considered as ‘super-cold’ gliomas with significantly deficient in TILs, suggesting that B7-H4 might inhibit T-cell trafficking into the central nervous system. This study demonstrated that PD-L1 and B7-H4 may serve as mutually compensatory immune checkpoint molecules in gliomas for immune targeted or active-specific immunotherapy. The distinct B7-H4 pathways modulating T-cell function and immune evasion in glioma patients deserved to be further explored in the future during immunotherapy.

Expression, methylation and prognostic feature of EMILIN2 in Low-Grade-Glioma

Low-grade gliomas (LGGs) are slow-growing brain cancer in central nervous system neoplasms. EMILIN2 is an extracellular matrix (ECM) protein which could influence the progress of some tumour which is unclear in LGG. In our study, the methylation, expression, prognosis and immune value of EMILIN2 in LGG were analysed through bioinformatics analysis. We analysed the LGG data from The Cancer Genome Atlas (TCGA) and discovered that the EMILIN2 expression, negatively correlated to the EMILIN2 methylation, could predict a poor prognosis and was associated with different clinical parameters. Moreover, univariate and multivariate Cox regression were performed in CGGA, which showed that the EMILIN2 could be an independent prognostic biomarker in LGG. Moreover, EMILIN2 expression showed a correlation with gene makers in some immune cells, which identified the significance of EMILIN2 in immune infiltration. Finally, we used RT-PCR to verify the EMILIN2 expression level in different grades which showed there were significantly different (P < 0.05). Similarly, high expression of EMILIN2 could predict a poor prognosis (P = 0.0078). In conclusion, EMILIN2 could act as an independent prognostic biomarker which might be associated with the malignancy and development of gliomas and play a crucial role in glioma in immune infiltration.

Integrated Analysis of the Clinical and Molecular Characteristics of IDH Wild-Type Gliomas in the Chinese Glioma Genome Atlas

Purpose

Current studies and guidelines suggest that the biobehavior of IDH-wild type (IDH-wt) lower-grade glioma (LGG, WHO II-III) is similar to IDH-wt glioblastoma (GBM). However, differences in their clinical and molecular characteristics have not been reported. This study aimed to analyze the clinical and genetic information of gliomas with IDH-wt.

Methods

389 patients with IDH-wt were enrolled in the study (LGG=165, GBM=224), and their clinical and genetic information was collected from the Chinese Glioma Genome Atlas (CGGA). We conducted an analysis of this information between the two groups of patients and drew conclusions thereof.

Results

The median age of the LGG patients was 42 (18–74) years, whereas that of the GBM patients was 51 (18–79) years (P < 0.010). GBM patients were more likely to undergo total resection (P = 0.018) and had fewer epileptic seizure symptoms (P < 0.001). The median overall survival (OS) was 55 months for the LGG patients and only 14.83 months for the GBM patients (P < 0.01). The median progression-free survival (PFS) was 44 months for the LGG patients and only 9.767 months for the GBM patients (P < 0.001). GBM patients were more prone to PETN mutations (P = 0.010). Transcriptome analysis showed that the differentially expressed genes in LGG patients were mainly enriched in metabolic pathways and pathways in cancer and in the function of signal transduction and positive regulation of GTPase activity, whereas in GBM patients, they were mainly enriched in the PI3K-Akt signaling pathway and in the functions of apoptotic process and oxidation-reduction process.

Conclusions

Our data indicate that these two groups of patients should be re-evaluated and treated differently, despite both having IDH wild type.

HGF/c-MET pathway in cancer: from molecular characterization to clinical evidence

This review provides a comprehensive landscape of HGF/c-MET (hepatocyte growth factor (HGF) /mesenchymal-epithelial transition factor (c-MET)) signaling pathway in cancers. First, we generalize the compelling influence of HGF/c-MET pathway on multiple cellular processes. Then, we present the genomic characterization of HGF/c-MET pathway in carcinogenesis. Furthermore, we extensively illustrate the malignant biological behaviors of HGF/c-MET pathway in cancers, in which hyperactive HGF/c-MET signaling is considered as a hallmark. In addition, we investigate the current clinical trials of HGF/c-MET-targeted therapy in cancers. We find that although HGF/c-MET-targeted therapy has led to breakthroughs in certain cancers, monotherapy of targeting HGF/c-MET has failed to demonstrate significant clinical efficacy in most cancers. With the advantage of the combinations of HGF/c-MET-targeted therapy, the exploration of more options of combinational targeted therapy in cancers may be the major challenge in the future.

Low expression of CDHR1 is an independent unfavorable prognostic factor in glioma

Background: Analysis of the differentially expressed genes between lower grade glioma (LGG) and glioblastoma (GBM) will identify genes involved in a more aggressive phenotype of glioma. Methods: Differentially expressed genes between GBM and LGG were identified using published datasets. Kaplan-Meier estimator was used to determine the overall survival of different groups of glioma patients. The biological functions of CDHR1 in glioma were tested using CCK-8 and trans-well assays. Results: CCDC109B, CD58, CLIC1, EFEMP2, EMP3, LAMC1, LGALS1, PDLIM1 and TNFRSF1A were over-expressed, while, CDHR1 was down-regulated in GBM in The Cancer Genome Atlas (TCGA), Chinese Glioma Genome Atlas (CGGA), GSE4412 and GSE43378 datasets. Compared with normal brain tissues, CDHR1 was down-regulated in glioma tissues. And low expression of CDHR1 was an unfavorable prognostic factor in glioma. Moreover, CDHR1 was lowly expressed in mesenchymal GBM subtype and lower expression of CDHR1 was associated with the worse clinical prognosis of GBM. Furthermore, CDHR1 was down-regulated in astrocytoma LGG subtype and low expression of CDHR1 was a bad prognosis of LGG. CDHR1 expression levels were also associated with IDH mutation. IDH mutant LGG or GBM patients were with higher CDHR1 expression. High expression of CDHR1 was a favorable prognosis in IDH mutant or IDH wild type LGG patients. CHDR1 expression was associated with MGMT methylation and CDHR1 was down-regulated in chemotherapy un-responsive LGG patients. CDHR1 was an independent prognostic factor and negatively associated with EMP3 expression. Glioma patients with low CDHR1 and high EMP3 expression had worse clinical outcomes. At last, we showed that over-expression of CDHR1 could inhibit glioma cell growth and invasion. Conclusion: Low expression of CDHR1 was an independent unfavorable prognostic factor in glioma.

Elevated TYROBP expression predicts poor prognosis and high tumor immune infiltration in patients with low-grade glioma

Background

Tyrosine protein tyrosine kinase binding protein (TYROBP) binds non-covalently to activated receptors on the surface of various immune cells, and mediates signal transduction and cellular activation. It is dysregulated in various malignancies, although little is known regarding its role in low-grade glioma. The aim of this study is to explore the clinicopathological significance, prognostic value and immune signature of TYROBP expression in low-grade glioma (LGG).

Methods

The differentially expressed genes (DEGs) between glioma samples and normal tissues were identified from two GEO microarray datasets using the limma package. The DEGs overlapping across both datasets were functionally annotated by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. STRING database was used to establish the protein-protein interaction (PPI) of the DEGs. The PPI network was visualized by Cytoscape and cytoHubba, and the core module and hub genes were identified. The expression profile of TYROBP and patient survival were validated in the Oncomine, GEPIA2 and CGGA databases. The correlation between TYROBP expression and the clinicopathologic characteristics were evaluated. Gene Set Enrichment Analysis (GSEA) and single-sample GSEA (ssGSEA) were performed by R based on the LGG data from TCGA. The TIMER2.0 database was used to determine the correlation between TYROBP expression and tumor immune infiltrating cells in the LGG patients. Univariate and multivariate Cox regression analyses were performed to determine the prognostic impact of clinicopathological factors via TCGA database.

Results

Sixty-two overlapping DEGs were identified in the 2 datasets, and were mainly enriched in the response to wounding, focal adhesion, GTPase activity and Parkinson disease pathways. TYROBP was identified through the PPI network and cytoHubba. TYROBP expression levels were significantly higher in the LGG tissues compared to the normal tissues, and was associated with worse prognosis and poor clinicopathological parameters. In addition, GSEA showed that TYROBP was positively correlated to neutrophil chemotaxis, macrophage activation, chemokine signaling pathway, JAK-STAT signaling pathway, and negatively associated with gamma aminobutyric acid signaling pathway, neurotransmitter transport, neuroactive ligand receptor intersection etc. TIMER2.0 and ssGSEA showed that TYROBP expression was significantly associated with the infiltration of neutrophils, macrophages, myeloid dendritic cells and monocytes. The infiltration of the M2 phenotype macrophages, cancer-associated fibroblasts and myeloid dendritic cells correlated to worse prognosis in LGG patients. Finally, multivariate analysis showed that elevated TYROBP expression is an independent risk factor for LGG.

Conclusion

TYROBP is dysregulated in LGG and correlates with immune infiltration. It is a potential therapeutic target and prognostic marker for LGG.

The CTSC-RAB38 Fusion Transcript Is Associated With the Risk of Hemorrhage in Brain Arteriovenous Malformations

Brain arteriovenous malformations (bAVMs) are congenital anomalies of blood vessels that cause intracranial hemorrhage in children and young adults. Chromosomal rearrangements and fusion genes play an important role in tumor pathogenesis, though the role of fusion genes in bAVM pathophysiological processes is unclear. The aim of this study was to identify fusion transcripts in bAVMs and analyze their effects. To identify fusion transcripts associated with bAVM, RNA sequencing was performed on 73 samples, including 66 bAVM and 7 normal cerebrovascular samples, followed by STAR-Fusion analysis. Reverse transcription polymerase chain reaction and Sanger sequencing were applied to verify fusion transcripts. Functional pathway analysis was performed to identify potential effects of different fusion types. A total of 21 fusion transcripts were detected. Cathepsin C (CTSC)-Ras-Related Protein Rab-38 (RAB38) was the most common fusion and was detected in 10 of 66 (15%) bAVM samples. In CTSC-RAB38 fusion-positive samples, CTSC and RAB38 expression was significantly increased and activated immune/inflammatory signaling. Clinically, CTSC-RAB38 fusion bAVM cases had a higher hemorrhage rate than non-CTSC-RAB38 bAVM cases (p < 0.05). Our study identified recurrent CTSC-RAB38 fusion transcripts in bAVMs, which may be associated with bAVM hemorrhage by promoting immune/inflammatory signaling.

The Role of N-myc Downstream-Regulated Gene Family in Glioma Based on Bioinformatics Analysis

Glioma is the most common type of primary tumor in the central nervous system, and the molecular mechanisms remain elusive. N-myc downstream-regulated gene (NDRG) family is reported to take part in the pathogenesis of various diseases, including some preliminary exploration in glioma. However, there has been no bioinformatics analysis of NDRG family in glioma yet. Herein, we focused on the expression changes of NDRGs with their value in predicting patients' prognoses, upstream regulatory mechanisms (DNA mutation, DNA methylation, transcription factors, and microRNA regulation) and gene enrichment analysis based on co-expressed genes with data from public databases. Furthermore, the expression pattern of NDRGs was verified by the paired glioma and peritumoral samples in our institute. It was suggested that NDRGs were differentially expressed genes in glioma. In particular, the lower expression of NDRG2 or NDRG4 could serve as a predictor of higher grade tumor and poorer prognosis. Also, NDRGs might play a crucial role in signal transduction, energy metabolism, and cross-talk among cells in glioma, under the control of a complex regulatory network. This study enables us to better understand the role of NDRGs in glioma and with further research, it may contribute to the development of glioma treatment.

Oncogenic FGFR Fusions Produce Centrosome and Cilia Defects by Ectopic Signaling

A single primary cilium projects from most vertebrate cells to guide cell fate decisions. A growing list of signaling molecules is found to function through cilia and control ciliogenesis, including the fibroblast growth factor receptors (FGFR). Aberrant FGFR activity produces abnormal cilia with deregulated signaling, which contributes to pathogenesis of the FGFR-mediated genetic disorders. FGFR lesions are also found in cancer, raising a possibility of cilia involvement in the neoplastic transformation and tumor progression. Here, we focus on FGFR gene fusions, and discuss the possible mechanisms by which they function as oncogenic drivers. We show that a substantial portion of the FGFR fusion partners are proteins associated with the centrosome cycle, including organization of the mitotic spindle and ciliogenesis. The functions of centrosome proteins are often lost with the gene fusion, leading to haploinsufficiency that induces cilia loss and deregulated cell division. We speculate that this complements the ectopic FGFR activity and drives the FGFR fusion cancers.

RUNX1 and REXO2 are associated with the heterogeneity and prognosis of IDH wild type lower grade glioma

Based on isocitrate dehydrogenase (IDH) alterations, lower grade glioma (LGG) is divided into IDH mutant and wild type subgroups. However, the further classification of IDH wild type LGG was unclear. Here, IDH wild type LGG patients in The Cancer Genome Atlas and Chinese Glioma Genome Atlas were divided into two sub-clusters using non-negative matrix factorization. IDH wild type LGG patients in sub-cluster2 had prolonged overall survival and low frequency of CDKN2A alterations and low immune infiltrations. Differentially expressed genes in sub-cluster1 were positively correlated with RUNX1 transcription factor. Moreover, IDH wild type LGG patients with higher stromal score or immune score were positively correlated with RUNX1 transcription factor. RUNX1 and its target gene REXO2 were up-regulated in sub-cluster1 and associated with the worse prognosis of IDH wild type LGG. RUNX1 and REXO2 were associated with the higher immune infiltrations. Furthermore, RUNX1 and REXO2 were correlated with the worse prognosis of LGG or glioma. IDH wild type LGG in sub-cluster2 was hyper-methylated. REXO2 hyper-methylation was associated with the favorable prognosis of LGG or glioma. At last, we showed that, age, tumor grade and REXO2 expression were independent prognostic factors in IDH wild type LGG.

Robustness of differential gene expression analysis of RNA-seq

RNA-sequencing (RNA-seq) is a relatively new technology that lacks standardisation. RNA-seq can be used for Differential Gene Expression (DGE) analysis, however, no consensus exists as to which methodology ensures robust and reproducible results. Indeed, it is broadly acknowledged that DGE methods provide disparate results. Despite obstacles, RNA-seq assays are in advanced development for clinical use but further optimisation will be needed. Herein, five DGE models (DESeq2, voom + limma, edgeR, EBSeq, NOISeq) for gene-level detection were investigated for robustness to sequencing alterations using a controlled analysis of fixed count matrices. Two breast cancer datasets were analysed with full and reduced sample sizes. DGE model robustness was compared between filtering regimes and for different expression levels (high, low) using unbiased metrics. Test sensitivity estimated as relative False Discovery Rate (FDR), concordance between model outputs and comparisons of a ’population’ of slopes of relative FDRs across different library sizes, generated using linear regressions, were examined. Patterns of relative DGE model robustness proved dataset-agnostic and reliable for drawing conclusions when sample sizes were sufficiently large. Overall, the non-parametric method NOISeq was the most robust followed by edgeR, voom, EBSeq and DESeq2. Our rigorous appraisal provides information for method selection for molecular diagnostics. Metrics may prove useful towards improving the standardisation of RNA-seq for precision medicine.

The relationship between PLOD1 expression level and glioma prognosis investigated using public databases

Background

Glioma is the most common type of intracranial tumor with high malignancy and poor prognosis despite the use of various aggressive treatments. Targeted therapy and immunotherapy are not effective and new biomarkers need to be explored. Some Procollagen-lysine 2-oxyglutarate 5-dioxygenase (PLOD) family members have been found to be involved in the metastasis and progression of tumors. Both PLOD2 and PLOD3 had been reported to be highly expressed in gliomas, while the prognostic value of PLOD1 remains to be further illustrated, so we want to investigate the PLOD1 expression in glioma and its clinical implication.

Methods

We collected gene expression and corresponding clinical data of glioma from the Chinese Glioma Genome Atlas (CGGA) database, The Cancer Genome Atlas (TCGA) database and the Gene Expression Omnibus (GEO) database. First, we analyzed the expression and mutation of PLOD1 in gliomas and its relationship with clinicopathologic characteristics. Then, we conducted survival analysis, prognostic analysis and nomogram construction of the PLOD1 gene. Finally, we conducted gene ontology (GO) enrichment analysis and gene set enrichment analysis (GSEA) to explore possible mechanisms and gene co-expression analysis was also be performed.

Results

The results showed that the expression level of PLOD1 was higher in gliomas than normal tissues, and high expression of PLOD1 was related to poor survival which can serve as an oncogenic factor and an independent prognostic indicator for glioma patients. Both the GO and GSEA analysis showed high expression of PLOD1 were enriched in Extracellular matrix (ECM) related pathways, the co-expression analysis revealed that PLOD1 was positively related to HSPG2, COL6A2, COL4A2, FN1, COL1A1, COL4A1, CD44, COL3A1, COL1A2 and SPP1, and high expression of these genes were also correlated to poor prognosis of glioma.

Conclusions

The results showed that high expression of PLOD1 leads to poor prognosis, and PLOD1 is an independent prognostic factor and a novel biomarker for the treatment of glioma. Furthermore, targeting PLOD1 is most likely a potential therapeutic strategy for glioma patients.

Ferroptosis-related gene signature predicts prognosis and immunotherapy in glioma

Aims

Glioma is a highly invasive brain tumor, which makes prognosis challenging and renders patients resistant to various treatments. Induction of cell death is promising in cancer therapy. Ferroptosis, a recently discovered regulated cell death, can be induced for killing glioma cells. However, the prognostic prediction of ferroptosis‐related genes (FRGs) in glioma remains elusive.

Methods

The mRNA expression profiles and gene variation and corresponding clinical data of glioma patients and NON‐TUMOR control were downloaded from public databases. Risk score based on a FRGs signature was constructed in REMBRANDT cohort and validated in other datasets including CGGA‐693, CGGA‐325, and TCGA.

Results

Our results demonstrated that the majority of FRGs was differentially expressed among GBM, LGG, and NON‐TUMOR groups (96.6%). Furthermore, the glioma patients with low‐risk score exhibited a more satisfactory clinical outcome. The better prognosis was also validated in the glioma patients with low‐risk score no matter to which grade they were affiliated. Functional analysis revealed that the high‐risk score group was positively correlated with the enrichment scores for immune checkpoint blockade‐related positive signatures, indicating the critical role of glioma immunotherapy via risk score.

Conclusion

A novel FRGs‐related risk score can predict prognosis and immunotherapy in glioma patients.

Integrated microenvironment-associated genomic profiles identify LRRC15 mediating recurrent glioblastoma-associated macrophages infiltration

Glioblastoma (GBM) is the most common malignant intracranial tumour with intrinsic infiltrative characteristics, which could lead to most patients eventually relapse. The prognosis of recurrent GBM patients remains unsatisfactory. Cancer cell infiltration and their interaction with the tumour microenvironment (TME) could promote tumour recurrence and treatment resistance. In our study, we aimed to identify potential tumour target correlated with rGBM microenvironment based on the gene expression profiles and clinical information of rGBM patients from The Cancer Genome Atlas (TCGA) database. LRRC15 gene with prognostic value was screened by univariate and multivariate analysis, and the correlation between macrophages and LRRC15 was identified as well. Furthermore, the prognosis correlation and immune characteristics of LRRC15 were validated using the Chinese Glioma Genome Atlas (CGGA) database and our clinical tissues by immunochemistry assay. Additionally, we utilized the transwell assay and carboxy fluorescein succinimidyl ester (CFSE) tracking to further confirm the effects of LRRC15 on attracting microglia/macrophages and tumour cell proliferation in the TME. Gene profiles‐based rGBM microenvironment identified that LRRC15 could act in collusion with microglia/macrophages in the rGBM microenvironment to promote the poor prognosis, especially in mesenchymal subtype, indicating the strategies of targeting LRRC15 to improve macrophages‐based immunosuppressive effects could be promising for rGBM treatments.

Identification of CRYAB+ KCNN3+ SOX9+ Astrocyte-Like and EGFR+ PDGFRA+ OLIG1+ Oligodendrocyte-Like Tumoral Cells in Diffuse IDH1-Mutant Gliomas and Implication of NOTCH1 Signalling in Their Genesis

Simple Summary

Diffuse grade II IDH-mutant gliomas are rare brain tumors mainly affecting young patients. These tumors are composed of different populations of tumoral cells. Little is known of these cells and how they are generated. These different cells may show different sensitivity to treatments, so our aim was to study them in detail by directly using patient resections. We identified two clearly distinct tumoral populations and defined reliable markers for them. We also uncovered part of the molecular mechanisms that generate them. Finally, we found that the two cell types have different electrical activity. This article provides unique data and new issues on these rare tumors, which need to be further investigated to develop innovative treatments.

Abstract

Diffuse grade II IDH-mutant gliomas are slow-growing brain tumors that progress into high-grade gliomas. They present intratumoral cell heterogeneity, and no reliable markers are available to distinguish the different cell subtypes. The molecular mechanisms underlying the formation of this cell diversity is also ill-defined. Here, we report that SOX9 and OLIG1 transcription factors, which specifically label astrocytes and oligodendrocytes in the normal brain, revealed the presence of two largely nonoverlapping tumoral populations in IDH1-mutant oligodendrogliomas and astrocytomas. Astrocyte-like SOX9⁺ cells additionally stained for APOE, CRYAB, ID4, KCNN3, while oligodendrocyte-like OLIG1⁺ cells stained for ASCL1, EGFR, IDH1, PDGFRA, PTPRZ1, SOX4, and SOX8. GPR17, an oligodendrocytic marker, was expressed by both cells. These two subpopulations appear to have distinct BMP, NOTCH1, and MAPK active pathways as stainings for BMP4, HEY1, HEY2, p-SMAD1/5 and p-ERK were higher in SOX9⁺ cells. We used primary cultures and a new cell line to explore the influence of NOTCH1 activation and BMP treatment on the IDH1-mutant glioma cell phenotype. This revealed that NOTCH1 globally reduced oligodendrocytic markers and IDH1 expression while upregulating APOE, CRYAB, HEY1/2, and an electrophysiologically-active Ca²⁺-activated apamin-sensitive K⁺ channel (KCNN3/SK3). This was accompanied by a reduction in proliferation. Similar effects of NOTCH1 activation were observed in nontumoral human oligodendrocytic cells, which additionally induced strong SOX9 expression. BMP treatment reduced OLIG1/2 expression and strongly upregulated CRYAB and NOGGIN, a negative regulator of BMP. The presence of astrocyte-like SOX9⁺ and oligodendrocyte-like OLIG1⁺ cells in grade II IDH1-mutant gliomas raises new questions about their role in the pathology.

Intratumor heterogeneity, microenvironment, and mechanisms of drug resistance in glioma recurrence and evolution

Glioma is the most common lethal tumor of the human brain. The median survival of patients with primary World Health Organization grade IV glioma is only 14.6 months. The World Health Organization classification of tumors of the central nervous system categorized gliomas into lower-grade gliomas and glioblastomas. Unlike primary glioblastoma that usually develop de novo in the elderly, secondary glioblastoma enriched with an isocitrate dehydrogenase mutant typically progresses from lower-grade glioma within 5-10 years from the time of diagnosis. Based on various evolutional trajectories brought on by clonal and subclonal alterations, the evolution patterns of glioma vary according to different theories. Some important features distinguish the normal brain from other tissues, e.g., the composition of the microenvironment around the tumor cells, the presence of the blood-brain barrier, and others. The underlying mechanism of glioma recurrence and evolution patterns of glioma are different from those of other types of cancer. Several studies correlated tumor recurrence with tumor heterogeneity and the immune microenvironment. However, the detailed reasons for the progression and recurrence of glioma remain controversial. In this review, we introduce the different mechanisms involved in glioma progression, including tumor heterogeneity, the tumor microenvironment and drug resistance, and their pre-clinical implements in clinical trials. This review aimed to provide new insights into further clinical strategies for the treatment of patients with recurrent and secondary glioma.

Comprehensive analysis of multi-omics data of recurrent gliomas identifies a recurrence-related signature as a novel prognostic marker

Tumor recurrence is a common clinical dilemma in diffuse gliomas. We aimed to identify a recurrence-related signature to predict the prognosis for glioma patients. In the public Chinese Glioma Genome Atlas dataset, we enrolled multi-omics data including genome, epigenome and transcriptome across primary and recurrent gliomas. We included RNA sequencing data from the batch 1 patients (325 patients) as the training set, while RNA sequencing data from the batch 2 patients (693 patients) were selected as the validation set. The R language was used for subsequent analysis. Compared with primary gliomas, more somatic mutations and copy number alterations were revealed in recurrent gliomas. In recurrent gliomas, we identified 113 genes whose methylation levels were significantly different from those of the primary glioma. Through differential expression analysis between primary and recurrent gliomas, we screened 121 recurrence-related genes. Based on these 121 gene expression profiles, consensus clustering of 325 patients yielded two robust groups with different molecular and prognostic features. We developed a recurrence-related risk signature with the lasso regression algorithm. High-risk group had shorter survival and earlier tumor recurrence than the low-risk group. Compared with traditional indicators, the signature showed better prognostic value. In addition, we constructed a nomogram model to predict glioma survival. Functional characteristics analysis found that the signature was associated with cell division and cell cycle. Immune analysis suggested that immunosuppressive status and macrophages might promote glioma recurrence. We demonstrated a novel 18-gene signature that could effectively predict recurrence and prognosis for glioma patients.

Oncogenic mechanism-based pharmaceutical validation of therapeutics targeting MET receptor tyrosine kinase

Aberrant expression and/or activation of the MET receptor tyrosine kinase is characterized by genomic recombination, gene amplification, activating mutation, alternative exon-splicing, increased transcription, and their different combinations. These dysregulations serve as oncogenic determinants contributing to cancerous initiation, progression, malignancy, and stemness. Moreover, integration of the MET pathway into the cellular signaling network as an addiction mechanism for survival has made this receptor an attractive pharmaceutical target for oncological intervention. For the last 20 years, MET-targeting small-molecule kinase inhibitors (SMKIs), conventional therapeutic monoclonal antibodies (TMABs), and antibody-based biotherapeutics such as bispecific antibodies, antibody–drug conjugates (ADC), and dual-targeting ADCs have been under intensive investigation. Outcomes from preclinical studies and clinical trials are mixed with certain successes but also various setbacks. Due to the complex nature of MET dysregulation with multiple facets and underlying mechanisms, mechanism-based validation of MET-targeting therapeutics is crucial for the selection and validation of lead candidates for clinical trials. In this review, we discuss the importance of various types of mechanism-based pharmaceutical models in evaluation of different types of MET-targeting therapeutics. The advantages and disadvantages of these mechanism-based strategies for SMKIs, conventional TMABs, and antibody-based biotherapeutics are analyzed. The demand for establishing new strategies suitable for validating novel biotherapeutics is also discussed. The information summarized should provide a pharmaceutical guideline for selection and validation of MET-targeting therapeutics for clinical application in the future.

A pan-cancer perspective analysis reveals the opposite prognostic significance of CD133 in lower grade glioma and papillary renal cell carcinoma

CD133 is a valuable prognostic marker in multiple types of cancer. However, the expression, methylation levels, and prognostic relevance of CD133 have not been evaluated in a pan-cancer perspective. The expression and methylation levels of CD133 across different types of cancer were determined using The Cancer Genome Atlas (TCGA) dataset. Univariate cox regression and Kaplan-Meier survival were used to determine the prognostic significance of CD133 expression and methylation. CD133 was highly expressed in papillary renal cell carcinoma (PRCC) or pancreatic adenocarcinoma (PAAD). Correspondingly, PAAD and PRCC had low CD133 methylation levels. Through pan-cancer perspective analysis, we found that CD133 high expression was a poor prognostic factor in lower grade glioma (LGG), while, CD133 high expression was a good prognostic factor in PRCC. Moreover, genes positively correlated with CD133 expression were associated with the poor clinical outcomes of LGG. In PRCC, genes negatively correlated with CD133 expression were correlated with the poor overall survival. Furthermore, CD133 expression levels were highly correlated with the CD133 methylation levels in LGG or PRCC. Correspondingly, CD133 hypermethylation was a good prognostic factor in LGG. On the contrary, CD133 hypomethylation was a good prognostic factor in PRCC. We also found that CD133 was highly expressed and hypomethylated in wild type IDH subgroup of LGG. CD133 was highly expressed and hypomethylated in low stages and type1 of PRCC. CD133 high expression and hypomethylation were bad prognostic factors in LGG, while, CD133 high expression and hypomethylation were good prognostic factors in PRCC.

Identification of hub genes related to prognosis in glioma

Glioma, a common malignant tumor of the central nervous system, has high invasiveness. The objective of the present study was to identify genes playing an important role in the development of glioma and to reveal their potential research value. Conjoint analysis on the GSE16011 dataset in the Gene Expression Omnibus (GEO) database and the 'Messenger RNA Expression Microarray of Diffuse Gliomas and Controls' dataset and 'RNA sequencing of Diffuse Gliomas' dataset in the Chinese Glioma Genome Atlas (CGGA) database is carried out in the study. The weighted correlation network analysis (WGCNA) was used to carry out co-expression analysis on the GSE16011. Finally, 24 genes significantly related to grade and prognosis were obtained. In addition, there is no report about CACNG2, JPH3, TUBB6 (tubulin β 6 class V), NRSN1, FAM19A2, NALCN, CDH18, GNAL on glioma.

CARD-Associated Risk Score Features the Immune Landscape and Predicts the Responsiveness to Anti-PD-1 Therapy in IDH Wild-Type Gliomas

Background

Proteins containing the caspase recruitment domain (CARD) play critical roles in cell apoptosis and immunity. However, the impact of CARD genes in tumor immune cell infiltration, responsiveness to checkpoint immunotherapy, and clinical outcomes of gliomas remains unclear. Here, we explore using CARD genes to depict the immune microenvironment and predict the responsiveness of gliomas to anti-PD-1 therapy.

Methods

The genome and transcriptome data of 231 patients with isocitrate dehydrogenase wild-type (IDH-wt) gliomas were retrieved from The Cancer Genome Atlas (TCGA) database to screen CARD genes associated with T lymphocyte infiltration in gliomas. Weighted co-expression network and LASSO penalized regression were employed to generate a CARD-associated risk score (CARS). Two independent and publicly available datasets were used to validate the effectiveness of CARS.

Results

The CARS divided the 231 glioma patients into high- and low-risk subgroups with distinct immune microenvironment and molecular features. The high-risk group had high CARS and was characterized by enrichment of dysfunctional T lymphocytes in a profound immunosuppressive microenvironment, whereas the low-risk group had low CARS and exhibited an immune exclusion genotype. Moreover, signaling aberrations including upregulation of PI3K/Akt/mTOR, NF-κB, and TGF-β were found in the high-risk group. In contrast, the activated WNT pathway was more evident in the low-risk group. Furthermore, we found that an elevated CARS indicated a decreased overall survival for IDH-wt gliomas under standard care but a clinical benefit from checkpoint immunotherapy.

Conclusion

This study developed an immune- and prognosis-relevant risk score, which could be used to enhance our understanding of the heterogeneity of immune microenvironment of gliomas and facilitate to identify patients who will benefit from checkpoint immunotherapy.

Prognostic Biomarker ZNF311 and Its Correlation With Tumor Progression and Immune Infiltrates in Glioma

BACKGROUND

Gliomas, particularly high-grade gliomas, are the most common primary brain tumors. From the Chinese Glioma Genome Atlas (CGGA) database, the relationships between the altered molecular pathways and gliomas could be easily observed. A close connection in the occurrence of the pathogenesis exists between the microenvironment, the glioma, and the associated genes.

METHODS

Validation of the role of ZNF311 oncogene was confirmed by data from the CGGA dataset on glioblastoma and low-grade glioma. Furthermore, we used CIBERSORT to analyze the correlation between ZNF311 and cancer immune infiltrates.

RESULTS

According to our analysis, ZNF311 was expressed higher in patients with grade-depended glioma with poor prognosis. In addition, we obtained valuable prognostic results between isocitrate dehydrogenase 1 (IDH1) and ZNF311 through the analysis of integrated correlations. Similarly, we simultaneously revealed the prognostic results between 1p/19q and ZNF311. In addition, we found that ZNF311 is correlated with a large number of tumor-infiltrating immune cells.

CONCLUSIONS

Based on the study findings, we conclude that ZNF311 is potentially a novel biomarker for assessing prognosis and immune infiltration in glioblastoma and diffuse glioma cases.

Aberrations of Genomic Imprinting in Glioblastoma Formation

In human glioblastoma (GBM), the presence of a small population of cells with stem cell characteristics, the glioma stem cells (GSCs), has been described. These cells have GBM potential and are responsible for the origin of the tumors. However, whether GSCs originate from normal neural stem cells (NSCs) as a consequence of genetic and epigenetic changes and/or dedifferentiation from somatic cells remains to be investigated. Genomic imprinting is an epigenetic marking process that causes genes to be expressed depending on their parental origin. The dysregulation of the imprinting pattern or the loss of genomic imprinting (LOI) have been described in different tumors including GBM, being one of the earliest and most common events that occurs in human cancers. Here we have gathered the current knowledge of the role of imprinted genes in normal NSCs function and how the imprinting process is altered in human GBM. We also review the changes at particular imprinted loci that might be involved in the development of the tumor. Understanding the mechanistic similarities in the regulation of genomic imprinting between normal NSCs and GBM cells will be helpful to identify molecular players that might be involved in the development of human GBM.

Developing a lncRNA Signature to Predict the Radiotherapy Response of Lower-Grade Gliomas Using Co-expression and ceRNA Network Analysis

BACKGROUND

Lower-grade glioma (LGG) is a type of central nervous system tumor that includes WHO grade II and grade III gliomas. Despite developments in medical science and technology and the availability of several treatment options, the management of LGG warrants further research. Surgical treatment for LGG treatment poses a challenge owing to its often inaccessible locations in the brain. Although radiation therapy (RT) is the most important approach in this condition and offers more advantages compared to surgery and chemotherapy, it is associated with certain limitations. Responses can vary from individual to individual based on genetic differences. The relationship between non-coding RNA and the response to radiation therapy, especially at the molecular level, is still undefined.

METHODS

In this study, using The Cancer Genome Atlas dataset and bioinformatics, the gene co-expression network that is involved in the response to radiation therapy in lower-grade gliomas was determined, and the ceRNA network of radiotherapy response was constructed based on three databases of RNA interaction. Next, survival analysis was performed for hub genes in the co-expression network, and the high-efficiency biomarkers that could predict the prognosis of patients with LGG undergoing radiotherapy was identified.

RESULTS

We found that some modules in the co-expression network were related to the radiotherapy responses in patients with LGG. Based on the genes in those modules and the three databases, we constructed a ceRNA network for the regulation of radiotherapy responses in LGG. We identified the hub genes and found that the long non-coding RNA, DRAIC, is a potential molecular biomarker to predict the prognosis of radiotherapy in LGG.

High-sensitive clinical diagnostic method for PTPRZ1-MET and the characteristic protein structure contributing to ligand-independent MET activation

Background

PTPRZ1‐MET (ZM) is a critical genetic alteration driving the progression of lower‐grade glioma. Glioma patients harboring ZM could benefit from MET inhibitors. According to the remarkable role of ZM as a driver of glioma progression and indicator of MET inhibitor sensitivity, it is necessary to detect this alteration even when it presents in glioma with relatively fewer copies.

Methods

Herein, we proposed that ZM could be detected with a high‐sensitive method of reverse transcriptase PCR with 50 amplification cycles. Via this newly proposed detection method, we depicted the incidence preference of ZM fusion in a cohort of 485 glioma patients. To further explore the oncogenic nature of ZM, we predicated the protein structure alteration of MET kinase brought by its fusion partner.

Results

The incidence of ZM fusions was much higher than previous report. ZM fusions exhibited a striking preference in lower‐grade glioma and secondary glioblastoma. By contrast, none of patients with primary glioblastoma was detected harboring ZM fusion. In each of the four variants of ZM, the fusion partner segment of MET contained a remarkable coiled‐coil motif. In glioma cells expressing ZM, MET kinase could be activated in a ligand‐independent manner, which might be contributed by the special coiled‐coil structure brought by the fusion partner. Corresponding to the 3D structural analysis and cell line experiment, the ZM positive clinical specimens showed hyperactivations of MET signaling.

Conclusions

ZM fusions are critical drivers of glioma progression and effective target of MET inhibitor. Early detection could be performed with a high‐sensitive method of reverse transcriptase PCR. The hyperactivations of MET signaling driving glioma progression might be contributed by a ligand‐independent activation enabled by the protein structure modification of extracellular domain of MET in ZM fusions.

Tumor hypoxia represses γδ T cell-mediated antitumor immunity against brain tumors

The anatomic location and immunologic characteristics of brain tumors result in strong lymphocyte suppression. Consequently, conventional immunotherapies targeting CD8 T cells are ineffective against brain tumors. Tumor cells escape immunosurveillance by various mechanisms and tumor cell metabolism can affect the metabolic states and functions of tumor-infiltrating lymphocytes. Here, we discovered that brain tumor cells had a particularly high demand for oxygen, which affected γδ T cell-mediated antitumor immune responses but not those of conventional T cells. Specifically, tumor hypoxia activated the γδ T cell protein kinase A pathway at a transcriptional level, resulting in repression of the activatory receptor NKG2D. Alleviating tumor hypoxia reinvigorated NKG2D expression and the antitumor function of γδ T cells. These results reveal a hypoxia-mediated mechanism through which brain tumors and γδ T cells interact and emphasize the importance of γδ T cells for antitumor immunity against brain tumors.

SH3BGRL2 functions as a crucial tumor suppressor in glioblastoma tumorigenesis

Glioblastoma is the most common and severe primary intrinsic tumor of the central nervous system. Glioblastoma harbors glioma stem cells (GSCs) as it not only possesses self-renewal and differentiation properties but also accounts for significant chemotherapy resistance and recurrence. Thus, targeting GSCs may be essential in overcoming the resistance and recurrence thereby improving GBM treatment. However, the underlying mechanism to sustain GSCs remains largely unknown. Here, we report that SH3 domain binding glutamate-rich protein like 2 (SH3BGRL2) is weakly expressed in glioblastoma multiforme (GBM) and isocitrate dehydrogenase1 (IDH1) wildtype GBM and correlated with glioma patients' poor prognosis. Moreover, ectopic expression of SH3BGRL2 significantly inhibited GBM cell growth, migration, and GSCs self-renewal in vitro as well as tumor growth in vivo. Additionally, we found that SH3BGRL2 suppressed SOX2 and CD133 expression, which are key regulators involved in GSCs self-renewal. Collectively, our findings shed additional light on SH3BGRL2 has potential to serve as a biomarker and a potent therapeutic target for patients with glioma.

Epigenetic modulator inhibition overcomes temozolomide chemoresistance and antagonizes tumor recurrence of glioblastoma

Glioblastoma multiforme (GBM) heterogeneity causes a greater number of deaths than any other brain tumor, despite the availability of alkylating chemotherapy. GBM stem-like cells (GSCs) contribute to GBM complexity and chemoresistance, but it remains challenging to identify and target GSCs or factors that control their activity. Here, we identified a specific GSC subset and show that activity of these cells is positively regulated by stabilization of methyl CpG binding domain 3 (MBD3) protein. MBD3 binds to CK1A and to BTRCP E3 ubiquitin ligase, triggering MBD3 degradation, suggesting that modulating this circuit could antagonize GBM recurrence. Accordingly, xenograft mice treated with the CK1A activator pyrvinium pamoate (Pyr-Pam) showed enhanced MBD3 degradation in cells expressing high levels of O6-methylguanine-DNA methyltransferase (MGMT) and in GSCs, overcoming temozolomide chemoresistance. Pyr-Pam blocked recruitment of MBD3 and the repressive nucleosome remodeling and deacetylase (NuRD) complex to neurogenesis-associated gene loci and increased acetyl-histone H3 activity and GSC differentiation. We conclude that CK1A/BTRCP/MBD3/NuRD signaling modulates GSC activation and malignancy, and that targeting this signaling could suppress GSC proliferation and GBM recurrence.

Novel roles of VAT1 expression in the immunosuppressive action of diffuse gliomas

Standard treatment regimen of gliomas has almost reached a bottleneck in terms of survival benefit. Immunotherapy has been explored and applied in glioma treatment. Immunosuppression, as a hallmark of glioma, could be alleviated by inhibiting certain abnormally expressed biomarkers. Here, transcriptome data of 325 whole grade gliomas were collected from the CGGA database. The TCGA RNA sequencing database was used for validation. Western blot was used to verify the expression level of VAT1 on cellular level. The results showed that the expression of VAT1 was positively correlated with the grades of glioma as classified by WHO. A higher expression level of VAT1 was observed in the mesenchymal subtype of gliomas. The area under the curve suggested that the expression level of VAT1 might be a potential prognostic marker of mesenchymal subtype. In survival analysis, we found that patients with high VAT1 expression level tended to have shorter overall survival, which indicated the prognostic value of VAT1 expression. The results of gene ontology analysis showed that most biological processes of VAT1-related genes were involved in immune and inflammatory responses. The results of GSEA analysis showed a negative correlation between VAT1 expression and immune cells. We also identified that the expression of immune checkpoints increased with VAT1 expression. Therefore, the high expression level of VAT1 in patients with glioma was a potential indicator of a lower survival rate for patients with gliomas. Remarkably, VAT1 contributed to glioma-induced immunosuppression and might be a novel target in glioma immunotherapy.

Potential Molecular Mechanism of TNF Superfamily-Related Genes in Glioblastoma Multiforme Based on Transcriptome and Epigenome

Objective: This study aimed to investigate the molecular mechanism of tumor necrosis factor (TNF) superfamily-related genes and potential therapeutic drugs for glioblastoma multiforme (GBM) patients based on transcriptome and epigenome.

Methods: Gene expression data, corresponding clinical data, and methylation data of GBM samples and normal samples in the TCGA-GBM and GTEx datasets were downloaded. The TNF-related genes were obtained, respectively, from two groups in the TCGA dataset. Then, the TNF-related differentially expressed genes (DEGs) were investigated between two groups, followed by enrichment analysis. Moreover, TNF superfamily-related gene expression and upstream methylation regulation were investigated to explore candidate genes and the prognostic model. Finally, the protein expression level of candidate genes was performed, followed by drug prediction analysis.

Results: A total of 41 DEGs including 4 ligands, 18 receptors, and 19 downstream signaling molecules were revealed between two groups. These DEGs were mainly enriched in pathways like TNF signaling and functions like response to TNF. A total of 5 methylation site-regulated prognosis-related genes including TNF Receptor Superfamily Member (TNFRSF) 12A, TNFRSF11B, and CD40 were explored. The prognosis model constructed by 5 genes showed a well-prediction effect on the current dataset and verification dataset. Finally, drug prediction analysis showed that zoledronic acid (ZA)-TNFRSF11B was the unique drug–gene relation in both two databases.

Conclusion: Methylation-driven gene TNFRSF12A might participate in the development of GBM via response to the TNF biological process and TNF signaling pathway and significantly associated with prognosis. ZA that targets TNFRSF11B expression might be a potential effective drug for clinical treatment of GBM.

Identification of a Glycolysis-Related LncRNA Signature to Predict Survival in Diffuse Glioma Patients

Glycolysis refers to one of the critical phenotypes of tumor cells, regulating tumor cell phenotypes and generating sufficient energy for glioma cells. A range of noticeable genes [such as isocitrate dehydrogenase (IDH), phosphatase, and tensin homolog (PTEN), or Ras] overall impact cell proliferation, invasion, cell cycle, and metastasis through glycolysis. Moreover, long non-coding RNAs (LncRNAs) are increasingly critical to disease progression. Accordingly, this study aimed to identify whether glycolysis-related LncRNAs have potential prognostic value for glioma patients. First, co-expression network between glycolysis-related protein-coding RNAs and LncRNAs was established according to Pearson correlation (Filter: |r| > 0.5 & P < 0.001). Furthermore, based on univariate Cox regression, the Least Absolute Shrinkage and Selection Operator (LASSO) analysis and multivariate Cox regression, a predictive model were built; vital glycolysis-related LncRNAs were identified; the risk score of every single patient was calculated. Moreover, receiver operating characteristic (ROC) curve analysis, gene set enrichment analysis (GSEA), GO and KEGG enrichment analysis were performed to assess the effect of risk score among glioma patients. 685 cases (including RNA sequences and clinical information) from two different cohorts of the Chinese Glioma Genome Atlas (CGGA) database were acquired. Based on the mentioned methods, the risk score calculation formula was yielded as follows: Risk score = (0.19 × EXP_FOXD2-AS1) + (−0.27 × EXP_AC062021.1) + (−0.16 × EXP_AF131216.5) + (−0.05 × EXP_LINC00844) + (0.11 × EXP_CRNDE) + (0.35 × EXP_LINC00665). The risk score was independently related to prognosis, and every single mentioned LncRNAs was significantly related to the overall survival of patients. Moreover, functional enrichment analysis indicated that the biologic process of the high-risk score was mainly involved in the cell cycle and DNA replication signaling pathway. This study confirmed that glycolysis-related LncRNAs significantly impact poor prognosis and short overall survival and may act as therapeutic targets in the future.

Systematic analysis of enhancer regulatory circuit perturbation driven by copy number variations in malignant glioma

Background: Enhancers are emerging regulatory regions controlling gene expression in diverse cancer types. However, the functions of enhancer regulatory circuit perturbations driven by copy number variations (CNVs) in malignant glioma are unclear. Therefore, we aimed to investigate the comprehensive enhancer regulatory perturbation and identify potential biomarkers in glioma.

Results: We performed a meta-analysis of the enhancer centered regulatory circuit perturbations in 683 gliomas by integrating CNVs, gene expression, and transcription factors (TFs) binding. We found widespread CNVs of enhancers during glioma progression, and CNVs were associated with the perturbations of enhancer activities. In particular, the degree of perturbations for amplified enhancers was much greater accompanied by the glioma malignant progression. In addition, CNVs and enhancers cooperatively regulated the expressions of cancer-related genes. Genome-wide TF binding profiles revealed that enhancers were pervasively regulated by TFs. A network-based analysis of TF-enhancer-gene regulatory circuits revealed a core TF-gene module (58 interactions including seven genes and 14 TFs) that was associated survival of patients with glioma (p < 0.001). Finally, we validated this prognosis-associated TF-gene regulatory module in an independent cohort. In summary, our analyses provided new molecular insights for enhancer-centered transcriptional perturbation in glioma therapy.

Conclusion: Integrative analysis revealed enhancer regulatory perturbations in glioma and also identified a network module that was associated with patient survival, thereby providing novel insights into enhancer-centered cancer therapy.

Recurrent Glioblastoma: From Molecular Landscape to New Treatment Perspectives

Glioblastoma is the most frequent and aggressive form among malignant central nervous system primary tumors in adults. Standard treatment for newly diagnosed glioblastoma consists in maximal safe resection, if feasible, followed by radiochemotherapy and adjuvant chemotherapy with temozolomide; despite this multimodal treatment, virtually all glioblastomas relapse. Once tumors progress after first-line therapy, treatment options are limited and management of recurrent glioblastoma remains challenging. Loco-regional therapy with re-surgery or re-irradiation may be evaluated in selected cases, while traditional systemic therapy with nitrosoureas and temozolomide rechallenge showed limited efficacy. In recent years, new clinical trials using, for example, regorafenib or a combination of tyrosine kinase inhibitors and immunotherapy were performed with promising results. In particular, molecular targeted therapy could show efficacy in selected patients with specific gene mutations. Nonetheless, some molecular characteristics and genetic alterations could change during tumor progression, thus affecting the efficacy of precision medicine. We therefore reviewed the molecular and genomic landscape of recurrent glioblastoma, the strategy for clinical management and the major phase I-III clinical trials analyzing recent drugs and combination regimens in these patients.

The lncRNA Toolkit: Databases and In Silico Tools for lncRNA Analysis

Long non-coding RNAs (lncRNAs) are a rapidly expanding field of research, with many new transcripts identified each year. However, only a small subset of lncRNAs has been characterized functionally thus far. To aid investigating the mechanisms of action by which new lncRNAs act, bioinformatic tools and databases are invaluable. Here, we review a selection of computational tools and databases for the in silico analysis of lncRNAs, including tissue-specific expression, protein coding potential, subcellular localization, structural conformation, and interaction partners. The assembled lncRNA toolkit is aimed primarily at experimental researchers as a useful starting point to guide wet-lab experiments, mainly containing multi-functional, user-friendly interfaces. With more and more new lncRNA analysis tools available, it will be essential to provide continuous updates and maintain the availability of key software in the future.

N6-methyladenine-related genes affect biological behavior and the prognosis of glioma

Background

Although aberrant expression of N6‐methyladenine (m⁶A) methylation‐related genes contribute to tumorigenesis in many solid tumors, the prognostic value of the m⁶A‐related genes and their correlation with clinicopathological features in gliomas need advanced study.

Methods

The clinical and sequencing data of 288 patients with glioma were extracted from Chinese Glioma Genome Atlas database. By univariate and multivariable Cox regression analysis, the m⁶A‐related prognostic genes were identified, and their correlation with clinicopathological features was further analysis. A nomogram was constructed by R software and the performance of it was assessed by calibration and time‐dependent receiver operating characteristic curve.

Results

Nine m⁶A‐related genes were identified as independent prognostic factors, which were mostly enriched in RNA splicing, regulation of immune response and vesicle‐mediated transport. By expression value and regression coefficient of these genes, we constructed risk score of each patient, which was highly associated with clinicopathological features. Kaplan–Meier curve showed that the prognosis of patients with high‐risk scores was significantly worse than that with low‐risk scores (HR = 4.30, 95% CI = 3.16–5.85, p < 0.0001). A nomogram was constructed based on the nine m⁶A‐related genes signature and clinicopathological features with well‐fitted calibration curves (c‐index = 0.82), showing high specificity and sensitivity (area under the curve for 1‐, 3‐, and 5‐years survival probability = 0.874, 0.918, and 0.934).

Conclusions

A nine m⁶A‐related genes signature was identified in gliomas. The m⁶A‐related risk score is a novel prognostic factor for patients with glioma, and is associated with clinicopathological features. Moreover, the nomogram based on the nine m⁶A‐related genes signature and clinicopathological features had good efficacy in predicting the survival probability.

FGFR3 overexpression is a useful detection tool for FGFR3 fusions and sequence variations in glioma

Background

Fibroblast growth factor receptor (FGFR) inhibitors are currently used in clinical development. A subset of glioblastomas carries gene fusion of FGFR3 and transforming acidic coiled-coil protein 3. The prevalence of other FGFR3 alterations in glioma is currently unclear.

Methods

We performed RT-PCR in 101 glioblastoma samples to detect FGFR3-TACC3 fusions ("RT-PCR cohort") and correlated results with FGFR3 immunohistochemistry (IHC). Further, we applied FGFR3 IHC in 552 tissue microarray glioma samples ("TMA cohort") and validated these results in two external cohorts with 319 patients. Gene panel sequencing was carried out in 88 samples ("NGS cohort") to identify other possible FGFR3 alterations. Molecular modeling was performed on newly detected mutations.

Results

In the "RT-PCR cohort," we identified FGFR3-TACC3 fusions in 2/101 glioblastomas. Positive IHC staining was observed in 73/1024 tumor samples of which 10 were strongly positive. In the "NGS cohort," we identified FGFR3 fusions in 9/88 cases, FGFR3 amplification in 2/88 cases, and FGFR3 gene mutations in 7/88 cases in targeted sequencing. All FGFR3 fusions and amplifications and a novel FGFR3 K649R missense mutation were associated with FGFR3 overexpression (sensitivity and specificity of 93% and 95%, respectively, at cutoff IHC score > 7). Modeling of these data indicated that Tyr647, a residue phosphorylated as a part of FGFR3 activation, is affected by the K649R mutation.

Conclusions

FGFR3 IHC is a useful screening tool for the detection of FGFR3 alterations and could be included in the workflow for isocitrate dehydrogenase (IDH) wild-type glioma diagnostics. Samples with positive FGFR3 staining could then be selected for NGS-based diagnostic tools.

Genetic and epigenetic landscape of IDH-wildtype glioblastomas with FGFR3-TACC3 fusions

A subset of glioblastomas (GBMs) harbors potentially druggable oncogenic FGFR3-TACC3 (F3T3) fusions. However, their associated molecular and clinical features are poorly understood. Here we analyze the frequency of F3T3-fusion positivity, its associated genetic and methylation profiles, and its impact on survival in 906 IDH-wildtype GBM patients. We establish an F3T3 prevalence of 4.1% and delineate its associations with cancer signaling pathway alterations. F3T3-positive GBMs had lower tumor mutational and copy-number alteration burdens than F3T3-wildtype GBMs. Although F3T3 fusions were predominantly mutually exclusive with other oncogenic RTK pathway alterations, they did rarely co-occur with EGFR amplification. They were less likely to harbor TP53 alterations. By methylation profiling, they were more likely to be assigned the mesenchymal or RTK II subclass. Despite being older at diagnosis and having similar frequencies of MGMT promoter hypermethylation, patients with F3T3-positive GBMs lived about 8 months longer than those with F3T3-wildtype tumors. While consistent with IDH-wildtype GBM, F3T3-positive GBMs exhibit distinct biological features, underscoring the importance of pursuing molecular studies prior to clinical trial enrollment and targeted treatment.

Weighted gene correlation network analysis identifies microenvironment-related genes signature as prognostic candidate for Grade II/III glioma

Glioma is the most common malignant tumor in the central nervous system. Evidence shows that clinical efficacy of immunotherapy is closely related to the tumor microenvironment. This study aims to establish a microenvironment-related genes (MRGs) model to predict the prognosis of patients with Grade II/III gliomas. Gene expression profile and clinical data of 459 patients with Grade II/III gliomas were extracted from The Cancer Genome Atlas. Then according to the immune/stromal scores generated by the ESTIMATE algorithm, the patients were scored one by one. Weighted gene co-expression network analysis (WGCNA) was used to construct a gene co-expression network to identify potential biomarkers for predicting the prognosis of patients. When adjusting clinical features including age, histology, grading, IDH status, we found that these features were independently associated with survival. The predicted value of the prognostic model was then verified in 440 samples in CGGA part B dataset and 182 samples in CGGA part C dataset by univariate and multivariate cox analysis. The clinical samples of 10 patients further confirmed our signature. Our findings suggested the eight-MRGs signature identified in this study are valuable prognostic predictors for patients with Grade II/III glioma.

Clinical practice guidelines for the management of adult diffuse gliomas

To follow the revision of the fourth edition of WHO classification and the recent progress on the management of diffuse gliomas, the joint guideline committee of Chinese Glioma Cooperative Group (CGCG), Society for Neuro-Oncology of China (SNO-China) and Chinese Brain Cancer Association (CBCA) updated the clinical practice guideline. It provides recommendations for diagnostic and management decisions, and for limiting unnecessary treatments and cost. The recommendations focus on molecular and pathological diagnostics, and the main treatment modalities of surgery, radiotherapy, and chemotherapy. In this guideline, we also integrated the results of some clinical trials of immune therapies and target therapies, which we think are ongoing future directions. The guideline should serve as an application for all professionals involved in the management of patients with adult diffuse glioma and also a source of knowledge for insurance companies and other institutions involved in the cost regulation of cancer care in China and other countries.

MicroRNA-1269a Promotes Proliferation and Arrest of Apoptosis of Glioma Cells by Directly Targeting ATRX

Glioma is one of the deadliest malignant brain tumors in adults worldwide. MicroRNA (miR) has been reported to be a pivotal regulator in human tumors. The aim of this study was to determine the expression, function, and mechanism of action of miR-1269a in glioma progression. The expression of miR-1269a was higher in both glioma cases reported in databases and glioma cell lines, and it was highly associated with poorer prognosis. Next, it was shown in vitro that mimic of miR-1269a could promote glioma progression and arrest apoptosis, whereas the inhibition of miR-1269a exhibited the opposite effects. In addition, miR-1269a was found to directly target ATRX chromatin remodeler by a dual-luciferase reporter assay. Moreover, ATRX overexpression could reverse the suppressive effects of miR-1269a on proliferation and apoptosis in vitro. In vivo subcutaneous xenograft tumor assay was also performed to confirm the phenotypes and molecular mechanism involved. Taking the findings together, our study implies that the miR-1269a/ATRX axis is a novel therapeutic target of glioma.

Molecular and Clinical Characterization of a Novel Prognostic and Immunologic Biomarker FAM111A in Diffuse Lower-Grade Glioma

Background

Family with sequence similarity 111 member A (FAM111A), as a replication factor required for proliferating cell nuclear antigen (PCNA) loading, has been demonstrated a possible association with carcinogenesis. However, the role of FAM111A in lower-grade glioma (LGG) remains unclear. We aim at investigating the expression and function of FAM111A in lower-grade glioma at the molecular and clinical levels.

Methods

In total, 711 lower-grade glioma samples were analyzed in our research, including 182 RNA-seq data from the Chinese Glioma Genome Atlas (CGGA) dataset and 529 RNA-seq data from The cancer Genome Atlas (TCGA) dataset. R language and the GraphPad software were used for the majority of statistical analysis and graphical work.

Results

FAM111A expression was overexpressed in WHO grade III and IDH-wildtype lower-grade glioma. FAM111A was significantly downregulated in the IDHmut-Codel molecular subtype. Univariate and multivariate Cox analysis demonstrated that FAM111A was an independent prognostic factor in LGG patients. Functional characterization of FAM111A revealed that it was associated with inflammatory response and immune response to tumor cells. FAM111A could also act as an indicator of the stromal and immune population, especially for monocytic lineage, myeloid dendritic cells and fibroblasts. It was positively correlated with macrophages, especially the M2 macrophage cells. Furthermore, FAM111A revealed predictive value for the immune subtypes and immune checkpoint blockade therapy.

Conclusion

FAM111A expression was closely related to the malignant phenotype, molecular pathology and immune response of lower-grade glioma. It might be a promising target for LGG immunotherapeutic strategies.

A machine learning analysis of a "normal-like" IDH-WT diffuse glioma transcriptomic subgroup associated with prolonged survival reveals novel immune and neurotransmitter-related actionable targets

BACKGROUND

Classification of primary central nervous system tumors according to the World Health Organization guidelines follows the integration of histologic interpretation with molecular information and aims at providing the most precise prognosis and optimal patient management. According to the cIMPACT-NOW update 3, diffuse isocitrate dehydrogenase-wild type (IDH-WT) gliomas should be graded as grade IV glioblastomas (GBM) if they possess one or more of the following molecular markers that predict aggressive clinical course: EGFR amplification, TERT promoter mutation, and whole-chromosome 7 gain combined with chromosome 10 loss.

METHODS

The Cancer Genome Atlas (TCGA) glioma expression datasets were reanalyzed in order to identify novel tumor subcategories which would be considered as GBM-equivalents with the current diagnostic algorithm. Unsupervised clustering allowed the identification of previously unrecognized transcriptomic subcategories. A supervised machine learning algorithm (k-nearest neighbor model) was also used to identify gene signatures specific to some of these subcategories.

RESULTS

We identified 14 IDH-WT infiltrating gliomas displaying a "normal-like" (NL) transcriptomic profile associated with a longer survival. Genes such as C5AR1 (complement receptor), SLC32A1 (vesicular gamma-aminobutyric acid transporter), MSR1 (or CD204, scavenger receptor A), and SYT5 (synaptotagmin 5) were differentially expressed and comprised in gene signatures specific to NL IDH-WT gliomas which were validated further using the Chinese Glioma Genome Atlas datasets. These gene signatures showed high discriminative power and correlation with survival.

CONCLUSION

NL IDH-WT gliomas represent an infiltrating glioma subcategory with a superior prognosis which can only be detected using genome-wide analysis. Differential expression of genes potentially involved in immune checkpoint and amino acid signaling pathways is providing insight into mechanisms of gliomagenesis and could pave the way to novel treatment targets for infiltrating gliomas.

Aging-related tumor associated fibroblasts changes could worsen the prognosis of GBM patients

Background

Glioblastoma multiforme (GBM) is the most malignant tumor in human brain, with highly heterogeneity among different patients. Age could function as an incidence and prognosis risk factor for many tumors.

Method

A series of bioinformatic experiments were conducted to evaluate the differences of incidence, differential expressed genes, enriched pathways with the data from Surveillance, Epidemiology, and End Results (SEER) program, the cancer genome atlas (TCGA) and Chinese glioma genome atlas (CGGA) project.

Results

We discovered in our present study that distinct difference of incidence and prognosis of different aged GBM patients. By a series of bioinformatic method, we found that the tumor associated fibroblasts (TAFs) was the most crucial tumor microenvironment (TME) component that led to this phenomenon. Epithelial-mesenchymal transition (EMT) could be the mechanism by which TAFs regulate the progression of GBM.

Conclusion

We have proposed a close correlation between age and GBM incidence and prognosis, and propose the underlying mechanism behind this correlation by mining different databases, which laid the foundation for future research.

Profiling pro-neural to mesenchymal transition identifies a lncRNA signature in glioma

Background

Molecular classification has laid the framework for exploring glioma biology and treatment strategies. Pro-neural to mesenchymal transition (PMT) of glioma is known to be associated with aggressive phenotypes, unfavorable prognosis, and treatment resistance. Recent studies have highlighted that long non-coding RNAs (lncRNAs) are key mediators in cancer mesenchymal transition. However, the relationship between lncRNAs and PMT in glioma has not been systematically investigated.

Methods

Gene expression profiles from The Cancer Genome Atlas (TCGA), the Chinese Glioma Genome Atlas (CGGA), GSE16011, and Rembrandt with available clinical and genomic information were used for analyses. Bioinformatics methods such as weighted gene co-expression network analysis (WGCNA), gene set enrichment analysis (GSEA), Cox analysis, and least absolute shrinkage and selection operator (LASSO) analysis were performed.

Results

According to PMT scores, we confirmed that PMT status was positively associated with risky behaviors and poor prognosis in glioma. The 149 PMT-related lncRNAs were identified by WGCNA analysis, among which 10 (LINC01057, TP73-AS1, AP000695.4, LINC01503, CRNDE, OSMR-AS1, SNHG18, AC145343.2, RP11-25K21.6, RP11-38L15.2) with significant prognostic value were further screened to construct a PMT-related lncRNA risk signature, which could divide cases into two groups with distinct prognoses. Multivariate Cox regression analyses indicated that the signature was an independent prognostic factor for high-grade glioma. High-risk cases were more likely to be classified as the mesenchymal subtype, which confers enhanced immunosuppressive status by recruiting macrophages, neutrophils, and regulatory T cells. Moreover, six lncRNAs of the signature could act as competing endogenous RNAs to promote PMT in glioblastoma.

Conclusions

We profiled PMT status in glioma and established a PMT-related 10-lncRNA signature for glioma that could independently predict glioma survival and trigger PMT, which enhanced immunosuppression.

Long-term efficacy of surgical resection with or without adjuvant therapy for treatment of secondary glioblastoma in adults

Background

There are limited studies on treatment strategies and associated clinical outcomes in patients with secondary glioblastoma (sGBM). We sought to investigate the prognostic factors and treatment decisions in a retrospective cohort of patients with sGBM.

Methods

One hundred and seventy-one patients with sGBM who met the screening criteria were included in this study. Kaplan–Meier survival analysis and Cox survival analysis were used to detect prognostic factors. R (v3.5.0) and SPSS software (v25.0, IBM) were used to perform statistical analyses.

Results

The median overall survival was 303 days (range 23–2237 days) and the median progression-free survival was 229 days (range 33–1964 days) in patients with sGBM. When assessing the relationship between adjuvant treatment outcome and extent of resection (EOR), the results showed that patients underwent gross total resection can benefit from postoperative radiotherapy and chemotherapy, but not in patients underwent subtotal resection. In addition, we also found that aggressive adjuvant therapy can significantly improve clinical outcomes of IDH1-mutated patients but no significant prognostic value for IDH1-wildtyped patients. The univariate Cox regression analyses demonstrated that EOR, adjuvant therapy, and postoperative Karnofsky Performance Scores were prognostic factors for patients with sGBM, and multivariate COX analysis confirmed that adjuvant therapy and EOR were independent prognostic factors.

Conclusions

For patients with sGBM, aggressive postoperative adjuvant therapy after gross total resection was recommended. However, we did not detect a benefit in IDH1-wildtype patients in our cohort.

Global RNA editome landscape discovers reduced RNA editing in glioma: loss of editing of gamma-amino butyric acid receptor alpha subunit 3 (GABRA3) favors glioma migration and invasion

Background

Gliomas are the most common and lethal type of intracranial tumors. With the current treatment regime, the median survival of patients with grade IV glioma (glioblastoma/GBM) remains at 14-16 months. RNA editing modifies the function and regulation of transcripts. The development of glial tumors may be caused by altered RNA editing events.

Methods

In this study, we uncover the global RNA editome landscape of glioma patients from RNA-seq data of control, lower grade glioma (LGG) and GBM samples (n = 1,083).

Results

A-to-I editing events were found to comprise 80% of the total editing events of which 96% were located in the Alu regions. The total RNA editing events were found to be reduced in glioma compared to control samples. More specifically, we found Gamma-aminobutyric acid type A receptor alpha3 (GABRA3) to be edited (c.1026 A-to-G; pI343M) in 73% (editing ratio 0.8) of control samples compared to LGG (28.96%; 0.47) and GBM (5.2%; 0.53) samples. GABRA3 transcript level was found to be downregulated in glioma compared to control in a grade-specific manner with GBMs having the lowest level of the transcript. Further, GABRA3 transcripts were observed to be higher in edited compared to unedited glioma samples. The transcript and protein levels of exogenously expressed gene were found to be higher for edited compared to unedited GABRA3 in glioma cells. Further, exogenously expressed edited GABRA3 inhibited migration and invasion of glioma cells efficiently but not the unedited GABRA3.

Conclusion

Collectively, our study discovered a reduction in RNA editing during glioma development. We further demonstrate that elevated RNA editing maintains a high level of GABRA3 RNA and protein in normal glial cells which provides a less migratory environment for the normal functioning of the brain. In contrast, the reduction in GABRA3 protein levels, due to lower stability of unedited RNA, results in the loss of function which confers an aggressive phenotype to GBM tumor.

A PRMT5-RNF168-SMURF2 Axis Controls H2AX Proteostasis

H2AX safeguards genomic stability in a dose-dependent manner; however, mechanisms governing its proteostasis are poorly understood. Here, we identify a PRMT5-RNF168-SMURF2 cascade that regulates H2AX proteostasis. We show that PRMT5 sustains the expression of RNF168, an E3 ubiquitin ligase essential for DNA damage response (DDR). Suppression of PRMT5 occurs in methylthioadenosine phosphorylase (MTAP)-deficient glioblastoma cells and attenuates the expression of RNF168, leading to destabilization of H2AX by E3 ubiquitin ligase SMURF2. RNF168 and SMURF2 serve as a stabilizer and destabilizer of H2AX, respectively, via their dynamic interactions with H2AX. In supporting an important role of this signaling cascade in regulating H2AX, MTAP-deficient glioblastoma cells display higher levels of DNA damage spontaneously or in response to genotoxic agents. These findings reveal a regulatory mechanism of H2AX proteostasis and define a signaling cascade that is essential to DDR and that is disrupted by the loss of a metabolic enzyme in tumor cells.