Investigation of Genetic Determinants of Glioma Immune Phenotype by Integrative Immunogenomic Scale Analysis

Glioma hexokinase 3 positively correlates with malignancy and macrophage infiltration

BACKGROUND

Glioma is the main subtype of primary central nervous system (CNS) tumor with high malignancy and poor prognosis under current therapeutic approaches. Glycolysis and suppressive tumor microenvironment (TME) are key markers of glioma with great importance for aggressive features of glioma and inferior clinical outcomes. Hexokinase 3 (HK3) is an important rate-limiting enzyme in glycolysis, but its function in glioma remains unknown.

METHODS

This study comprehensively assessed the expression distribution and immunological effect of HK3 via pan-cancer analysis based on datasets from Genotype Tissue Expression (GTEx), Cancer Cell Line Encyclopedia (CCLE), and The Cancer Genome Atlas (TCGA). Furthermore, it explored the malignant phenotype and genomic landscape between low-HK3 and high-HK3 expression groups in gliomas from Chinese Glioma Genome Atlas (CGGA) and TCGA. Moreover, data from the TIMER website predicted the relationship between macrophage infiltration and HK3 expression. Also, single-cell sequencing data were used to validate the relationship.

RESULTS

For pan-cancer patients, HK3 was expressed in various cancers. The results showed that HK3 was highly expressed in gliomas and positively correlated with tumor-infiltrating immune cells (TIICs), immune checkpoints, immunomodulators, and chemokines. Meanwhile, HK3 expression was highest in normal immune cells and tissues. In gliomas, the expression of HK3 was found to be closely correlated with the malignant clinical characteristics and the infiltration of macrophages. Also, HK3 was proven to be positively associated with macrophage through single-cell sequencing data and immunohistochemistry techniques. Finally, it is predicted that samples with high HK3 expression are often malignant entities and also significant genomic aberrations of driver oncogenes.

CONCLUSIONS

This is the first comprehensive research to figure out the relationship between HK3 and TME characteristics in gliomas. HK3 is positively associated with macrophage infiltration and can induce the immunosuppressive TME and malignant phenotype of gliomas.

Spatial distance between tumor and lymphocyte can predict the survival of patients with resectable lung adenocarcinoma

Background and objective

Spatial interaction between tumor-infiltrating lymphocytes (TILs) and tumor cells is valuable in predicting the effectiveness of immune response and prognosis amongst patients with lung adenocarcinoma (LUAD). Recent evidence suggests that the spatial distance between tumor cells and lymphocytes also influences the immune responses, but the distance analysis based on Hematoxylin and Eosin (H&E) -stained whole-slide images (WSIs) remains insufficient. To address this issue, we aim to explore the relationship between distance and prognosis prediction of patients with LUAD in this study.

Methods

We recruited patients with resectable LUAD from three independent cohorts in this multi-center study. We proposed a simple but effective deep learning-driven workflow to automatically segment different cell types in the tumor region using the HoVer-Net model, and quantified the spatial distance (DIST) between tumor cells and lymphocytes based on H&E-stained WSIs. The association of DIST with disease-free survival (DFS) was explored in the discovery set (D1, n = 276) and the two validation sets (V1, n = 139; V2, n = 115).

Results

In multivariable analysis, the low DIST group was associated with significantly better DFS in the discovery set (D1, HR, 0.61; 95 % CI, 0.40-0.94; p = 0.027) and the two validation sets (V1, HR, 0.54; 95 % CI, 0.32-0.91; p = 0.022; V2, HR, 0.44; 95 % CI, 0.24-0.81; p = 0.009). By integrating the DIST with clinicopathological factors, the integrated model (full model) had better discrimination for DFS in the discovery set (C-index, D1, 0.745 vs. 0.723) and the two validation sets (V1, 0.621 vs. 0.596; V2, 0.671 vs. 0.650). Furthermore, the computerized DIST was associated with immune phenotypes such as immune-desert and inflamed phenotypes.

Conclusions

The integration of DIST with clinicopathological factors could improve the stratification performance of patients with resectable LUAD, was beneficial for the prognosis prediction of LUAD patients, and was also expected to assist physicians in individualized treatment.

Integrating transcriptomics and machine learning for immunotherapy assessment in colorectal cancer

Colorectal cancer (CRC) is a highly prevalent and lethal malignancy worldwide. Although immunotherapy has substantially improved CRC outcomes, intolerance remains a major concern among most patients. Considering the pivotal role of the tumor microenvironment (TME) in tumor progression and treatment outcomes, profiling the TME at the transcriptomic level can provide novel insights for developing CRC treatment strategies. Seventy-seven TME-associated signatures were acquired from previous studies. To elucidate variations in prognosis, clinical features, genomic alterations, and responses to immunotherapy in CRC, we employed a non-negative matrix factorization algorithm to categorize 2595 CRC samples of 27 microarrays from the Gene Expression Omnibus database. Three machine learning techniques were employed to identify a signature specific to immunotherapy. Subsequently, the mechanisms by which this signature interacts with TME subtypes and immunotherapy were investigated. Our findings revealed five distinct TME subtypes (TMESs; TMES1-TMES5) in CRC, each exhibiting a unique pattern of immunotherapy response. TMES1, TMES4, and TMES5 had relatively inferior outcomes, TMES2 was associated with the poorest prognosis, and TMES3 had a superior outcome. Subsequent investigations revealed that activated dendritic cells could enhance the immunotherapy response rate, with their augmentation effect closely associated with the activation of CD8+T cells. We successfully classified CRC into five TMESs, each demonstrating varying response rates to immunotherapy. Notably, the application of machine learning to identify activated dendritic cells helped elucidate the underlying mechanisms contributing to these differences. We posit that these TMESs hold promising clinical implications for prognostic evaluation and guidance of immunotherapy strategies, thereby providing valuable insights to inform clinical decision-making.

The Crosstalk between the EGFR and IFN-γ Pathways and Synergistic Roles in Survival Prediction and Immune Escape in Gliomas

Glioma is the most common primary malignant brain tumor. The poor prognosis of gliomas, especially glioblastoma (GBM), is associated with their unique molecular landscape and tumor microenvironment (TME) features. The epidermal growth factor receptor (EGFR) gene is one of the frequently altered loci in gliomas, leading to the activation of the EGFR signaling pathway and thus, promoting the genesis of gliomas. Whether there exist factors within the TME that can lead to EGFR activation in the context of gliomas is currently unexplored. In total, 702 samples from The Cancer Genome Atlas (TCGA) and 325 samples from The Chinese Glioma Genome Atlas (CGGA) were enrolled in this study. Gene signatures related to EGFR signaling and interferon-γ (IFN-γ) response were established via the LASSO-COX algorithm. Gene Set Enrichment Analysis (GSEA) and Gene Ontology (GO) analysis were applied for function exploration. Kaplan-Meier (KM) curves and single sample GSEA (ssGSEA) of immune cell subpopulations were performed to analyze the prognosis and TME characteristics of different subgroups. Moreover, Western blotting (WB) and flow cytometry (FCM) demonstrated the correlation between IFN-γ and EGFR signaling activation and the subsequent induction of programmed death ligand 1 (PD-L1) expression. An EGFR signaling-related risk score was established, and a higher score was correlated with poorer prognosis and a more malignant phenotype in gliomas. Biological function analysis revealed that a higher EGFR-related score was significantly associated with various cytokine response pathways, especially IFN-γ. Long-term (7 days) exposure to IFN-γ (400 ng/mL) induced the activation of EGFR signaling in the u87 cell line. Next, an IFN-γ response-related risk score was established; the combination of these two scores could be used to further reclassify gliomas into subtypes with different clinical features and TME features. Double high-risk samples tended to have a poorer prognosis and more immunosuppressive TME. Additionally, FCM discovered that the activation of EGFR signaling via EGF (100 ng/mL) could trigger PD-L1 protein expression. This research indicates that IFN-γ, an inflammatory cytokine, can activate the EGFR pathway. The combination of EGFR signaling and IFN-γ response pathway can establish a more precise classification of gliomas.

Radiomics for characterization of the glioma immune microenvironment

Increasing evidence suggests that besides mutational and molecular alterations, the immune component of the tumor microenvironment also substantially impacts tumor behavior and complicates treatment response, particularly to immunotherapies. Although the standard method for characterizing tumor immune profile is through performing integrated genomic analysis on tissue biopsies, the dynamic change in the immune composition of the tumor microenvironment makes this approach not feasible, especially for brain tumors. Radiomics is a rapidly growing field that uses advanced imaging techniques and computational algorithms to extract numerous quantitative features from medical images. Recent advances in machine learning methods are facilitating biological validation of radiomic signatures and allowing them to "mine" for a variety of significant correlates, including genetic, immunologic, and histologic data. Radiomics has the potential to be used as a non-invasive approach to predict the presence and density of immune cells within the microenvironment, as well as to assess the expression of immune-related genes and pathways. This information can be essential for patient stratification, informing treatment decisions and predicting patients' response to immunotherapies. This is particularly important for tumors with difficult surgical access such as gliomas. In this review, we provide an overview of the glioma microenvironment, describe novel approaches for clustering patients based on their tumor immune profile, and discuss the latest progress on utilization of radiomics for immune profiling of glioma based on current literature.

The clinical relevance and prediction efficacy from therapy of tumor microenvironment related signature score in colorectal cancer

Introduction

As the top 3 cancer in terms of incidence and mortality, the first-line treatment for CRC includes FOLFOX, FOLFIRI, Cetuximab or immunotherapy. However, the drug sensitivity of patients to regimens is different. There has been increasing evidence that immune components of TME can affect the sensitivity of patients to drugs. Therefore, it is necessary to define novo molecular subtypes of CRC based on TME immune components, and screen patients who are sensitive to the treatments, to make personalized therapy possible.

Methods

We analyzed the expression profiles and 197 TME-related signatures of 1775 patients using ssGSEA, univariate Cox proportional risk model and LASSO-Cox regression model, and defined a novo molecular subtype (TMERSS) of CRC. Simultaneously, we compared the clinicopathological factors, antitumor immune activity, immune cell abundance and differences of cell states in different TMERSS subtypes. In addition, patients sensitive to the therapy were screened out by correlation analysis between TMERSS subtypes and drug responses.

Results

Compared with low TMERSS subtype, high TMERSS subtype has a better outcome, which may be associated to higher abundance of antitumor immune cell in high TMERSS subtype. Our findings suggested that the high TMERSS subtype may have a higher proportion of respondents to Cetuximab agent and immunotherapy, while the low TMERSS subtype may be more suitable for treatment with FOLFOX and FOLFIRI regimens.

Discussion

In conclusion, the TMERSS model may provide a partial reference for the prognosis evaluation of patients, the prediction of drug sensitivity, and the implementation of clinical decision-making.

Immune landscape-based machine-learning-assisted subclassification, prognosis, and immunotherapy prediction for glioblastoma

Introduction

As a malignant brain tumor, glioblastoma (GBM) is characterized by intratumor heterogeneity, a worse prognosis, and highly invasive, lethal, and refractory natures. Immunotherapy has been becoming a promising strategy to treat diverse cancers. It has been known that there are highly heterogeneous immunosuppressive microenvironments among different GBM molecular subtypes that mainly include classical (CL), mesenchymal (MES), and proneural (PN), respectively. Therefore, an in-depth understanding of immune landscapes among them is essential for identifying novel immune markers of GBM.

Methods and results

In the present study, based on collecting the largest number of 109 immune signatures, we aim to achieve a precise diagnosis, prognosis, and immunotherapy prediction for GBM by performing a comprehensive immunogenomic analysis. Firstly, machine-learning (ML) methods were proposed to evaluate the diagnostic values of these immune signatures, and the optimal classifier was constructed for accurate recognition of three GBM subtypes with robust and promising performance. The prognostic values of these signatures were then confirmed, and a risk score was established to divide all GBM patients into high-, medium-, and low-risk groups with a high predictive accuracy for overall survival (OS). Therefore, complete differential analysis across GBM subtypes was performed in terms of the immune characteristics along with clinicopathological and molecular features, which indicates that MES shows much higher immune heterogeneity compared to CL and PN but has significantly better immunotherapy responses, although MES patients may have an immunosuppressive microenvironment and be more proinflammatory and invasive. Finally, the MES subtype is proved to be more sensitive to 17-AAG, docetaxel, and erlotinib using drug sensitivity analysis and three compounds of AS-703026, PD-0325901, and MEK1-2-inhibitor might be potential therapeutic agents.

Conclusion

Overall, the findings of this research could help enhance our understanding of the tumor immune microenvironment and provide new insights for improving the prognosis and immunotherapy of GBM patients.

Radiogenomics to characterize the immune-related prognostic signature associated with biological functions in glioblastoma

OBJECTIVES

The tumor microenvironment and immune cell infiltration (ICI) associated with glioblastoma (GBM) play a vital role in cancer development, progression, and prognosis. This study aimed to establish an ICI-related prognostic biomarker and explore the associations between ICI signatures and radiomic features in patients with GBM.

METHODS

The gene expression and survival data of patients with GBM were obtained from three databases. Based on the ICI pattern, an individualized ICI score for each GBM patient was developed in the discovery set (n = 400) and independently verified in the validation set (n = 374). A total of 5915 radiomic features were extracted from the intratumoral and peritumoral regions. Recursive feature elimination and support vector machine methods were performed to select the key features and generate a model predictive of low- or high- ICI scores. The prognostic value of the identified radio genomic model was examined in an independent dataset (n = 149) using imaging and survival data.

RESULTS

We found that higher ICI scores often indicated worse patient prognosis (multivariable hazard ratio: 0.48 and 0.63 in discovery and validation set, respectively) and higher expression levels of immune checkpoint-related genes. A model that combined 11 radiomic features could well distinguish tumors with different ICI scores (AUC = 0.96, accuracy = 94%). This model was proven to be helpful for noninvasive prognostic stratification in an independent validation cohort.

CONCLUSIONS

ICI scores may serve as an effective prognostic biomarker to characterize potential biological processes in patients with GBM. This ICI signature can be evaluated noninvasively through radiogenomic analysis.

KEY POINTS

• Immune cell infiltration (ICI) scores can serve as an effective prognostic biomarker in patients with glioblastoma. • The ICI signature can be evaluated noninvasively through radiomic features derived from the intratumoral and peritumoral regions. • The prognostic value of the radiogenomic model can be verified by independent survival and MRI data.

PD-L2 Serves as a Potential Prognostic Biomarker That Correlates With Immune Infiltration and May Predict Therapeutic Sensitivity in Lower-Grade Gliomas

Although patients with lower-grade gliomas (LGGs; grades II and III) have a relatively favorable prognosis, patients frequently relapse and tend to progress to higher-grade gliomas, leading to treatment resistance, poor survival, and ultimately treatment failure. However, until now, thorough research has not yet been reported on the relationship between PD-L2 and immune infiltration and therapeutic sensitivity to immunotherapy and TMZ-based chemotherapy of LGGs. In this study, we found that the expression of PD-L2 is upregulated in glioma, with high PD-L2 expression predicting a worse prognosis. Univariate and multivariate Cox regression analysis both indicated that PD-L2 represented an independent prognostic factor with high accuracy in survival prediction for LGGs. A nomogram comprising of age, grade, IDH mutation, and PD-L2 was established for predicting OS. Additionally, PD-L2 was found to be remarkably correlated with immune infiltration and some anti-tumor immune functions. The degree of PD-L2 expression was also found to be strongly related to the prediction of therapeutic sensitivity to immunotherapy and TMZ-based chemotherapy. Furthermore, immunohistochemistry demonstrated that PD-L2 and the macrophage biomarker CD68 were both increased in glioma, with PD-L2 expression having a strong positive connection with CD68 expression. Taken together, PD-L2 is a prognostic biomarker for LGGs patients that may provide novel insights into glioma individualized therapeutic strategies and guide effective immunotherapy and chemotherapy.

A Novel Classification Model for Lower-Grade Glioma Patients Based on Pyroptosis-Related Genes

Recent studies demonstrated that pyroptosis plays a crucial role in shaping the tumor-immune microenvironment. However, the influence of pyroptosis on lower-grade glioma regarding immunotherapy and targeted therapy is still unknown. This study analyzed the variations of 33 pyroptosis-related genes in lower-grade glioma and normal tissues. Our study found considerable genetic and expression alterations in heterogeneity among lower-grade gliomas and normal brain tissues. There are two pyroptosis phenotypes in lower-grade glioma, and they exhibited differences in cell infiltration characteristics and clinical characters. Then, a PyroScore model using the lasso-cox method was constructed to measure the level of pyroptosis in each patient. PyroScore can refine the lower-grade glioma patients with a stratified prognosis and a distinct tumor immune microenvironment. Pyscore may also be an effective factor in predicting potential therapeutic benefits. In silico analysis showed that patients with a lower PyroScore are expected to be more sensitive to targeted therapy and immunotherapy. These findings may enhance our understanding of pyroptosis in lower-grade glioma and might help optimize risk stratification for the survival and personalized management of lower-grade glioma patients.

Prognostic Features of the Tumor Immune Microenvironment in Glioma and Their Clinical Applications: Analysis of Multiple Cohorts

Glioma is the most common malignant tumor of the central nervous system. Tumor purity is a source of important prognostic factor for glioma patients, showing the key roles of the microenvironment in glioma prognosis. In this study, we systematically screened functional characterization related to the tumor immune microenvironment and constructed a risk model named Glioma MicroEnvironment Functional Signature (GMEFS) based on eight cohorts. The prognostic value of the GMEFS model was also verified in another two glioma cohorts, glioblastoma (GBM) and low-grade glioma (LGG) cohorts, from The Cancer Genome Atlas (TCGA). Nomograms were established in the training and testing cohorts to validate the clinical use of this model. Furthermore, the relationships between the risk score, intrinsic molecular subtypes, tumor purity, and tumor-infiltrating immune cell abundance were also evaluated. Meanwhile, the performance of the GMEFS model in glioma formation and glioma recurrence was systematically analyzed based on 16 glioma cohorts from the Gene Expression Omnibus (GEO) database. Based on multiple-cohort integrated analysis, risk subpathway signatures were identified, and a drug–subpathway association network was further constructed to explore candidate therapy target regions. Three subpathways derived from Focal adhesion (path: 04510) were identified and contained known targets including platelet derived growth factor receptor alpha (PDGFRA), epidermal growth factor receptor (EGFR), and erb-b2 receptor tyrosine kinase 2 (ERBB2). In conclusion, the novel functional signatures identified in this study could serve as a robust prognostic biomarker, and this study provided a framework to identify candidate therapeutic target regions, which further guide glioma patients’ clinical decision.

Molecular landscape of IDH-mutant astrocytoma and oligodendroglioma grade 2 indicate tumor purity as an underlying genomic factor

Background

IDH-mutant astrocytoma and oligodendroglioma have an indolent natural history and are recognized as distinct entities of neoplasms. There is little knowledge on the molecular differences between IDH-mutant astrocytoma and oligodendroglioma grade 2. Therefore, we investigated the multiomics and clinical data regarding these two types of tumors.

Method

In silico analyses were performed around mRNA, somatic mutations, copy number alternations (CNAs), DNA methylation, microRNA (miRNA), epigenetics, immune microenvironment characterization and clinical features of the two types of gliomas. A diagnostic model incorporating tumor purity was further established using machine learning algorithms, and the predictive value was evaluated by receiver operative characteristic curves.

Results

Both types of gliomas shared chromosomal instability, and astrocytomas exhibited increased total CNAs compared to oligodendrogliomas. Oligodendrogliomas displayed distinct chromosome 4 (chr 4) loss, and subtyping of chr 7 gain/chr 4 loss (+ 7/− 4) presented the worst survival (P = 0.004) and progression-free interval (PFI) (P < 0.001). In DNA damage signatures, oligodendroglioma had a higher subclonal genome fraction (P < 0.001) and tumor purity (P = 0.001), and astrocytoma had a higher aneuploidy score (P < 0.001). Furthermore, astrocytomas exhibited inflamed immune cell infiltration, activated T cells and a potential response to immune checkpoint inhibitors (ICIs), while oligodendrogliomas were more homogeneous with increased tumor purity and decreased aggression. The tumor purity-involved diagnostic model exhibited great accuracy in identifying astrocytoma and oligodendroglioma.

Conclusion

This study addresses the similarities and differences between IDH-mutant astrocytoma and oligodendroglioma grade 2 and facilitates a deeper understanding of their molecular features, immune microenvironment, tumor purity and prognosis. The diagnostic tool developed using machine learning may offer support for clinical decisions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s10020-022-00454-z.