The hypoxic tumor microenvironment was reported to be involved in different tumorigenesis mechanisms of triple-negative breast cancer (TNBC), such as invasion, immune evasion, chemoresistance, and metastasis. However, a systematic analysis of the prognostic prediction models based on multiple hypoxia-related genes (HRGs) has not been established in TNBC before in the literature. We aimed to develop and verify a hypoxia gene signature for prognostic prediction in TNBC patients.
The RNA sequencing profiles and clinical data of TNBC patients were generated from the TCGA, GSE103091, and METABRIC databases. The TNBC-specific differential HRGs (dHRGs) were obtained from differential expression analysis of hypoxia cultured TNBC cell lines compared with normoxic cell lines from the GEO database. Non-negative matrix factorization (NMF) method was then performed on the TNBC patients using the dHRGs to explore a novel molecular classification on the basis of the dHRG expression patterns. Prognosis-associated dHRGs were identified by univariate and multivariate Cox regression analysis to establish the prognostic risk score model.
Based on the expressions of 205 dHRGs, all the patients in the TCGA training cohort were categorized into two subgroups, and the patients in Cluster 1 demonstrated worse OS than those in Cluster 2, which was validated in two independent cohorts. Additionally, the effects of somatic copy number variation (SCNV), somatic single nucleotide variation (SSNV), and methylation level on the expressions of dHRGs were also analyzed. Then, we performed Cox regression analyses to construct an HRG-based risk score model (3-gene dHRG signature), which could reliably discriminate the overall survival (OS) of high-risk and low-risk patients in TCGA, GSE103091, METABRIC, and BMCHH (qRT-PCR) cohorts.
In this study, a robust predictive signature was developed for patients with TNBC, indicating that the 3-gene dHRG model might serve as a potential prognostic biomarker for TNBC.