Tumor Immune Microenvironment Landscape in Glioma Identifies a Prognostic and Immunotherapeutic Signature

Chronic Stress Exacerbates the Immunosuppressive Microenvironment and Progression of Gliomas by Reducing Secretion of CCL3

As understanding of cancer has deepened, increasing attention has been turned to the roles of psychological factors, especially chronic stress-induced depression, in the occurrence and development of tumors. However, whether and how depression affects the progression of gliomas are still unclear. In this study, we have revealed that chronic stress inhibited the recruitment of tumor-associated macrophages (TAM) and other immune cells, especially M1-type TAMs and CD8+ T cells, and decreased the level of proinflammatory cytokines in gliomas, leading to an immunosuppressive microenvironment and glioma progression. Mechanistically, by promoting the secretion of stress hormones, chronic stress inhibited the secretion of the chemokine CCL3 and the recruitment of M1-type TAMs in gliomas. Intratumoral administration of CCL3 reprogrammed the immune microenvironment of gliomas and abolished the progression of gliomas induced by chronic stress. Moreover, levels of CCL3 and M1-type TAMs were decreased in the tumor tissues of glioma patients with depression, and CCL3 administration enhanced the antitumor effect of anti-PD-1 therapy in orthotopic models of gliomas undergoing chronic stress. In conclusion, our study has revealed that chronic stress exacerbates the immunosuppressive microenvironment and progression of gliomas by reducing the secretion of CCL3. CCL3 alone or in combination with an anti-PD-1 may be an effective immunotherapy for the treatment of gliomas with depression. See related Spotlight by Cui and Kang, p. 514.

Construction of a pyroptosis-related lncRNAs signature for predicting prognosis and immunotherapy response in glioma

Recent studies have proved that pyroptosis-related long non-coding RNAs (PRlncRNAs) are closely linked to tumor progression, prognosis, and immunity. Here, we systematically evaluated the correlation of PRlncRNAs with glioma prognosis. This study included 3 glioma cohorts (The Cancer Genome Atlas, Chinese Glioma Genome Atlas, and Gravendeel). Through Pearson correlation analysis, PRlncRNAs were screened from these 3 cohorts. Univariate Cox regression analysis was then carried out to determine the prognostic PRlncRNAs. A pyroptosis-related lncRNAs signature (PRLS) was then built by least absolute shrinkage and selection operator and multivariate Cox analyses. We systematically evaluated the correlation of the PRLS with the prognosis, immune features, and tumor mutation burden in glioma. A total of 14 prognostic PRlncRNAs overlapped in all cohorts and were selected as candidate lncRNAs. Based on The Cancer Genome Atlas cohort, a PRLS containing 7 PRlncRNAs was built. In all cohorts, the PRLS was proved to be a good predictor of glioma prognosis, with a higher risk score related to a poorer prognosis. We observed obvious differences in the immune microenvironment, immune cell infiltration level, and immune checkpoint expression in low- and high-risk subgroups. Compared with low-risk cases, high-risk cases had lower Tumor Immune Dysfunction and Exclusion scores and greater tumor mutation burden, indicating that high-risk cases can be more sensitive to immunotherapy. A nomogram combining PRLS and clinical parameters was constructed, which showed more robust and accurate predictive power. In conclusion, the PRLS is a potentially useful indicator for predicting prognosis and response to immunotherapy in glioma. Our findings may provide a useful insight into clinically individualized treatment strategies for patients.

Identification of two robust subclasses of sepsis with both prognostic and therapeutic values based on machine learning analysis

Background

Sepsis is a heterogeneous syndrome with high morbidity and mortality. Optimal and effective classifications are in urgent need and to be developed.

Methods and results

A total of 1,936 patients (sepsis samples, n=1,692; normal samples, n=244) in 7 discovery datasets were included to conduct weighted gene co-expression network analysis (WGCNA) to filter out candidate genes related to sepsis. Then, two subtypes of sepsis were classified in the training sepsis set (n=1,692), the Adaptive and Inflammatory, using K-means clustering analysis on 90 sepsis-related features. We validated these subtypes using 617 samples in 5 independent datasets and the merged 5 sets. Cibersort method revealed the Adaptive subtype was related to high infiltration levels of T cells and natural killer (NK) cells and a better clinical outcome. Immune features were validated by single-cell RNA sequencing (scRNA-seq) analysis. The Inflammatory subtype was associated with high infiltration of macrophages and a disadvantageous prognosis. Based on functional analysis, upregulation of the Toll-like receptor signaling pathway was obtained in Inflammatory subtype and NK cell-mediated cytotoxicity and T cell receptor signaling pathway were upregulated in Adaptive group. To quantify the cluster findings, a scoring system, called, risk score, was established using four datasets (n=980) in the discovery cohorts based on least absolute shrinkage and selection operator (LASSO) and logistic regression and validated in external sets (n=760). Multivariate logistic regression analysis revealed the risk score was an independent predictor of outcomes of sepsis patients (OR [odds ratio], 2.752, 95% confidence interval [CI], 2.234-3.389, P<0.001), when adjusted by age and gender. In addition, the validation sets confirmed the performance (OR, 1.638, 95% CI, 1.309-2.048, P<0.001). Finally, nomograms demonstrated great discriminatory potential than that of risk score, age and gender (training set: AUC=0.682, 95% CI, 0.643-0.719; validation set: AUC=0.624, 95% CI, 0.576-0.664). Decision curve analysis (DCA) demonstrated that the nomograms were clinically useful and had better discriminative performance to recognize patients at high risk than the age, gender and risk score, respectively.

Conclusions

In-depth analysis of a comprehensive landscape of the transcriptome characteristics of sepsis might contribute to personalized treatments and prediction of clinical outcomes.

Construction and validation of an angiogenesis-related gene expression signature associated with clinical outcome and tumor immune microenvironment in glioma

Background: Glioma is the most prevalent malignant intracranial tumor. Many studies have shown that angiogenesis plays a crucial role in glioma tumorigenesis, metastasis, and prognosis. In this study, we conducted a comprehensive analysis of angiogenesis-related genes (ARGs) in glioma.

Methods: RNA-sequencing data of glioma patients were obtained from TCGA and CGGA databases. Via consensus clustering analysis, ARGs in the sequencing data were distinctly classified into two subgroups. We performed univariate Cox regression analysis to determine prognostic differentially expressed ARGs and least absolute shrinkage and selection operator Cox regression to construct a 14-ARG risk signature. The CIBERSORT algorithm was used to explore immune cell infiltration, and the ESTIMATE algorithm was applied to calculate immune and stromal scores.

Results: We found that the 14-ARG signature reflected the infiltration characteristics of different immune cells in the tumor immune microenvironment. Additionally, total tumor mutational burden increased significantly in the high-risk group. We combined the 14-ARG signature with patient clinicopathological data to construct a nomogram for predicting 1-, 3-, and 5-year overall survival with good accuracy. The predictive value of the prognostic model was verified in the CGGA cohort. SPP1 was a potential biomarker of glioma risk and was involved in the proliferation, invasion, and angiogenesis of glioma cells.

Conclusion: In conclusion, we established and validated a novel ARG risk signature that independently predicted the clinical outcomes of glioma patients and was associated with the tumor immune microenvironment.

Global Characterization of Metabolic Genes Regulating Survival and Immune Infiltration in Osteosarcoma

Background: The alterations in metabolic profile of tumors have been identified as one of the prognostic hallmarks of cancers, including osteosarcoma. These alterations are majorly controlled by groups of metabolically active genes. However, the regulation of metabolic gene signatures in tumor microenvironment of osteosarcoma has not been well explained.

Objectives: Thus, we investigated the sets of previously published metabolic genes in osteosarcoma patients and normal samples.

Methods: We applied computational techniques to identify metabolic genes involved in the immune function of tumor microenvironment (TME) and survival and prognosis of the osteosarcoma patients. Potential candidate gene PAICS (phosphoribosyl aminoimidazole carboxylase, phosphoribosyl aminoimidazole succino carboxamide synthetase) was chosen for further studies in osteosarcoma cell lines for its role in cell proliferation, migration and apoptosis.

Results: Our analyses identified a list of metabolic genes differentially expressed in osteosarcoma tissues. Next, we scrutinized the list of genes correlated with survival and immune cells, followed by clustering osteosarcoma patients into three categories: C1, C2, and C3. These analyses led us to choose PAICS as potential candidate gene as its expression showed association with poor survival and negative correlation with the immune cells. Furthermore, we established that loss of PAICS induced apoptosis and inhibited proliferation, migration, and wound healing in HOS and MG-63 cell lines. Finally, the results were supported by constructing and validating a prediction model for prognosis of the osteosarcoma patients.

Conclusion: Here, we conclude that metabolic genes specifically PAICS play an integral role in the immune cell infiltration in osteosarcoma TME, as well as cancer development and metastasis.