Gene Expression Profiling Stratifies IDH-Wildtype Glioblastoma With Distinct Prognoses

The Three Pillars of Glioblastoma: A Systematic Review and Novel Analysis of Multi-Omics and Clinical Data

Glioblastoma is the most fatal and common malignant brain tumor, excluding metastasis and with a median survival of approximately one year. While solid tumors benefit from newly approved drugs, immunotherapy, and prevention, none of these scenarios are opening for glioblastoma. The key to unlocking the peculiar features of glioblastoma is observing its molecular and anatomical features tightly entangled with the host's central nervous system (CNS). In June 2024, we searched the PUBMED electronic database. Data collection and analysis were conducted independently by two reviewers. Results: A total of 215 articles were identified, and 192 were excluded based on inclusion and exclusion criteria. The remaining 23 were used for collecting divergent molecular pathways and anatomical features of glioblastoma. The analysis of the selected papers revealed a multifaced tumor with extreme variability and cellular reprogramming that are observable within the same patient. All the variability of glioblastoma could be clustered into three pillars to dissect the physiology of the tumor: 1. necrotic core; 2. vascular proliferation; 3. CNS infiltration. These three pillars support glioblastoma survival, with a pivotal role of the neurovascular unit, as supported by the most recent paper published by experts in the field.

Identification of MDK as a Hypoxia- and Epithelial-Mesenchymal Transition-Related Gene Biomarker of Glioblastoma Based on a Novel Risk Model and In Vitro Experiments

BACKGROUND

Tumor cells are commonly exposed to a hypoxic environment, which can easily induce the epithelial-mesenchymal transition (EMT) of tumor cells, further affecting tumor proliferation, invasion, metastasis, and drug resistance. However, the predictive role of hypoxia and EMT-related genes in glioblastoma (GBM) has not been investigated.

METHODS

Intersection genes were identified by weighted correlation network analysis (WGCNA) and differential expression analyses, and a risk model was further constructed by LASSO and Cox analyses. Clinical, immune infiltration, tumor mutation, drug treatment, and enrichment profiles were analyzed based on the risk model. The expression level of the MDK gene was tested using RT-PCR, immunohistochemistry, and immunofluorescence. CCK8 and EdU were employed to determine the GBM cells' capacity for proliferation while the migration and invasion ability were detected by a wound healing assay and transwell assay, respectively.

RESULTS

Based on the GBM data of the TCGA and GTEx databases, 58 intersection genes were identified, and a risk model was constructed. The model was verified in the CGGA cohort, and its accuracy was confirmed by the ROC curve (AUC = 0.807). After combining clinical subgroups, univariate and multivariate Cox regression analyses showed that risk score and age were independent risk factors for GBM patients. Furthermore, our subsequent analysis of immune infiltration, tumor mutation, and drug treatment showed that risk score and high- and low-risk groups were associated with multiple immune cells, mutated genes, and drugs. Enrichment analysis indicated that the differences between high- and low-risk groups were manifested in tumor-related pathways, including the PI3K-AKT and JAK-STAT pathways. Finally, in vivo experiments proved that the hypoxia environment promoted the expression of MDK, and MDK knockdown reduced the proliferation, migration, and EMT of GBM cells induced by hypoxia.

CONCLUSIONS

Our novel prognostic correlation model provided more potential treatment strategies for GBM patients.

Multi-Modal Learning for Predicting the Genotype of Glioma

The isocitrate dehydrogenase (IDH) gene mutation is an essential biomarker for the diagnosis and prognosis of glioma. It is promising to better predict glioma genotype by integrating focal tumor image and geometric features with brain network features derived from MRI. Convolutional neural networks show reasonable performance in predicting IDH mutation, which, however, cannot learn from non-Euclidean data, e.g., geometric and network data. In this study, we propose a multi-modal learning framework using three separate encoders to extract features of focal tumor image, tumor geometrics and global brain networks. To mitigate the limited availability of diffusion MRI, we develop a self-supervised approach to generate brain networks from anatomical multi-sequence MRI. Moreover, to extract tumor-related features from the brain network, we design a hierarchical attention module for the brain network encoder. Further, we design a bi-level multi-modal contrastive loss to align the multi-modal features and tackle the domain gap at the focal tumor and global brain. Finally, we propose a weighted population graph to integrate the multi-modal features for genotype prediction. Experimental results on the testing set show that the proposed model outperforms the baseline deep learning models. The ablation experiments validate the performance of different components of the framework. The visualized interpretation corresponds to clinical knowledge with further validation. In conclusion, the proposed learning framework provides a novel approach for predicting the genotype of glioma.

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.

Feature selection combined with top-down and bottom-up strategies for survival analysis: A case of prognostic prediction in glioblastoma

Over the last decades, molecular signatures have attracted extensive attention in cancer research. However, most of the reported biomarkers show a weak distinguishing ability in predicting the survival risks of patients. Actually, univariate analysis is generally considered in regression analysis, which makes the existing statistical methods ineffective. Furthermore, there is too much human involvement in the ways of classifying patients with high and low risk. Last but not least, the participation of therapy after conservative surgery also makes the survival analysis more complex. In order to solve these problems, we propose a solid method of feature selection which combines top-down and bottom-up strategies. The top-down strategy is to randomly extract some genes each time and select candidate genes through cumulative voting. The bottom-up strategy is to fully enumerate the selected genes and to use a clustering algorithm to classify samples. We analyzed glioblastoma data from the Cancer Genome Atlas (TCGA) and got candidate signatures. The results of simulation data, as well as an independent test set the Chinese Glioma Genome Atlas (CGGA), verified the reliability of the method and validity of the selected features.

Expression of hub genes of endothelial cells in glioblastoma-A prognostic model for GBM patients integrating single-cell RNA sequencing and bulk RNA sequencing

BACKGROUND

This study aimed to use single-cell RNA-seq (scRNA-seq) to discover marker genes in endothelial cells (ECs) and construct a prognostic model for glioblastoma multiforme (GBM) patients in combination with traditional high-throughput RNA sequencing (bulk RNA-seq).

METHODS

Bulk RNA-seq data was downloaded from The Cancer Genome Atlas (TCGA) and The China Glioma Genome Atlas (CGGA) databases. 10x scRNA-seq data for GBM were obtained from the Gene Expression Omnibus (GEO) database. The uniform manifold approximation and projection (UMAP) were used for downscaling and cluster identification. Key modules and differentially expressed genes (DEGs) were identified by weighted gene correlation network analysis (WGCNA). A non-negative matrix decomposition (NMF) algorithm was used to identify the different subtypes based on DEGs, and multivariate cox regression analysis to model the prognosis. Finally, differences in mutational landscape, immune cell abundance, immune checkpoint inhibitors (ICIs)-associated genes, immunotherapy effects, and enriched pathways were investigated between different risk groups.

RESULTS

The analysis of scRNA-seq data from eight samples revealed 13 clusters and four cell types. After applying Fisher's exact test, ECs were identified as the most important cell type. The NMF algorithm identified two clusters with different prognostic and immunological features based on DEGs. We finally built a prognostic model based on the expression levels of four key genes. Higher risk scores were significantly associated with poorer survival outcomes, low mutation rates in IDH genes, and upregulation of immune checkpoints such as PD-L1 and CD276.

CONCLUSION

We built and validated a 4-gene signature for GBM using 10 scRNA-seq and bulk RNA-seq data in this work.

Aggressiveness of Grade 4 Gliomas of Adults

Grade 4 adult gliomas are IDH-mutant astrocytomas and IDH-wildtype glioblastomas. They have a very high mortality rate, with survival at 5 years not exceeding 5%. We aimed to conduct a clinical imaging and morphogenetic characterization of them, as well as to identify the main negative prognostic factors that give them such aggressiveness. We conducted a ten-year retrospective study. We followed the clinical, imaging, and morphogenetic aspects of the cases. We analyzed immunohistochemical markers (IDH1, Ki-67, and nestin) and FISH tests based on the CDKN2A gene. The obtained results were analyzed using SPSS Statistics with the appropriate parameters. The clinical aspects representing negative prognostic factors were represented by patients’ comorbidities: hypertension (HR = 1.776) and diabetes mellitus/hyperglycemia (HR = 2.159). The lesions were mostly supratentorial, and the temporal lobe was the most affected. The mean volume was 88.05 cm3 and produced a midline shift with an average of 8.52 mm. Subtotal surgical resection was a negative prognostic factor (HR = 1.877). The proliferative index did not influence survival rate, whereas CDKN2A gene mutations were shown to have a major impact on survival. We identified the main negative prognostic factors that support the aggressiveness of grade 4 gliomas: patient comorbidities, type of surgical resection, degree of cell differentiation, and CDKN2A gene mutations.

An Immune-Related Prognostic Signature for Predicting Clinical Outcomes and Immune Landscape in IDH-Mutant Lower-Grade Gliomas

BACKGROUND

IDH mutation is the most common in diffuse LGGs, correlated with a favorable prognosis. However, the IDH-mutant LGGs patients with poor prognoses need to be identified, and the potential mechanism leading to a worse outcome and treatment options needs to be investigated.

METHODS

A six-gene immune-related prognostic signature in IDH-mutant LGGs was constructed based on two public datasets and univariate, multivariate, and LASSO Cox regression analysis. Patients were divided into low- and high-risk groups based on the median risk score in the training and validation sets. We analyzed enriched pathways and immune cell infiltration, applying the GSEA and the immune evaluation algorithms.

RESULTS

Stratification and multivariate Cox analysis unveiled that the six-gene signature was an independent prognostic factor. The signature (0.806/0.795/0.822) showed a remarkable prognostic performance, with 1-, 3-, and 5-year time-dependent AUC, higher than for grade (0.612/0.638/0.649) and 1p19q codeletion status (0.606/0.658/0.676). High-risk patients had higher infiltrating immune cells. However, the specific immune escape was observed in the high-risk group after immune activation, owing to increasing immunosuppressive cells, inhibitory cytokines, and immune checkpoint molecules. Moreover, a novel nomogram model was developed to evaluate the survival in IDH-mutant LGGs patients.

CONCLUSION

The six-gene signature could be a promising prognostic biomarker, which is promising to promote individual therapy and improve the clinical outcomes of IDH-mutant gliomas. The study also refined the current classification system of IDH-mutant gliomas, classifying patients into two subtypes with distinct immunophenotypes and overall survival.

Prognostic biomarkers for predicting papillary thyroid carcinoma patients at high risk using nine genes of apoptotic pathway

Aberrant expressions of apoptotic genes have been associated with papillary thyroid carcinoma (PTC) in the past, however, their prognostic role and utility as biomarkers remains poorly understood. In this study, we analysed 505 PTC patients by employing Cox-PH regression techniques, prognostic index models and machine learning methods to elucidate the relationship between overall survival (OS) of PTC patients and 165 apoptosis related genes. It was observed that nine genes (ANXA1, TGFBR3, CLU, PSEN1, TNFRSF12A, GPX4, TIMP3, LEF1, BNIP3L) showed significant association with OS of PTC patients. Five out of nine genes were found to be positively correlated with OS of the patients, while the remaining four genes were negatively correlated. These genes were used for developing risk prediction models, which can be utilized to classify patients with a higher risk of death from the patients which have a good prognosis. Our voting-based model achieved highest performance (HR = 41.59, p = 3.36x10^-4, C = 0.84, logrank-p = 3.8x10^-8). The performance of voting-based model improved significantly when we used the age of patients with prognostic biomarker genes and achieved HR = 57.04 with p = 10⁻⁴ (C = 0.88, logrank-p = 1.44x10^-9). We also developed classification models that can classify high risk patients (survival ≤ 6 years) and low risk patients (survival > 6 years). Our best model achieved AUROC of 0.92. Further, the expression pattern of the prognostic genes was verified at mRNA level, which showed their differential expression between normal and PTC samples. Also, the immunostaining results from HPA validated these findings. Since these genes can also be used as potential therapeutic targets in PTC, we also identified potential drug molecules which could modulate their expression profile. The study briefly revealed the key prognostic biomarker genes in the apoptotic pathway whose altered expression is associated with PTC progression and aggressiveness. In addition to this, risk assessment models proposed here can help in efficient management of PTC patients.

Universal and cross-cancer prognostic biomarkers for predicting survival risk of cancer patients from expression profile of apoptotic pathway genes

Numerous cancer-specific prognostic models have been developed in the past, wherein one model is applicable for only one type of cancer. In this study, an attempt has been made to identify universal or multi-cancer prognostic biomarkers and develop models for predicting survival risk across different types of cancer patients. In order to accomplish this, we gauged the prognostic role of mRNA expression of 165 apoptosis-related genes across 33 cancers in the context of patient survival. Firstly, we identified specific prognostic biomarker genes for 30 cancers. The cancer-specific prognostic models achieved a minimum Hazard Ratio, HR_SKCM = 1.99 and maximum HR_THCA = 41.59. Secondly, a comprehensive analysis was performed to identify universal biomarkers across many cancers. Our best prognostic model consisted of 11 genes (TOP2A, ISG20, CD44, LEF1, CASP2, PSEN1, PTK2, SATB1, SLC20A1, EREG, and CD2) and stratified risk groups across 27 cancers (HR_OV = 1.53-HR_UVM = 11.74). The model was validated on eight independent cancer cohorts and exhibited a comparable performance. Further, we clustered cancer-types on the basis of shared survival related apoptosis genes. This approach proved helpful in development of cross-cancer prognostic models. To show its efficacy, a prognostic model consisting of 15 genes was thereby developed for LGG-KIRC pair (HR_KIRC = 3.27, HR_LGG = 4.23). Additionally, we predicted potential therapeutic candidates for LGG-KIRC high risk patients.

The Influence of Gene Aberrations on Survival in Resected IDH Wildtype Glioblastoma Patients: A Single-Institution Study

This prospective population-based study on a group of 132 resected IDH-wildtype (IDH-wt) glioblastoma (GBM) patients assesses the prognostic and predictive value of selected genetic biomarkers and clinical factors for GBM as well as the dependence of these values on the applied therapeutic modalities. The patients were treated in our hospital between June 2006 and June 2015. Clinical data and tumor samples were analyzed to determine the frequencies of TP53, MDM2, EGFR, RB1, BCR, and CCND1 gene aberrations and the duplication/deletion statuses of the 9p21.3, 1p36.3, 19q13.32, and 10p11.1 chromosome regions. Cut-off values distinguishing low (LCN) and high (HCN) copy number status for each marker were defined. Additionally, MGMT promoter methylation and IDH1/2 mutation status were investigated retrospectively. Young age, female gender, Karnofsky scores (KS) above 80, chemoradiotherapy, TP53 HCN, and CCND1 HCN were identified as positive prognostic factors, and smoking was identified as a negative prognostic factor. Cox proportional regression models of the chemoradiotherapy patient group revealed TP53 HCN and CCND1 HCN to be positive prognostic factors for both progression-free survival and overall survival. These results confirmed the influence of key clinical factors (age, KS, adjuvant oncotherapy, and smoking) on survival in GBM IDH-wt patients and demonstrated the prognostic and/or predictive importance of CCND1, MDM2, and 22q12.2 aberrations.

Hub gene identification and prognostic model construction for isocitrate dehydrogenase mutation in glioma

•

We identified ten hub genes which were driving IDH status in GBM and LGG.

•

We constructed a prognostic model for IDH-mutant patients.

•

Our findings have important clinical implications for accurate treatment in glioma.

Our study attempted to identify hub genes related to isocitrate dehydrogenase (IDH) mutation in glioma and develop a prognostic model for IDH-mutant glioma patients. In a first step, ten hub genes significantly associated with the IDH status were identified by weighted gene coexpression analysis (WGCNA). The functional enrichment analysis demonstrated that the most enriched terms of these hub genes were cadherin binding and glutathione metabolism. Three of these hub genes were significantly linked with the survival of glioma patients. 328 samples of IDH-mutant glioma were separated into two datasets: a training set (N = 228) and a test set (N = 100). Based on the training set, we identified two IDH-mutant subtypes with significantly different pathological features by using consensus clustering. A 31 gene-signature was identified by the least absolute shrinkage and selection operator (LASSO) algorithm and used for establishing a differential prognostic model for IDH-mutant patients. In addition, the test set was employed for validating the prognostic model, and the model was proven to be of high value in classifying prognostic information of samples. The functional annotation revealed that the genes related to the model were mainly enriched in nuclear division, DNA replication, and cell cycle. Collectively, this study provided novel insights into the molecular mechanism of IDH mutation in glioma, and constructed a prognostic model which can be effective for predicting prognosis of glioma patients with IDH-mutation, which might promote the development of IDH target agents in glioma therapies and contribute to accurate prognostication and management in IDH-mutant glioma patients.

Clinical implication of cellular vaccine in glioma: current advances and future prospects

Gliomas, especially glioblastomas, represent one of the most aggressive and difficult-to-treat human brain tumors. In the last few decades, clinical immunotherapy has been developed and has provided exceptional achievements in checkpoint inhibitors and vaccines for cancer treatment. Immunization with cellular vaccines has the advantage of containing specific antigens and acceptable safety to potentially improve cancer therapy. Based on T cells, dendritic cells (DC), tumor cells and natural killer cells, the safety and feasibility of cellular vaccines have been validated in clinical trials for glioma treatment. For TAA engineered T cells, therapy mainly uses chimeric antigen receptors (IL13Rα2, EGFRvIII and HER2) and DNA methylation-induced technology (CT antigen) to activate the immune response. Autologous dendritic cells/tumor antigen vaccine (ADCTA) pulsed with tumor lysate and peptides elicit antigen-specific and cytotoxic T cell responses in patients with malignant gliomas, while its pro-survival effect is biased. Vaccinations using autologous tumor cells modified with TAAs or fusion with fibroblast cells are characterized by both effective humoral and cell-mediated immunity. Even though few therapeutic effects have been observed, most of this therapy showed safety and feasibility, asking for larger cohort studies and better guidelines to optimize cellular vaccine efficiency in anti-glioma therapy.

Analyses of metastasis-associated genes in IDH wild-type glioma

Background

Glioma is the most common malignant tumor of the brain. The existence of metastatic tumor cells is an important cause of recurrence even after radical glioma resection.

Methods

Single-cell sequencing data and high-throughput data were downloaded from GEO database and TCGA/CGGA database. By means of PCA and tSNE clustering methods, metastasis-associated genes in glioma were identified. GSEA explored possible biological functions that these metastasis-associated genes may participate in. Univariate and multivariate Cox regression were used to construct a prognostic model.

Results

Glioma metastatic cells and metastasis-associated genes were identified. The prognostic model based on metastasis-associated genes had good sensitivity and specificity for the prognosis of glioma. These genes may be involved in signal pathways such as cellular protein catabolic process, p53 signaling pathway, transcriptional misregulation in cancer and JAK-STAT signaling pathway.

Conclusion

This study explored glioma metastasis-associated genes through single-cell sequencing data mining, and aimed to identify prognostic metastasis-associated signatures for glioma and may provide potential targets for further cancer research.

Identification of the prognostic value of immune gene signature and infiltrating immune cells for esophageal cancer patients

BACKGROUND

Esophageal cancer (ESCA) is one of the deadliest solid malignancies with worse survival rate worldwide. Here, we aimed to establish an immune-gene prognostic signature for predicting patients' survival and providing accurate targets for personalized therapy or immunotherapy.

METHODS

Gene expression profile of patients with ESCA were download from The Cancer Genome Atlas (TCGA) database (dataset 1: n = 159) and immune-related genes from the ImmPORT database. Dataset 1 was subdivided into two groups (dataset 2: n = 80; dataset 3: n = 79). Kaplan-Meier and receiver operating characteristic (ROC) curves were plotted to validate the predictive effect of the prognostic signature on the three datasets. TIMER and CIBERSORT analysis were used to evaluate the correlation between the prognostic signature and infiltrating immune cells.

RESULTS

We constructed a prognostic signature composed of six immune genes (HSPA6, S100A12, FABP3, DKK1, OSM and NR2F2). Kaplan-Meier curves validated the good predictive ability of the prognostic signature in datasets 1, 2 and 3 (P = 0.0034, P = 0.0081, and P = 0.0363, respectively). The area under the curve (AUC) of the ROC curves validated the predictive accuracy of the immune signature (AUCs = 0.757, 0.800, and 0.701, respectively). We also revealed the good prognostic value of the immune cells, including activated memory CD4 T cells, T follicular helper cells and monocytes. Potential target drugs, including Olopatadine and Amlexanox, were identified for clinical therapies to improve patients' survival outcomes.

CONCLUSION

Our study indicated that the immune-related prognostic signature could serve as a novel biomarker for predicting patients' prognosis and providing new immunotherapy targets in ESCA.