An iron metabolism and immune related gene signature for the prediction of clinical outcome and molecular characteristics of triple-negative breast cancer
BMC Cancer volume 22, Article number: 619 (2022)
An imbalance of intracellular iron metabolism can lead to the occurrence of ferroptosis. Ferroptosis can be a factor in the remodeling of the immune microenvironment and can affect the efficacy of cancer immunotherapy. How to combine ferroptosis-promoting modalities with immunotherapy to suppress triple-negative breast cancer (TNBC) has become an issue of great interest in cancer therapy. However, potential biomarkers related to iron metabolism and immune regulation in TNBC remain poorly understand.
We constructed an optimal prognostic TNBC-IMRGs (iron metabolism and immune-related genes) signature using least absolute shrinkage and selection operator (LASSO) cox regression. Survival analysis and ROC curves were analyzed to identify the predictive value in a training cohort and external validation cohorts. The correlations of gene signature with ferroptosis regulators and immune infiltration are also discussed. Finally, we combined the gene signature with the clinical model to construct a combined model, which was further evaluated using a calibration curve and decision curve analysis (DCA).
Compared with the high-risk group, TNBC patients with low-risk scores had a remarkably better prognosis in both the training set and external validation sets. Both the IMRGs signature and combined model had a high predictive capacity, 1/3/5- year AUC: 0.866, 0.869, 0.754, and 1/3/5-yaer AUC: 0.942, 0.934, 0.846, respectively. The calibration curve and DCA also indicate a good predictive performance of the combined model. Gene set enrichment analysis (GSEA) suggests that the high-risk group is mainly enriched in metabolic processes, while the low-risk group is mostly clustered in immune related pathways. Multiple algorithms and single sample GSEA further show that the low-risk score is associated with a high tumor immune infiltration level. Differences in expression of ferroptosis regulators are also observed among different risk groups.
The IMRGs signature based on a combination of iron metabolism and immune factors may contribute to evaluating prognosis, understanding molecular characteristics and selecting treatment options in TNBC.
Breast cancer, an extremely heterogeneous disease, is the most common cancer for women worldwide . Triple-negative breast cancer (TNBC) is the most aggressive subtype of breast cancer with a high risk of distant relapse in the first 3 to 5 years following diagnosis . More importantly, since there is no powerful targeted therapy currently available, the overall survival of TNBC patients is worse than patients with other subtypes of breast cancer [3, 4]. Therefore, identification of new sensitive biomarkers to determine prognosis or develop novel therapies for TNBC patients remains an urgent clinical need.
Iron is an essential trace element for the human body, and either an excess or shortage can influence many physiological and pathological processes including energy metabolism, mitochondrial function and DNA synthesis and repair . It is reported that iron metabolism has dual effects in tumor growth and development. On the one hand, cancer cells are more dependent on iron intake for proliferation and more susceptible to iron deficiency than non-cancer cells, which is termed “iron addition” . On the other hand, highly increased iron concentrations lead to cell death driven by excessive reactive oxygen species (ROS) and lipid peroxidation, known as ferroptosis . Ferroptosis is a novel form of programmed cell death that is different from apoptosis, pyroptosis, necroptosis, and autophagy . TNBC cells are more prone to ferroptosis due to their complex metabolic characteristics and cellular signaling pathways, which make it a promising choice to overcome the refractory disease .
An increasing amount of evidence indicates that immunotherapy can be an effective treatment strategy for TNBC patients . Immune cells like macrophages and lymphocytes can play major roles in iron homeostasis through iron recycling of senescent erythrocytes  and non-transferrin bound iron uptake , respectively. In addition, a report has implicated a new mechanism by which CD8 + T cells repress tumor development through inducing ferroptosis . Based on the crosstalk between iron metabolism and the tumor microenvironment (TME), it would be of significance to identify any potential applications of iron metabolism and immune related genes in targeted therapy and immunotherapy for TNBC patients.
Compared with a single biomarker, multigene signatures have been proven to produce a higher accuracy prognosis . Previous studies have primarily focused on the prognostic value of the ferroptosis-related signature  and immune-related signature  in breast cancer. However, the identification of iron metabolism and immune-related signature in TNBC remains unclear.
The main purpose of this study was to identify the iron metabolism and immune-related genes (IMRGs) associated with the prognosis of TNBC through bioinformatics analysis, and then to build a prognosis model of TNBC patients according to LASSO regression analysis. The molecular characteristics of this model were identified by iron metabolism and tumor infiltrating lymphocytes. Finally, we combined the gene signature with the clinical model to construct a combined model, which had excellent predictive performance in terms of prognosis.
Materials and methods
Patient data sets
The gene expression profiles and clinical information of 119 normal samples and 123 TNBC samples were extracted from The Cancer Genome Atlas (TCGA) database (https://portal.gdc.cancer.gov/). The validation cohorts GSE2603 and GSE21653 (the available TNBC patient datasets with integral microarray data and survival data) were downloaded from the NCBI-GEOdatabase (https://www.ncbi.nlm.nih.gov/gds). There were no ethical conflicts because all data came from the public database.
Identification of IMRGs
To meet the high demand of iron, cancer cells have remodeled iron metabolism pathways, including acquisition, storage, and efflux, affecting iron homeostasis . Therefore, iron metabolism related genes were extracted from the gene sets “M15748: iron ion homeostasis”, “M11074: cellular iron ion homeostasis”, “M24450: cellular response to iron ion”, “M37743: abnormality of iron homeostasis”, “M18915: iron ion binding” and “M962: iron uptake and transport” from the Molecular Signatures Database (MSigDB). Ferroptosis related genes were obtained from the FerrDb database, “ferroptosis pathway (map04216)” from the KEGG PATHWAY Database  and “M39768: ferroptosis” from MSigDB. After integration of intersecting genes and elimination of unrelated genes, 499 iron metabolism and ferroptosis related genes were included in the follow-up studies. We also found 2483 immune related genes from the ImmPort database.
Differential expressed analysis
Differential expressed genes (DEGs) of TNBC were identified between 119 normal samples and 123 TNBC samples with the cut-off values: |log2 fold change|> 1 and P < 0.05. Differently expressed genes were obtained between high-risk and low-risk groups with |log2 fold change|> 1 and P < 0.05 as the threshold. The results were calculated using the ‘limma’ package.
Gene signature construction
The least absolute shrinkage and selection operator (LASSO) was used to construct the TNBC-IMRGs signature by TCGA transcriptome data. Risk score = -0.05076 × ExpEGR3-0.22386 × ExpCCND2-0.06754 × ExpSOCS3-0.06622 × ExpJunD-0.52219 × ExpSLC27A6.
Functional enrichment analysis
KEGG/GO and GSEA were analyzed by the ‘clusterProfiler’ package . The annotated gene sets of GSEA were selected, c2.cp.v7.2.symbols.gmt and c5.all.v7.2.symbols.gmt, from the MSigDB Collections (https://www.gsea-msigdb.org/gsea/msigdb/index.jsp), the number of permutations was set to 1000 times. The significance threshold was placed at FDR < 0.25, and P—value < 0.05.
Immune infiltration analysis
The ESTIMATE algorithm was determined by the ‘estimate’ package (version 1.0.13) . The association between risk score and immune cell infiltration levels were analyzed via TIMER , QUANTISEQ , MCPCOUNTER , XCELL , EPIC , and CIBERSPRT  algorithms from TIMER2.0 (https://cistrome.shinyapps.io/timer/). Single sample GSEA (ssGSEA) was performed with the ‘GSVA’ package (version 1.34.0) .
All R packages were executed using R Studio software (version 4.0.5). The ‘ggplot2’ R package (version 3.3.3) was used to visualize the volcano plot and heatmap. A Pearson correlation was used to describe the correlation analysis. Univariate and multivariate Cox regression analysis were used to identify the prognostic factors. Kaplan–Meier survival curves were compared using the log-rank test. The ROC analysis was performed with the ‘pROC’ package. Significance was statistically considered at ***P < 0.001, **P < 0.01, *P < 0.05, and ns, P ≥ 0.05.
Identification and validation of TNBC-IMRGs signature-based prognostic model
This study was conducted according to the flow chart presented in Fig. 1. The detailed clinicopathological features of TNBC patients from different cohorts are summarized in Table S1. A total of 2087 differently expressed genes (TNBC-DEGs) were identified between 119 normal samples and 123 TNBC samples from the TCGA dataset (Fig. 2a). We crossed over 499 iron metabolism and ferroptosis-related genes with 2483 immune-related genes to yielded 56 IMRGs. A Pearson correlation analysis was performed between TNBC-DEGs and IMRGs (Correlation Coefficient > 0.4 or < 0.4 and P < 0.01) and subsequently we acquired 1244 TNBC-IMRGs that had potential functions in both iron metabolism and tumor immunity. This dataset was integrated with the prognostic gene dataset from TCGA-BRCA (P.cox < 0.01) to obtain 30 candidate prognostic TNBC-IMRGs. The heatmap in Fig. 2b shows the expression profiles of these genes in TNBC and normal samples. Some TNBC-IMRGs show interrelationships in Fig. 2c. To further identify the best candidate genes, we performed LASSO cox regression to establish a TNBC-IMRGs signature-based prognostic model. When the model reached the minimum of lambda (λ), an optimal prognostic model with five non-zero coefficient genes (EGR3, CCND2, SOCS3, JunD, and SLC27A6) was constructed (Fig. 3a, 3b). A median risk score was calculated as follows: risk score = -0.05076 × ExpEGR3-0.22386 × ExpCCND2-0.06754 × ExpSOCS3-0.06622 × ExpJunD-0.52219 × ExpSLC27A6, and then we categorized patients into high-risk and low-risk groups. The distributions of the risk scores, survival time, survival status, and the expression patterns of five TNBC-IMRGs are displayed in Fig. 3c. As expected, all of these five protective genes were highly expressed in the low-risk group. Kaplan–Meier survival analysis showed that patients with low-risk scores exhibited a significantly better outcome than those with high-risk scores in the TCGA training set (P = 0.001, Fig. 3d). We used GSE2603 (Fig. 3f) and GSE21653 (Fig. 3i) as external validation sets to test the robustness of the model based on the TCGA-TNBC cohort. We consistently found that high-risk groups which were stratified using the same calculation formula as the training set were markedly related to worse prognosis in the validation sets (Fig. 3g, 3j). Next, through ROC curves, we evaluated the predictive accuracy of the prognostic model. The AUC of the TCGA cohort (Fig. 3e) at 1-, 3- and 5-years were 0.866, 0.869, and 0.754. The AUC at 1-, 3-, and 5-years of the GSE2603 cohort were 0.444, 0.652, and 0.772 (Fig. 3h) and for GSE21653 cohort were 0.643, 0.610, and 0.649 (Fig. 3k), respectively. These results indicate a very good prognostic value of this TNBC-IMRGs signature.
Stratified prognostic analysis of the TNBC-IMRGs signature
Stratified prognostic analysis suggests that the high-risk score is significantly associated with a worse OS in TNBC patients of an older age (> 55) (Figs. 4a-b) or smaller tumor size (T1-2) (Figs. 4e-f). Additionally, no matter the lymph node score (Figs. 4c-d) and stage status (Figs. 4g-h), patients with a low-risk score shown a better outcome than the high-risk group. This TNBC-IMRGs prognostic signature may act as an effective biomarker in assessing the prognosis of TNBC patients.
Functional enrichment analysis and immune infiltration analysis
To explore the potential molecular characteristics in different risk groups, differently expressed genes between the high-risk group and low-risk group were screened out. The volcano map is depicted to visualize 44 DEGs (Fig. 5a). GO enrichment analysis revealed that the DEGs were significantly clustered in extracellular matrix organization, response to metal ions, positive regulation lipid metabolic process, establishment of T cell polarity, dendritic cell chemotaxis, and others (Fig. 5b). The KEGG pathway analysis demonstrated that most DEGs were clustered in physiological signaling pathways such as the PPAR signaling pathway, estrogen signaling pathway, NF-kappa B signaling pathway, etc. (Fig. 5c). Additionally, using GSEA analysis we found that multiple metabolism processes including cellular response to hypoxia, glucose metabolism, nucleotides and fatty acid metabolism, and oxidative phosphorylation were highly clustered in the high-risk group (Fig. 5d). While in the low-risk group, immune-related responses were enriched; for example, T cell, B cell, and lymphocyte activation, as well as the response to interferon (Fig. 5e, 5f).
Recent studies have already investigated a series of key regulatory genes that affect ferroptosis through abnormal accumulation of iron, free radicals, and ROS. Our results show that a total of four ferroptosis drivers (e.g., ACSL4 ) were upregulated in the low-risk group, while five anti-ferroptosis regulators (e.g., GPX4 , and FTH1 ) were overexpressed in the high-risk group (Fig. 5g).
Infiltrating stromal and immune cells, which make up the majority of normal cells in tumor tissue, not only interfere with tumor signaling but also regulate cancer biology. To further determine the relation of biological function with the TNBC-IMRGs signature in the immune response, we made use of various algorithms to investigate the immune microenvironment landscape. The ESTIMATE algorithm suggests that there might be a weak trend toward a negative correlation between risk score and immune-related score at the macro-levels (Figs. 6a-c). Therefore, we applied other algorithms, including TIMER, QUANTISEQ, MCPCOUNTER, EPIC, and CIBERSORT to further explore specific infiltrating immune cells (Fig. 6d). The results indicate that immunoreactive cells (like B cells, CD4 + T cells, and T follicular helper cells) were more abundant in the low-risk group, while immunosuppressive cells (e.g., M2 macrophage) were enriched in the high-risk group. There was a moderate negative correlation between CD4 + T cells and risk score (r = -0.501, P < 0.001) (Fig. 6e). Subsequently, differences were distinguished between the high- and low- risk groups in the quantities of 16 types of immune cells through ssGSEA (Fig. 6f). Similarly, the enrichment scores of immunoreactive cells (e.g., B cells, DC, iDC, Mast cells, NK cells, TFH, and Tgd) were higher in the low-risk group compared with the high-risk group (P < 0.05).
Immune checkpoint inhibitors are antitumor immunotherapies that are increasing used in clinical practice. Differences in immune checkpoint gene expression between high- and low-risk groups may confer differential susceptibility to immune checkpoint inhibitors. In our observations, we found that expression of some checkpoint molecules such as VTCN1, CD160, HLA-DOB, and IL6 in the low-risk group were all higher than that in the high-risk group (Fig. 6g). Taken together, these results suggest that the TNBC-IMRGs prognostic signature may be related to some extent to iron metabolism and immune cell infiltration.
Construction and validation of the combined gene and clinical model
The forest plots are presented to illustrate the summary of the univariate and multivariate Cox regression analyses of some clinical features including risk score based on the TNBC-IMRGs signature (Fig. 7a). Node stage (P = 0.03) and risk score (P < 0.001) were determined to be robust independent prognostic predictors in the TCGA cohort. To further compare the superiority of gene-based signature, we constructed a clinical prognostic model that included Tumor stage and Node stage. The formula of the clinical risk score was as follows: clinical risk score = -1.18194 + 0.33509*Tumor stage + 1.01152*Node stage (Fig. 7b). The Fig. 7c depicts the AUC of the clinical model for 1-, 3-, and 5-year survival at 0.805, 0.791, and 0.657. It is obvious that the predictive power of the clinical model is not as good as the gene signature. We also established a combined model integrating the genomic risk score and clinical variables. The combined risk score was calculated as follows: combined risk score = 4.10626 + 0.32847*Tumor stage + 1.33369*Node stage + 2.94864*Risk Score (gene signature). The prognostic nomogram in the TCGA dataset is shown in Fig. 7d, and the C-index value for this combined model is 0.901. The combined model shows better prognostic performance than the gene signature and clinical model. There is a marked difference in prognosis between the high- and low-risk groups (Fig. 7f. P < 0.001). Likewise, the same trend is also observed in the validation sets; GSE2603 and GSE21653 (Fig. 8a, 8c). In addition, the AUC of the combined model for 1-, 3- and 5-year survival are 0.942, 0.934 and 0.846 (Fig. 7g). The calibration analysis of the 1-, 3- and 5- year outcome prediction is illustrated in Fig. 7e. The blue and red lines, not the green one, have a closer fit to the dotted gray line, revealing that the nomogram may do better predicting short-term prognosis than long-term prognosis. Additionally, decision curve analysis (DCA) indicates that the combined model has the best performance in prognosis prediction compared with the clinical model and gene signature (Fig. 7h). We also used two validation sets to evaluate this combined model. From Fig. 8, we can observe that the 1-, 3-, and 5-year AUC of the combined model in the GSE2603 set and GSE21653 set are higher than for the gene signature (0.444–0.779 vs. 0.444–0.772 and 0.657–0.735 vs. 0.610–0.649), which indicates the better prognostic performance of the combined model.
The imbalance of iron metabolism leads to excessive intracellular iron storage and the accumulation of lipid reactive oxygen species, which in turn induces ferroptosis. Intriguingly, some reports reveal that ferroptosis has been involved in the remodeling of the immune microenvironment and it can modify the efficacy of cancer immunotherapy. How to combine ferroptosis-promoting modalities with immunotherapy to produce a synergetic effect of tumor suppression has become a hot research topic. Recent evidence suggests that TNBC is particularly susceptible to ferroptosis, making it an attractive drug target. However, potential biomarkers related to both iron metabolism and immune regulation in TNBC have not been sufficiently examined.
There are already a variety of multigene signatures applied in breast cancer, for example Oncotype Dx® (21-genes) , Mammaprint (70-genes) , Breast Cancer Index (BCI)  and EndoPredict® (EP) . However, the multigene prognostic model is poorly designed and validated in TNBC. As far as we know, our study is the first one to synthetically identify prognostic IMRGs in TNBC. We constructed a five-IMRGs signature and a combined prognostic model integrating the IMRGs signature with Tumor stage and Node stage to predict the likelihood of adverse events in TNBC patients. Compared with the high-risk group, TNBC patients with low-risk scores had remarkably better outcomes in both the training set and external validation sets. Both the IMRGs signature and the combined model had higher AUC values for 1-, 3- and 5-year survival than the clinical model.
TNBC cells are prone to ferroptosis, which is determined by intracellular metabolic processes that include lipid metabolism, iron metabolism, and amino acid metabolism [35, 36]. In our study, we found the high-risk group mostly enriched in various metabolic pathways (Fig. 5d). Apart from the common crucial regulatory genes of ferroptosis presented in Fig. 5g, five IMGRs (EGR3, SOCS3, JunD, SLC27A6, and CCND2) which were used to construct the gene signature were also shown to be partially involved in the process of ferroptosis. SOCS3, a member of the suppressor of cytokine signaling (SOCS) family, has been found to negatively regulate the upregulation of hepcidin which is a key regulator of iron metabolism . More importantly, it has been shown to participate in tumor suppression and regulate ferroptosis by interacting with p53 . Ferritin is a major cytosolic storage protein for iron that chelates excessive free iron, thereby reducing the occurrence of iron-mediated oxidative stress. Tsuji, Y reported that JunD was able to bind to the ferritin heavy chain and activate its transcription, which then protected cells from ROS in hepatocarcinoma cells . Ni et al. found that cyclin D2(CCND2) downregulation promotes cardiomyocyte ferroptosis . Whether EGR3 or SLC27A6 is involved in iron metabolism and ferroptosis has not been reported, our finding may provide clues for an underlying mechanism.
Since immune-related genes (EGR3 , CCND2 , SOCS3 , and JunD ) were also included in our IMRGs signature, we hypothesized that patients in different risk groups would have different proportion of infiltrating immune cells. Through various algorithms (Fig. 6), we found that the IMRGs risk scores tended to be negatively linked to immune cell infiltration levels, suggesting that increased IMRGs risk was associated with fewer infiltrating stromal and immune cells and higher tumor cells component . In line with the tumor immunoediting hypothesis, some cancers may be eliminated by antitumor immune responses, whereas others develop due to evasion or partial control of immune surveillance . Clearly, the high-risk group had a higher immunosuppression signature (e.g., M2 macrophage) and lower antitumor immune cell populations (e.g., B cells, NK cells, and T follicular helper cells) in the tumor microenvironment (Figs. 6d-f), which coincided with tumor progression and a worse prognosis for the high-risk group. Interestingly, most algorithms revealed that the B cell infiltration level was upregulated in the low-risk group (Fig. 6d). Lglesia et al. reported that the improved metastasis-free/progression-free survival was related to the B cell gene signature, mainly observed in basal-like and HER2-enriched breast cancer subtypes , which is in accordance with our results. We found that some immune checkpoint molecules expressed in immune cells (e.g., VTCN1 and CD160) were upregulated in the low-risk group, which may be driven by the increased infiltration of immune cells such as T cells and NK cells (Fig. 6g). However, due to the lower proportion of immunoreactive cells and lower expression of immune checkpoint genes, the high-risk group with a poorer prognosis may benefit less from immune checkpoint inhibitors than the low-risk group. A report has indicated that IFN-γ induces cell ferroptosis by activating the JAK1-2/STAT1/SLC7A11 signaling pathway. Meanwhile, treatment with IFN-γ reduces the level of GPX4, making cells more sensitive to ferroptosis . In our data, we also found that IFN-γ related responses were enriched in the low-risk group (Fig. 5f). Maybe combination therapies with IFN-γ will be a basis for treatment of TNBC patients.
However, our study is not without limitations. Firstly, our study was mainly based on retrospective data of TCGA and GEO datasets, and prospective studies of the gene signature and combined model in multi-center cohorts will be necessary. Secondly, further biological experiments in vitro and vivo should be performed to reveal the potential mechanisms of IMRGs in TNBC.
In conclusion, our study established a predictive signature associated with iron metabolism and tumor immunity that can accurately distinguish TNBC patients with different clinical outcomes. Moreover, this study provides evidence to support the combination of ferroptosis-induced cell death and antitumor immune in future clinical therapy in TNBC patients.
Availability of data and materials
The data supporting this study's findings are openly available in the TCGA dataset and GEO database.
Triple-negative breast cancer
Iron metabolism and immune-related genes
Reactive oxygen species
The Cancer Genome Atlas
Gene expression omnibus
Differential expressed genes
Least absolute shrinkage and selection operator
Early growth response 3
Suppressor of cytokine signaling 3
Solute carrier family 27 member 6
Receiver operating characteristic
Area under the curve
Decision curve analysis
Kyoto Encyclopedia of Genes and Genomes
Gene set enrichment analysis
Single sample GSEA
Acyl-CoA synthetase long chain family member 4
Activin A receptor type 1B
Ferritin heavy chain 1
Glutathione peroxidase 4
Interferon − gamma
- TIMER 2.0:
Tumor immune estimation resource, version 2
Siegel R, Miller K, Jemal A. Cancer statistics, 2020. CA: a cancer journal for clinicians. 2020;70(1):7–30.
Foulkes W, Smith I, Reis-Filho J. Triple-negative breast cancer. N Engl J Med. 2010;363(20):1938–48.
Liedtke C, Mazouni C, Hess K, André F, Tordai A, Mejia J, Symmans W, Gonzalez-Angulo A, Hennessy B, Green M, et al. Response to neoadjuvant therapy and long-term survival in patients with triple-negative breast cancer. Journal of clinical oncology : official journal of the American Society of Clinical Oncology. 2008;26(8):1275–81.
Tan D, Marchió C, Jones R, Savage K, Smith I, Dowsett M, Reis-Filho J. Triple negative breast cancer: molecular profiling and prognostic impact in adjuvant anthracycline-treated patients. Breast Cancer Res Treat. 2008;111(1):27–44.
Muckenthaler M, Rivella S, Hentze M, Galy B. A Red Carpet for Iron Metabolism. Cell. 2017;168(3):344–61.
Greten F. The Irony of Tumor-Induced Inflammation. Cell Metab. 2016;24(3):368–9.
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.
Galluzzi L, Vitale I, Aaronson S, Abrams J, Adam D, Agostinis P, Alnemri E, Altucci L, Amelio I, Andrews D, et al. Molecular mechanisms of cell death: recommendations of the Nomenclature Committee on Cell Death 2018. Cell Death Differ. 2018;25(3):486–541.
Li Z, Chen L, Chen C, Zhou Y, Hu D, Yang J, Chen Y, Zhuo W, Mao M, Zhang X, et al. Targeting ferroptosis in breast cancer. Biomarker research. 2020;8(1):58.
Keenan TE, Tolaney SM. Role of Immunotherapy in Triple-Negative Breast Cancer. J Natl Compr Canc Netw. 2020;18(4):479–89.
Andrews N, Schmidt P. Iron homeostasis. Annu Rev Physiol. 2007;69:69–85.
Arezes J, Costa M, Vieira I, Dias V, Kong X, Fernandes R, Vos M, Carlsson A, Rikers Y, Porto G, et al. Non-transferrin-bound iron (NTBI) uptake by T lymphocytes: evidence for the selective acquisition of oligomeric ferric citrate species. PLoS ONE. 2013;8(11): e79870.
Tang R, Xu J, Zhang B, Liu J, Liang C, Hua J, Meng Q, Yu X, Shi S. Ferroptosis, necroptosis, and pyroptosis in anticancer immunity. J Hematol Oncol. 2020;13(1):110.
Puppe J, Seifert T, Eichler C, Pilch H, Mallmann P, Malter W. Genomic Signatures in Luminal Breast Cancer. Breast care (Basel, Switzerland). 2020;15(4):355–65.
Wang D, Wei G, Ma J, Cheng S, Jia L, Song X, Zhang M, Ju M, Wang L, Zhao L, et al. Identification of the prognostic value of ferroptosis-related gene signature in breast cancer patients. BMC Cancer. 2021;21(1):645.
Wang X, Su W, Tang D, Jing J, Xiong J, Deng Y, Liu H, Ma W, Liu Z, Zhang Q. An Immune-Related Gene Prognostic Index for Triple-Negative Breast Cancer Integrates Multiple Aspects of Tumor-Immune Microenvironment. Cancers. 2021;13(21):5342.
Wang Y, Yu L, Ding J, Chen Y. Iron Metabolism in Cancer. Int J Mol Sci. 2018;20(1):95.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics-a Journal of Integrative Biology. 2012;16(5):284–7.
Yoshihara K, Shahmoradgoli M, Martínez E, Vegesna R, Kim H, Torres-Garcia W, Treviño V, Shen H, Laird P, Levine D, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612.
Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, Li B, Liu XS. TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res. 2017;77(21):e108–10.
Finotello F, Mayer C, Plattner C, Laschober G, Rieder D, Hackl H, Krogsdam A, Loncova Z, Posch W, Wilflingseder D, et al. Molecular and pharmacological modulators of the tumor immune contexture revealed by deconvolution of RNA-seq data. Genome medicine. 2019;11(1):34.
Becht E, Giraldo N, Lacroix L, Buttard B, Elarouci N, Petitprez F, Selves J, Laurent-Puig P, Sautès-Fridman C, Fridman W, et al. Estimating the population abundance of tissue-infiltrating immune and stromal cell populations using gene expression. Genome Biol. 2016;17(1):218.
Aran D, Hu Z, Butte A. xCell: digitally portraying the tissue cellular heterogeneity landscape. Genome Biol. 2017;18(1):220.
Racle J, Gfeller D. EPIC: A Tool to Estimate the Proportions of Different Cell Types from Bulk Gene Expression Data. Methods in molecular biology (Clifton, NJ). 2020;2120:233–48.
Chen B, Khodadoust M, Liu C, Newman A, Alizadeh A. Profiling Tumor Infiltrating Immune Cells with CIBERSORT. Methods in molecular biology (Clifton, NJ). 2018;1711:243–59.
Bindea G, Mlecnik B, Tosolini M, Kirilovsky A, Waldner M, Obenauf AC, Angell H, Fredriksen T, Lafontaine L, Berger A, et al. Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer. Immunity. 2013;39(4):782–95.
Yuan H, Li X, Zhang X, Kang R, Tang D. Identification of ACSL4 as a biomarker and contributor of ferroptosis. Biochem Biophys Res Commun. 2016;478(3):1338–43.
Yang W, SriRamaratnam R, Welsch M, Shimada K, Skouta R, Viswanathan V, Cheah J, Clemons P, Shamji A, Clish C, et al. Regulation of ferroptotic cancer cell death by GPX4. Cell. 2014;156:317–31.
Tian Y, Lu J, Hao X, Li H, Zhang G, Liu X, Li X, Zhao C, Kuang W, Chen D, et al. FTH1 Inhibits Ferroptosis Through Ferritinophagy in the 6-OHDA Model of Parkinson’s Disease. Neurotherapeutics : the journal of the American Society for Experimental NeuroTherapeutics. 2020;17(4):1796–812.
Siow Z, De Boer R, Lindeman G, Mann G. Spotlight on the utility of the Oncotype DX breast cancer assay. Int J Women’s Health. 2018;10:89–100.
Buyse M, Loi S, van't Veer L, Viale G, Delorenzi M, Glas AM, d'Assignies MS, Bergh J, Lidereau R, Ellis P, et al. Validation and clinical utility of a 70-gene prognostic signature for women with node-negative breast cancer. J Natl Cancer Inst. 2006;98(17):1183–1192.
Noordhoek I, Treuner K, Putter H, Zhang Y, Wong J, Meershoek-Klein Kranenbarg E, Duijm-de Carpentier M, van de Velde C, Schnabel C, Liefers G. Breast Cancer Index Predicts Extended Endocrine Benefit to Individualize Selection of Patients with HR Early-stage Breast Cancer for 10 Years of Endocrine Therapy. Clinical cancer research : an official journal of the American Association for Cancer Research. 2021;27(1):311–9.
Filipits M, Rudas M, Jakesz R, Dubsky P, Fitzal F, Singer C, Dietze O, Greil R, Jelen A, Sevelda P, et al. A new molecular predictor of distant recurrence in ER-positive, HER2-negative breast cancer adds independent information to conventional clinical risk factors. Clinical cancer research : an official journal of the American Association for Cancer Research. 2011;17(18):6012–20.
Stockwell B, Jiang X. The Chemistry and Biology of Ferroptosis. Cell Chem Biol. 2020;27(4):365–75.
Zheng J, Conrad M. The Metabolic Underpinnings of Ferroptosis. Cell Metab. 2020;32(6):920–37.
Miyachi H, Kobayashi Y, Relja B, Fujita N, Iwasa M, Gabazza E, Takei Y. Effect of suppressor of cytokine signaling on hepcidin production in hepatitis C virus replicon cells. Hepatology research : the official journal of the Japan Society of Hepatology. 2011;41(4):364–74.
Lessard F, Saint-Germain E, Mignacca L, Ferbeyre G. SOCS1: phosphorylation, dimerization and tumor suppression. Oncoscience. 2019;6:386–9.
Tsuji Y. JunD activates transcription of the human ferritin H gene through an antioxidant response element during oxidative stress. Oncogene. 2005;24(51):7567–78.
Ni T, Huang X, Pan S, Lu Z. Inhibition of the long non-coding RNA ZFAS1 attenuates ferroptosis by sponging miR-150-5p and activates CCND2 against diabetic cardiomyopathy. J Cell Mol Med. 2021;25(21):9995–10007.
Taefehshokr S, Key Y, Khakpour M, Dadebighlu P, Oveisi A. Early growth response 2 and Egr3 are unique regulators in immune system. Central-European journal of immunology. 2017;42(2):205–9.
Chiles T. Regulation and function of cyclin D2 in B lymphocyte subsets. J Immunol (Baltimore, Md : 1950). 2004;173(5):2901–2907.
Yoshimura A, Ito M, Chikuma S, Akanuma T, Nakatsukasa H. Negative Regulation of Cytokine Signaling in Immunity. Cold Spring Harb Perspect Biol. 2018;10(7):a028571.
Davidson C, Cameron L, Burshtyn D. The AP-1 transcription factor JunD activates the leukocyte immunoglobulin-like receptor 1 distal promoter. Int Immunol. 2014;26(1):21–33.
Dunn G, Bruce A, Ikeda H, Old L, Schreiber R. Cancer immunoediting: from immunosurveillance to tumor escape. Nat Immunol. 2002;3(11):991–8.
Iglesia M, Vincent B, Parker J, Hoadley K, Carey L, Perou C, Serody J. Prognostic B-cell signatures using mRNA-seq in patients with subtype-specific breast and ovarian cancer. Clinical cancer research : an official journal of the American Association for Cancer Research. 2014;20(14):3818–29.
Wei T, Zhang M, Zheng X, Xie T, Wang W, Zou J, Li Y, Li H, Cai J, Wang X, et al. Interferon-γ induces retinal pigment epithelial cell Ferroptosis by a JAK1–2/STAT1/SLC7A11 signaling pathway in Age-related Macular Degeneration. The FEBS journal. 2022;289(7):1968-1983.
This study was supported by the fund of the Provincial Natural Science Foundation of Fujian, China (2021J01751). The authors acknowledge the comments of the reviewers that helped to improve the study.
The Provincial Natural Science Foundation of Fujian, China (2021J01751).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Li, XF., Fu, WF., Zhang, J. et al. An iron metabolism and immune related gene signature for the prediction of clinical outcome and molecular characteristics of triple-negative breast cancer. BMC Cancer 22, 619 (2022). https://doi.org/10.1186/s12885-022-09679-x
- Iron metabolism
- Immune infiltration
- Gene signature
- Triple-negative breast cancer