Establishment and Validation of a 5 m6A RNA Methylation Regulatory Gene Prognostic Model in Low-Grade Glioma

Development of stemness-related signature to optimize prognosis prediction and identify XMD8-85 as a novel therapeutic compound for glioma

Glioma is a highly invasive and aggressive type of brain cancer with poor treatment response. Stemness-related transcription factors form a regulatory network that sustains the malignant phenotype of gliomas. We conducted an integrated analysis of stemness-related transcription factors using The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) datasets, established the characteristics of stemness-related transcription factors, including Octamer-Binding Protein 4 (OCT4), Meis Homeobox 1 (MEIS1), E2F Transcription Factor 1 (E2F1), Transcription Factor CP2 Like 1 (TFCP2L1), and RUNX Family Transcription Factor 1 (RUNX1). The characteristic of stemness-related transcription factors was identified as an independent prognostic factor for glioma patients. Patients in the high-risk group have a worse prognosis than those in the low-risk group. The glioma microenvironment in the high-risk group exhibited a more active immune status. Single-cell level analysis revealed that stem cell-like cells exhibited stronger intercellular communication than glioma cells. Meanwhile, patients in different risk stratification exhibited varying sensitivities to immunotherapy and small molecule drug therapy. XMD8-85 was more effective in the high-risk group, and its antitumor effects were validated both in vivo and in vitro. Our results indicate that this prognostic feature will assist clinicians in predicting the prognosis of glioma patients, guiding immunotherapy and personalized treatment, as well as the potential clinical application of XMD8-85 in glioma treatment, and helping to develop effective treatment strategies.

Systematic integration of m6A regulators and autophagy-related genes in combination with long non-coding RNAs predicts survival in glioblastoma multiforme

Glioblastoma multiforme (GBM) is probably the only tumor in which a unique epigenetic alteration, namely methylation of the MGMT gene, possesses direct clinical relevance. Now with the emergence of aberrant N6 methyladenosine (m6A) modifications (the most common epigenetic modification of mRNA, closely linked to the autophagy process) in cancer, the epi-transcriptomic landscape of GBM pathobiology has been expanded. Considering this, herein, we systematically analyzed m6A regulators, assessed their correlation with autophagy-related genes (ATG), and established a long non-coding RNAs (lncRNA)-dependent prognostic signature (m6A-autophagy-lncRNAs) for GBM. Our analysis identified a novel signature of five long non-coding RNAs (lncRNAs: ITGA6-AS1, AC124248.1, NFYC-AS1, AC025171.1, and AC005229.3) associated with survival of GBM patients, and four among them clearly showed cancer-associated potential. We further validated and confirmed the altered expression of two lncRNAs (AC124248.1, AC005229.3) in GBM associated clinical samples using RT-PCR. Concerning the prognostic ability, the obtained signature determined high-/low-risk groups in GBM patients and showed sensitivity to anticancer drugs. Collectively, the m6A-autophagy-lncRNAs signature presented in the study is clinically relevant and is the first attempt to systematically predict the potential interaction between the three key determinants (m6A, autophagy, lncRNA) in cancer, particularly in GBM.

Understanding the Epitranscriptome for Avant-Garde Brain Tumour Diagnostics

RNA modifications are diverse, dynamic, and reversible transcript alterations rapidly gaining attention due to their newly defined RNA regulatory roles in cellular pathways and pathogenic mechanisms. The exciting emerging field of 'epitranscriptomics' is predominantly centred on studying the most abundant mRNA modification, N6-methyladenine (m⁶A). The m⁶A mark, similar to many other RNA modifications, is strictly regulated by so-called 'writer', 'reader', and 'eraser' protein species. The abundance of genes coding for the expression of these regulator proteins and m⁶A levels shows great potential as diagnostic and predictive tools across several cancer fields. This review explores our current understanding of RNA modifications in glioma biology and the potential of epitranscriptomics to develop new diagnostic and predictive classification tools that can stratify these highly complex and heterogeneous brain tumours.

Establishment of three heterogeneous subtypes and a risk model of low-grade gliomas based on cell senescence-related genes

Background

Cellular senescence is a key element in the occurrence and progression of a variety of tumors. As a result, cellular senescence-related markers can be categorized based on the prognosis status of patients. Due to the heterogeneity and the complexity of the tumor microenvironment (TME), the long-term effectiveness of low-grade glioma (LGG) treatment remains a clinical challenge. Consequently, developing and refining effective treatment approaches to aid with LGG management is critical.

Methods

Based on the expressions of cell senescence-related genes (CSRGs) acquired from the cellAge database, consensus clustering was utilized to identify stable molecular subtypes. Clinical features, immune infiltration, route modifications, and genetic changes of various subtypes were also assessed. Following that, the least absolute shrinkage and selection operator (LASSO) regression and univariate Cox regression analysis were used for developing the cell senescence-related risk score (CSRS) model. Finally, a correlation study of the CSRS model with molecular, immunological, and immunotherapy parameters was performed.

Results

C1, C2, and C3, are the three senescence-related subtypes that were identified. Patients belonging to the C1 subtype had poor prognoses and a substantial proportion of them was in the grade G3. The differentially expressed genes (DEGs) among the three subtypes were used to develop the CSRS model. In both the training and independent validation cohort, the model had a high area under the receiver operating characteristic (ROC) curve in predicting the overall survival (OS) of patients. As a result, this model can predict clinical features and responses to immunotherapy in a variety of patients and it is a potential independent prognostic factor for LGG.

Conclusion

This research discovered three LGG subtypes related to cell senescence and created a CSRS model for six genes. Cell senescence was highly associated with unfavorable prognosis in LGG. The CSRS model can be used to predict the prognosis of patients and identify patients who would benefit from immunotherapy.

Identification of a novel cuproptosis-related gene signature and integrative analyses in patients with lower-grade gliomas

Background

Cuproptosis is a newly discovered unique non-apoptotic programmed cell death distinguished from known death mechanisms like ferroptosis, pyroptosis, and necroptosis. However, the prognostic value of cuproptosis and the correlation between cuproptosis and the tumor microenvironment (TME) in lower-grade gliomas (LGGs) remain unknown.

Methods

In this study, we systematically investigated the genetic and transcriptional variation, prognostic value, and expression patterns of cuproptosis-related genes (CRGs). The CRG score was applied to quantify the cuproptosis subtypes. We then evaluated their values in the TME, prognostic prediction, and therapeutic responses in LGG. Lastly, we collected five paired LGG and matched normal adjacent tissue samples from Sun Yat-sen University Cancer Center (SYSUCC) to verify the expression of signature genes by quantitative real-time PCR (qRT-PCR) and Western blotting (WB).

Results

Two distinct cuproptosis-related clusters were identified using consensus unsupervised clustering analysis. The correlation between multilayer CRG alterations with clinical characteristics, prognosis, and TME cell infiltration were observed. Then, a well-performed cuproptosis-related risk model (CRG score) was developed to predict LGG patients' prognosis, which was evaluated and validated in two external cohorts. We classified patients into high- and low-risk groups according to the CRG score and found that patients in the low-risk group showed significantly higher survival possibilities than those in the high-risk group (P<0.001). A high CRG score implies higher TME scores, more significant TME cell infiltration, and increased mutation burden. Meanwhile, the CRG score was significantly correlated with the cancer stem cell index, chemoradiotherapy sensitivity-related genes and immune checkpoint genes, and chemotherapeutic sensitivity, indicating the association with CRGs and treatment responses. Univariate and multivariate Cox regression analyses revealed that the CRG score was an independent prognostic predictor for LGG patients. Subsequently, a highly accurate predictive model was established for facilitating the clinical application of the CRG score, showing good predictive ability and calibration. Additionally, crucial CRGs were further validated by qRT-PCR and WB.

Conclusion

Collectively, we demonstrated a comprehensive overview of CRG profiles in LGG and established a novel risk model for LGG patients' therapy status and prognosis. Our findings highlight the potential clinical implications of CRGs, suggesting that cuproptosis may be the potential therapeutic target for patients with LGG.

Recurrence- and Malignant Progression-Associated Biomarkers in Low-Grade Gliomas and Their Roles in Immunotherapy

Despite a generally better prognosis than high-grade glioma (HGG), recurrence and malignant progression are the main causes for the poor prognosis and difficulties in the treatment of low-grade glioma (LGG). It is of great importance to learn about the risk factors and underlying mechanisms of LGG recurrence and progression. In this study, the transcriptome characteristics of four groups, namely, normal brain tissue and recurrent LGG (rLGG), normal brain tissue and secondary glioblastoma (sGBM), primary LGG (pLGG) and rLGG, and pLGG and sGBM, were compared using Chinese Glioma Genome Atlas (CGGA) and Genotype-Tissue Expression Project (GTEx) databases. In this study, 296 downregulated and 396 upregulated differentially expressed genes (DEGs) with high consensus were screened out. Univariate Cox regression analysis of data from The Cancer Genome Atlas (TCGA) yielded 86 prognostically relevant DEGs; a prognostic prediction model based on five key genes (HOXA1, KIF18A, FAM133A, HGF, and MN1) was established using the least absolute shrinkage and selection operator (LASSO) regression dimensionality reduction and multivariate Cox regression analysis. LGG was divided into high- and low-risk groups using this prediction model. Gene Set Enrichment Analysis (GSEA) revealed that signaling pathway differences in the high- and low-risk groups were mainly seen in tumor immune regulation and DNA damage-related cell cycle checkpoints. Furthermore, the infiltration of immune cells in the high- and low-risk groups was analyzed, which indicated a stronger infiltration of immune cells in the high-risk group than that in the low-risk group, suggesting that an immune microenvironment more conducive to tumor growth emerged due to the interaction between tumor and immune cells. The tumor mutational burden and tumor methylation burden in the high- and low-risk groups were also analyzed, which indicated higher gene mutation burden and lower DNA methylation level in the high-risk group, suggesting that with the accumulation of genomic mutations and epigenetic changes, tumor cells continued to evolve and led to the progression of LGG to HGG. Finally, the value of potential therapeutic targets for the five key genes was analyzed, and findings demonstrated that KIF18A was the gene most likely to be a potential therapeutic target. In conclusion, the prediction model based on these five key genes can better identify the high- and low-risk groups of LGG and lay a solid foundation for evaluating the risk of LGG recurrence and malignant progression.

A Novel CRISPR/Cas9 Screening Potential Index for Prognostic and Immunological Prediction in Low-Grade Glioma

Glioma is a malignancy with the highest mortality in central nervous system disorders. Here, we implemented the computational tools based on CRISPR/Cas9 to predict the clinical outcomes and biological characteristics of low-grade glioma (LGG). The transcriptional expression profiles and clinical phenotypes of LGG patients were retrieved from The Cancer Genome Atlas and Chinese Glioma Genome Atlas. The CERES algorithm was used to screen for LGG-lethal genes. Cox regression and random survival forest were adopted for survival-related gene selection. Nonnegative matrix factorization distinguished patients into different clusters. Single-sample gene set enrichment analysis was employed to create a novel CRISPR/Cas9 screening potential index (CCSPI), and patients were stratified into low- and high-CCSPI groups. Survival analysis, area under the curve values (AUCs), nomogram, and tumor microenvironment exploration were included for the model validation. A total of 20 essential genes in LGG were used to classify patients into two clusters and construct the CCSPI system. High-CCSPI patients were associated with a worse prognosis of both training and validation set (p < 0.0001) and higher immune fractions than low-CCSPI individuals. The CCSPI system had a promising performance with 1-, 3-, and 5-year AUCs of 0.816, 0.779, 0.724, respectively, and the C-index of the nomogram model reached 0.743 (95% CI = 0.725–0.760). Immune-infiltrating cells and immune checkpoints such as PD-1/PD-L1 and POLD3 were positively associated with CCSPI. In conclusion, the CCSPI had prognostic value in LGG, and the model will deepen our cognition of the interaction between the CNS and immune system in different LGG subtypes.