A novel hypoxia- and lactate metabolism-related signature to predict prognosis and immunotherapy responses for breast cancer by integrating machine learning and bioinformatic analyses.

Breast cancer is the most common cancer worldwide. Hypoxia and lactate metabolism are hallmarks of cancer. This study aimed to construct a novel hypoxia- and lactate metabolism-related gene signature to predict the survival, immune microenvironment, and treatment response of breast cancer patients. RNA-seq and clinical data of breast cancer from The Cancer Genome Atlas database and Gene Expression Omnibus were downloaded. Hypoxia- and lactate metabolism-related genes were collected from publicly available data sources. The differentially expressed genes were identified using the "edgeR" R package. Univariate Cox regression, random survival forest (RSF), and stepwise multivariate Cox regression analyses were performed to construct the hypoxia-lactate metabolism-related prognostic model (HLMRPM). Further analyses, including functional enrichment, ESTIMATE, CIBERSORTx, Immune Cell Abundance Identifier (ImmuCellAI), TIDE, immunophenoscore (IPS), pRRophetic, and CellMiner, were performed to analyze immune status and treatment responses. We identified 181 differentially expressed hypoxia-lactate metabolism-related genes (HLMRGs), 24 of which were valuable prognostic genes. Using RSF and stepwise multivariate Cox regression analysis, five HLMRGs were included to establish the HLMRPM. According to the medium-risk score, patients were divided into high- and low-risk groups. Patients in the high-risk group had a worse prognosis than those in the low-risk group (P < 0.05). A nomogram was further built to predict overall survival (OS). Functional enrichment analyses showed that the low-risk group was enriched with immune-related pathways, such as antigen processing and presentation and cytokine-cytokine receptor interaction, whereas the high-risk group was enriched in mTOR and Wnt signaling pathways. CIBERSORTx and ImmuCellAI showed that the low-risk group had abundant anti-tumor immune cells, whereas in the high-risk group, immunosuppressive cells were dominant. Independent immunotherapy datasets (IMvigor210 and GSE78220), TIDE, IPS and pRRophetic analyses revealed that the low-risk group responded better to common immunotherapy and chemotherapy drugs. We constructed a novel prognostic signature combining lactate metabolism and hypoxia to predict OS, immune status, and treatment response of patients with breast cancer, providing a viewpoint for individualized treatment.

Li J ，Qiao H ，Wu F ，Sun S ，Feng C ，Li C ，Yan W ，Lv W ，Wu H ，Liu M ，Chen X ，Liu X ，Wang W ，Cai Y ，Zhang Y ，Zhou Z ，Zhang Y ，Zhang S ... -  《Frontiers in Immunology》

Comprehensive analysis of nicotinamide metabolism-related signature for predicting prognosis and immunotherapy response in breast cancer.

Breast cancer (BC) is the most common malignancy among women. Nicotinamide (NAM) metabolism regulates the development of multiple tumors. Herein, we sought to develop a NAM metabolism-related signature (NMRS) to make predictions of survival, tumor microenvironment (TME) and treatment efficacy in BC patients. Transcriptional profiles and clinical data from The Cancer Genome Atlas (TCGA) were analyzed. NAM metabolism-related genes (NMRGs) were retrieved from the Molecular Signatures Database. Consensus clustering was performed on the NMRGs and the differentially expressed genes between different clusters were identified. Univariate Cox, Lasso, and multivariate Cox regression analyses were sequentially conducted to develop the NAM metabolism-related signature (NMRS), which was then validated in the International Cancer Genome Consortium (ICGC) database and Gene Expression Omnibus (GEO) single-cell RNA-seq data. Further studies, such as gene set enrichment analysis (GSEA), ESTIMATE, CIBERSORT, SubMap, and Immunophenoscore (IPS) algorithm, cancer-immunity cycle (CIC), tumor mutation burden (TMB), and drug sensitivity were performed to assess the TME and treatment response. We identified a 6-gene NMRS that was significantly associated with BC prognosis as an independent indicator. We performed risk stratification according to the NMRS and the low-risk group showed preferable clinical outcomes (P < 0.001). A comprehensive nomogram was developed and showed excellent predictive value for prognosis. GSEA demonstrated that the low-risk group was predominantly enriched in immune-associated pathways, whereas the high-risk group was enriched in cancer-related pathways. The ESTIMATE and CIBERSORT algorithms revealed that the low-risk group had a higher abundance of anti-tumor immunocyte infiltration (P < 0.05). Results of Submap, IPS, CIC, TMB, and external immunotherapy cohort (iMvigor210) analyses showed that the low-risk group were indicative of better immunotherapy response (P < 0.05). The novel signature offers a promising way to evaluate the prognosis and treatment efficacy in BC patients, which may facilitate clinical practice and management.

Cui H ，Ren X ，Dai L ，Chang L ，Liu D ，Zhai Z ，Kang H ，Ma X ... -  《Frontiers in Immunology》

Identification of Hypoxia-Related Subtypes, Establishment of Prognostic Models, and Characteristics of Tumor Microenvironment Infiltration in Colon Cancer.

Background: Immunotherapy is a treatment that can significantly improve the prognosis of patients with colon cancer, but the response to immunotherapy is different in patients with colon cancer because of the heterogeneity of colon carcinoma and the complex nature of the tumor microenvironment (TME). In the precision therapy mode, finding predictive biomarkers that can accurately identify immunotherapy-sensitive types of colon cancer is essential. Hypoxia plays an important role in tumor proliferation, apoptosis, angiogenesis, invasion and metastasis, energy metabolism, and chemotherapy and immunotherapy resistance. Thus, understanding the mechanism of hypoxia-related genes (HRGs) in colon cancer progression and constructing hypoxia-related signatures will help enrich our treatment strategies and improve patient prognosis. Methods: We obtained the gene expression data and corresponding clinical information of 1,025 colon carcinoma patients from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) databases, respectively. We identified two distinct hypoxia subtypes (subtype A and subtype B) according to unsupervised clustering analysis and assessed the clinical parameters, prognosis, and TME cell-infiltrating characteristics of patients in the two subtypes. We identified 1,132 differentially expressed genes (DEGs) between the two hypoxia subtypes, and all patients were randomly divided into the training group (n = 513) and testing groups (n = 512). Following univariate Cox regression with DEGs, we construct the prognostic model (HRG-score) including six genes (S1PR3, ETV5, CD36, FOXC1, CXCL10, and MMP12) through the LASSO-multivariate cox method in the training group. We comprehensively evaluated the sensitivity and applicability of the HRG-score model from the training group and the testing group, respectively. We explored the correlation between HRG-score and clinical parameters, tumor microenvironment, cancer stem cells (CSCs), and MMR status. In order to evaluate the value of the risk model in clinical application, we further analyzed the sensitivity of chemotherapeutics and immunotherapy between the low-risk group and high-risk group and constructed a nomogram for improving the clinical application of the HRG-score. Result: Subtype A was significantly enriched in metabolism-related pathways, and subtype B was significantly enriched in immune activation and several tumor-associated pathways. The level of immune cell infiltration and immune checkpoint-related genes, stromal score, estimate score, and immune dysfunction and exclusion (TIDE) prediction score was significantly different in subtype A and subtype B. The level of immune checkpoint-related genes and TIDE score was significantly lower in subtype A than that in subtype B, indicating that subtype A might benefit from immune checkpoint inhibitors. Finally, an HRG-score signature for predicting prognosis was constructed through the training group, and the predictive capability was validated through the testing group. The survival analysis and correlation analysis of clinical parameters revealed that the prognosis of patients in the high-risk group was significantly worse than that in the low-risk group. There were also significant differences in immune status, mismatch repair status (MMR), and cancer stem cell index (CSC), between the two risk groups. The correlation analysis of risk scores with IC50 and IPS showed that patients in the low-risk group had a higher benefit from chemotherapy and immunotherapy than those in the high-risk group, and the external validation IMvigor210 demonstrated that patients with low risk were more sensitive to immunotherapy. Conclusion: We identified two novel molecular subgroups based on HRGs and constructed an HRG-score model consisting of six genes, which can help us to better understand the mechanisms of hypoxia-related genes in the progression of colon cancer and identify patients susceptible to chemotherapy or immunotherapy, so as to achieve precision therapy for colon cancer.

Wang C ，Tang Y ，Ma H ，Wei S ，Hu X ，Zhao L ，Wang G ... -  《Frontiers in Genetics》

Novel hypoxia- and lactate metabolism-related molecular subtyping and prognostic signature for colorectal cancer.

Colorectal cancer (CRC) is a serious global health burden because of its high morbidity and mortality rates. Hypoxia and massive lactate production are hallmarks of the CRC microenvironment. However, the effects of hypoxia and lactate metabolism on CRC have not been fully elucidated. This study aimed to develop a novel molecular subtyping based on hypoxia-related genes (HRGs) and lactate metabolism-related genes (LMRGs) and construct a signature to predict the prognosis of patients with CRC and treatment efficacy. Bulk and single-cell RNA-sequencing and clinical data of CRC were downloaded from the TCGA and GEO databases. HRGs and LMRGs were obtained from the Molecular Signatures Database. The R software package DESeq2 was used to perform differential expression analysis. Molecular subtyping was performed using unsupervised clustering. A predictive signature was developed using univariate Cox regression, random forest model, LASSO, and multivariate Cox regression analyses. Finally, the sensitivity of tumor cells to chemotherapeutic agents before and after hypoxia was verified using in vitro experiments. We classified 575 patients with CRC into three molecular subtypes and were able to distinguish their prognoses clearly. The C1 subtype, which exhibits high levels of hypoxia, has a low proportion of CD8 + T cells and a high proportion of macrophages. The expression of immune checkpoint genes is generally elevated in C1 patients with severe immune dysfunction. Subsequently, we constructed a predictive model, the HLM score, which effectively predicts the prognosis of patients with CRC and the efficacy of immunotherapy. The HLM score was validated in GSE39582, GSE106584, GSE17536, and IMvigor210 datasets. Patients with high HLM scores exhibit high infiltration of CD8 + exhausted T cells (Tex), especially terminal Tex, and oxidative phosphorylation (OXPHOS)-Tex in the immune microenvironment. Finally, in vitro experiments confirmed that CRC cell lines were less sensitive to 5-fluorouracil, oxaliplatin, and irinotecan under hypoxic conditions. We constructed novel hypoxia- and lactate metabolism-related molecular subtypes and revealed their immunological and genetic characteristics. We also developed an HLM scoring system that could be used to predict the prognosis and efficacy of immunotherapy in patients with CRC.

Huang A ，Sun Z ，Hong H ，Yang Y ，Chen J ，Gao Z ，Gu J ... -  《Journal of Translational Medicine》

Tumor immune microenvironment-based clusters in predicting prognosis and guiding immunotherapy in breast cancer.

Liu Y ，He X ，Yang YI 《-》