A combined hypoxia and immune gene signature for predicting survival and risk stratification in triple-negative breast cancer

Background: Increasing evidence showed that the clinical significance of the interaction between hypoxia and immune status in tumor microenvironment. However, reliable biomarkers based on the hypoxia and immune status in triple-negative breast cancer (TNBC) have not been well established. This study aimed to explore a gene signature based on the hypoxia and immune status for predicting prognosis, risk stratification, and individual treatment in TNBC. Methods: Hypoxia-related genes (HRGs) and Immune-related genes (IRGs) were identified using the weighted gene co-expression network analysis (WGCNA) method and the single-sample gene set enrichment analysis (ssGSEA Z-score) with the transcriptomic profiles from Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) cohort. Then, prognostic hypoxia and immune based genes were identified in TNBC patients from the METABRIC (N = 221), The Cancer Genome Atlas (TCGA) (N = 142), and GSE58812 (N = 107) using univariate cox regression model. A robust hypoxia-immune based gene signature for prognosis was constructed using the least absolute shrinkage and selection operator (LASSO) method. Based on the cross-cohort prognostic hypoxia–immune related gene signature, a comprehensive index of hypoxia and immune was developed and two risk groups with distinct hypoxia–immune status were identified. The prognosis value, hypoxia and immune status, and therapeutic response in different risk groups were analyzed. Furthermore, a nomogram was constructed to predict the prognosis for individual patients, and an independent cohort from the gene expression omnibus (GEO) database was used for external validation. Results: Six cross-cohort prognostic hypoxia–immune related genes were identified to establish the comprehensive index of hypoxia and immune. Then, patients were clustered into high- and low-risk groups based on the hypoxia–immune status. Patients in the high-risk group showed poorer prognoses to their low-risk counterparts, and the nomogram we constructed yielded favorable performance to predict survival and risk stratification. Besides, the high-risk group had a higher expression of hypoxia-related genes and correlated with hypoxia status in tumor microenvironment. The high-risk group had lower fractions of activated immune cells, and exhibited lower expression of immune checkpoint markers. Furthermore, the ratio of complete response (CR) was greatly declined, and the ratio of breast cancer related events were significantly elevated in the high-risk group. Conclusion: The hypoxia–immune based gene signature we constructed for predicting prognosis was developed and validated, which may contribute to the optimization of risk stratification for prognosis and personalized treatment in TNBC patients.