Redox Regulator GLRX Is Associated With Tumor Immunity in Glioma

The Significance of the Redox Gene in the Prognosis and Therapeutic Response of Glioma

OBJECTIVES

Glioblastoma (GBM) is a fatal adult central nervous system tumor. Due to its high heterogeneity, the survival rate and prognosis of patients are poor. Thousands of people die of this disease every year all over the world. At present, the treatment of GBM is mainly through surgical resection and the combination of later drugs, radiotherapy, and chemotherapy. An abnormal redox system is involved in the malignant progression and treatment tolerance of glioma, which is the main reason for poor survival and prognosis. The construction of a GBM redox-related prognostic model may be helpful in improving the redox immunotherapy and prognosis of GBM.

METHODS

Based on glioma transcriptome data and clinical data from The Cancer Genome Atlas, databases, a risk model of redox genes was constructed by univariate and multivariate Cox analysis. The good prediction performance of the model was verified by the internal validation set of The Cancer Genome Atlas, and the external data of Chinese Glioma Genome Atlas.

RESULTS

The results confirmed that the higher the risk score, the worse the survival of patients. Age and isocitrate dehydrogenase status were significantly correlated with risk scores. The analysis of immune infiltration and immunotherapy found that there were significant differences in the immune score, matrix score, and ESTIMATE score between high and low-risk groups. reverse transcription polymerase chain reaction and immunohistochemical staining of glioma samples confirmed the expression of the hub gene.

CONCLUSION

Our study suggests that the 5 oxidative-related genes nitricoxidesynthase3 , NCF2 , VASN , FKBP1B , and TXNDC2 are hub genes, which may provide a reliable prognostic tool for glioma clinical treatment.

Glutaredoxin 2 Protein (Grx2) as an Independent Prognostic Factor Associated with the Survival of Colon Adenocarcinoma Patients

Glutaredoxin 2 (Grx2; Glrx2) is a glutathione-dependent oxidoreductase located in mitochondria, which is central to the regulation of glutathione homeostasis and mitochondrial redox, and plays a crucial role in highly metabolic tissues. In response to mitochondrial redox signals and oxidative stress, Grx2 can catalyze the oxidation and S-glutathionylation of membrane-bound thiol proteins in mitochondria. Therefore, it can have a significant impact on cancer development. To investigate this further, we performed an immunohistochemical analysis of Grx2 protein expression in colon adenocarcinoma samples collected from patients with primary colon adenocarcinoma (stage I and II) and patients with metastasis to regional lymph nodes (stage III). The results of our study revealed a significant relationship between the immunohistochemical expression of Grx2 and tumor histological grade, depth of invasion, regional lymph node involvement, angioinvasion, staging, and PCNA immunohistochemical expression. It was found that 87% of patients with stage I had high levels of Grx2 expression. In contrast, only 33% of patients with stage II and 1% of patients with stage III had high levels of Grx2 expression. Moreover, the multivariate analysis revealed that the immunohistochemical expression of Grx2 protein apart from the grade of tumor differentiation was an independent prognostic factors for the survival of patients with colon adenocarcinoma. Studies analyzing Grx2 levels in patients' blood confirmed that the highest levels of serum Grx2 protein was also found in stage I patients, which was reflected in the survival curves. A higher level of Grx2 in the serum has been associated with a more favorable outcome. These results were supported by in vitro analysis conducted on colorectal cancer cell lines that corresponded to stages I, II, and III of colorectal cancer, using qRT-PCR and Western Blot.

Comprehensive analysis of oxidative stress-related lncRNA signatures in glioma reveals the discrepancy of prognostic and immune infiltration

Oxidative stress refers to the process of reactive oxide species (ROS) increase in human body due to various factors, which leads to oxidative damage in human tissues. Current studies have confirmed that sustained oxidative stress is one of the distinctive features throughout the development of tumors. Numerous reports have shown that lncRNAs can regulate the process of oxidative stress through multiple pathways. However, the relationship between glioma-associated oxidative stress and lncRNAs is not clearly investigated. RNA sequencing data of GBM (glioblastoma) and LGG (low grade glioma) and corresponding clinical data were retrieved from the TCGA database. Oxidative stress related lncRNAs (ORLs) were identified by Pearson correlation analysis. Prognostic models for 6-ORLs were structured in the training cohort by univariate Cox regression analysis, multivariate Cox regression analysis and LASSO regression analysis. We constructed the nomogram and verified its predictive efficacy by Calibration curves and DCA decision curves. The biological functions and pathways of 6-ORLs-related mRNAs were inferred by Gene Set Enrichment Analysis. Immune cell abundance and immune function associated with risk score (RS) were estimated by ssGSEA, CIBERSORT and MCPcounter synthetically. External validation of the signature was completed using the CGGA-325 and CGGA-693 datasets. 6-ORLs signature-AC083864.2, AC107294.1, AL035446.1, CRNDE, LINC02600, and SNAI3-AS1-were identified through our analysis as being predictive of glioma prognosis. Kaplan-Meier and ROC curves indicated that the signature has a dependable predictive efficacy in the TCGA training cohort, validation cohort and CGGA-325/CGGA-693 test cohort. The 6-ORLs signature were verified to be independent prognostic predictors by multivariate cox regression and stratified survival analysis. Nomogram built with risk scores had strong predictive efficacy for patients' overall survival (OS). The outcomes of the functional enrichment analysis revealing potential molecular regulatory mechanisms for the 6-ORLs. Patients in the high-risk subgroup presented a significant immune microenvironment of macrophage M0 and cancer-associated fibroblast infiltration which was associated with a poorer prognosis. Finally, the expression levels of 6-ORLs in U87/U251/T98/U138 and HA1800 cell lines were verified by RT-qPCR. The nomogram in this study has been made available as a web version for clinicians. This 6-ORLs risk signature has the capabilities to predict the prognosis of glioma patients, assist in evaluating immune infiltration, and assess the efficacy of various anti-tumor systemic therapy regimens.

Impact of redox-related genes on tumor microenvironment immune characteristics and prognosis of high-grade gliomas

Introduction

High-grade glioma (HGG) defines a group of brain gliomas characterized by contrast enhancement, high tumor heterogeneity, and poor clinical outcome. Disturbed reduction-oxidation (redox) balance has been frequently associated with the development of tumor cells and their microenvironment (TME).

Methods

To study the influence of redox balance on HGGs and their microenvironment, we collected mRNA-sequencing and clinical data of HGG patients from TCGA and CGGA databases and our own cohort. Redox-related genes (ROGs) were defined as genes in the MSigDB pathways with keyword "redox" that were differentially expressed between HGGs and normal brain samples. Unsupervised clustering analysis was used to discover ROG expression clusters. Over-representation analysis (ORA), gene set enrichment analysis (GSEA) and gene set variation analysis (GSVA) were also employed to understand the biological implication of differentially expressed genes between HGG clusters. CIBERSORTx and ESTIMATE were used to profile the immune TME landscapes of tumors, and TIDE was used to evaluated the potential response to immune checkpoint inhibitors. Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression was used to construct HGG-ROG expression risk signature (GRORS).

Results

Seventy-five ROGs were found and consensus clustering using the expression profile of ROGs divided the both IDH-mutant (IDHmut) and IDH-wildtype (IDHwt) HGGs into subclusters with different prognosis. Functional enrichment analysis revealed that the differential aggressiveness between redox subclusters in IDHmut HGGs were significantly associated with cell cycle regulation pathways, while IDHwt HGG redox subclusters showed differentially activated immune-related pathways. In silico TME analysis on immune landscapes in the TME showed that the more aggressive redox subclusters in both IDHmut and IDHwt HGGs may harbor a more diverse composition of tumor-infiltrating immune cells, expressed a higher level of immune checkpoints and were more likely to respond to immune checkpoint blockade. Next, we established a GRORS which showed AUCs of 0.787, 0.884, and 0.917 in predicting 1-3-year survival of HGG patients in the held-out validation datasets, and the C-index of a nomogram combining the GRORS and other prognostic information reached 0.835.

Conclusion

Briefly, our results suggest that the expression pattern of ROGs was closely associated with the prognosis as well as the TME immune profile of HGGs, and may serve as a potential indicator for their response to immunotherapies.

BCAT2 Shapes a Noninflamed Tumor Microenvironment and Induces Resistance to Anti-PD-1/PD-L1 Immunotherapy by Negatively Regulating Proinflammatory Chemokines and Anticancer Immunity

To improve response rate of monotherapy of immune checkpoint blockade (ICB), it is necessary to find an emerging target in combination therapy. Through analyzing tumor microenvironment (TME)-related indicators, it is validated that BCAT2 shapes a noninflamed TME in bladder cancer. The outcomes of multiomics indicate that BCAT2 has an inhibitory effect on cytotoxic lymphocyte recruitment by restraining activities of proinflammatory cytokine/chemokine-related pathways and T-cell-chemotaxis pathway. Immunoassays reveal that secretion of CD8+ T-cell-related chemokines keeps a robust negative correlation with BCAT2, generating a decreasing tendency of CD8+ T cells around BCAT2+ tumor cells from far to near. Cotreatment of BCAT2 deficiency and anti-PD-1 antibody has a synergistic effect in vivo, implying the potential of BCAT2 in combination therapy. Moreover, the value of BCAT2 in predicting efficacy of immunotherapy is validated in multiple immunotherapy cohorts. Together, as a key molecule in TME, BCAT2 is an emerging target in combination with ICB and a biomarker of guiding precision therapy.

MCM10 as a novel prognostic biomarker and its relevance to immune infiltration in gliomas

BACKGROUND

Gliomas are one of the most common malignancies in the central nervous system (CNS). Members of the minichromosomal maintenance protein (MCM) family play an essential role in diagnosing and prognosis of malignant tumors. MCM10 is found in gliomas, but the prognosis and immune infiltration of gliomas has not been elucidated.

OBJECTIVE

To explore the biological function and immune infiltration of MCM10 in gliomas and provide a reference for the diagnosis, treatment, and prognostic evaluation.

METHODS

The MCM10 expression profile and the clinical information database of glioma patients were obtained from the China Glioma Genome Atlas (CGGA) and Cancer Genome Atlas (TCGA) glioma data. We analyzed the MCM10 expression levels in various cancers from The TCGA.RNA sequencing data were analyzed using the R packages to determine differentially expressed genes (DEGs) between high- and low MCM10 expressing GBM tissues from the TCGA-GBM database. The Wilcoxon rank sum test was used to compare MCM10 expression levels in glioma and normal brain tissue. To evaluate the value of MCM10 expressions in the prognosis of glioma patients by the Kaplan-Meier survival analysis, a univariate Cox analysis, multivariate Cox analysis, and a ROC curve analysis were used to analyze the correlation of MCM10 expression and the clinicopathological features of glioma patients using the TCGA database data. Subsequently, a functional enrichment analysis was performed to explore its potential signaling pathways and biological functions. Moreover, a single-sample gene set enrichment analysis was used to assess the extent of immune cell infiltration. Lastly, the authors constructed a nomogram to predict the overall survival rate (OS) of gliomas at 1, 3 and 5 years after diagnosis.

RESULTS

MCM10 is highly expressed in 20 cancer types including gliomas, and MCM10 expression was an independent adverse prognostic factor in glioma patients. Similarly, high expression of MCM10 was associated with advanced age (60 years), increased tumor grade, tumor recurrence or development of a secondary tumor, IDH wild-type, and non-codeletion of 1p19q (p< 0.01). The OS nomogram generated a consistency index of 0.821. The results of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment and Gene Ontology (GO) functional analysis showed that the cell-cycle-related and tumor-related signaling pathways were significantly enriched in the MCM10 high expression phenotype. Moreover, signaling pathways were significantly enriched in Gene Set Enrichment Analysis (GSEA), including Rho GTPases, M phase, DNA repair, extracellular matrix organization, and nuclear receptors. Furthermore, MCM10 over expression was negatively correlated with the level of immune cell infiltration in natural killer CD56 bright cells, follicular helper T cells, plasmacytoma dendritic cells, and dendritic cells.

CONCLUSION

MCM10 is an independent prognostic index of glioma patients, and the high expression of MCM10 suggests a poor prognosis; MCM10 expression is closely related to the immune cell infiltration of gliomas, and MCM10 may be related to drug resistance and development of gliomas.

TREM2 is associated with tumor immunity and implies poor prognosis in glioma

Triggering receptor expressed on myeloid cells 2 (TREM2) is expressed in myeloid cells of the central nervous system (CNS), which mediate the immunological response in a variety of diseases. Uncertain is the function of TREM2 in glioma and tumor immune responses. In this research, the expression patterns of TREM2 in glioma were analyzed, along with its prognostic value and functional roles. TREM2 expression is increased in glioblastomas, gliomas with a mesenchymal subtype, gliomas with wild-type isocitrate dehydrogenase, and gliomas without 1p/19q deletion, all of which suggest the aggressiveness and poor prognosis of gliomas. Gene ontology, KEGG, and Gene set variation analyses indicated that TREM2 may serve as an immune response mediator. However, the function of T cells against tumor cells was negatively correlated with TREM2, suggesting that TREM2 may suppress tumor immunity. Further investigation demonstrated a correlation between TREM2 expression and immune checkpoint expression. CIBERSORT research revealed a link between a higher TREM2 expression level and the enrichment of tumor-associated macrophages, especially M2 subtype. Single-cell analysis and multiple immunohistochemical staining results showed that microglia and macrophage cells expressed TREM2. Immunofluorescent staining indicated that knocking down the expression of TREM2 would result in a decrease in M2 polarization. TREM2 was discovered to be an independent prognostic factor in glioma. In conclusion, our findings revealed that TREM2 was significantly expressed in microglia and macrophage cells and was intimately associated with the tumor immune microenvironment. Thus, it is expected that small-molecule medications targeting TREM2 or monoclonal antibodies would enhance the efficacy of glioma immunotherapy.

New insights from the single-cell level: Tumor associated macrophages heterogeneity and personalized therapy

Tumor-associated macrophages (TAMs) are important immune cells in the tumor microenvironment, and their invasion in tumors is closely related to poor prognosis. Although TAMs are recognized as therapeutic targets, their heterogeneity makes studying tumor mechanism and developing drugs targeting TAMs difficult. The study of TAMs heterogeneity can be used to analyze the mechanism of tumor progression and drug resistance, and may provide possible treatment strategies for cancer patients. Single-cell RNA sequencing (scRNA-seq) can reveal the RNA expression profile for each TAM to distinguish heterogeneity, thereby providing a more efficient detection method and more accurate information for TAM-related studies. In this review, by summarizing the research progress in macrophage heterogeneity and other aspects of scRNA-seq over the past five years, we introduced the development of scRNA-seq technology and its application status in solid tumors, analyzed the advantages and selections of scRNA-seq in TAMs, and summarized the detailed specific research fields. To explore the mechanism of tumor progression and drug intervention from single cell level will provide new perspective for personalized treatment strategies targeting macrophages.

Integrative Single-Cell RNA-Seq and ATAC-Seq Analysis of Mesenchymal Stem/Stromal Cells Derived from Human Placenta

Mesenchymal stem/stromal cells derived from placenta (PMSCs) are an attractive source for regenerative medicine because of their multidifferentiation potential and immunomodulatory capabilities. However, the cellular and molecular heterogeneity of PMSCs has not been fully characterized. Here, we applied single-cell RNA sequencing (scRNA-seq) and assay for transposase-accessible chromatin sequencing (scATAC-seq) techniques to cultured PMSCs from human full-term placenta. Based on the inferred characteristics of cell clusters, we identify several distinct subsets of PMSCs with specific characteristics, including immunomodulatory-potential and highly proliferative cell states. Furthermore, integrative analysis of gene expression and chromatin accessibility showed a clearer chromatin accessibility signature than those at the transcriptional level on immunomodulatory-related genes. Cell cycle gene-related heterogeneity can be more easily distinguished at the transcriptional than the chromatin accessibility level in PMSCs. We further reveal putative subset-specific cis-regulatory elements regulating the expression of immunomodulatory- and proliferation-related genes in the immunomodulatory-potential and proliferative subpopulations, respectively. Moreover, we infer a novel transcription factor PRDM1, which might play a crucial role in maintaining immunomodulatory capability by activating PRDM1-regulon loop. Collectively, our study first provides a comprehensive and integrative view of the transcriptomic and epigenomic features of PMSCs, which paves the way for a deeper understanding of cellular heterogeneity and offers fundamental biological insight of PMSC subset-based cell therapy.

The MicroRNA-106a/20b Strongly Enhances the Antitumour Immune Responses of Dendritic Cells Pulsed with Glioma Stem Cells by Targeting STAT3

BACKGROUND

Evaluate the effect of the miRNA-106a/20b on the efficacy of DCs pulsed with GSCs in activating GSC-specific T cell responses.

METHODS

We cultured GSCs and prepared GSC antigen lysates by apoptosis. Then, immature DCs were pulsed with GSC antigen lysates in vitro. STAT3 levels in DCs were assessed by Western blotting, and the expression of CD80, CD86, and MHC-II was tested by fluorescence-activated cell sorting. The production and secretion of the cytokines IL-6, IL-12, TNF-α, and IL-10 in DCs induced by GSCs were determined by enzyme-linked immunosorbent assay. Finally, the cytotoxic functions of T cells stimulated by GSC-DC fusion cells transfected with a miR-106a/20b mimic in vitro and the antitumour activity in vivo were detected.

RESULTS

We found that the levels of miR-106a/20b were downregulated, but the expression of STAT3 was significantly upregulated. Simultaneously, the inhibition of STAT3 in the fusion cells by STAT3-specific siRNA caused significant upregulation of the expression of CD80, CD86, and MHC-II, and the secretion of the cytokines IL-6 and IL-12 was substantially increased, IL-10 was markedly decreased. These findings revealed that STAT3 is an important regulator of DC maturation. Furthermore, the interactional binding sites between the 3'-untranslated region (3'-UTR) of STAT3 mRNA and miR-106a/20b were predicted by bioinformatics and verified by a dual-luciferase assay. Moreover, the reduction in STAT3 levels in GSC-DCs enhanced the generation of CD8+ T cells and reduced the generation of Foxp3+ regulatory T cells. Meanwhile, the secretion of the T cell cytokine IFN-γ was significantly increased. Further research showed that DCs after miR-106a/20b-mimics transfection could promote the inhibition of GSC proliferation by T cells in vitro and suppress tumour growth in vivo.

CONCLUSIONS

This study indicted that the miR-106a/20b activation could be one of the important molecular mechanisms leading to enhance antitumour immune responses of GSC-mediated DCs, which downregulated the expression of STAT3 to alleviate its the inhibitory effect.

Clinical characterization and immunosuppressive regulation of CD161 (KLRB1) in glioma through 916 samples

Background

Glioblastoma is a paradigm of cancer‐associated immunosuppression, limiting the effects of immunotherapeutic strategies. Thus, identifying the molecular mechanisms underlying immune surveillance evasion is critical. Recently, the preferential expression of inhibitory natural killer (NK) cell receptor CD161 on glioma‐infiltrating cytotoxic T cells was identified. Focusing on the molecularly annotated, large‐scale clinical samples from different ethnic origins, the data presented here provide evidence of this immune modulator's essential roles in brain tumor biology.

Methods

Retrospective RNA‐seq data analysis was conducted in a cohort of 313 patients with glioma in the Chinese Glioma Genome Atlas (CGGA) database and 603 patients in The Cancer Genome Atlas (TCGA) database. In addition, single‐cell sequencing data from seven surgical specimens of glioblastoma patients and a model in which patient‐derived glioma stem cells were cocultured with peripheral lymphocytes, were used to analyze the molecular evolution process during gliomagenesis.

Results

CD161 was enriched in high‐grade gliomas and isocitrate dehydrogenase (IDH)‐wildtype glioma. CD161 acted as a potential biomarker for the mesenchymal subtype of glioma and an independent prognostic factor for the overall survival (OS) of patients with glioma. In addition, CD161 played an essential role in inhibiting the cytotoxicity of T cells in glioma patients. During the process of gliomagenesis, the expression of CD161 on different lymphocytes dynamically evolved.

Conclusion

The expression of CD161 was closely related to the pathology and molecular pathology of glioma. Meanwhile, CD161 promoted the progression and evolution of gliomas through its unique effect on T cell dysfunction. Thus, CD161 is a promising novel target for immunotherapeutic strategies in glioma treatment.

Applications of single-cell sequencing in cancer research: progress and perspectives

Single-cell sequencing, including genomics, transcriptomics, epigenomics, proteomics and metabolomics sequencing, is a powerful tool to decipher the cellular and molecular landscape at a single-cell resolution, unlike bulk sequencing, which provides averaged data. The use of single-cell sequencing in cancer research has revolutionized our understanding of the biological characteristics and dynamics within cancer lesions. In this review, we summarize emerging single-cell sequencing technologies and recent cancer research progress obtained by single-cell sequencing, including information related to the landscapes of malignant cells and immune cells, tumor heterogeneity, circulating tumor cells and the underlying mechanisms of tumor biological behaviors. Overall, the prospects of single-cell sequencing in facilitating diagnosis, targeted therapy and prognostic prediction among a spectrum of tumors are bright. In the near future, advances in single-cell sequencing will undoubtedly improve our understanding of the biological characteristics of tumors and highlight potential precise therapeutic targets for patients.

Characteristics of hypoxic tumor microenvironment in non-small cell lung cancer, involving molecular patterns and prognostic signature

Background

The mechanisms of hypoxia or immune microenvironment in cancer have been studied respectively, but the role of hypoxia immune microenvironment in non-small cell lung cancer (NSCLC) still needs further exploration.

Methods

By applying the K-means algorithm, 1,121 patients with NSCLC were divided into three categories. We evaluated the constructed signature in order to link it with the prognosis, which was constructed by univariate and least absolute shrinkage operator (LASSO) Cox regression analysis.

Results

A total of three clusters were obtained by clustering five Gene Expression Omnibus (GEO) data sets. Gene Set Variation Analysis (GSVA) and immune infiltration analysis were performed to explore the biological behavior. Cluster one presented an activated state of oncogenic pathways, and compared with the other two clusters, the median risk score was the highest, which was the reason for its poor survival. Cluster three showed that the immune pathway was active and the median risk score was the lowest, so the survival was the best. However, cluster two presented a state in which both immune and matrix pathways were activate. This was manifested as mutual antagonism, and its risk score was in the middle. Its survival was in the middle.

Conclusions

This work revealed the role of hypoxia related genes (HRGs) modification in tumor microenvironment, which was conducive to our comprehensive analysis of the prognosis of NSCLC, and provided direction and guidance for clinical immunotherapy.

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.