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》

A New Genetic Signature of Lactate Metabolism-Associated Genes Predicting Clinically Distinctive Features and Tumor Microenvironment in Colorectal Cancer.

Wang K ，Lou Y ，Tao Z 《-》

Development of a novel hypoxia-immune-related LncRNA risk signature for predicting the prognosis and immunotherapy response of colorectal cancer.

Colorectal cancer (CRC) is one of the most common digestive system tumors worldwide. Hypoxia and immunity are closely related in CRC; however, the role of hypoxia-immune-related lncRNAs in CRC prognosis is unknown. Data used in the current study were sourced from the Gene Expression Omnibus and The Cancer Genome Atlas (TCGA) databases. CRC patients were divided into low- and high-hypoxia groups using the single-sample gene set enrichment analysis (ssGSEA) algorithm and into low- and high-immune groups using the Estimation of STromal and Immune cells in MAlignant Tumours using Expression data (ESTIMATE) algorithm. Differentially expressed lncRNAs (DElncRNAs) between low- and high-hypoxia groups, low- and high-immune groups, and tumor and control samples were identified using the limma package. Hypoxia-immune-related lncRNAs were obtained by intersecting these DElncRNAs. A hypoxia-immune-related lncRNA risk signature was developed using univariate Cox regression and least absolute shrinkage and selection operator (LASSO) analyses. The tumor microenvironments in the low- and high-risk groups were evaluated using ssGSEA, ESTIMATE, and the expression of immune checkpoints. The therapeutic response in the two groups was assessed using TIDE, IPS, and IC50. A ceRNA network based on signature lncRNAs was constructed. Finally, we used RT-qPCR to verify the expression of hypoxia-immune-related lncRNA signatures in normal and cancer tissues. Using differential expression analysis, and univariate Cox and LASSO regression analyses, ZNF667-AS1, LINC01354, LINC00996, DANCR, CECR7, and LINC01116 were selected to construct a hypoxia-immune-related lncRNA signature. The performance of the risk signature in predicting CRC prognosis was validated in internal and external datasets, as evidenced by receiver operating characteristic curves. In addition, we observed significant differences in the tumor microenvironment and immunotherapy response between low- and high-risk groups and constructed a CECR7-miRNA-mRNA regulatory network in CRC. Furthermore, RT-qPCR results confirmed that the expression patterns of the six lncRNA signatures were consistent with those in TCGA-CRC cohort. Our study identified six hypoxia-immune-related lncRNAs for predicting CRC survival and sensitivity to immunotherapy. These findings may enrich our understanding of CRC and help improve CRC treatment. However, large-scale long-term follow-up studies are required for verification.

Luan L ，Dai Y ，Shen T ，Yang C ，Chen Z ，Liu S ，Jia J ，Li Z ，Fang S ，Qiu H ，Cheng X ，Yang Z ... -  《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》

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》