Molecular and clinical characterization of IDH associated immune signature in lower-grade gliomas

Isocitrate Dehydrogenase Inhibitors in Glioma: From Bench to Bedside

Isocitrate dehydrogenase (IDH) mutant gliomas are a primary malignancy of the central nervous system (CNS) malignancies, most commonly affecting adults under the age of 55. Standard of care therapy for IDH-mutant gliomas involves maximal safe resection, radiotherapy, and chemotherapy. However, despite good initial responses to multimodality treatment, recurrence is virtually universal. IDH-mutant gliomas represent a life-limiting prognosis. For this reason, there is a great need for novel treatments that can prolong survival. Uniquely for IDH-mutant gliomas, the IDH mutation is the direct driver of oncogenesis through its oncometabolite 2-hydroxygluterate. Inhibition of this mutated IDH with a corresponding reduction in 2-hydroxygluterate offers an attractive treatment target. Researchers have tested several IDH inhibitors in glioma through preclinical and early clinical trials. A phase III clinical trial of an IDH1 and IDH2 inhibitor vorasidenib yielded promising results among patients with low-grade IDH-mutant gliomas who had undergone initial surgery and no radiation or chemotherapy. However, many questions remain regarding optimal use of IDH inhibitors in clinical practice. In this review, we discuss the importance of IDH mutations in oncogenesis of adult-type diffuse gliomas and current evidence supporting the use of IDH inhibitors as therapeutic agents for glioma treatment. We also examine unresolved questions and propose potential directions for future research.

Transcriptomic characterization and construction of M2 macrophage-related prognostic and immunotherapeutic signature in ovarian metastasis of gastric cancer

BACKGROUND

Ovarian metastasis (OM) poses a major threat to the outcome of gastric cancer (GC) patients. Recently, immunotherapy emerged as a novel promising therapeutic strategy to treat late-stage GC, whereas its efficacy is influenced by tumor immune microenvironment (TIME). M2 macrophage, a key subset within TIME, plays dual immunosuppressive and pro-tumorigenic roles in cancer progression and is recognized as a potential therapeutic target. However, molecular mechanisms underlying OM remain elusive and the TIME-related prognostic and immunotherapeutic index for these patients is yet to establish.

METHODS

Differential expressed genes (DEGs) between paired normal mucosa, primary GC and OM of patients from Fudan University Shanghai Cancer Center (FUSCC) cohort (n = 6) were identified by transcriptome sequencing, followed by the functional annotation of enriched hallmark pathways of DEGs between them. CIBERSORT was used to profile the relative expression level of 22 immune cell subsets in normal tissues, primary and metastatic tumors, followed by weighted gene coexpression network analysis (WGCNA) uncovering immune cell-correlated gene sets. The intersected genes between DEGs and M2 macrophage-related genes were processed by least absolute shrinkage and selection operator (LASSO) regression analysis to construct a predictive signature, M2GO, which was further validated by training set and test set of The Cancer Genome Atlas-Stomach Adenocarcinoma (TCGA-STAD), GSE62254 and GSE84437 cohorts. GC patients were divided into M2GO-high and -low subgroup according to the optimal cutoff value of the M2GO score. Furthermore, the clinical, molecular and immune features between M2GO-high and -low subgroups were analyzed. Clinical cohorts of immunotherapy were used to validate the predictive value of M2GO in regard to immunotherapy effectiveness.

RESULTS

Transcriptomic sequencing and follow-up analyses of triple-matched normal tissues, primary and ovarian metastatic tumors identified distinctive sets of DEGs and enriched immune-, cancer- and metastasis-related pathways between them. Of note, M2 macrophage, a major immunosuppressive and pro-tumorigenic component within TIME, was significantly up-regulated in OMs. WGCNA and LASSO regression were applied to establish a novel OM- and M2 macrophage-related predictive signature, M2GO, based on M2 macrophage-related prognostic genes including GJA1, MAGED1 and SERPINE1. M2GO served as an independent prognostic factor of GC patients. Comprehensive molecular and immune characterization of M2GO-based subgroups uncovered their distinctive features in terms of enriched functional pathways, tumor mutation burden, key immune checkpoints, major regulators of natural immune cGAS-STING pathway, infiltrated subsets of immune cells and tumor immune exclusion/dysfunction (TIDE) score. Notably, the M2GO score was significantly lower in responsive group than non-responsive group (P < 0.05) in clinical cohort of metastatic GC patients undergoing immunotherapy.

CONCLUSION

Transcriptomic characterization of paired normal mucosae, primary and ovarian metastatic tumors revealed their unique molecular and immune features. Follow-up analyses established a novel OM- and M2 macrophage-related signature, M2GO, which served as a promising prognostic and immunotherapeutic biomarker to distinguish the clinical outcome, molecular and immune features of GC patients and predict their differential responses to immunotherapy.

An MRI radiomics approach to predict survival and tumour-infiltrating macrophages in gliomas

Preoperative MRI is one of the most important clinical results for the diagnosis and treatment of glioma patients. The objective of this study was to construct a stable and validatable preoperative T₂-weighted MRI-based radiomics model for predicting the survival of gliomas.

A total of 652 glioma patients across three independent cohorts were covered in this study including their preoperative T₂-weighted MRI images, RNA-seq and clinical data. Radiomic features (1731) were extracted from preoperative T₂-weighted MRI images of 167 gliomas (discovery cohort) collected from Beijing Tiantan Hospital and then used to develop a radiomics prediction model through a machine learning-based method. The performance of the radiomics prediction model was validated in two independent cohorts including 261 gliomas from the The Cancer Genomae Atlas database (external validation cohort) and 224 gliomas collected in the prospective study from Beijing Tiantan Hospital (prospective validation cohort). RNA-seq data of gliomas from discovery and external validation cohorts were applied to establish the relationship between biological function and the key radiomics features, which were further validated by single-cell sequencing and immunohistochemical staining.

The 14 radiomic features-based prediction model was constructed from preoperative T₂-weighted MRI images in the discovery cohort, and showed highly robust predictive power for overall survival of gliomas in external and prospective validation cohorts. The radiomic features in the prediction model were associated with immune response, especially tumour macrophage infiltration.

The preoperative T₂-weighted MRI radiomics prediction model can stably predict the survival of glioma patients and assist in preoperatively assessing the extent of macrophage infiltration in glioma tumours.

Gasdermin D Is a Novel Prognostic Biomarker and Relates to TMZ Response in Glioblastoma

The gasdermin (GSDM) family act as executioners during pyroptosis. However, its expression and biological role in glioma remain to be determined. This study carried out gene expression from six public datasets. Westerns blots and immunohistochemistry (IHC) staining were employed to examine GSDM expression in glioma in an in-house cohort. Kaplan-Meier and Cox regression analyses were performed to evaluate the prognostic role of GSDMs in glioma. Association between gene expression and immune infiltration was assessed by IHC and immunofluorescence (IF) staining of tissue sections. TMZ-induced pyroptosis was assessed by observation of morphological changes, WB and ELISA detection. Only GSDMD expression was upregulated in glioma compared with nontumor brain tissues both in the public datasets and in-house cohort. High GSDMD expression was significantly associated with WHO grade IV, IDH 1/2 wild-type and mesenchymal subtypes. Besides, high GSDMD expression was associated with shorter overall survival and could be used as an independent risk factor for poor outcomes in LGG and GBM. GO enrichment analysis and IHC validation revealed that GSDMD expression might participate in regulating macrophage infiltration and polarization. TMZ treatment induced the pyroptosis in GBM cells and GSDMD expression increased with after treating with TMZ in a time-dependent manner. Moreover, knocking down GSDMD obviously decreased IL-1β expression and reduced TMZ-induced pyroptosis in in vitro. GSDMD was a novel prognostic biomarker, as well as TMZ-treatment response marker in glioma.

Effects of 1p/19q Codeletion on Immune Phenotype in Low Grade Glioma

Background: Chromosome 1p/19q codeletion is one of the most important genetic alterations for low grade gliomas (LGGs), and patients with 1p/19q codeletion have significantly prolonged survival compared to those without the codeletion. And the tumor immune microenvironment also plays a vital role in the tumor progression and prognosis. However, the effect of 1p/19q codeletion on the tumor immune microenvironment in LGGs is unclear. Methods: Immune cell infiltration of 281 LGGs from The Cancer Genome Atlas (TCGA) and 543 LGGs from the Chinese Glioma Genome Atlas (CGGA) were analyzed for immune cell infiltration through three bioinformatics tools: ESTIMATE algorithm, TIMER, and xCell. The infiltrating level of immune cells and expression of immune checkpoint genes were compared between different groups classified by 1p/19q codeletion and IDH (isocitrate dehydrogenase) mutation status. The differential biological processes and signaling pathways were evaluated through Gene Set Enrichment Analysis (GSEA). Correlations were analyzed using Spearman correlation. Results: 1p/19q codeletion was associated with immune-related biological processes in LGGs. The infiltrating level of multiple kinds of immune cells and expression of immune checkpoint genes were significantly lower in 1p/19q codeletion LGGs compared to 1p/19q non-codeletion cohorts. There are 127 immune-related genes on chromosome 1p or 19q, such as TGFB1, JAK1, and CSF1. The mRNA expression of these genes was positively correlated with their DNA copy number. These genes are distributed in multiple immune categories, such as chemokines/cytokines, TGF-β family members, and TNF family members, regulating immune cell infiltration and expression of the immune checkpoint genes in tumors. Conclusion: Our results indicated that 1p/19q codeletion status is closely associated with the immunosuppressive microenvironment in LGGs. LGGs with 1p/19q codeletion display less immune cell infiltration and lower expression of immune checkpoint genes than 1p/19q non-codeletion cases. Mechanistically, this may be, at least in part, due to the deletion of copy number of immune-related genes in LGGs with 1p/19q codeletion. Our findings may be relevant to investigate immune evasion in LGGs and contribute to the design of immunotherapeutic strategies for patients with LGGs.

SAMD9 Is Relating With M2 Macrophage and Remarkable Malignancy Characters in Low-Grade Glioma

Immunoreactions regulated by TAMs (Tumor-associated macrophages) play a pivotal role in tumorigenesis and metastasis. In recent decades, treatments based on immune regulation have achieved revolutionary breakthroughs in cancer targeted therapies. The phenotypes of TAMs in gliomas are more heterogeneous and inherently complex than can be simply defined by classification into the M1 and M2 polarized states. The detailed mechanisms surrounding infiltrating macrophage phenotype and glioma characteristics remain undefined. SAMD9 (Sterile Alpha Motif Domain-Containing Protein 9) was found to be highly expressed in glioma and closely related to histological and genetic features in CGGA and TCGA databases. Simultaneously, we present evidence to show that there was a positive association between SAMD9 and malignancy characters in LGG. Univariable and Multivariate proportional hazard Cox analysis showed that SAMD9 was an independent prognostic factor for LGG. Surprisingly, Gene Ontology (GO) analysis showed SAMD9 expression level was remarkably well correlated with immunological responses and the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis supported the connection with immune responses and tumorigenesis. Immune infiltration analysis demonstrated that high SAMD9 expression resulted in an accumulation of macrophages by CIBERSORT and TIMER databases, especially positively related to macrophage total marker gene AIF1 and Macrophage M2 marker gene CD163. IHC staining further indicated a high correlation of SAMD9 with those specific macrophage markers in the immune response. Human THP-1 cells were induced into M2 macrophages, which were then co-cultured with LN229 cells. Silencing of SAMD9 by shRNA in LN229 cells attenuated the infiltration abilities of M2 macrophage. SAMD9 explored immune response via relating of M2 macrophage in vitro. Our results revealed SAMD9 acted as the malignancy characters in LGG, enrichment with M2 macrophage.

TEAD4 is a novel independent predictor of prognosis in LGG patients with IDH mutation

TEA domain family members (TEADs) play important roles in tumor progression. Till now, the genomic status of TEADs in patients with glioma has not been well investigated. To confirm whether the genomic status of TEADs could affect the prognosis of patients with glioma, the copy number variation (CNV), mutation and expression data of glioma cohorts in The Cancer Genome Atlas, Gene Expression Omnibus and Chinese Glioma Genome Atlas were comprehensively analyzed. Results showed that TEAD CNV frequency in lower grade gliomas (LGGs) was higher than in glioblastoma multiforme (GBM). Multivariate cox regression analysis showed that TEAD4 CNV increase was significantly associated with overall survival (OS) and disease-free survival (DFS) in LGGs (OS p = 0.022, HR = 1.444, 95% CI: 1.054–1.978; DFS p = 0.005, HR = 1.485, 95% CI: 1.124–1.962), while not in GBM. Patients with TEAD4 CNV increase showed higher expression level of TEAD4 gene. In LGG patients with IDH mutation, those with higher TEAD4 expression levels had shorter OS and DFS. Integrating TEAD4 CNV increase, IDH mutations, TP53 mutation, ATRX mutation and 1p19q co-deletion would separate patients with LGG into four groups with significant differences in prognosis. These study results suggested that TEAD4 variations were independent predictive biomarkers for the prognosis in patients with LGG with IDH mutation.

Molecular and Clinicopathological Characterization of a Prognostic Immune Gene Signature Associated With MGMT Methylation in Glioblastoma

Background: O6-methylguanine-DNA methyltransferase (MGMT) methylation status affects tumor chemo-resistance and the prognosis of glioblastoma (GBM) patients. We aimed to investigate the role of MGMT methylation in the regulation of GBM immunophenotype and discover an effective biomarker to improve prognosis prediction of GBM patients. Methods: A total of 769 GBM patients with clinical information from five independent cohorts were enrolled in the present study. Samples from the Cancer Genome Atlas (TCGA) dataset were used as the training set, whereas transcriptome data from the Chinese Glioma Genome Atlas (CGGA) RNA-seq, CGGA microarray, GSE16011, and the Repository for Molecular Brain Neoplasia (REMBRANDT) cohort were used for validation. A series of bioinformatics approaches were carried out to construct a prognostic signature based on immune-related genes, which were tightly related to the MGMT methylation status. In silico analyses were performed to investigate the influence of the signature on immunosuppression and remodeling of the tumor microenvironment. Then, the utility of this immune gene signature was analyzed by the development and evaluation of a nomogram. In vitro experiments were further used to verify the immunologic function of the genes in the signature. Results: We found that MGMT unmethylation was closely associated with immune-related biological processes in GBM. Sixty-five immune genes were more highly expressed in the MGMT unmethylated than the MGMT-methylated group. An immune gene-based risk model was further established to divide patients into high and low-risk groups, and the prognostic value of this signature was validated in several GBM cohorts. Functional analyses manifested a universal up-regulation of immune-related pathways in the high-risk group. Furthermore, the risk score was highly correlated to the immune cell infiltration, immunosuppression, inflammatory activities, as well as the expression levels of immune checkpoints. A nomogram was developed for clinical application. Knockdown of the five genes in the signature remodeled the immunosuppressive microenvironment by restraining M2 macrophage polarization and suppressing immunosuppressive cytokines production. Conclusions: MGMT methylation is strongly related to the immune responses in GBM. The immune gene-based signature we identified may have potential implications in predicting the prognosis of GBM patients and mechanisms underlying the role of MGMT methylation.

Magnetic resonance imaging-based radiomic features for extrapolating infiltration levels of immune cells in lower-grade gliomas

PURPOSE

To extrapolate the infiltration levels of immune cells in patients with lower-grade gliomas (LGGs) using magnetic resonance imaging (MRI)-based radiomic features.

METHODS

A retrospective dataset of 516 patients with LGGs from The Cancer Genome Atlas (TCGA) database was analysed for the infiltration levels of six types of immune cells using Tumor IMmune Estimation Resource (TIMER) based on RNA sequencing data. Radiomic features were extracted from 107 patients whose pre-operative MRI data are available in The Cancer Imaging Archive; 85 and 22 of these patients were assigned to the training and testing cohort, respectively. The least absolute shrinkage and selection operator (LASSO) was applied to select optimal radiomic features to build the radiomic signatures for extrapolating the infiltration levels of immune cells in the training cohort. The developed radiomic signatures were examined in the testing cohort using Pearson's correlation.

RESULTS

The infiltration levels of B cells, CD4+ T cells, CD8+ T cells, macrophages, neutrophils and dendritic cells negatively correlated with overall survival in the 516 patient cohort when using univariate Cox's regression. Age, Karnofsky Performance Scale, WHO grade, isocitrate dehydrogenase mutant status and the infiltration of neutrophils correlated with survival using multivariate Cox's regression analysis. The infiltration levels of the 6 cell types could be estimated by radiomic features in the training cohort, and their corresponding radiomic signatures were built. The infiltration levels of B cells, CD8+ T cells, neutrophils and macrophages estimated by radiomics correlated with those estimated by TIMER in the testing cohort. Combining clinical/genomic features with the radiomic signatures only slightly improved the prediction of immune cell infiltrations.

CONCLUSION

We developed MRI-based radiomic models for extrapolating the infiltration levels of immune cells in LGGs. Our results may have implications for treatment planning.

YTHDF2 correlates with tumor immune infiltrates in lower-grade glioma

Immunotherapy is an effective treatment for many cancer types. However, YTHDF2 effects on the prognosis of different tumors and correlation with tumor immune infiltration are unclear. Here, we analyzed The Cancer Genome Atlas and Gene Expression Omnibus data obtained through various web-based platforms. The analyses showed that YTHDF2 expression and associated prognoses may depend on cancer type. High YTHDF2 expression was associated with poor overall survival in lower-grade glioma (LGG). In addition, YTHDF2 expression positively correlated with expression of several immune cell markers, including PD-1, TIM-3, and CTLA-4, as well as tumor-associated macrophage gene markers, and isocitrate dehydrogenase 1 in LGG. These findings suggest that YTHDF2 is a potential prognostic biomarker that correlates with LGG tumor-infiltrating immune cells.

A novel signature based on microvascular invasion predicts the recurrence of HCC

Background and objectives

In hepatocellular carcinoma (HCC) patients, microvascular invasion (MVI) is associated with worse outcomes regardless of treatment. No single reliable preoperative factor exists to predict MVI. The aim of the work described here was to develop a new MVI− based mRNA biomarker to differentiate between high and low risk patients.

Methods

Using The Cancer Genome Atlas (TCGA) database, we collected data from 315 HCC patients, including mRNA expression and complete clinical data. We generated a seven-mRNA signature to predict patient outcomes. The mRNA signature was validated using the GSE36376 cohort. Finally, we tested the formula in our own 53 HCC patients using qPCR for the seven mRNAs and analyzing the computed tomography (CT) features.

Results

This seven‐mRNA signature significantly correlated with length of recurrence-free survival (RFS) and overall survival (OS) for both the training and validation groups. RFS and OS were briefer in high risk versus low risk patients. A Kaplan–Meier analysis also indicated that survival time was significantly shortened in the high risk group versus the low risk group. Time-dependent receiver operating characteristic analysis demonstrated good predictive performance for the seven-mRNA signature. The mRNA signature also acts as an independent factor according to a Multivariate analysis. Our results are consistent with the seven-mRNA formula risk score.

Conclusion

Our research showed a novel seven-mRNA biomarker based on MVI predicting RFS and OS in HCC patients. This mRNA signature can stratify patients into subgroups based on their risk of recurrence to help guide individualized treatment and precision management in HCC.

Systematic identification, development, and validation of prognostic biomarkers involving the tumor-immune microenvironment for glioblastoma

Gliomas are infiltrative neoplasms with a highly invasive nature. Due to its distinct genomic, genetic and epigenetic features, the immune prognostic signature (IPS) and immune microenvironment of glioblastoma (GBM) merit further research. We aimed to explore prognosis-related immune genes and develop an IPS model for predicting prognosis in GBM. RNA-sequencing data, as well as clinical information, from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) public cohorts were analyzed. To develop the IPS, least absolute shrinkage and selection operator (LASSO) Cox analysis was performed for immune-related genes that were differentially expressed between GBM and normal tissues. Then, interaction effects of the IPS on the immune microenvironment were systematically analyzed; the precise prognostic model was developed based on the IPS and clinical data and was then further validated. A total of 21 immune prognostic genes were identified based on GBM microenvironment status. An 8-gene IPS was established, and the GBM patients were effectively stratified into low- and high-risk groups in the TCGA cohort as a training set. Univariate and multivariate Cox analyses revealed that IPS was an independent prognostic factor, and the prognostic performance of individual IPS genes was systematically illustrated. In addition, a comprehensive and novel nomogram model was initially established to estimate overall survival in TCGA-GBM patients, and high-risk patients had higher levels of dendritic cell and neutrophil infiltration. Furthermore, the nomogram model was developed and validated in the CGGA validation set. The low-risk IPS was linked to a stronger response to anti-PD-L1 immunotherapy and clinical advantages in the IMvigor210 cohort. This novel IPS with promising biomarkers classifies GBM patients into subgroups with distinct clinical outcomes and immunophenotypes. Our findings and this resource may help to characterize the immune microenvironment, inform cancer immunotherapy and facilitate the development of precision immuno-oncology.

Inflammation and lymphocyte infiltration are associated with shorter survival in patients with high-grade glioma

Glioma represents a serious health burden in terms of morbidity and mortality. The prognostic significance of the lymphoid and myeloid infiltrates in glioma is not clearly determined. Moreover, the characterization of different leukocyte subsets in the tumor microenvironment relies mainly on immunohistochemistry observations, and data about their association with prognosis are contradictory. Here, we performed acomprehensive study of both the tumor-infiltrating and circulating immune compartments of patients with high-grade glioma. Nineteen tumor biopsies and 30 PBMC samples were analyzed by RNA sequencing. Validation was performed on The Cancer Genome Atlas (TCGA) RNA sequencing data from glioma and on additional 39 tumor biopsies analyzed by flow cytometry. We identified prognostic tumor and peripheral immune signatures, which associate increased inflammation, immune infiltration and activation with shorter overall survival in high-grade glioma patients. Importantly, we confirmed our observations by flow cytometry analysis and validated the tumor-signature using the TCGA dataset. In addition, both tumor genotype and grade associated with the degree of glioma immune infiltration. Unlike in the majority of cancers, lymphocyte infiltration at the tumor site is anegative prognostic factor in glioma, suggesting the ambivalent pro-tumorigenic role of immune responses in glioma.

Novel Immune-Related Gene Signature for Risk Stratification and Prognosis of Survival in Lower-Grade Glioma

Objective

Despite several clinicopathological factors being integrated as prognostic biomarkers, the individual variants and risk stratification have not been fully elucidated in lower grade glioma (LGG). With the prevalence of gene expression profiling in LGG, and based on the critical role of the immune microenvironment, the aim of our study was to develop an immune-related signature for risk stratification and prognosis prediction in LGG.

Methods

RNA-sequencing data from The Cancer Genome Atlas (TCGA), Genome Tissue Expression (GTEx), and Chinese Glioma Genome Atlas (CGGA) were used. Immune-related genes were obtained from the Immunology Database and Analysis Portal (ImmPort). Univariate, multivariate cox regression, and Lasso regression were employed to identify differentially expressed immune-related genes (DEGs) and establish the signature. A nomogram was constructed, and its performance was evaluated by Harrell’s concordance index (C-index), receiver operating characteristic (ROC), and calibration curves. Relationships between the risk score and tumor-infiltrating immune cell abundances were evaluated using CIBERSORTx and TIMER.

Results

Noted, 277 immune-related DEGs were identified. Consecutively, 6 immune genes (CANX, HSPA1B, KLRC2, PSMC6, RFXAP, and TAP1) were identified as risk signature and Kaplan–Meier curve, ROC curve, and risk plot verified its performance in TCGA and CGGA datasets. Univariate and multivariate Cox regression indicated that the risk group was an independent predictor in primary LGG. The prognostic signature showed fair accuracy for 3- and 5-year overall survival in both internal (TCGA) and external (CGGA) validation cohorts. However, predictive performance was poor in the recurrent LGG cohort. The CIBERSORTx algorithm revealed that naïve CD4⁺ T cells were significant higher in low-risk group. Conversely, the infiltration levels of M1-type macrophages, M2-type macrophages, and CD8⁺T cells were significant higher in high-risk group in both TCGA and CGGA cohorts.

Conclusion

The present study constructed a robust six immune-related gene signature and established a prognostic nomogram effective in risk stratification and prediction of overall survival in primary LGG.

Prognostic and predictive value of an immune infiltration signature in diffuse lower-grade gliomas

BACKGROUNDLower-grade gliomas (LGGs) vary widely in terms of the patient's overall survival (OS). There is no current, valid method that could exactly predict the survival. The effects of intratumoral immune infiltration on clinical outcome have been widely reported. Thus, we aim to develop an immune infiltration signature to predict the survival of LGG patients.METHODSWe analyzed 1216 LGGs from 5 public data sets, including 2 RNA sequencing data sets and 3 microarray data sets. Least absolute shrinkage and selection operator (LASSO) Cox regression was used to select an immune infiltration signature and build a risk score. The performance of the risk score was assessed in the training set (329 patients), internal validation set (140 patients), and 4 external validation sets (405, 118, 88, and 136 patients).RESULTSAn immune infiltration signature consisting of 20 immune metagenes was used to generate a risk score. The performance of the risk score was thoroughly verified in the training and validation sets. Additionally, we found that the risk score was positively correlated with the expression levels of TGF-β and PD-L1, which were important targets of combination immunotherapy. Furthermore, a nomogram incorporating the risk score, patient's age, and tumor grade was developed to predict the OS, and it performed well in all the training and validation sets (C-index: 0.873, 0.881, 0.781, 0.765, 0.721, and 0.753).CONCLUSIONThe risk score based on the immune infiltration signature has reliable prognostic and predictive value for patients with LGGs and is a potential biomarker for the cotargeting immunotherapy.FUNDINGThis work was supported by The National Natural Science Foundation of China (grant nos. 81472370 and 81672506), the Natural Science Foundation of Beijing (grant no. J180005), the National High Technology Research and Development Program of China (863 Program, grant no. 2014AA020610), and the National Basic Research Program of China (973 Program, grant no. 2014CB542006).

Gene Expression Profiling Stratifies IDH-Wildtype Glioblastoma With Distinct Prognoses

Objectives: In the present study, we aimed to determine the candidate genes that may function as biomarkers to further distinguish patients with isocitrate dehydrogenase (IDH)-wildtype glioblastoma (GBM), which are heterogeneous with respect to clinical outcomes.

Materials and Methods: We selected 41 candidate genes associated with overall survival (OS) using univariate Cox regression from IDH-wildtype GBM patients based on RNA sequencing (RNAseq) expression data from the Chinese Glioma Genome Atlas (CGGA, n = 105) and The Cancer Genome Atlas (TCGA, n = 139) cohorts. Next, a seven-gene-based risk signature was formulated according to Least Absolute Shrinkage and Selection Operator (LASSO) regression algorithm in the CGGA RNAseq database as a training set, while another 525 IDH-wildtype GBM patient TCGA datasets, consisting of RNA sequencing and microarray data, were used for validation. Patient survival in the low- and high-risk groups was calculated using Kaplan-Meier survival curve analysis and the log-rank test. Uni-and multivariate Cox regression analysis was used to assess the prognosis value. Gene oncology (GO) and gene set enrichment analysis (GSEA) were performed for the functional analysis of the seven-gene-based risk signature.

Results: We developed a seven-gene-based signature, which allocated each patient to a risk group (low or high). Patients in the high-risk group had dramatically shorter overall survival than their low-risk counterparts in three independent cohorts. Univariate and multivariate analysis showed that the seven-gene signature remained an independent prognostic factor. Moreover, the seven-gene risk signature exhibited a striking prognostic validity, with AUC of 78.4 and 73.9%, which was higher than for traditional “age” (53.7%, 62.4%) and “GBM sub-type” (57.7%, 52.9%) in the CGGA- and TCGA-RNAseq databases, respectively. Subsequent bioinformatics analysis predicted that the seven-gene signature was involved in the inflammatory response, immune response, cell adhesion, and apoptotic process.

Conclusions: Our findings indicate that the seven-gene signature could be a potential prognostic biomarker. This study refined the current classification system of IDH-wildtype GBM and may provide a novel perspective for the research and individual therapy of IDH-wildtype GBM.

Genetic Programs Driving Oncogenic Transformation: Lessons from in Vitro Models

Cancer complexity relies on the intracellular pleiotropy of oncogenes/tumor suppressors and in the strong interplay between tumors and micro- and macro-environments. Here we followed a reductionist approach, by analyzing the transcriptional adaptations induced by three oncogenes (RAS, MYC, and HDAC4) in an isogenic transformation process. Common pathways, in place of common genes became dysregulated. From our analysis it emerges that, during the process of transformation, tumor cells cultured in vitro prime some signaling pathways suitable for coping with the blood supply restriction, metabolic adaptations, infiltration of immune cells, and for acquiring the morphological plasticity needed during the metastatic phase. Finally, we identified two signatures of genes commonly regulated by the three oncogenes that successfully predict the outcome of patients affected by different cancer types. These results emphasize that, in spite of the heterogeneous mutational burden among different cancers and even within the same tumor, some common hubs do exist. Their location, at the intersection of the various signaling pathways, makes a therapeutic approach exploitable.

Identification of a Tumor Microenvironment-Related Eight-Gene Signature for Predicting Prognosis in Lower-Grade Gliomas

Lower-grade gliomas (LrGG), characterized by invasiveness and heterogeneity, remain incurable with current therapies. Clinicopathological features provide insufficient information to guide optimal individual treatment and cannot predict prognosis completely. Recently, an increasing amount of studies indicate that the tumor microenvironment plays a pivotal role in tumor malignancy and treatment responses. However, studies focusing on the tumor microenvironment (TME) of LrGG are still limited. In this study, taking advantage of the currently popular computational methods for estimating the infiltration of tumor-associated normal cells in tumor samples and using weighted gene co-expression network analysis, we screened the co-expressed gene modules associated with the TME and further identified the prognostic hub genes in these modules. Furthermore, eight prognostic hub genes (ARHGDIB, CLIC1, OAS3, PDIA4, PARP9, STAT1, TAP2, and TAGLN2) were selected to construct a prognostic risk score model using the least absolute shrinkage and selection operator method. Univariate and multivariate Cox regression analysis demonstrated that the risk score was an independent prognostic factor for LrGG. Moreover, time-dependent ROC curves indicated that our model had favorable efficiency in predicting both short- and long-term prognosis in LrGG patients, and the stratified survival analysis demonstrated that our model had prognostic value for different subgroups of LrGG patients. Additionally, our model had potential value for predicting the sensitivity of LrGG patients to radio- and chemotherapy. Besides, differential expression analysis showed that the eight genes were aberrantly expressed in LrGG compared to normal brain tissue. Correlation analysis revealed that the expression of the eight genes was significantly associated with the infiltration levels of six types of immune cells in LrGG. In summary, the TME-related eight-gene signature was significantly associated with the prognosis of LrGG patients. They might act as potential prognostic biomarkers for LrGG patients, offer better stratification for future clinical trials, and be candidate targets for immunotherapy, thus deserving further investigation.

Prognostic value of preoperative hematological markers combined with molecular pathology in patients with diffuse gliomas

The prediction of clinical outcome for patients with infiltrative gliomas is challenging. Although preoperative hematological markers have been proposed as predictors of survival in glioma and other cancers, systematic investigations that combine these data with other relevant clinical variables are needed to improve prognostic accuracy and patient outcomes. We investigated the prognostic value of preoperative hematological markers, alone and in combination with molecular pathology, for the survival of 592 patients with Grade II-IV diffuse gliomas. On univariate analysis, increased neutrophil-to-lymphocyte ratio (NLR), platelet-to-lymphocyte ratio (PLR), and monocyte-to-lymphocyte ratio (MLR), and decreased albumin-to-globulin ratio (AGR), all predicted poor prognosis in Grade II/III gliomas. Multivariate analysis incorporating tumor status based on the presence of IDH mutations, TERT promoter mutations, and 1p/19q codeletion showed that in lower-grade gliomas, high NLR predicted poorer survival for the triple-negative, IDH mutation only, TERT mutation only, and IDH and TERT mutation groups. NLR was an independent prognostic factor in Grade IV glioma. We therefore propose a prognostic model for diffuse gliomas based on the presence of IDH and TERT promoter mutations, 1p/19q codeletion, and NLR. This model classifies lower-grade gliomas into nine subgroups that can be combined into four main risk groups based on survival projections.

Identification of a Five-Pseudogene Signature for Predicting Survival and Its ceRNA Network in Glioma

Background: Glioma is the most common primary brain tumor with a dismal prognosis. It is urgent to develop novel molecular biomarkers and conform to individualized schemes. Methods: Differentially expressed pseudogenes between low grade glioma (LGG) and glioblastoma multiforme (GBM) were identified in the training cohort. Least absolute shrinkage and selection operator (LASSO) regression and multivariate Cox proportional hazards regression analyses were used to select pseudogenes associated with prognosis of glioma. A risk signature was constructed based on the selected pseudogenes for predicting the survival of glioma patients. A pseudogene-miRNA-mRNA regulatory network was established and visualized using Cytoscape 3.5.1. Gene Oncology (GO) and signaling pathway analyses were performed on the targeted genes to investigate functional roles of the risk signature. Results: Five pseudogenes (ANXA2P2, EEF1A1P9, FER1L4, HILS1, and RAET1K) correlating with glioma survival were selected and used to establish a risk signature. Time-dependent receiver operating characteristic (ROC) curves revealed that the risk signature could accurately predict the 1, 3, and 5-year survival of glioma patients. GO and signaling pathway analyses showed that the risk signature was involved in regulation of proliferation, migration, angiogenesis, and apoptosis in glioma. Conclusions: In this study, a risk signature with five pseudogenes was constructed and shown to accurately predict 1-, 3-, and 5-year survival for glioma patient. The risk signature may serve as a potential target against glioma.

A Risk Classification System With Five-Gene for Survival Prediction of Glioblastoma Patients

Objective: Glioblastoma (GBM) is the most common and fatal primary brain tumor in adults. It is necessary to identify novel and effective biomarkers or risk signatures for GBM patients. Methods: Differentially expressed genes (DEGs) between GBM and low-grade glioma (LGG) in TCGA samples were screened out and weight correlation network analysis (WGCNA) was performed to confirm WHO grade-related genes. Five genes were selected via multivariate Cox proportional hazards regression analysis and were used to construct a risk signature. A nomogram composed of the risk signature and clinical characters (age, radiotherapy, and chemotherapy experience) was established to predict 1, 3, 5-year survival rate for GBM patients. Results: One hundred ninety-four DEGs in blue gene module were found to be positively related to WHO grade via WGCNA. Five genes (DES, RANBP17, CLEC5A, HOXC11, POSTN) were selected to construct a risk signature for GBM via R language. This risk signature was identified to independently predict the outcome of GBM patients, as well as stratified by IDH1 status, MGMT promoter status, and radio-chemotherapy. The nomogram was established which combined the risk signature with clinical factors. The results of c-index, ROC curve and calibration plot revealed the nomogram showing a good accuracy for predicting 1, 3, or 5-year survival of GBM patients. Conclusion: The risk signature with five genes could serve as an independent factor for predicting the prognosis of patients with GBM. Moreover, the nomogram with the risk signature and clinical traits proved to perform better for predicting 1, 3, 5-year survival rate.

Amino acid metabolism-related gene expression-based risk signature can better predict overall survival for glioma

Metabolic reprogramming has been proposed to be a hallmark of cancer. Aside from the glycolytic pathway, the metabolic changes of cancer cells primarily involve amino acid metabolism. However, in glioma, the characteristics of the amino acid metabolism‐related gene set have not been systematically profiled. In the present study, RNA sequencing expression data from 309 patients in the Chinese Glioma Genome Atlas database were included as a training set, while another 550 patients within The Cancer Genome Atlas database were used to validate. Consensus clustering of the 309 samples yielded two robust groups. Compared with Cluster1, Cluster2 correlated with a better clinical outcome. We then developed an amino acid metabolism‐related risk signature for glioma. Our results showed that patients in the high‐risk group had dramatically shorter overall survival than low‐risk counterparts in any subgroup, stratified by isocitrate dehydrogenase and 1p/19q status based on the 2016 World Health Organization classification guidelines. The 30‐gene signature showed better prognostic value than the traditional factors “age” and “grade” by analyzing the receiver operating characteristic curve with areas under curve of 0.966, 0.692, 0.898 and 0.975, 0.677, 0.885 for 3‐ and 5‐year survival, respectively. Moreover, univariate and multivariate analysis showed that the 30‐gene signature was an independent prognostic factor for glioma. Furthermore, Gene Ontology analysis and Gene Set Enrichment Analysis showed that tumors with a high risk score correlated with various aspects of the malignancy of glioma. In summary, we demonstrated a novel amino acid metabolism‐related risk signature for predicting prognosis for glioma.

Bridging Cancer Biology with the Clinic: Comprehending and Exploiting IDH Gene Mutations in Gliomas

Isocitrate dehydrogenases 1 and 2 (IDH1/2) are enzymes that play a major role in the Krebs cycle. Mutations in these enzymes are found in the majority of lower gliomas and secondary glioblastomas, but also in myeloid malignancies and other cancers. IDH1 and IDH2 mutations are restricted to specific arginine residues in the active site of the enzymes and are gain-of-function, i.e. they confer a neomorphic enzyme activity resulting in the accumulation of D-2-hydroxyglutarate (2-HG). 2-HG is an oncometabolite causing profound metabolic dysregulation which, among others, results in methylator phenotypes and in defects in homologous recombination repair. In this review, we summarize current knowledge regarding the function of normal and mutated IDH, explain the possible mechanisms through which these mutations might drive malignant transformation of progenitor cells in the central nervous system, and provide a comprehensive review of potential treatment strategies for IDH-mutated malignancies, focusing on gliomas.

Preoperative Changes in Hematological Markers and Predictors of Glioma Grade and Survival

Background: Preoperative hematological markers that indicate nutritional, coagulation, and inflammation statuses have prognostic value for gliomas. This study aimed to investigate hematological markers with regard to tumor grades, isocitrate dehydrogenase mutations (IDH), age, and sex in patients with gliomas. Methods: From 2008 to 2017, patients with a pathological diagnosis of glioma who underwent surgery were retrospectively enrolled in this study. Information from clinical records, including age, sex, preoperative experiment tests (routine blood tests, biochemistry, and coagulation examinations), pathological results, and IDH status, was collected. A univariable survival analysis was performed. Hematological factors such as the neutrophil-to-lymphocyte ratio (NLR), platelet-to-lymphocyte-ratio (PLR), and albumin-to-globulin (AGR) were calculated. The prognostic nutrition index (PNI) was calculated as 10 × serum albumin value (g/dl) + 0.005 × peripheral lymphocyte count (per mm³). Results: Our study included 706 patients. The univariate analysis showed that age, IDH-1, and hematological factors were all significantly associated with overall survival (OS) in patients with gliomas. Our results showed that inflammation markers (NLR, PLR, and fibrinogen) were positively associated with age, whereas AGR was negatively associated with age. The PLR was significantly increased, whereas the AGR and PNI were decreased in women with gliomas, as compared with men. We found that inflammation markers increased and nutrition markers decreased with gliomas grade. However, these hematological markers did not significantly differ with IDH status. NLR was the best single hematological marker for distinguishing glioblastoma (GBM) [0.684 (0.645-0.723)], IDH-wt GBM [0.672 (0.631-0.71)] from other gliomas subtypes. Combinations of age with PNI and age with AGR were the best predictors of GBM [0.750 (0.713-0.786)] and IDH-wt GBM [0.759 (0.719-0.798)], respectively. Conclusion: Preoperative hematological marker levels vary among glioma grades and have high predictive values for GBM.