Jianxin Zhong
Xi Shen
Junjie Zhou3†- 1Key Laboratory of Carcinogenesis and Translational Research (Ministry of Education), Department of Breast Oncology, Peking University Cancer Hospital and Institute, Beijing, China
- 2Department of Head and Neck Oncology and Department of Radiation Oncology, Cancer Center, West China Hospital, Sichuan University, Chengdu, Sichuan, China
- 3Department of Pathology, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
- 4Department of Thyroid and Breast Surgery, Wuhan Fourth Hospital, Wuhan, China
- 5Department of Thoracic Surgery, Wuhan Fourth Hospital, Wuhan, China
Background: Hypoxia is involved in tumor biological processes and disease progression. Ferroptosis, as a newly discovered programmed cell death process, is closely related to breast cancer (BC) occurrence and development. However, reliable prognostic signatures based on a combination of hypoxia and ferroptosis in BC have not been developed.
Method: We set The Cancer Genome Atlas (TCGA) breast cancer cohort as training set and the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) BC cohort as the validation set. Least Absolute Shrinkage and Selection Operator (LASSO) and COX regression approaches were used to construct ferroptosis-related genes (FRGs) and hypoxia-related genes (HRGs) prognostic signature (HFRS). The CIBERSORT algorithm and ESTIMATE score were used to explore the relationship between HFRS and tumor immune microenvironment. Immunohistochemical staining was used to detect protein expression in tissue samples. A nomogram was developed to advance the clinical application of HFRS signature.
Results: Ten ferroptosis-related genes and hypoxia-related genes were screened to construct the HFRS prognostic signature in TCGA BC cohort, and the predictive capacity was verified in METABRIC BC cohort. BC patients with high-HFRS had shorter survival time, higher tumor stage, and a higher rate of positive lymph node. Moreover, high HFRS was associated with high hypoxia, ferroptosis, and immunosuppression status. A nomogram that was constructed with age, stage, and HFRS signature showed a strong prognostic capability to predict overall survival (OS) for BC patients.
Conclusion: We developed a novel prognostic model with hypoxia and ferroptosis-related genes to predict OS, and characterize the immune microenvironment of BC patients, which might provide new cures for clinical decision-making and individual treatment of BC patients.
1 Introduction
Breast cancer has been the most prevalent tumor in women worldwide and is the leading cause of cancer-related death among women with malignant diseases (1, 2). With the advancement of diagnosis and treatment, the overall survival of primary breast cancer has been greatly improved; however, advanced breast cancer is still refractory, and some patients who were diagnosed with distant metastasis that lose the chance of surgery therapy were required more efficient target drugs to improve their prognosis. Breast cancer is a heterogeneous tumor; individual treatment and biological feature depiction of each patient is a field that calls for exploration (3). Although age at diagnosis, tumor stage, and histological grade are considered as prognostic factors, there are few reliable biomarkers based on personal gene expression pattern to facilitate clinical assessment (4–7). Therefore, it is important to discover novel prognostic factors and potential therapeutic targets for individual treatment.
Hypoxia is a feature of solid tumors generated since the supply could not meet the consumption of oxygen under rapid tumor proliferation, and form a tumor microenvironment (8, 9). Hypoxia could induce tumor angiogenesis, cell proliferation, metastasis, and invasion and promote tumor immune suppression and escape, while reducing apoptosis, differentiation, and ferroptosis to accelerate tumor progression (8, 10–12). Ferroptosis, as a newly discovered regulatory cell death, was closely related to tumor development (13, 14). Increasing evidence revealed that targeting ferroptosis induced treatment response in BC (15–17). More importantly, hypoxia has been proven to participate in the regulation of ferroptosis (18, 19). Some studies demonstrated that hypoxia blocks ferroptosis in hepatocellular carcinoma, and HIF-1α-induced lncRNA PMAN promoted gastric cancer peritoneal dissemination by inhibiting ferroptosis (20, 21). However, the association between hypoxia and ferroptosis in BC has not been reported yet.
Immune microenvironment regulation was critical in tumor progression, which has been wildly verified to be associated with hypoxia status (22, 23). Some studies have reported that hypoxia inducible factor-1α (HIF-1α) increased PD-L1 expression and antigen non-specific T-cell suppression, and promote the differentiation of MDSC to immune suppressive TAM in various kinds of tumors including breast cancer (24–27). HIF-1α could negatively regulate the functions of CD4+ and CD8+ T lymphocytes, and depletion of HIF-1α enhanced T cell response (28, 29). Interestingly, many evidences indicated that ferroptosis was also involved in the regulation of the immune microenvironment and immunotherapy resistance in cancers (30–33). Thus, there was potential interaction between hypoxia and ferroptosis, and either of them was associated with immune microenvironment regulation in cancers.
Given that hypoxia and ferroptosis are related to breast cancer prognosis, there were few studies that reported the crosstalk between hypoxia and ferroptosis, and no prognostic signature has been established in BC for risk stratification and immune microenvironment profiling. This study firstly combined ferroptosis-related genes (FRGs) with hypoxia related genes (HRGs) to construct a prognostic signature HFRS to predict BC prognosis and immune status.
2 Materials and methods
2.1 Data acquisition
The mRNA expression data and corresponding clinicopathological information of BC patients were obtained from the TCGA and METABRIC websites. A total of 1075 BC patients from the TCGA database were enrolled in the training cohort, and 1399 patients with completed clinical information from the METABRIC database were included as a validation cohort, after excluding patients who lacked tumor stage and survival information. By intersecting the ferroptosis-related genes in the FerrDb database and Molecular Signatures Database3, 47 FRGs were retrieved; 243 HRGs were downloaded from the hypoxia-related gene set “winter_hypoxia_metagenes” in Molecular Signatures Database 3 (MSigDB: https://www.gsea-msigdb.org/gsea/msigdb). Gene expression data from these databases were normalized by the R package “limma”. Supplementary Table 1 shows the clinicopathological information of TCGA and METABRIC cohort in this study. Identification of PAM50 subtypes of all the patients was performed by the ‘genefu’ R package based on gene expression profiles.
2.2 Development of the HFRS
Univariate COX regression analysis was used to screen prognostic genes among 47 FRGs and 243 HRGs in the TCGA cohort. Then, 15 FRGs and HRGs significantly associated with prognosis identified in univariate regression analysis (p< 0.001) of BC patients were input into the Least Absolute Shrinkage and Selection Operator (LASSO) The COX regression model was used to identify the critical genes and the corresponding regression coefficient by using the R package “glmnet” (Friedman et al., 2010). We constructed a hypoxia and ferroptosis prognostic signature (HFRS) for the BC patients with 10 FRGs and HRGs selected by LASSO COX analysis. HFRS scores were calculated for all patients according to the formula: lambda.min = 0.0027
Where xi is the expression level of each FRG or HRG and Coefi is the coefficient.
Then the R package “survminer” was used to calculate the optimal cut-off value (this is an outcome-oriented method providing a value of a cut-off point that corresponds to the most significant relation with survival) and the patients were divided into two subgroups (low-HFPS and high-HFRS group) according to the optimal cut-off value.
2.3 Functional analysis
Gene Set Enrichment Analysis (GSEA) was used to investigate the pathways enriched in the low-HFRS subgroup and high-HFRS subgroup and identified significant enrichment pathways with normalized enrichment score >1, nominal p< 0.05, and false discovery rate q< 0.25. Differentially expressed genes (DEGs) between the high HFRS and low-HFRS groups were obtained using the R package “Deseq2” (| log2(Fold change) |>1 and adjust p<0.05) and were input into “ClusterProfiler” R package for functional enrichment and pathway analysis, including the Kyoto Encyclopedia of Genes and Genomes Pathway (KEGG pathway) and Gene Ontology (GO) analysis. The FRGs and HRGs significantly associated with the prognosis (p< 0.001) of BC patients were subjected to construct a protein–protein interaction (PPI) network by MetaScape (https://metascape.org/).
2.4 Analysis of immune cell infiltration
To investigate the difference in immune infiltration status between patients in the high- and low-HFRS group, the CIBERSORT algorithm was used to analyze the immune cell type-specific gene expression profiles of BC patients with the LM22 signatures. Moreover, we also used the ESTIMATE method to calculate immune cell characteristics for BC patients. We downloaded the Immunophenoscores (IPS) of BC patients of the TCGA cohort from the TCIA database (https://tcia.at/) to predict the sensitivity of immune therapy of the high- and low-HFRS groups.
2.5 Analysis of genetic alteration in BC patients
The R package “Maftools” was used to visualize the single nucleotide variation (SNV) profile of TCGA BC patients with mutation data, and compared the different mutation patterns between the high- and low-HFRS groups. The copy number variation (CNV) of HFRS genes and their correlation with mRNA expression were analyzed in the GSCAlite (http://bioinfo.life.hust.edu.cn/web/GSCALite/) website.
2.6 Immunohistochemistry staining of HFRS gene protein expression in tissues
We collected 20 pairs of BC tissues and adjacent normal breast tissues from Wuhan Pu-Ai Hospital, which was approved by the ethics committees of Pu-Ai Hospital (No. KY2022-050-02). The BC tissues and adjacent normal breast tissues were fixed with 10% formalin, embedded by paraffin, and sectioned; then we selected the optimal tissue sections for degreasing and immunohistochemistry staining. Protein expression levels were evaluated semi-quantitatively following the Allred scoring system guidelines and scored separately by two qualified pathologists (34). Then, the sections were scanned to obtain high-resolution digital images using a 3DHISTECH scanner (Pannoramic, TaiBei). Antibodies used in this study are as follows: BTG1 (Proteintech, Cat No. 14879-1-AP). SLC16A2 (Abcam, ab192828). In addition, immunohistochemical staining images of the remaining eight HFRS genes were obtained from the Human Protein Atlas (HPA).
2.7 Statistical analysis
In this study, all statistical analysis was conducted by R 4.1.1. Univariate COX regression was used to identify independent prognostic risk factors. Multivariate COX regression analysis was used to construct a nomogram to predict OS for BC patients. The predictive efficiency of the nomogram was verified in METABRIC cohorts. The R package ‘rms’ was used in the construction and validation of the nomogram. In addition, The ROC curve and AUC were used to analyze the prognosis predictive accuracy of nomogram and other prognostic factors via R package “timeROC”. For descriptive statistics, mean ± standard deviation (SD) or median (range) was used for continuous variables; Student’s t-test and Mann–Whitney U test were used to analyze the difference between two groups of continuous variables. Fisher exact test or Wilcoxon’s test was used to compare the difference of clinical features of categorical variables between two groups when appropriate. Two-tailed p< 0.05 was considered as statistically significant.
3 Results
3.1 Construction of the HFRS in the TCGA cohort
The study design is shown in the flow chart (Figure 1), The GO pathways analysis conducted in the Metascape website showed that these genes were enriched in hypoxia, metabolites, energy, and oxidative stress-related pathways (Figure 2A). Univariate COX regression was used to screen for hypoxia-related prognostic genes (HRGs) and ferroptosis-related prognostic genes (FRGs) in the TCGA cohort. In the condition of p< 0.001, there were 15 prognostic HRGs and FRGs that were significantly associated with the prognosis of BC patients (Supplementary Table 2). Subsequently, 15 prognostic FRGs and HRGs were subjected to the LASSO-Cox regression analysis, and we screened 10 genes (BTG1, CCT6A, KRT14, P4HA2, PGK1, SLC16A2, TPD52 and STC2 as HRGs, CISD1 as FRGs, and TF as both HRG and FRG) to construct a hypoxia and ferroptosis prognostic combined signature (HFRS) (Figure 2B). HFRS scored the BC patients in TCGA cohort and the patients were further divided into the high-HFRS group (n = 414) and low-HFRS group (n = 661) according to the optimal cut-off value. Kaplan–Meier curves and log-rank test showed that patients in the high-HFRS group had a significantly worse prognosis than the low-HFRS group (p< 0.001) (Figure 2C). The distributions of the survival status and HFRS score are shown in Figure 2D. The ROC curves indicated an efficient prognostic predictive capacity of HFRS for the overall survival of BC patients; the AUC of 1, 3, and 5 years were 0.72, 0.73, and 0.72, respectively (Figure 2E).
Figure 1 Study flowchart.
Figure 2 Construction and validation of hypoxia and ferroptosis-related gene risk signature(HFRS) in BC patients. (A) The barplot shows the enrichment of Go pathways significantly associated with prognostic HRGs and FRGs. (B) LASSO regression analysis identified 10 hub genes to construct HFRS signature. (C) Kaplan–Meier curves show the significant difference in overall survival between high- and low-HFRS groups in TCGA cohort. (D) The ranked dot plot indicates the HFRS_score distribution; scatter plot present the patients’ survival status in TCGA cohort. (E) The ROC curves of HFRS for predicting 1, 3, and 5 years overall survival in TGCA cohort. (F) The K-M curves show the significant difference in overall survival between high- and low-HFRS groups in the METABRIC cohort. (G) The ranked dot plot indicates the HFRS_score distribution; scatter plot presents the patients’ survival status in the METABRIC cohort (H) The ROC curves of HFRS for predicting 1, 3, and 5 years overall survival in METABRIC cohort. (I) The PCA plot based on HFRS gene expression show the distinct subgroups of TCGA BC cohort. (J) The PCA plot based on HFRS gene expression divided the METABRIC BC cohort to two subgroups.
3.2 Validation of the HFRS in the METABRIC cohort
HFRS scores of BC patients in the METABRIC cohort were calculated by the same signature model, and the patients were divided into the low-HFRS group (n = 639) and high-HFRS group (n = 760) according to the optimal cut-off value. The results of the METABRIC cohort are generally consistent with those of the TCGA cohort; patients in the high-HFRS group had significantly poorer prognosis (Figure 2F). The distribution of survival status and HFRS score also indicated that patients with higher HFRS scores had shorter overall survival time and higher mortality (Figure 2G). The ROC curve showed that HFRS score also had strong predictive power in the METABRIC cohort. The AUCs were: 0.73 (1 year), 0.68 (2 years), and 0.65 (3 years) (Figure 2H). In addition, the PCA suggested that the BC patients could be distinctively clustered by PCA according to HFRS as well (Figures 2I, J).
3.3 Prognostic analysis and genetic alteration of the 10 FHRS genes
Univariate COX regression analysis of 10 HFRS genes showed that BTG1, KRT14, STC2, and TF were protective factors in BC (0< Hazard Ratio (HR)< 1; p< 0.001), while CCT6A, P4HA2, PGK1, TPD52, SLC16A2, and CISD1 were risk factors (HR > 1; p< 0.001) for the overall survival of BC patients (Figure 3A). In addition, the heat map shows the differential expressions of 10 HFRS genes in TCGA breast cancer samples. The expressions of PGK1, CCT6A, P4HA2, TPD52, SLC16A2, and CISD1 increased with the HFRS scores while the expression of BTG1, KRT14, STC2, and TF decreased with the HFRS scores. In addition, the distribution of HFRS gene expressions was also associated with the tumor stage of BC patients (Figure 3B). We further explored the genetic alteration of HFRS genes in cancers. We investigated the single nucleotide variations (SNVs) of HFRS genes in different cancers and observed that some genes (CCT6A, TF, KRT14, P4HA2, PGK1, SLC16A2, and STC2) were frequently mutated in COAD, STAD, UCEC, and SKCM (Figure 3C). In BRCA, BLCA, HNSC, and LUAD, the copy number variations (CNVs) of some genes were positively correlated with mRNA level (Figure 3D). Moreover, we found that TPD52, TF, and CCT6A were more frequently heterozygous amplified while STC2, CISD1, and P4HA2 were more likely to occur in heterozygous deletion in cancers (Figure 3E). In contrast, homozygous amplification and deletion of HFRS genes were very rare in cancers (Figure 3F).
Figure 3 Prognosis value and expression of HFRS hub genes. (A) Forest plot shows the prognostic value of 10 prognostic genes in signature. (B) Heat map shows the relationship between mRNA expression levels of 10 HFRS genes, HFRS score, and tumor stage features in the TCGA cohort. (C) The mutation frequencies of 10 HFRS genes in pan-cancer. (D) The correlation between CNV and mRNA expression of HFRS genes in pan-cancer using Pearson analysis. The size of bubble indicated the -log10 (FDR) value. (E, F) The profiles of heterozygous (E) and homozygous (F) amplification/deletion of 10 HFRS genes in cancers.
3.4 Clinical effects of HFRS on breast cancer patients
To investigate whether the HFRS score was associated with clinicopathological characteristics of BC patients, we compared the HFRS score of BC patients according to different clinical subgroups such as BC patients of age >65; positive lymph node status and TNM stage III/IV had significantly higher HFRS scores. The result suggested that HFRS related to clinical features of BC patients, and might reflect tumor burden (Figures 4A, D, G). In addition, we identified overall survival in different clinical subgroups using Kaplan–Meier curves. The result suggested that BC patients with low HFRS had better OS than the patients with high HFRS in both age >65 and ≤65 subgroups (Figures 4B, C), and the same results were observed in negative and positive lymph node subgroups (Figures 4E, F), stage I/II, and III/IV subgroups (Figures 4H, I). Moreover, except for the age subgroup (p = 0.26), the similar HFRS distribution in these subgroups was observed in the subgroups of the METABRIC BC cohort (Figures 4J–L). In addition, we found that patients with higher tumor grade had significantly higher HFRS (Figure 4M). These results indicated that the HFRS was an effective signature to predict prognosis and was associated with BC clinical characteristics.
Figure 4 The relationship between HFRS score and clinicopathological features in BC patients. (A, D, G) The boxplots shows the comparison of HFRS risk score of BC subgroups stratified with different clinicopathological features (age>65 or<=65; positive or negative lymph node status; III/IV or I/II tumor stage) in TCGA cohort. (B, C) Comparison of the overall survival of patients with high- and low- HFRS risk score in age>65 (B) and age<=65 (C) subgroups. (E, F) KM curves to show the different overall survival of patients with high or low HFRS risk scores in lymph node positive (E) and lymph node negative (F) subgroups. (H, I) Comparison the survival of patients with high or low HFRS risk score in tumor stage III/IV (H) and tumor stage I/II (I) subgroups. (J–M) The boxplots shows the significant difference HFRS risk score levels of patients with different clinicopathological features(age, tumor stage, lymph node status and primary tumor grade) in METABRIC cohort.
3.5 Analysis of ferroptosis and hypoxia status
To explore whether HFRS could assess the ferroptosis status of BC patients, we first compared the expression of ferroptosis suppressors and drivers in high and low-HFRS groups. As shown in Figure 5A, in the TCGA cohort, except for STAT3 and HSPB1, the expression of the rest ferroptosis suppressors (ACSL3, ATF4, CA9, CD44, FTH1, GPX4, HELLS, HMOX1, HSF1, HSPA5, HSPB1, NQO1, SCD, SLC7A11) was significantly higher in the high-HFRS group (Figure 5A). To validate the result, we also analyzed the expression of the above genes in the METABRIC cohort, and the result was similar to those of the TCGA cohort (Figure 5B). In addition, we also compared the expression of ferroptosis drivers between the two groups in the TCGA cohort and METABRIC cohort. The results showed that in the low-HFRS group, the expression of more than half of the ferroptosis drivers (ALOX12 ANO6, ATF3, ATG5, ATG7, EGFR, CHAC1, EGLN2, ELAVL1, IREB2, KEAP1, NCOA4, and VDAC2 in TCGA cohort; ALOX12, ALOX15, ANO6, ATM, BAP1, DPP4, EGLN2, ELAVL1, IDH1, IREB2, KEAP1, NCOA4, SAT1, and VDAC1 in the METABRIC cohort) were significantly higher than in the high-HFRS group (Supplementary Figures 1A, B). These results suggested that ferroptosis might be induced in patients of the low-HFRS group.
Figure 5 The different expression of ferroptosis and hypoxia regulations between high and low HFRS risk groups. (A, B) The boxplots shows the difference in ferroptosis suppressors mRNA expression between the high- and low-HFRS groups of the TCGA (A) and METABRIC (B) cohorts. (C) The heat map shows the association between Reactome hypoxia-related pathways and HFRS score in TCGA BC samples. (D) The boxplot shows the significant difference in the enrichment score of HALLMARK_CANCER_HYPOXIA between high- and low-HFRS risk groups. (*p< 0.05; **p< 0.01; ***p< 0.001; ****p< 0.0001, n.s., not significant).
We further explored the hypoxia status of BC patients in the TCGA cohort by using GSEA analysis to estimate the enrichment scores of hypoxia and hypoxia-induced factors (HIFs) signal-related gene sets from REACTOME website. As shown in Figure 5C, these gene sets were enriched in the high-HFRS group, indicating that hypoxia status may be induced in the high-HFRS group. Furthermore, the box gram shows that the HALLMARK_CANCER_HYPOXIA enrichment score was significantly higher in high-HFRS group than low-HFRS group (p< 0.05) (Figure 5D). These results also implied that BC patients with high HFRS exhibit high hypoxia status.
3.6 Analysis of tumor immune cell infiltration
To investigate whether HFRS was associated with tumor immune microenvironment, the GSVA was used to analyze the enrichment of KEGG and GO pathways in high-HFRS and low-HFRS groups. KEGG analysis revealed that the following pathways were significantly activated in the high-HFRS group: fatty acid metabolism, HEDGEHOG signaling pathway, MAPK signaling pathway, TP53 signaling pathway, and others. The T-cell receptor signaling pathway, antigen processing, and presentation were significantly enriched in the low-HFRS group (Figure 6A). In addition, gene ontology (GO) biopathway analysis revealed that hypoxia and ferroptosis genes in low-HFRS group were significantly enriched in the immune-related functional sets such as natural killer cell-mediated immunity, T cell activation involved in immune response, and immune response regulating signaling pathway (Figure 6B). Thus, besides reflecting hypoxia and ferroptosis status, the HFRS might also be related to the tumor immune microenvironment.
Figure 6 The relationship between tumor immune cell infiltration and HFRS in BC patients. (A, B) KEGG analyses and GO analyses for hypoxia and ferroptosis-related genes of the high- and low-HFRS groups. (C) Comparison of the immune cells infiltration between the high- and low-HFRS groups of the TCGA cohort by the CIBERSORT algorithm. (D) The violin plots show significant difference in stromal, immune ESTIMATE scores and tumor purity of between high- and low-HFRS risk groups in TCGA cohort. (*p< 0.05; **p< 0.01; ****p< 0.0001, n.s., not significant).
We further investigated the difference in tumor immune cell infiltration between low- and high-HFRS patients. In CIBERSORT analysis, the fraction of B_cells_naive, Monocytes, NK_cells_activated, T_cells_CD4_memory_resting, T_cells_CD8, and T_cells_gamma_delta was significantly higher in the low-HFRS group, while Macrophages_M0/M2 and Tregs were significantly lower in the TCGA BC cohort (Figure 6C). The different infiltration fraction of 28 immune cells in the high- and low-HFRS groups was compared by ssGSEA as well (Supplementary Figure 2). Results of ssGSEA were generally consistent with those of CIBERSORT analysis. For instance, Activated_B_cell, Natural_Killer_cell, Activated_CD8_T_cells, and Monocyte were significantly highly infiltrated in samples of the low-HFRS group, while the Regulatory_T_cell was highly infiltrated in high-HFRS BC samples. Moreover, high infiltration of myeliod-derived suppressor cell and Immature_dendritic_cell and a lower Mast_cell were detected in the high-HFRS group by ssGSEA. These data revealed that the HFRS score was associated with immune cell infiltration in breast cancer.
In addition, the results of the ESTIMATE analysis showed that the immune score and tumor purity of the low-HFRS group were significantly higher than in the high-HFRS group, while the ESTIMATE score was lower than in the high-HFRS group (Figure 6D). Moreover, the IPS scores of the low-HFRS group were significantly higher than in the high-HFRS group in all the four subgroups (ips_ctla4_neg_pd1_neg, ips_ctla4_neg_pd1_pos, ips_ctla4_pos_pd1_neg, ips_ctla4_pos_pd1_pos) (Supplementary Figure 3). Thus, these results revealed that high HFRS might be associated with reduced anti-tumor immunity and decreased tumor purity.
3.7 Analysis of pathway and process enrichment
We identified a total of 272 DEGs between the high-HFRS group and the low-HFRS group with the criteria | log2(fold change) | > 1 and p< 0.05 (Supplementary Figure 4A). Then, DEGs were subjected to GSEA analysis based on REACTOME gene sets. The results showed that the DEGs were significantly enriched in the following terms: cell cycle and cellular response to hypoxia (Supplementary Figures 4B, C). Moreover, several tumor-related pathways and metabolism-related gene sets were enriched in the high-HFRS subgroup, such as CHOLESTEROL_HOMEOSTASIS, MTORC1_SIGNALING, TP53-PTEN related gene sets, MYC targets V1, E2F_TARGETS, PIK3-AKT-MTOR signaling (Supplementary Figures 4D, E–H). These results suggested that HFRS might be related to multiple tumor biology processes via communicating with cell cycle regulation, hypoxia microenvironment, energy metabolism, and oncogenic signal pathways, and may provide a new perspective and help us to find the potential therapeutic targets from cancer-related pathways.
3.8 Analysis of the gene mutation profile of different HFRS groups
To investigate the difference of gene mutation between the high- and low-HFRS groups, we analyzed the simple nucleoside variation profile of two groups in the TCGA cohort. As shown in Figures 7A, B, in the low-HFRS group, the top five genes with mutation frequency were PIK3CA (39%), TP53 (26%), CDH1 (18%), TTN (13%), and GATA3 (12%), while those in the high-HFRS group were TP53 (47%), PIK3CA (27%), TTN (22%), GATA3 (12%), and MUC16 (11%). TP53 is one of the most important tumor suppressor genes, whose mutation could lead to tumor occurrence and progression, and might be associated with suppressed ferroptosis and anti-tumor immunity (35, 36). TP53 mutation indicated worse prognosis in breast cancer (37). In this study, we found that the mutation frequency of TP53 in high HFRS was higher than in low HFRS (47% vs. 26%), which also suggested that high HFRS might be associated with TP53 mutation-induced prognostic risk.
Figure 7 Mutation landscape and prognostic factors of BC patients. (A, B) Oncoplots show the mutated genes in the high-HFRS (A) and low-HFRS (B) groups of the TCGA cohort. (C) Forrest plots of univariate and multivariate analyses show the independent prognostic predictors in the TCGA and METABRIC cohorts.
3.9 HFRS is an independent prognostic factor for BC
To identify the clinical factors to predict the prognosis in BC, we used univariate and multivariate Cox regression analysis to estimate the hazard ratio with HFRS score and other clinicopathological features in both cohorts. The results of univariate Cox regression analysis showed that HFRS was a strong risk factor for OS in BC patients (in the TCGA cohort, HR: 2.941, 95% confidence interval (CI): 2.214–3.906, p< 0.001; in the METABRIC cohort, HR: 1.429, 95% CI:1.263–1.617, p< 0.001; Figures 7C, E). The multivariate regression analysis showed that HFRS was an independent prognostic factor for BC patients (in TCGA, HR: 5.512, 95% CI: 2.756–8.932, p< 0.001; In METABRIC, HR: 1.251, 95% CI:1.101–1.420, p< 0.001; Figures 7D, F). Then, the survival analysis of DFS and RFS of BC patients showed that patients with higher HFRS score have significantly shorter DFS and RFS (Supplementary Figure 5). The above results indicated that HFRS was an independent prognostic factor for BC patients.
3.10 Construction and validation of nomogram base on HFRS
We developed a nomogram based on HFRS and other independent prognostic factors (TNM stage) in the TCGA cohort to affiliate the application of HFRS in clinical practice (Figure 8A), which was validated in the METABRIC cohort. Calibration curves show that the predicted rates were highly concordant with the actual rates for 1-, 3-, and 5-year survival in the TCGA cohort (Figures 8B–D), and 1-, 3-, 5-, 8- and 10-year survival in the METABRIC cohort (Figures 8E–I). Moreover, ROC curves show that the prognostic predictive ability of the nomogram model in BC patients was better than other factors (including HFRS score, age, and TNM stage). The AUCs of 1, 3, and 5 years reached 0.81, 0.81, and 0.79 in TCGA cohort (Figures 8J–L) and 0.74, 0.70, and 0.67 (1, 3, and 5 years) in the METABRIC cohort (Figures 8M–O). These results indicated that the nomogram, based on HFRS score and TNM stage, has a strong and stable ability to predict the OS of BC patients.
Figure 8 Construction and validation of nomogram. (A) The nomogram based on HFRS risk score, age and tumor stage for predicting overall survival of BC patients. (B–D) Calibration plots of the nomogram for predicting the probability of OS at 1, 3, and 5 years in the TCGA cohort. (E-I) The calibration plots of the nomogram for predicting 1-, 3-, 5-, 8-, and 10-year OS in the METABRIC cohort. (J–L) ROC curves of nomogram, risk score, tumor stage, age for predicting 1-, 3-, 5-year OS in the TCGA cohort. (M–O) ROC curves of nomogram, risk score, tumor stage, age for predicting 1-, 3-, 5-year OS in the METABRIC cohort.
3.11 Protein expression of 10 HFRS genes in normal breast tissues and BC tissues
To explore the protein expression of HFRS genes in BC tumor tissues and normal breast tissues, we first collected the immunohistochemical staining images of several HFRS genes from HPA; the protein expression of CCT6A, CISD1, P4HA2, PGK1, TPD52 were higher in tumor tissues than in normal tissues, while KRT14, TF, STC2 were more highly expressed in normal breast tissues (Figure 9A). Then, we found that the protein expression of SLC16A2 was higher in BC tumor tissues, whereas that of BTG1 was higher in adjacent normal tissues by immunohistochemical staining assay (Figure 9B, Supplementary Figure 6).
Figure 9 The protein expression of HFRS genes in normal and tumor samples. (A) Representative immunohistochemistry images of CCT6A, CISD, KRT14, P4HA2, PGK1, STC2, TF, and TPD52 expression between BC tissues and adjacent normal tissues. (B) IHC staining shows the protein expression of SLC16A2 and BTG1 in BC tissues and adjacent normal tissues. Scale bar: 100μm.
4 Discussion
Breast cancer is a highly heterogeneous solid tumor, calling for individualized treatment for BC patients to improve prognosis. Though current surgery, endocrine, chemotherapy, and target therapy could improve the overall survival of BC patients, a large proportion of patients recur or progress, which leads to poor prognosis. Thus, investigating the differentially expressed genes and their roles in tumor malignant biological processes might help to analyze clinicopathological features of individual patients and offer precise therapeutic regimens and estimate outcomes for BC patients.
Hypoxia is one of the most impactful hallmarks of solid tumors that could influence tumor progression such as promoting tumor cell proliferation, invasion, and regulating cell cycle, energy metabolism, and immune escape (38–40). Recently, many studies have demonstrated that hypoxia status was considered as an important characteristic of the tumor microenvironment that has a close relationship with immune therapy sensitivity (41–43). Regarding breast cancer, hypoxia could induce cell growth by activating the glycogen metabolic program, improve migration, angiogenesis, and regulate apoptosis (44, 45). Moreover, hypoxia suppresses immune effector gene expression in immune cells, leading to immune effector cell dysfunction and resistance to anti-PD-1 therapy in triple-negative breast cancer (46).
Ferroptosis, a newly identified programmed cell death, has been found to have a relationship with tumor occurrence and development (47). Repression of ferroptosis could promote tumor progression, using anti-tumor drugs such as PI3K-AKT-mTOR pathway inhibitors (GDC-0941, MK-2206) which could promote sensitivity to ferroptosis in breast cancer cells (48). Some studies revealed that metformin could induce ferroptosis in breast cancer, which suggested that the patients simultaneously suffering from BC and type 2 diabetes could prolong their survival by metformin; however, further clinical trials were needed to provide more convincing evidence (49–51).
Despite the progress that has been made, the identification of effective prognostic biomarkers and the development of drugs targeting hypoxia and ferroptosis remain scarce. Antitumor therapeutic strategies directly targeting the hypoxic microenvironment are mostly focused on developing nanoparticles, and drugs targeting ferroptosis were far from clinical application (52). Reassuringly, recent studies have attracted our attention. It has been found that hypoxia can inhibit ferroptosis in hepatocellular carcinoma and breast cancer (53). Alternatively, ferroptosis could enhance the radiosensitivity of hypoxic tumor cells by amplifying oxidative stress or inhibiting antioxidant regulation (54). The nanoplatform-based tumor reoxygenation, which could generate the active superoxide radical (), together with H2O2-participated iron-involved Fenton reactions of ferroptosis, plays a synergistic role in overcoming hypoxia-induced chemotherapy resistance of osteosarcoma in vivo (55). These inspiring findings suggest that simultaneously suppressing hypoxia and inducing ferroptosis sensitivity of tumors may produce a potently synergistic antitumor effect.
Herein, we hypothesize that we could identify hypoxia and ferroptosis-related genes that associated with breast cancer prognosis, clinical characteristics, and immune microenvironment, to provide potential biomarkers and molecular targets for anti-tumor drug development and construct efficient gene signatures that could predict prognosis and simultaneously reflect tumor microenvironment characteristics as well as the ferroptosis status of breast cancer.
In this study, we included BC patients from the TCGA dataset and screened prognostic genes related to hypoxia and ferroptosis from known gene sets (MsigDB and ferroDB). The Lasso–Cox method was used to construct a predictive model (HFPS) based on these genes and divided BC patients into high- and low-risk groups. Notably, the high-risk group showed significantly worse overall survival, higher TNM stage, higher rate of lymph node invasion, and a lower rate of ER positive than the low risk group which could help us to implement individual strategies in clinical practice. The results were validated in the METABRIC BC cohort. Moreover, to predict the prognosis of patients, we constructed a nomogram model based on HFRS and prognostic clinical factors including age and TNM stage to predict 1-, 3-, 5-year overall survival of BC patients. The high predictive ability was validated by calibration curves, and ROC curves in both TCGA and METABRIC cohorts, which could help to make individual clinical decisions for patients.
Then we identified that there were more suppressive ferroptosis and higher hypoxia status in high-HFRS patients than in low-HFRS patients. In GSEA analysis, the key gene sets that assess tumor hypoxia status such as WINTER_HYPOXIA_METAGENE, WINTER_HYPOXIA_UP, HARRIS_HYPOXIA, and REACTOME_CELLULAR_RESPONSE_TO_HYPOXIA were highly enriched in patients with high HFRS. These results implied that BC patients with high HFRS were more likely to form hypoxia microenvironment. The mRNA expression of ferroptosis drivers were reduced in high HFRS group, while ferroptosis suppressors were highly expressed. Some evidence demonstrated that targeting ferroptosis-related genes and promoting ferroptosis sensitivity of tumor cells were promising approaches for reducing cancer progression (56). Thus, patients with a higher HFRS which also present a repressed ferroptosis status might profit from therapeutic drugs that could induce ferroptosis such as Lapatinib and Cisplatin (57–59).
By analyzing the immune microenvironment phenotype of the high- and low-HFRS groups, we found that the immune cell infiltration of B_cell_naive, T_cell_CD4, T_cell_CD8, and Tregs, macrohpage_M0, M2, monocytes, NK_cell_activated, and NK_cell_resting were significantly different between the high-and low-HFRS group by CIBERSORT analysis. In further ssGSEA analysis, we confirmed that HFRS score was negatively correlated with activated CD8 T cell and activated B cell, while it was positively correlated with MDSCs and Tregs. Studies have reported that repressed CD4, CD8 and T cell and activated B cells infiltration indicated an immune suppression microenvironment in cancer. MDSCs mediated immune suppression via expansion and differentiation of Tregs and limiting NKs, DCs, and the polarization of macrophages to M2-phenotype, and were associated with clinical outcome of BC (60). Tregs has been wildly reported to promote cancer immune escape and contribute to BC progression (61, 62). Studies have reported that tumor immune cell infiltration could be regulated by hypoxia status and related pathways. In breast cancer, hypoxia boosted CD8+ T cell infiltration in tumor tissue and increased sensitivity to immune checkpoint blockade (63). The CD4+ T cell has been wildly demonstrated to possess cytotoxic programs and can directly kill cancer cells (64). Hypoxia or HIF-1α signal pathway could influence CD4+ T cell function, metabolism, differentiation, and infiltration to enhance immunosuppression in tumors (65, 66). Suthen et al. reported that Tregs and immunosuppressive myeloid subsets were found to be significantly enriched in the hypoxia tumor tissue regions (67). Furthermore, patients in high the HFRS group showed significantly higher Stromal score and Immune score, while they have lower tumor purity in ESTIMATE analysis, which indicated that HFRS was associated with immune microenvironment of BC. In conclusion, high HFRS was positively correlated with immunosuppression in BC cancers.
Some of the HFRS genes have been demonstrated to serve as tumor oncogenes or tumor suppressors and might be correlated with tumor biological behavior and prognosis in diverse kinds of cancers. BTG1 has been reported as a tumor suppressor inhibiting tumor proliferation and migration and increasing anti-tumor therapy sensitivity in some kinds of tumors (68–70), including breast cancer (71–74). However, some studies reported that its overexpression promoted tumor malignancy in colorectal cancer (75). Researches demonstrated that STC2 could impair breast cancer cell growth, migration, and cell viability, which was consistent with our results (76, 77). In colorectal cancer, the upregulated STC2 was associated with a poorer prognosis (78). Lin et al. also reported that STC2 promoted pancreatic cancer migration, invasion, and EMT (79). A pan-cancer research found that STC2 was closely related to tumor immune microenvironment including immune cell infiltration, ICGs, MMRs, TMB, and MSI (80). One bioinformatic analysis reported that overexpression of CCT6A in tumor tissue was associated with poor breast cancer prognosis (81). Jie Jiang et al. reported that upregulated CCT6A in Ewing sarcoma was correlated with a worse prognosis (82). A similar result was observed in hepatocellular carcinoma (83). Studies of Bilandzic et al. implicated the basal epithelial marker KRT14 as an absolute determinant for ovarian cancer cells’ spheroid integrity, mesothelial attachment, invasive potential, and chemotherapy resistance, which could provide some in vitro evidences to explain the role KRT14 plays in cancer (84, 85). Thus, exploring the detail function of KRT14 in breast cancer is required in further studies. Recently, P4HA2 has been demonstrated to play important roles in tumor, but its function in cancers might be different. For example, P4HA2 induced EMT and promote tumor growth, migration, and invasion in cervical cancer and glioma (86, 87), while in prostate and pancreatic cancer, it served as a tumor suppressor (88, 89). Consistent with our findings, studies demonstrated that high P4HA2 expression was associated with poor survival in breast cancer (90, 91). PGK1 is a glycolytic enzyme that catalyzes the conversion of 1,3-diphosphoglycerate to 3-phosphoglycerate and participates in tumor angiogenesis (92, 93). Many studies reported PGK1 as a prognostic gene in cancers, and it has been demonstrated to promote EMT and the progression of breast cancer (94–99). SLC16A2 is a member of SLC16 gene family, that encodes monocarboxylate transporters, but its function in cancer has not been identified yet, which required further investigation (100). TF is also known as Transferrin which is essential for ferric iron transporting into cells and could influence iron metabolism in human, and might be involved in ferroptosis regulation in tumor cells indirectly (101). In addition, the knockdown of transferrin leads to decreased lapatinib-related BC cell death, but further in vivo experiments were absent (16). TPD52 is an oncogene and closely associated with prostate, breast cancer, and other cancers (102–104), which was consistent with our results. CISD1 reduces ferroptosis via iron-sulfur cluster biogenesis and was identified as prognostic ferroptosis-related genes in bladder cancer, lung cancer, and hepatocellular carcinoma (53, 105, 106). Although these researches provided some evidence to demonstrate its relation with tumor disease, its function in breast cancer has not been investigated yet. Thus, our study identified HFRS genes that might provide potential targets for the development of clinical therapeutic regimens.
Previous studies mainly focus on ferroptosis or hypoxia-related genes to develop prognostic models, but there is no study to consider the cell death regulation and microenvironment heterogeneity of breast cancer together. Our study firstly explored the effect of combining hypoxia and ferroptosis on breast cancer prognosis by constructing a novel predictive signature (HFRS) with hypoxia and ferroptosis-related genes. Additionally, HFRS could distinguish ferroptosis, hypoxia status, immune cell infiltration, and clinical characteristics of BC patients, which might help to make individual therapeutic strategies. Meanwhile, to improve the sensitivity and specificity of HFRS, we established a nomogram based on HFRS and clinical prognostic factors, which could also facilitate the clinical application of HFRS.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Ethics statement
The study was reviewed and approved by Ethics committees of Wuhan Fourth Hospital. The patients/participants provided their written informed consent to participate in this study.
Author contributions
JXZ and FL designed and constructed this study. XS, HY, and BW performed the data analysis, figure plotting, and writing. JJZ did the immunohistochemistry staining. FL, JW, and JS were responsible for the data acquisition and critical reading of the manuscript. All authors contributed to the article and approved the submitted version.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2023.1077342/full#supplementary-material
References
1. DeSantis CE, Ma J, Gaudet MM, Newman LA, Miller KD, Goding Sauer A, et al. Breast cancer statistics, 2019. CA Cancer J Clin (2019) 69(6):438–51. doi: 10.3322/caac.21583
2. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin (2021) 71(3):209–49. doi: 10.3322/caac.21660
3. Goldhirsch A, Winer EP, Coates AS, Gelber RD, Piccart-Gebhart M, Thurlimann B, et al. Personalizing the treatment of women with early breast cancer: highlights of the St gallen international expert consensus on the primary therapy of early breast cancer 2013. Ann Oncol (2013) 24(9):2206–23. doi: 10.1093/annonc/mdt303
4. Amin MB, Greene FL, Edge SB, Compton CC, Gershenwald JE, Brookland RK, et al. The eighth edition AJCC cancer staging manual: Continuing to build a bridge from a population-based to a more "personalized" approach to cancer staging. CA Cancer J Clin (2017) 67(2):93–9. doi: 10.3322/caac.21388
5. Hashmi AA, Hashmi KA, Irfan M, Khan SM, Edhi MM, Ali JP, et al. Ki67 index in intrinsic breast cancer subtypes and its association with prognostic parameters. BMC Res Notes. (2019) 12(1):605. doi: 10.1186/s13104-019-4653-x
6. Leong AS, Raymond WA. Prognostic parameters in breast cancer. Pathology (1989) 21(3):169–75. doi: 10.3109/00313028909061053
7. Rakha EA, Reis-Filho JS, Baehner F, Dabbs DJ, Decker T, Eusebi V, et al. Breast cancer prognostic classification in the molecular era: the role of histological grade. Breast Cancer Res (2010) 12(4):207. doi: 10.1186/bcr2607
8. Harris AL. Hypoxia–a key regulatory factor in tumour growth. Nat Rev Cancer. (2002) 2(1):38–47. doi: 10.1038/nrc704
9. Vaupel P, Kelleher DK, Thews O. Modulation of tumor oxygenation. Int J Radiat Oncol Biol Phys (1998) 42(4):843–8. doi: 10.1016/S0360-3016(98)00324-1
10. Gilkes DM, Semenza GL, Wirtz D. Hypoxia and the extracellular matrix: drivers of tumour metastasis. Nat Rev Cancer. (2014) 14(6):430–9. doi: 10.1038/nrc3726
11. Barsoum IB, Smallwood CA, Siemens DR, Graham CH. A mechanism of hypoxia-mediated escape from adaptive immunity in cancer cells. Cancer Res (2014) 74(3):665–74. doi: 10.1158/0008-5472.CAN-13-0992
12. Riera-Domingo C, Audige A, Granja S, Cheng WC, Ho PC, Baltazar F, et al. Immunity, hypoxia, and metabolism-the menage a trois of cancer: Implications for immunotherapy. Physiol Rev (2020) 100(1):1–102. doi: 10.1152/physrev.00018.2019
13. Lu B, Chen XB, Ying MD, He QJ, Cao J, Yang B. The role of ferroptosis in cancer development and treatment response. Front Pharmacol (2017) 8:992. doi: 10.3389/fphar.2017.00992
14. Mou Y, Wang J, Wu J, He D, Zhang C, Duan C, et al. Ferroptosis, a new form of cell death: opportunities and challenges in cancer. J Hematol Oncol (2019) 12(1):34. doi: 10.1186/s13045-019-0720-y
15. Li Z, Chen L, Chen C, Zhou Y, Hu D, Yang J, et al. Targeting ferroptosis in breast cancer. biomark Res (2020) 8(1):58. doi: 10.1186/s40364-020-00230-3
16. Ma S, Henson ES, Chen Y, Gibson SB. Ferroptosis is induced following siramesine and lapatinib treatment of breast cancer cells. Cell Death Dis (2016) 7:e2307. doi: 10.1038/cddis.2016.208
17. Song X, Wang X, Liu Z, Yu Z. Role of GPX4-mediated ferroptosis in the sensitivity of triple negative breast cancer cells to gefitinib. Front Oncol (2020) 10:597434. doi: 10.3389/fonc.2020.597434
18. Fuhrmann DC, Mondorf A, Beifuss J, Jung M, Brune B. Hypoxia inhibits ferritinophagy, increases mitochondrial ferritin, and protects from ferroptosis. Redox Biol (2020) 36:101670. doi: 10.1016/j.redox.2020.101670
19. Wang Y, Zhang L, Zhou X. Activation of Nrf2 signaling protects hypoxia-induced HTR-8/SVneo cells against ferroptosis. J Obstet Gynaecol Res (2021) 47(11):3797–806. doi: 10.1111/jog.15009
20. Fan Z, Yang G, Zhang W, Liu Q, Liu G, Liu P, et al. Hypoxia blocks ferroptosis of hepatocellular carcinoma via suppression of METTL14 triggered YTHDF2-dependent silencing of SLC7A11. J Cell Mol Med (2021) 25(21):10197–212. doi: 10.1111/jcmm.16957
21. Lin Z, Song J, Gao Y, Huang S, Dou R, Zhong P, et al. Hypoxia-induced HIF-1alpha/lncRNA-PMAN inhibits ferroptosis by promoting the cytoplasmic translocation of ELAVL1 in peritoneal dissemination from gastric cancer. Redox Biol (2022) 52:102312. doi: 10.1016/j.redox.2022.102312
22. Emens LA. Breast cancer immunotherapy: Facts and hopes. Clin Cancer Res (2018) 24(3):511–20. doi: 10.1158/1078-0432.CCR-16-3001
23. Venetis K, Invernizzi M, Sajjadi E, Curigliano G, Fusco N. Cellular immunotherapy in breast cancer: The quest for consistent biomarkers. Cancer Treat Rev (2020) 90:102089. doi: 10.1016/j.ctrv.2020.102089
24. Ding XC, Wang LL, Zhang XD, Xu JL, Li PF, Liang H, et al. The relationship between expression of PD-L1 and HIF-1alpha in glioma cells under hypoxia. J Hematol Oncol (2021) 14(1):92. doi: 10.1186/s13045-021-01102-5
25. Kumar V, Gabrilovich DI. Hypoxia-inducible factors in regulation of immune responses in tumour microenvironment. Immunology (2014) 143(4):512–9. doi: 10.1111/imm.12380
26. Liu G, Bi Y, Shen B, Yang H, Zhang Y, Wang X, et al. SIRT1 limits the function and fate of myeloid-derived suppressor cells in tumors by orchestrating HIF-1alpha-dependent glycolysis. Cancer Res (2014) 74(3):727–37. doi: 10.1158/0008-5472.CAN-13-2584
27. Zhang H, Lu H, Xiang L, Bullen JW, Zhang C, Samanta D, et al. HIF-1 regulates CD47 expression in breast cancer cells to promote evasion of phagocytosis and maintenance of cancer stem cells. Proc Natl Acad Sci U S A. (2015) 112(45):E6215–23. doi: 10.1073/pnas.1520032112
28. Lukashev D, Klebanov B, Kojima H, Grinberg A, Ohta A, Berenfeld L, et al. Cutting edge: hypoxia-inducible factor 1alpha and its activation-inducible short isoform I.1 negatively regulate functions of CD4+ and CD8+ T lymphocytes. J Immunol (2006) 177(8):4962–5. doi: 10.4049/jimmunol.177.8.4962
29. Thiel M, Caldwell CC, Kreth S, Kuboki S, Chen P, Smith P, et al. Targeted deletion of HIF-1alpha gene in T cells prevents their inhibition in hypoxic inflamed tissues and improves septic mice survival. PloS One (2007) 2(9):e853. doi: 10.1371/journal.pone.0000853
30. Friedmann Angeli JP, Krysko DV, Conrad M. Ferroptosis at the crossroads of cancer-acquired drug resistance and immune evasion. Nat Rev Cancer. (2019) 19(7):405–14. doi: 10.1038/s41568-019-0149-1
31. Jiang Z, Lim SO, Yan M, Hsu JL, Yao J, Wei Y, et al. TYRO3 induces anti-PD-1/PD-L1 therapy resistance by limiting innate immunity and tumoral ferroptosis. J Clin Invest (2021) 131(8):e139434. doi: 10.1172/JCI139434
32. Li H, Yang P, Wang J, Zhang J, Ma Q, Jiang Y, et al. HLF regulates ferroptosis, development and chemoresistance of triple-negative breast cancer by activating tumor cell-macrophage crosstalk. J Hematol Oncol (2022) 15(1):2. doi: 10.1186/s13045-021-01223-x
33. Wang W, Green M, Choi JE, Gijon M, Kennedy PD, Johnson JK, et al. CD8(+) T cells regulate tumour ferroptosis during cancer immunotherapy. Nature (2019) 569(7755):270–4. doi: 10.1038/s41586-019-1170-y
34. Ilic IR, Stojanovic NM, Radulovic NS, Zivkovic VV, Randjelovic PJ, Petrovic AS, et al. The quantitative ER immunohistochemical analysis in breast cancer: Detecting the 3 + 0, 4 + 0, and 5 + 0 allred score cases. Medicina (Kaunas) (2019) 55(8):461. doi: 10.3390/medicina55080461
35. Mogi A, Kuwano H. TP53 mutations in nonsmall cell lung cancer. J BioMed Biotechnol (2011) 2011:583929. doi: 10.1155/2011/583929
36. Silwal-Pandit L, Langerod A, Borresen-Dale AL. TP53 mutations in breast and ovarian cancer. Cold Spring Harb Perspect Med (2017) 7(1):a026252. doi: 10.1101/cshperspect.a026252
37. Kaur RP, Vasudeva K, Kumar R, Munshi A. Role of p53 gene in breast cancer: Focus on mutation spectrum and therapeutic strategies. Curr Pharm Des (2018) 24(30):3566–75. doi: 10.2174/1381612824666180926095709
38. Bao MH, Wong CC. Hypoxia, metabolic reprogramming, and drug resistance in liver cancer. Cells (2021) 10(7):1715. doi: 10.3390/cells10071715
39. Infantino V, Santarsiero A, Convertini P, Todisco S, Iacobazzi V. Cancer cell metabolism in hypoxia: Role of HIF-1 as key regulator and therapeutic target. Int J Mol Sci (2021) 22(11):5703. doi: 10.3390/ijms22115703
40. Luo W, Wang Y. Hypoxia mediates tumor malignancy and therapy resistance. Adv Exp Med Biol (2019) 1136:1–18. doi: 10.1007/978-3-030-12734-3_1
41. Jiang M, Qin B, Luo L, Li X, Shi Y, Zhang J, et al. A clinically acceptable strategy for sensitizing anti-PD-1 treatment by hypoxia relief. J Control Release. (2021) 335:408–19. doi: 10.1016/j.jconrel.2021.06.001
42. Reeves KM, Song PN, Angermeier A, Della Manna D, Li Y, Wang J, et al. (18)F-FMISO PET imaging identifies hypoxia and immunosuppressive tumor microenvironments and guides targeted evofosfamide therapy in tumors refractory to PD-1 and CTLA-4 inhibition. Clin Cancer Res (2022) 28(2):327–37. doi: 10.1158/1078-0432.CCR-21-2394
43. Tang K, Zhu L, Chen J, Wang D, Zeng L, Chen C, et al. Hypoxia promotes breast cancer cell growth by activating a glycogen metabolic program. Cancer Res (2021) 81(19):4949–63. doi: 10.1158/0008-5472.CAN-21-0753
44. Rana NK, Singh P, Koch B. CoCl2 simulated hypoxia induce cell proliferation and alter the expression pattern of hypoxia associated genes involved in angiogenesis and apoptosis. Biol Res (2019) 52(1):12. doi: 10.1186/s40659-019-0221-z
45. Wang R, Godet I, Yang Y, Salman S, Lu H, Lyu Y, et al. Hypoxia-inducible factor-dependent ADAM12 expression mediates breast cancer invasion and metastasis. Proc Natl Acad Sci U.S.A. (2021) 118(19):e2020490118. doi: 10.1073/pnas.2020490118
46. Ma S, Zhao Y, Lee WC, Ong LT, Lee PL, Jiang Z, et al. Hypoxia induces HIF1alpha-dependent epigenetic vulnerability in triple negative breast cancer to confer immune effector dysfunction and resistance to anti-PD-1 immunotherapy. Nat Commun (2022) 13(1):4118. doi: 10.1038/s41467-022-31764-9
47. Lei G, Zhuang L, Gan B. Targeting ferroptosis as a vulnerability in cancer. Nat Rev Cancer. (2022) 22(7):381–96. doi: 10.1038/s41568-022-00459-0
48. Yi J, Zhu J, Wu J, Thompson CB, Jiang X. Oncogenic activation of PI3K-AKT-mTOR signaling suppresses ferroptosis via SREBP-mediated lipogenesis. Proc Natl Acad Sci U S A. (2020) 117(49):31189–97. doi: 10.1073/pnas.2017152117
49. Yang J, Zhou Y, Xie S, Wang J, Li Z, Chen L, et al. Metformin induces ferroptosis by inhibiting UFMylation of SLC7A11 in breast cancer. J Exp Clin Cancer Res (2021) 40(1):206. doi: 10.1186/s13046-021-02012-7
50. Hou Y, Cai S, Yu S, Lin H. Metformin induces ferroptosis by targeting miR-324-3p/GPX4 axis in breast cancer. Acta Biochim Biophys Sin (Shanghai). (2021) 53(3):333–41. doi: 10.1093/abbs/gmaa180
51. Dias Lopes NM, Marinello PC, Sanches LJ, da Silva Brito WA, Lovo-Martins MI, Pinge-Filho P, et al. Patterns of cell death induced by metformin in human MCF-7 breast cancer cells. Pathol Res Pract (2020) 216(11):153199. doi: 10.1016/j.prp.2020.153199
52. Wang H, Li J, Wang Y, Gong X, Xu X, Wang J, et al. Nanoparticles-mediated reoxygenation strategy relieves tumor hypoxia for enhanced cancer therapy. J Control Release. (2020) 319:25–45. doi: 10.1016/j.jconrel.2019.12.028
53. Chen X, Kang R, Kroemer G, Tang D. Broadening horizons: the role of ferroptosis in cancer. Nat Rev Clin Oncol. (2021) 18(5):280–96. doi: 10.1038/s41571-020-00462-0
54. Su J, Zhao Q, Zheng Z, Wang H, Bian C, Meng L, et al. Prospective application of ferroptosis in hypoxic cells for tumor radiotherapy. Antioxidants (Basel) (2022) 11(5):921. doi: 10.3390/antiox11050921
55. Fu J, Li T, Yang Y, Jiang L, Wang W, Fu L, et al. Activatable nanomedicine for overcoming hypoxia-induced resistance to chemotherapy and inhibiting tumor growth by inducing collaborative apoptosis and ferroptosis in solid tumors. Biomaterials (2021) 268:120537. doi: 10.1016/j.biomaterials.2020.120537
56. Xu G, Wang H, Li X, Huang R, Luo L. Recent progress on targeting ferroptosis for cancer therapy. Biochem Pharmacol (2021) 190:114584. doi: 10.1016/j.bcp.2021.114584
57. Su Y, Zhao B, Zhou L, Zhang Z, Shen Y, Lv H, et al. Ferroptosis, a novel pharmacological mechanism of anti-cancer drugs. Cancer Lett (2020) 483:127–36. doi: 10.1016/j.canlet.2020.02.015
58. Guo J, Xu B, Han Q, Zhou H, Xia Y, Gong C, et al. Ferroptosis: A novel anti-tumor action for cisplatin. Cancer Res Treat (2018) 50(2):445–60. doi: 10.4143/crt.2016.572
59. Ni J, Chen K, Zhang J, Zhang X. Inhibition of GPX4 or mTOR overcomes resistance to lapatinib via promoting ferroptosis in NSCLC cells. Biochem Biophys Res Commun (2021) 567:154–60. doi: 10.1016/j.bbrc.2021.06.051
60. Weber R, Groth C, Lasser S, Arkhypov I, Petrova V, Altevogt P, et al. IL-6 as a major regulator of MDSC activity and possible target for cancer immunotherapy. Cell Immunol (2021) 359:104254. doi: 10.1016/j.cellimm.2020.104254
61. Hashemi V, Maleki LA, Esmaily M, Masjedi A, Ghalamfarsa G, Namdar A, et al. Regulatory T cells in breast cancer as a potent anti-cancer therapeutic target. Int Immunopharmacol. (2020) 78:106087. doi: 10.1016/j.intimp.2019.106087
62. Laine A, Labiad O, Hernandez-Vargas H, This S, Sanlaville A, Leon S, et al. Regulatory T cells promote cancer immune-escape through integrin alphavbeta8-mediated TGF-beta activation. Nat Commun (2021) 12(1):6228. doi: 10.1371/journal.pone.0000853
63. Gomes-Santos IL, Amoozgar Z, Kumar AS, Ho WW, Roh K, Talele NP, et al. Exercise training improves tumor control by increasing CD8(+) T-cell infiltration via CXCR3 signaling and sensitizes breast cancer to immune checkpoint blockade. Cancer Immunol Res (2021) 9(7):765–78. doi: 10.1158/2326-6066.CIR-20-0499
64. Oh DY, Fong L. Cytotoxic CD4(+) T cells in cancer: Expanding the immune effector toolbox. Immunity (2021) 54(12):2701–11. doi: 10.1016/j.immuni.2021.11.015
65. Westendorf AM, Skibbe K, Adamczyk A, Buer J, Geffers R, Hansen W, et al. Hypoxia enhances immunosuppression by inhibiting CD4+ effector T cell function and promoting treg activity. Cell Physiol Biochem (2017) 41(4):1271–84. doi: 10.1159/000464429
66. Yu Q, Dong L, Li Y, Liu G. SIRT1 and HIF1alpha signaling in metabolism and immune responses. Cancer Lett (2018) 418:20–6. doi: 10.1016/j.canlet.2017.12.035
67. Suthen S, Lim CJ, Nguyen PHD, Dutertre CA, Lai HLH, Wasser M, et al. Hypoxia-driven immunosuppression by treg and type-2 conventional dendritic cells in HCC. Hepatology (2022) 76(5):1329–44. doi: 10.1002/hep.32419
68. Li Y, Huo J, He J, Zhang Y, Ma X. BTG1 inhibits malignancy as a novel prognosis signature in endometrial carcinoma. Cancer Cell Int (2020) 20:490. doi: 10.1186/s12935-020-01591-3
69. Sun GG, Wang YD, Cheng YJ, Hu WN. BTG1 underexpression is an independent prognostic marker in esophageal squamous cell carcinoma. Tumour Biol (2014) 35(10):9707–16. doi: 10.1007/s13277-014-2245-x
70. Yuniati L, Scheijen B, van der Meer LT, van Leeuwen FN. Tumor suppressors BTG1 and BTG2: Beyond growth control. J Cell Physiol (2019) 234(5):5379–89. doi: 10.1002/jcp.27407
71. Li W, Zou ST, Zhu R, Wan JM, Xu Y, Wu HR. Bcell translocation 1 gene inhibits cellular metastasisassociated behavior in breast cancer. Mol Med Rep (2014) 9(6):2374–80. doi: 10.3892/mmr.2014.2118
72. Sheng SH, Zhao CM, Sun GG. BTG1 expression correlates with the pathogenesis and progression of breast carcinomas. Tumour Biol (2014) 35(4):3317–26. doi: 10.1007/s13277-013-1437-0
73. Vadgama JV, Scuric Z, Chakrabarti R, Marzo E, Shen D, Wu Y. Insulin-like growth factor I differentially regulates the expression of HIRF1/hCAF1 and BTG1 genes in human MCF-7 breast cancer cells. Int J Mol Med (2006) 18(1):129–39. doi: 10.3892/ijmm.18.1.129
74. Zhu R, Li W, Xu Y, Wan J, Zhang Z. Upregulation of BTG1 enhances the radiation sensitivity of human breast cancer in vitro and in vivo. Oncol Rep (2015) 34(6):3017–24. doi: 10.3892/or.2015.4311
75. Zhao S, Xue H, Hao CL, Jiang HM, Zheng HC. BTG1 overexpression might promote invasion and metastasis of colorectal cancer via decreasing adhesion and inducing epithelial-mesenchymal transition. Front Oncol (2020) 10:598192. doi: 10.3389/fonc.2020.598192
76. Jiang ST, Wang HQ, Yang TC, Wang DW, Yang LJ, Xi YQ, et al. Expression of stanniocalcin 2 in breast cancer and its clinical significance. Curr Med Sci (2019) 39(6):978–83. doi: 10.1007/s11596-019-2131-2
77. Raulic S, Ramos-Valdes Y, DiMattia GE. Stanniocalcin 2 expression is regulated by hormone signalling and negatively affects breast cancer cell viability in vitro. J Endocrinol (2008) 197(3):517–29. doi: 10.1677/JOE-08-0043
78. Zhang C, Chen S, Ma X, Yang Q, Su F, Shu X, et al. Upregulation of STC2 in colorectal cancer and its clinicopathological significance. Onco Targets Ther (2019) 12:1249–58. doi: 10.2147/OTT.S191609
79. Lin C, Sun L, Huang S, Weng X, Wu Z. STC2 is a potential prognostic biomarker for pancreatic cancer and promotes migration and invasion by inducing epithelial-mesenchymal transition. BioMed Res Int (2019) 2019:8042489. doi: 10.1155/2019/8042489
80. Jiang ZH, Shen X, Wei Y, Chen Y, Chai H, Xia L, et al. A pan-cancer analysis reveals the prognostic and immunotherapeutic value of stanniocalcin-2 (STC2). Front Genet (2022) 13:927046. doi: 10.3389/fgene.2022.927046
81. Huang K, Zeng Y, Xie Y, Huang L, Wu Y. Bioinformatics analysis of the prognostic value of CCT6A and associated signalling pathways in breast cancer. Mol Med Rep (2019) 19(5):4344–52. doi: 10.3892/mmr.2019.10100
82. Jiang J, Liu C, Xu G, Liang T, Yu C, Liao S, et al. CCT6A, a novel prognostic biomarker for Ewing sarcoma. Med (Baltimore). (2021) 100(4):e24484. doi: 10.1097/MD.0000000000024484
83. Zeng G, Wang J, Huang Y, Lian Y, Chen D, Wei H, et al. Overexpressing CCT6A contributes to cancer cell growth by affecting the G1-To-S phase transition and predicts a negative prognosis in hepatocellular carcinoma. Onco Targets Ther (2019) 12:10427–39. doi: 10.2147/OTT.S229231
84. Karimnia N, Wilson AL, Green E, Matthews A, Jobling TW, Plebanski M, et al. Chemoresistance is mediated by ovarian cancer leader cells in vitro. J Exp Clin Cancer Res (2021) 40(1):276. doi: 10.1186/s13046-021-02086-3
85. Bilandzic M, Rainczuk A, Green E, Fairweather N, Jobling TW, Plebanski M, et al. Keratin-14 (KRT14) positive leader cells mediate mesothelial clearance and invasion by ovarian cancer cells. Cancers (Basel) (2019) 11(9):1228. doi: 10.3390/cancers11091228
86. Cao Y, Han Q, Li J, Jia Y, Zhang R, Shi H. P4HA2 contributes to cervical cancer progression via inducing epithelial-mesenchymal transition. J Cancer. (2020) 11(10):2788–99. doi: 10.7150/jca.38401
87. Lin J, Jiang L, Wang X, Wei W, Song C, Cui Y, et al. P4HA2 promotes epithelial-to-Mesenchymal transition and glioma malignancy through the collagen-dependent PI3K/AKT pathway. J Oncol (2021) 2021:1406853. doi: 10.1155/2021/1406853
88. Hu D, Ansari D, Zhou Q, Sasor A, Said Hilmersson K, Andersson R. Low P4HA2 and high PRTN3 expression predicts poor survival in patients with pancreatic cancer. Scand J Gastroenterol (2019) 54(2):246–51. doi: 10.1080/00365521.2019.1574360
89. Zhu M, Peng R, Liang X, Lan Z, Tang M, Hou P, et al. P4HA2-induced prolyl hydroxylation suppresses YAP1-mediated prostate cancer cell migration, invasion, and metastasis. Oncogene (2021) 40(41):6049–56. doi: 10.1038/s41388-021-02000-3
90. Li M, Wang Q, Zheng Q, Wu L, Zhao B, Wu Y. Prognostic and diagnostic roles of prolyl 4-hydroxylase subunit alpha members in breast cancer. biomark Med (2021) 15(13):1085–95. doi: 10.2217/bmm-2020-0323
91. Xiong G, Deng L, Zhu J, Rychahou PG, Xu R. Prolyl-4-hydroxylase alpha subunit 2 promotes breast cancer progression and metastasis by regulating collagen deposition. BMC Cancer. (2014) 14:1. doi: 10.1186/1471-2407-14-1
92. Duncan L, Shay C, Teng Y. PGK1 : An essential player in modulating tumor metabolism. Methods Mol Biol (2022) 2343:57–70. doi: 10.1007/978-1-0716-1558-4_4
93. He Y, Luo Y, Zhang D, Wang X, Zhang P, Li H, et al. PGK1-mediated cancer progression and drug resistance. Am J Cancer Res (2019) 9(11):2280–302.
94. Hu H, Zhu W, Qin J, Chen M, Gong L, Li L, et al. Acetylation of PGK1 promotes liver cancer cell proliferation and tumorigenesis. Hepatology (2017) 65(2):515–28. doi: 10.1002/hep.28887
95. Li Y, Wang S, Zhang X, Yang R, Wei X, Yan R, et al. Expression characteristics and significant prognostic values of PGK1 in breast cancer. Front Mol Biosci (2021) 8:695420. doi: 10.3389/fmolb.2021.695420
96. Ye T, Liang Y, Zhang D, Zhang X. MicroRNA-16-1-3p represses breast tumor growth and metastasis by inhibiting PGK1-mediated warburg effect. Front Cell Dev Biol (2020) 8:615154. doi: 10.3389/fcell.2020.615154
97. Fu D, He C, Wei J, Zhang Z, Luo Y, Tan H, et al. PGK1 is a potential survival biomarker and invasion promoter by regulating the HIF-1alpha-Mediated epithelial-mesenchymal transition process in breast cancer. Cell Physiol Biochem (2018) 51(5):2434–44. doi: 10.1159/000495900
98. Tang W, Wu Y, Qi X, Yu R, Lu Z, Chen A, et al. PGK1-coupled HSP90 stabilizes GSK3beta expression to regulate the stemness of breast cancer stem cells. Cancer Biol Med (2021) 19(4):486–503. doi: 10.20892/j.issn.2095-3941.2020.0362
99. Sun S, Wu H, Wu X, You Z, Jiang Y, Liang X, et al. Silencing of PGK1 promotes sensitivity to paclitaxel treatment by upregulating XAF1-mediated apoptosis in triple-negative breast cancer. Front Oncol (2021) 11:535230. doi: 10.3389/fonc.2021.535230
100. Halestrap AP. The SLC16 gene family - structure, role and regulation in health and disease. Mol Aspects Med (2013) 34(2-3):337–49. doi: 10.1016/j.mam.2012.05.003
101. Gao M, Monian P, Quadri N, Ramasamy R, Jiang X. Glutaminolysis and transferrin regulate ferroptosis. Mol Cell (2015) 59(2):298–308. doi: 10.1016/j.molcel.2015.06.011
102. Wang Z, Li Y, Fan L, Zhao Q, Tan B, Liu R, et al. Silencing of TPD52 inhibits proliferation, migration, invasion but induces apoptosis of pancreatic cancer cells by deactivating akt pathway. Neoplasma (2020) 67(2):277–85. doi: 10.4149/neo_2019_190404N295
103. Roslan N, Bieche I, Bright RK, Lidereau R, Chen Y, Byrne JA. TPD52 represents a survival factor in ERBB2-amplified breast cancer cells. Mol Carcinog. (2014) 53(10):807–19. doi: 10.1002/mc.22038
104. Li G, Yao L, Zhang J, Li X, Dang S, Zeng K, et al. Tumor-suppressive microRNA-34a inhibits breast cancer cell migration and invasion via targeting oncogenic TPD52. Tumour Biol (2016) 37(6):7481–91. doi: 10.1007/s13277-015-4623-4
105. Yang Y, Hu H, Chen L, Zhang H, Yang J. A new survival model based on ferroptosis-related genes (FRGS) for prognostic prediction in bladder cancer. Actas Urol Esp (Engl Ed). (2022) 46(8):494–503. doi: 10.1016/j.acuro.2021.12.006
Keywords: hypoxia, ferroptosis, prognostic model, breast cancer, tumor microenvironment
Citation: Zhong J, Shen X, Zhou J, Yu H, Wang B, Sun J, Wang J and Liu F (2023) Development and validation of a combined hypoxia and ferroptosis prognostic signature for breast cancer. Front. Oncol. 13:1077342. doi: 10.3389/fonc.2023.1077342
Received: 22 October 2022; Accepted: 23 February 2023;
Published: 14 March 2023.
Edited by:
Marta Truffi, Scientific Clinical Institute Maugeri (ICS Maugeri), ItalyReviewed by:
Zeyan Li, Shandong University, ChinaXiongjian Rao, University of Kentucky, United States
Copyright © 2023 Zhong, Shen, Zhou, Yu, Wang, Sun, Wang and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jing Wang, 151680848@qq.com; Feng Liu, Lf5806@sina.com
†These authors have contributed equally to this work