Combined Therapy Sensitivity Index Based on a 13-Gene Signature Predicts Prognosis for IDH Wild-type and MGMT Promoter Unmethylated Glioblastoma Patients

mRNA markers for survival prediction in glioblastoma multiforme patients: a systematic review with bioinformatic analyses

BACKGROUND

Glioblastoma multiforme (GBM) is a type of fast-growing brain glioma associated with a very poor prognosis. This study aims to identify key genes whose expression is associated with the overall survival (OS) in patients with GBM.

METHODS

A systematic review was performed using PubMed, Scopus, Cochrane, and Web of Science up to Journey 2024. Two researchers independently extracted the data and assessed the study quality according to the New Castle Ottawa scale (NOS). The genes whose expression was found to be associated with survival were identified and considered in a subsequent bioinformatic study. The products of these genes were also analyzed considering protein-protein interaction (PPI) relationship analysis using STRING. Additionally, the most important genes associated with GBM patients' survival were also identified using the Cytoscape 3.9.0 software. For final validation, GEPIA and CGGA (mRNAseq_325 and mRNAseq_693) databases were used to conduct OS analyses. Gene set enrichment analysis was performed with GO Biological Process 2023.

RESULTS

From an initial search of 4104 articles, 255 studies were included from 24 countries. Studies described 613 unique genes whose mRNAs were significantly associated with OS in GBM patients, of which 107 were described in 2 or more studies. Based on the NOS, 131 studies were of high quality, while 124 were considered as low-quality studies. According to the PPI network, 31 key target genes were identified. Pathway analysis revealed five hub genes (IL6, NOTCH1, TGFB1, EGFR, and KDR). However, in the validation study, only, the FN1 gene was significant in three cohorts.

CONCLUSION

We successfully identified the most important 31 genes whose products may be considered as potential prognosis biomarkers as well as candidate target genes for innovative therapy of GBM tumors.

Genome-wide DNA methylation analysis identifies potent CpG signature for temzolomide response in non-G-CIMP glioblastomas with unmethylated MGMT promoter: MGMT-dependent roles of GPR81

PURPOSES

To identify potent DNA methylation candidates that could predict response to temozolomide (TMZ) in glioblastomas (GBMs) that do not have glioma-CpGs island methylator phenotype (G-CIMP) but have an unmethylated promoter of O-6-methylguanine-DNA methyltransferase (unMGMT).

METHODS

The discovery-validation approach was planned incorporating a series of G-CIMP-/unMGMT GBM cohorts with DNA methylation microarray data and clinical information, to construct multi-CpG prediction models. Different bioinformatic and experimental analyses were performed for biological exploration.

RESULTS

By analyzing discovery sets with radiotherapy (RT) plus TMZ versus RT alone, we identified a panel of 64 TMZ efficacy-related CpGs, from which a 10-CpG risk signature was further constructed. Both the 64-CpG panel and the 10-CpG risk signature were validated showing significant correlations with overall survival of G-CIMP-/unMGMT GBMs when treated with RT/TMZ, rather than RT alone. The 10-CpG risk signature was further observed for aiding TMZ choice by distinguishing differential outcomes to RT/TMZ versus RT within each risk subgroup. Functional studies on GPR81, the gene harboring one of the 10 CpGs, indicated its distinct impacts on TMZ resistance in GBM cells, which may be dependent on the status of MGMT expression.

CONCLUSIONS

The 64 TMZ efficacy-related CpGs and in particular the 10-CpG risk signature may serve as promising predictive biomarker candidates for guiding optimal usage of TMZ in G-CIMP-/unMGMT GBMs.

COVPRIG robustly predicts the overall survival of IDH wild-type glioblastoma and highlights METTL1+ neural-progenitor-like tumor cell in driving unfavorable outcome

BACKGROUND

Accurately predicting the outcome of isocitrate dehydrogenase (IDH) wild-type glioblastoma (GBM) remains hitherto challenging. This study aims to Construct and Validate a Robust Prognostic Model for IDH wild-type GBM (COVPRIG) for the prediction of overall survival using a novel metric, gene-gene (G × G) interaction, and explore molecular and cellular underpinnings.

METHODS

Univariate and multivariate Cox regression of four independent trans-ethnic cohorts containing a total of 800 samples. Prediction efficacy was comprehensively evaluated and compared with previous models by a systematic literature review. The molecular underpinnings of COVPRIG were elucidated by integrated analysis of bulk-tumor and single-cell based datasets.

RESULTS

Using a Cox-ph model-based method, six of the 93,961 G × G interactions were screened to form an optimal combination which, together with age, comprised the COVPRIG model. COVPRIG was designed for RNA-seq and microarray, respectively, and effectively identified patients at high risk of mortality. The predictive performance of COVPRIG was satisfactory, with area under the curve (AUC) ranging from 0.56 (CGGA693, RNA-seq, 6-month survival) to 0.79 (TCGA RNAseq, 18-month survival), which can be further validated by decision curves. Nomograms were constructed for individual risk prediction for RNA-seq and microarray-based cohorts, respectively. Besides, the prognostic significance of COVPRIG was also validated in GBM including the IDH mutant samples. Notably, COVPRIG was comprehensively evaluated and externally validated, and a systemic review disclosed that COVPRIG outperformed current validated models with an integrated discrimination improvement (IDI) of 6-16%. Moreover, integrative bioinformatics analysis predicted an essential role of METTL1+ neural-progenitor-like (NPC-like) malignant cell in driving unfavorable outcome.

CONCLUSION

This study provided a powerful tool for the outcome prediction for IDH wild-type GBM, and preliminary molecular underpinnings for future research.

Identification and Potential Mechanisms of a 7-MicroRNA Signature That Predicts Prognosis in Patients with Lower-Grade Glioma

Background

Lower-grade glioma is an intracranial cancer that may develop into glioblastoma with high mortality. The main objective of our study is to develop microRNA for LGG patients which will provide novel prognostic biomarkers along with therapeutic targets.

Methods

Clinicopathological data of LGG patients and their RNA expression profile were downloaded through The Cancer Genome Atlas Relevant expression profiles of RNA, and clinicopathological data of the LGG patients had been extracted from the database of "The Cancer Genome Atlas." Differential expression analysis had been conducted for identification of the differentially expressed microRNAs as well as mRNAs in LGG samples and normal ones. ROC curves and K-M plots were plotted to confirm performance and for predictive accuracy. For the confirmation of microRNAs as an independent prognostic factor, an independent prognosis analysis was conducted. Moreover, target differentially expressed genes of these identified prognostic microRNAs that were extracted and protein-protein interaction networks were developed. Moreover, the biological functions of signature were determined through Genome Ontology analysis, genome pathway analysis, and Kyoto Encyclopedia of Genes.

Results

7-microRNA signature was identified that has the ability of categorization of individuals with LGG into high- and low-risk groups on the basis of significant difference in survival during training and testing cohorts (P < 0.001). The 7-microRNA signature had appeared to be robust in predictive accuracy (all AUC> 0.65). It was also approved with multivariate Cox regression along with some traditional clinical practices that we can use 7-microRNA signature for therapeutic purposes as a self-regulating predictive OS factor (P < 0.001). KEGG and Gene Ontology (GO) analyses reported that 7-microRNAs had mainly developed in important pathways related with glioma, e.g., the "cAMP signaling pathway," "glutamatergic synapses," and "calcium signaling pathway".

Conclusion

A newly discovered 7-microRNA signature could be a potential target for the diagnosis and treatment for LGG patients.

Promoting Prognostic Model Application: A Review Based on Gliomas

Malignant neoplasms are characterized by poor therapeutic efficacy, high recurrence rate, and extensive metastasis, leading to short survival. Previous methods for grouping prognostic risks are based on anatomic, clinical, and pathological features that exhibit lower distinguishing capability compared with genetic signatures. The update of sequencing techniques and machine learning promotes the genetic panels-based prognostic model development, especially the RNA-panel models. Gliomas harbor the most malignant features and the poorest survival among all tumors. Currently, numerous glioma prognostic models have been reported. We systematically reviewed all 138 machine-learning-based genetic models and proposed novel criteria in assessing their quality. Besides, the biological and clinical significance of some highly overlapped glioma markers in these models were discussed. This study screened out markers with strong prognostic potential and 27 models presenting high quality. Conclusively, we comprehensively reviewed 138 prognostic models combined with glioma genetic panels and presented novel criteria for the development and assessment of clinically important prognostic models. This will guide the genetic models in cancers from laboratory-based research studies to clinical applications and improve glioma patient prognostic management.

Promoting Prognostic Model Application: A Review Based on Gliomas

Malignant neoplasms are characterized by poor therapeutic efficacy, high recurrence rate, and extensive metastasis, leading to short survival. Previous methods for grouping prognostic risks are based on anatomic, clinical, and pathological features that exhibit lower distinguishing capability compared with genetic signatures. The update of sequencing techniques and machine learning promotes the genetic panels-based prognostic model development, especially the RNA-panel models. Gliomas harbor the most malignant features and the poorest survival among all tumors. Currently, numerous glioma prognostic models have been reported. We systematically reviewed all 138 machine-learning-based genetic models and proposed novel criteria in assessing their quality. Besides, the biological and clinical significance of some highly overlapped glioma markers in these models were discussed. This study screened out markers with strong prognostic potential and 27 models presenting high quality. Conclusively, we comprehensively reviewed 138 prognostic models combined with glioma genetic panels and presented novel criteria for the development and assessment of clinically important prognostic models. This will guide the genetic models in cancers from laboratory-based research studies to clinical applications and improve glioma patient prognostic management.

Prognostic value and immune cell infiltration of hypoxic phenotype-related gene signatures in glioblastoma microenvironment

Glioblastoma (GBM) is a malignant intracranial tumour with the highest proportion and lethality. It is characterized by invasiveness and heterogeneity. However, the currently available therapies are not curative. As an essential environmental cue that maintains glioma stem cells, hypoxia is considered the cause of tumour resistance to chemotherapy and radiation. Growing evidence shows that immunotherapy focusing on the tumour microenvironment is an effective treatment for GBM; however, the current clinicopathological features cannot predict the response to immunotherapy and provide accurate guidance for immunotherapy. Based on the ESTIMATE algorithm, GBM cases of The Cancer Genome Atlas (TCGA) data set were classified into high- and low-immune/stromal score groups, and a four-gene tumour environment-related model was constructed. This model exhibited good efficiency at forecasting short- and long-term prognosis and could also act as an independent prognostic biomarker. Additionally, this model and four of its genes (CLECL5A, SERPING1, CHI3L1 and C1R) were found to be associated with immune cell infiltration, and further study demonstrated that these four genes might drive the hypoxic phenotype of perinecrotic GBM, which affects hypoxia-induced glioma stemness. Therefore, these might be important candidates for immunotherapy of GBM and deserve further exploration.

Expression profile analysis to predict potential biomarkers for glaucoma: BMP1, DMD and GEM

Purpose

Glaucoma is the second commonest cause of blindness. We assessed the gene expression profile of astrocytes in the optic nerve head to identify possible prognostic biomarkers for glaucoma.

Method

A total of 20 patient and nine normal control subject samples were derived from the GSE9944 (six normal samples and 13 patient samples) and GSE2378 (three normal samples and seven patient samples) datasets, screened by microarray-tested optic nerve head tissues, were obtained from the Gene Expression Omnibus (GEO) database. We used a weighted gene coexpression network analysis (WGCNA) to identify coexpressed gene modules. We also performed a functional enrichment analysis and least absolute shrinkage and selection operator (LASSO) regression analysis. Genes expression was represented by boxplots, functional geneset enrichment analyses (GSEA) were used to profile the expression patterns of all the key genes. Then the key genes were validated by the external dataset.

Results

A total 8,606 genes and 19 human optic nerve head samples taken from glaucoma patients in the GSE9944 were compared with normal control samples to construct the co-expression gene modules. After selecting the most common clinical traits of glaucoma, their association with gene expression was established, which sorted two modules showing greatest correlations. One with the correlation coefficient is 0.56 (P = 0.01) and the other with the correlation coefficient is −0.56 (P = 0.01). Hub genes of these modules were identified using scatterplots of gene significance versus module membership. A functional enrichment analysis showed that the former module was mainly enriched in genes involved in cellular inflammation and injury, whereas the latter was mainly enriched in genes involved in tissue homeostasis and physiological processes. This suggests that genes in the green–yellow module may play critical roles in the onset and development of glaucoma. A LASSO regression analysis identified three hub genes: Recombinant Bone Morphogenetic Protein 1 gene (BMP1), Duchenne muscular dystrophy gene (DMD) and mitogens induced GTP-binding protein gene (GEM). The expression levels of the three genes in the glaucoma group were significantly lower than those in the normal group. GSEA further illuminated that BMP1, DMD and GEM participated in the occurrence and development of some important metabolic progresses. Using the GSE2378 dataset, we confirmed the high validity of the model, with an area under the receiver operator characteristic curve of 85%.

Conclusion

We identified several key genes, including BMP1, DMD and GEM, that may be involved in the pathogenesis of glaucoma. Our results may help to determine the prognosis of glaucoma and/or to design gene- or molecule-targeted drugs.

The interaction between TERT promoter mutation and MGMT promoter methylation on overall survival of glioma patients: a meta-analysis

Background

There are controversial results concerning the prognostic implication of TERT promoter mutation in glioma patients concerning MGMT status. In this meta-analysis, we investigated whether there are any interactions of these two genetic markers on the overall survival (OS) of glioma patients.

Methods

Electronic databases including PubMed and Web of Science were searched for relevant studies. Hazard ratio (HR) and its 95% confidence interval (CI) for OS adjusted for selected covariates were calculated from the individual patient data (IPD), Kaplan-Meier curve (KMC), or directly obtained from the included studies.

Results

A total of nine studies comprising 2819 glioma patients were included for meta-analysis. Our results showed that TERT promoter mutation was associated with a superior outcome in MGMT-methylated gliomas (HR = 0.73; 95% CI = 0.55–0.98; p-value = 0.04), whereas this mutation was associated with poorer survival in gliomas without MGMT methylation (HR = 1.86; 95% CI = 1.54–2.26; p-value < 0.001). TERT-mutated glioblastoma (GBM) patients with MGMT methylation benefited from temozolomide (TMZ) treatment (HR = 0.33; 95% CI = 0.23–0.47; p-value < 0.001). MGMT methylation was not related with any improvement in OS in TERT-wild type GBMs (HR = 0.80; 95% CI = 0.56–1.15; p-value = 0.23).

Conclusions

The prognostic value of TERT promoter mutation may be modulated by MGMT methylation status. Not all MGMT-methylated GBM patients may benefit from TMZ; it is possible that only TERT-mutated GBM with MGMT methylation, in particular, may respond.

Integrating Genomic Data with Transcriptomic Data for Improved Survival Prediction for Adult Diffuse Glioma

Background: Glioma is the most common type of primary central nervous system tumors. However, the relationship between gene mutations and transcriptome is unclear in diffuse glioma, and there are no systemic analyses with regard to the genotype-phenotype association currently.

Methods: We performed the multi-omics analysis in large glioblastoma multiforme (GBM, n=126) and low-grade glioma (LGG, n=481) cohorts obtained from The Cancer Genome Atlas (TCGA) database. We used multivariate linear models to evaluate associations between driver gene mutations and global gene expression. We developed generalized linear models to evaluate associations between genetic/expression factors with clinicopathologic features. Multivariate Cox proportional hazards models were used to predict the overall survival.

Results: The potential relationship between genotype and genetics, clinical as well as pathologic features, on diffused glioma was observed. At least one driver mutation correlated with expression changes of about 10% of genes in GBMs while about 80% of genes in LGGs. The strongest association between mutations and expression changes was observed for DRG2 and LRCC41 gene in GBMs and LGGs, respectively. Additionally, the association between genomics features and clinicopathologic features suggested the different underlying molecular mechanisms in molecular subtypes or histology subtypes. For predicting survival, among genetics, transcriptome and clinical variables, transcriptome features made the largest contribution. By combining all the available data, the accuracy in predicting the prognosis of diffuse glioma in patients was also improved.

Conclusion: Our study results revealed the influences of driver gene mutations on global gene expression in diffuse glioma patients. A more accurate model in predicting the prognosis of patients was achieved when combining with all the available data than just transcriptomic data.

Identification and validation of a three-gene signature as a candidate prognostic biomarker for lower grade glioma

Background

Lower grade glioma (LGG) are a heterogeneous tumor that may develop into high-grade malignant glioma seriously shortens patient survival time. The clinical prognostic biomarker of lower-grade glioma is still lacking. The aim of our study is to explore novel biomarkers for LGG that contribute to distinguish potential malignancy in low-grade glioma, to guide clinical adoption of more rational and effective treatments.

Methods

The RNA-seq data for LGG was downloaded from UCSC Xena and the Chinese Glioma Genome Atlas (CGGA). By a robust likelihood-based survival model, least absolute shrinkage and selection operator regression and multivariate Cox regression analysis, we developed a three-gene signature and established a risk score to predict the prognosis of patient with LGG. The three-gene signature was an independent survival predictor compared to other clinical parameters. Based on the signature related risk score system, stratified survival analysis was performed in patients with different age group, gender and pathologic grade. The prognostic signature was validated in the CGGA dataset. Finally, weighted gene co-expression network analysis (WGCNA) was carried out to find the co-expression genes related to the member of the signature and enrichment analysis of the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway were conducted for those co-expression network. To prove the efficiency of the model, time-dependent receiver operating characteristic curves of our model and other models are constructed.

Results

In this study, a three-gene signature (WEE1, CRTAC1, SEMA4G) was constructed. Based on the model, the risk score of each patient was calculated with LGG (low-risk vs. high-risk, hazard ratio (HR) = 0.198 (95% CI [0.120-0.325])) and patients in the high-risk group had significantly poorer survival results than those in the low-risk group. Furthermore, the model was validated in the CGGA dataset. Lastly, by WGCNA, we constructed the co-expression network of the three genes and conducted the enrichment of GO and KEGG. Our study identified a three-gene model that showed satisfactory performance in predicting the 1-, 3- and 5-year survival of LGG patients compared to other models and may be a promising independent biomarker of LGG.