Identification of the prognostic value of ferroptosis-related gene signature in breast cancer patients
BMC Cancer volume 21, Article number: 645 (2021)
Breast cancer (BRCA) is a malignant tumor with high morbidity and mortality, which is a threat to women’s health worldwide. Ferroptosis is closely related to the occurrence and development of breast cancer. Here, we aimed to establish a ferroptosis-related prognostic gene signature for predicting patients’ survival.
Gene expression profile and corresponding clinical information of patients from The Cancer Genome Atlas (TCGA) database and Gene Expression Omnibus (GEO) database. The Least absolute shrinkage and selection operator (LASSO)-penalized Cox regression analysis model was utilized to construct a multigene signature. The Kaplan-Meier (K-M) and Receiver Operating Characteristic (ROC) curves were plotted to validate the predictive effect of the prognostic signature. Gene Ontology (GO) and Kyoto Encyclopedia of Genes, Genomes (KEGG) pathway and single-sample gene set enrichment analysis (ssGSEA) were performed for patients between the high-risk and low-risk groups divided by the median value of risk score.
We constructed a prognostic signature consisted of nine ferroptosis-related genes (ALOX15, CISD1, CS, GCLC, GPX4, SLC7A11, EMC2, G6PD and ACSF2). The Kaplan-Meier curves validated the fine predictive accuracy of the prognostic signature (p < 0.001). The area under the curve (AUC) of the ROC curves manifested that the ferroptosis-related signature had moderate predictive power. GO and KEGG functional analysis revealed that immune-related responses were largely enriched, and immune cells, including activated dendritic cells (aDCs), dendritic cells (DCs), T-helper 1 (Th1), were higher in high-risk groups (p < 0.001). Oppositely, type I IFN response and type II IFN response were lower in high-risk groups (p < 0.001).
Our study indicated that the ferroptosis-related prognostic signature gene could serve as a novel biomarker for predicting breast cancer patients’ prognosis. Furthermore, we found that immunotherapy might play a vital role in therapeutic schedule based on the level and difference of immune-related cells and pathways in different risk groups for breast cancer patients.
Breast cancer (BRCA) is the most common cancer among women worldwide, with high morbidity and mortality rates . According to breast cancer statistics, there are expected to be about 276,480 new cases of invasive breast cancer and 42,170 deaths for women in the United States in 2020 . Despite major advances in therapies, including surgery, radiation therapy, chemotherapy, hormonal therapy, targeted and immunotherapy, the mortality rate for breast cancer still remains high . It has been reported that the immune microenvironment was crucial to the development of breast cancer, especially the infiltration of immune cells. These immune cells either expressed different immune antigens for themselves or influenced other immune cells to help tumor cells escape immunity, but specific mechanism still isn’t very clear . Therefore, to determine the molecular mechanism of breast cancer occurrence and development is crucial to advance cancer therapies.
Ferroptosis is a novel form of regulated cell death characterized by destruction of intracellular redox balance and non-apoptosis . Ferroptosis has become a promising therapeutic option for triggering cancer cell death . It is reported that siramesine and lapatinib are effective ferroptosis inducers in breast cancer [7, 8]. In addition, Dihydroisonone I (DT), a pure compound present in Salvia miltiorrhiza, can inhibit GPX4 protein expression and induce ferroptosis through lipid peroxidation to improve the prognosis of breast cancer . In terms of mechanism, elastin induces ferroptosis in breast cancer via the Glycogen synthase kinase-3β (GSK3β) / (nuclear factor erythroid 2-related factor 2) Nrf2 signaling pathway . Besides, acyl-CoA synthetase long-chain family member 4 (ACSL4) can increase the content of intracellular lipids to promote ferroptosis . Another study shows that the down-regulation of activator of transcription factor 2 (ATF2) promote ferroptosis by a negative feedback manner . However, whether these ferroptosis-related genes are associated with the prognosis of BRCA patients remains non-statistical. And immunological therapy is also not ignored for breast cancer patients.
Immune checkpoint blockade therapy has been used for a variety of cancers, including breast cancer, most of which are specific to CTLA-4 and PD-1/ PD-L1 . And related drugs have been developed for clinical application, including Tremelimumab, Ipilimumab, Avelumab, Atezolizumab and Pembrolizumab. Meanwhile, now a number of new immune checkpoint inhibitors have emerged, such as LAG-3, TIM-3 and TIGIT. And targeted drugs were also in clinical trials . However, sometimes a single treatment never met our expectations for the desired effect. Thus, we have developed a combination of therapies to combat the onset and progression of tumors, such as chemotherapy and immunotherapy or radiotherapy and immune checkpoint blocking therapy .
In this study, we first acquired mRNA expression data and corresponding clinical profiles of BRCA patients from The Cancer Genome Atlas (TCGA). Then, we constructed a prognostic multigene signature with ferroptosis-related differentially expressed genes (DEGs), and validated the predictive power of prognostic signature. Meanwhile, we performed GEO database to verify the credibility. Finally, we further conducted functional enrichment analysis to explore the potential immune-related mechanisms.
Gene expression data (count) and corresponding clinical information of 1223 breast cancer patients were obtained from The Cancer Genome Atlas (TCGA) up to September 08, 2020 (https://portal.gdc.cancer.gov/repository). This cohort has 1097 breast cancer patients with the associated gene expression profiles and clinical characteristics. Then, 21 patients were removed due to transcriptomic and clinical data was incomplete. Thus, the remaining data (n = 1076) with complete follow-up information was included in our training data set for further analyses. The testing data set for validation was downloaded from the Gene Expression Omnibus database (GEO, https://www.ncbi.nlm.nih.gov/geo/). GSE42568 was conducted by GPL570 (Affymetrix Human Genome U133 Plus 2.0 Array), including 104 tumor samples with breast cancer and 17 non-tumor samples as a normal control. The ICGC database (https://icgc.org/) was also downloaded to verify the reliability of the model. Furthermore, 60 ferroptosis-related genes were retrieved from the previous literature and are provided in Supplementary Table S1 .
Construction and validation of a prognostic ferroptosis-related gene signature
The “limma” R package was performed to ascertain DEGs related to ferroptosis between tumor tissues and non-tumor tissues with a false discovery rate (FDR) < 0.05 in the TCGA cohort. Univariate Cox analysis of overall survival (OS) was used to screen ferroptosis-related genes with prognostic values and was visualized by Forest plots. The intersection of ferroptosis-related DEGs and prognostic genes was demonstrated by Venn diagram and was visualized by heatmap. An interaction network of prognostic DEGs was generated by Search Tool for the Retrieval of Interacting Genes (STRING) database (https://string-db.org). We performed Human Protein Atlas (HPA) database (http://www.proteinatlas.org/) to evaluated the expression of DEGs. The Least absolute shrinkage and selection operator (LASSO)-penalized Cox regression analysis was used to construct a prognostic model for minimizing the risk of overfitting by performing the function “glmnet” of R package. Subsequently, the risk score of patients was calculated based on gene expression and corresponding Cox regression coefficient as follows: score = esum (each gene’s expression × corresponding coefficient) (Table S2). Then, patients were divided into high-risk and low-risk groups according to the median risk score values. Based on the expression of genes signature, PCA was performed with the “prcomp” function of the “stats” R package. Furthermore, t-SNE were applied to explore the distribution of different groups using the “Rtsne” R package. Survival analysis between high-risk and low-risk groups was carried out by the “surv_cutpoint” function of the “survminer” R package. Time-dependent ROC curve analyze was performed by “survivalROC” R package to evaluate the predictive accuracy of the genes signature.
Functional enrichment analysis
Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed for patients between the high-risk and low-risk groups by using the “clusterProfiler” R package. GO terms and KEGG pathways with P values < 0.05 were statistically significant. The infiltrating score of 16 immune cells and the activity of 13 immune-related pathways were determined by the “single-sample gene set enrichment analysis (ssGSEA) “ function of the “gsva” R package (Table S3).
Statistical analyses were performed by using R (version 4.0.2) software packages. Perl language was used for data matrix and all data processing. A two-tailed P < 0.05 was considered statistically significant.
The detailed flowchart is shown in Fig. 1. In this study, a total of 1076 BRCA patients from the TCGA cohort were finally enrolled. The general clinical information of those patients was provided in Table 1.
Identification of prognostic ferroptosis-related DEGs in TCGA
Most of the ferroptosis-related genes (51/60, 85%) were differentially expressed between tumor tissues and adjacent nontumorous tissues, and all of them were associated with OS in the univariate Cox regression analysis (Fig. 2a). Upregulated genes, including 15-lipoxygenase (ALOX15), CDGSH iron sulfur domain 1 (CISD1), citrate synthase (CS), glutamate-cysteine ligase catalytic subunit (GCLC), cystine/glutamate antiporter solute carrier family 7 member11 (SLC7A11), estrogen receptor membrane complex 2 (EMC2), squalene epoxidase (SQLE) and glucose-6-phosphate dehydrogenase (G6PD), manifested an excellent prognosis in heatmap (Fig. 2b) and univariate Cox analysis (Fig. 2c). The protein-protein interaction network among these genes illuminates that GPX4 and G6PD were the hub gene (Fig. 2d).
The expression of six-gene is higher in BRCA tissue compared with normal tissue in HPA database
According to the immunohistochemical analyses in HPA database, the high staining intensity of these six genes (ALOX15, CS, GCLC, EMC2, SQLE, G6PD) in BRCA tissues contrasted starkly with the low intensity or lack of staining in normal tissues (Fig. 3a-f), while GPX4 and ACSF2 (Fig. 3g-h) didn’t show striking difference. Regrettably, CISD1 and SLC7A11 didn’t be founded in HPA database.
Construction and evaluation of nine-gene prognostic model in TCGA
The LASSO Cox regression analysis was applied to establish a prognostic model using the expression profile of the 10 genes mentioned above. Finally, a nine-gene (15-lipoxygenase (ALOX15), CDGSH iron sulfur domain 1 (CISD1), citrate synthase (CS), glutamate-cysteine ligase catalytic subunit (GCLC), selenoenzyme glutathione peroxidase (GPX4), cystine/glutamate antiporter solute carrier family 7 member11 (SLC7A11), estrogen receptor membrane complex 2 (EMC2), glucose-6-phosphate dehydrogenase (G6PD) and acyl-CoA synthetase family member2 (ACSF2)) prognostic model was constructed to predict prognosis based on the risk score = e (0.167 * expression level of ALOX15 + 0.212 * expression level of CISD1 + 0.228 * expression level of Cs + 0.057 * expression level of GCLC + (− 0.116) * expression level of GPX4 + 0.018 * expression level of SLC7A11 + 0.211 * expression level of EMC2 + 0.134 * expression level of G6PD + (− 0.098) * expression level of ACSF2). The patients were divided into a high-risk group (n = 535) or a low-risk group (n = 535) according to the median value of risk score in TCGA cohort (Fig. 4a), Principal component analysis (PCA) and t-distributed stochastic neighbor embedding (t-SNE) analysis suggested the patients in different risk groups were distributed in two directions (Fig. 4b-c). The patients in high-risk group had a higher probability of death earlier than those in low-risk group (Fig. 4d). Analogously, the Kaplan-Meier curve showed the prognostic signature clearly distinguished patients with high and low survival rate (Fig. 4e). The area under the curve (AUC) confirmed that the identified prognostic signature had a robust efficiency for predicting the OS for BRCA patients (AUC = 0.618, 0.653 and 0.663; at 1, 2 and 3 year, respectively, Fig. 4f).
Validation of the nine-gene signature in GEO and ICGC dataset
The breast cancer patients from the GEO and ICGC cohort were categorized into high-risk or low-risk groups by the median value of risk score (Fig. 5a). The results of PCA and t-SNE analysis in GEO and ICGC was similar to TCGA (Fig. 5b-c). The high-risk group had a poor prognosis compared to the low-risk group (Fig. 5d). the KM curve showed the prognostic signature clearly distinguished patients with different survival rate (Fig. 5e). The AUC of ROC indicated nine-gene signature had a moderate predictive ability (AUC = 0.621, 0.644 and 0.572; at 1, 2 and 3 year, respectively in GEO. AUC = 0.871, 0.843 and 0.505; at 1, 2 and 3 year, respectively in ICGC, Fig. 5f).
Independent prognostic value of the nine-gene signature in TCGA
In order to investigate the independence of prognostic signature for OS, we performed univariate and multivariate Cox regression tests to determine the relationship between the risk model and different clinicopathological parameters (Fig. 6a-b). The risk score was significantly associated with OS in univariate Cox regression analyses (HR = 3.584, 95% CI = 2.353–5.457, P < 0.001). Consistently, the risk score still indicated to be an independent predictor for OS in the multivariate Cox regression analysis (HR = 3.145, 95% CI = 2.087–4.738, P < 0.001) (Table 2).
Immune-related functions and pathways were enriched in GO and KEGG
GO enrichment and KEGG pathway analyses were performed to elucidate the biological functions and pathways that were associated with the risk score. GO analysis results showed that several immune -related molecular functions were significantly enriched (Fig. 7a-b). Likewise, KEGG analyses still found that genes were significantly enriched in IL-17 signaling and cytokine-cytokine receptor interaction pathway (Fig. 7c-d).
The distribution of immune-associated cells and processes in high- and low-risk group
To further explore the relevancy between the risk score and immune status, we quantified the enrichment scores of diverse immune cell subpopulations, associated functions or pathways with ssGSEA. Contents of the antigen presentation process, including the score of activated dendritic cells (aDCs), dendritic cells (DCs), antigen presenting cell (APC) co-stimulation, inflammation-promoting and MHC class-I (Fig. 7a-b), were significantly different between the high-risk and low-risk group. Moreover, the score of macrophages, natural killer cell (NK), T follicular helper cell (Tfh), T-helper 1 (Th1), T-cell-co − inhibition, type I IFN response, and type II IFN response were higher in the high-risk group, while the scores of mast cells was just the opposite (Fig. 8a-b). It is reported that the occurrence of breast cancer is closely related to the infiltration of immune cells. In the early stage, immune cells played a role in eliminating tumors and gradually were depleted, eventually leading to the weakening of the immune response [15, 16]. Our results might explain the occurrence of early breast cancer, the greater the threat of tumors to human body, the more serious the immune response.
Breast cancer is an extremely common malignancy in women and has a poor prognosis. In recent years, new biomarkers have emerged for the prognosis of cancer patients , but this is the first time to discuss the establishment of a prognostic model related to ferroptosis in breast cancer patients.
The prognostic model presented in this study consist of nine ferroptosis-related genes (ALOX15, CISD1, CS, GCLC, GPX4, SLC7A11, EMC2, G6PD and ACSF2). Many studies reported that these genes played an significant role in the pathogenesis of breast cancer. ALOX15 mediates lymphatic vessel invasion and lymph node metastasis in human breast cancer xenograft mouse . CISD1 promotes human breast cancer proliferation and confers autophagy resistance [19, 20]. CS restrains aggressive triple-negative breast cancer cells by targeting glycolysis and the cancer stem cell phenotype . 15-deoxy-Δ (12,14)-prostaglandin J2-induced GCLC is mediated by Multidrug resistance-associated protein 1 via Nrf2 signaling in human breast cancer cells . GPX4, as an oncogene, inhibits the ferroptosis effect of cancer cells, while GPX4 inhibition can enhance the anticancer effect of cisplatin . SLC7A11 over-expression promotes lipo-ROS accumulation in MCF-7 breast cancer cells . The blockade of G6PD by autophagy enhanced the inhibitory effect of tyrosine kinase on breast cancer cells . Although these genes are closely related to the pathogenesis of breast cancer, for the first time we have combined them as a marker of prognosis for BRCA patients.
In this study, we established nine genetic biomarkers as a new prognostic model and analyzed their ability to predict the prognosis of high-risk and low-risk groups. The prognostic performance of the model was verified by KM curve and ROC curve, and the results showed that the model had good predictive performance. Meanwhile, we also drew forest plot with different clinical parameters, and the results showed that there were significant correlations between them. Prognostic indicators were also good independent indicators of survival after adjustment for clinical parameters, including age (≤60 years or >60 years), cancer stage (stage I-II or stage II-IV), risk Score (high or low). The above confirmed that the model can predict the prognostic characteristics of breast cancer patients.
Based on DEGs between different risk groups, we performed GO and KEGG functional analyses, and found unexpectedly that many immune-related biological processes and pathways were enriched [16, 26]. Therefore, we further explored the function of immune cell function between the high-risk and low-risk groups and different subgroups. The results showed that most of immune-related cells and functions possessed a higher risk score. This might be related to the overactivation of the immune system in the early stages of breast cancer, for specific performance that the level of antigen presentation and Th1 cells were higher in high-risk groups, which suggested that relative genes were also significantly altered. Therefore, these targets might provide a possibility for immunotherapy of breast cancer patients.
In conclusion, this study established a new prognostic model associated with nine ferroptosis-related genes and the good prediction ability of the model was verified by three databases, including TCGA、GEO、ICGC database. Besides, we found immune-related cells and pathways were significant differences in high- and low-risk group, which might be helpful for illustrating the application of immunotherapy for breast cancer patients.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
False discovery rate
The Cancer Genome Atlas
Gene Expression Omnibus
International Cancer Genome Consortium
Human Protein Atlas
Search Tool for the Retrieval of Interacting Gene
The Least absolute shrinkage and selection operator
Receiver operating characteristic
Area under the curve
Kyoto Encyclopedia of Genes, Genomes
Single-sample gene set enrichment analysis
CDGSH iron sulfur domain 1
Glutamate-cysteine ligase catalytic subunit
Solute carrier family 7 member11
Estrogen receptor membrane complex 2
Antigen presenting cell
Natural killer cell
Activated dendritic cells
T follicular helper cell
Glycogen synthase kinase-3β
Nuclear factor erythroid 2-related factor 2
Acyl-CoA synthetase long-chain family member 4
Activator of transcription factor 2
Anti-cytotoxic T lymphocyte-associated antigen-4
Anti-programmed cell death protein 1/programmed cell death ligand 1
Lymphocyte activation Gene-3
T cell immunoglobulin and mucin-domain containing-3
T cell immunoglobulin and ITIM domain
Han L, Yan Y, Zhao L, Liu Y, Lv X, Zhang L, et al. LncRNA HOTTIP facilitates the stemness of breast cancer via regulation of miR-148a-3p/WNT1 pathway. J Cell Mol Med. 2020;24(11):6242–52. https://doi.org/10.1111/jcmm.15261.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin. 2020;70(1):7–30. https://doi.org/10.3322/caac.21590.
Zuo S, Yu J, Pan H, Lu L. Novel insights on targeting ferroptosis in cancer therapy. Biomark Res. 2020;8(1):50. https://doi.org/10.1186/s40364-020-00229-w.
Bates JP, Derakhshandeh R, Jones L, Webb TJ. Mechanisms of immune evasion in breast cancer. BMC Cancer. 2018;18(1):556. https://doi.org/10.1186/s12885-018-4441-3.
Liang JY, Wang DS, Lin HC, Chen XX, Yang H, Zheng Y, et al. A novel ferroptosis-related gene signature for overall survival prediction in patients with hepatocellular carcinoma. Int J Biol Sci. 2020;16(13):2430–41. https://doi.org/10.7150/ijbs.45050.
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. https://doi.org/10.1186/s40364-020-00230-3.
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(7):e2307. https://doi.org/10.1038/cddis.2016.208.
Ma S, Dielschneider RF, Henson ES, Xiao W, Choquette TR, Blankstein AR, et al. Ferroptosis and autophagy induced cell death occur independently after siramesine and lapatinib treatment in breast cancer cells. PLoS One. 2017;12(8):e0182921. https://doi.org/10.1371/journal.pone.0182921.
Lin YS, Shen YC, Wu CY, Tsai YY, Yang YH, Lin YY, et al. Danshen improves survival of patients with breast cancer and dihydroisotanshinone I induces ferroptosis and apoptosis of breast cancer cells. Front Pharmacol. 2019;10:1226. https://doi.org/10.3389/fphar.2019.01226.
Wu X, Liu C, Li Z, Gai C, Ding D, Chen W, et al. Regulation of GSK3beta/Nrf2 signaling pathway modulated erastin-induced ferroptosis in breast cancer. Mol Cell Biochem. 2020;473(1–2):217–28. https://doi.org/10.1007/s11010-020-03821-8.
Doll S, Proneth B, Tyurina YY, Panzilius E, Kobayashi S, Ingold I, et al. ACSL4 dictates ferroptosis sensitivity by shaping cellular lipid composition. Nat Chem Biol. 2017;13(1):91–8. https://doi.org/10.1038/nchembio.2239.
Wang L, Chen Y, Mi Y, Qiao J, Jin H, Li J, Lu Z, Wang Q, Zou Z. ATF2 inhibits ani-tumor effects of BET inhibitor in a negative feedback manner by attenuating ferroptosis. Biochem Biophys Res Commun. 2021;558:216-223.
Darvin P, Toor SM, Sasidharan Nair V, Elkord E. Immune checkpoint inhibitors: recent progress and potential biomarkers. Exp Mol Med. 2018;50(12):1–11. https://doi.org/10.1038/s12276-018-0191-1.
Qin S, Xu L, Yi M, Yu S, Wu K, Luo S. Novel immune checkpoint targets: moving beyond PD-1 and CTLA-4. Mol Cancer. 2019;18(1):155. https://doi.org/10.1186/s12943-019-1091-2.
Ribas A. Adaptive immune resistance: how cancer protects from immune attack. Cancer Discov. 2015;5(9):915–9. https://doi.org/10.1158/2159-8290.CD-15-0563.
Emens LA. Breast cancer immunotherapy: facts and hopes. Clin Cancer Res. 2018;24(3):511–20. https://doi.org/10.1158/1078-0432.CCR-16-3001.
Avgeris M, Panoutsopoulou K, Papadimitriou MA, Scorilas A. Circulating exosomal miRNAs: clinical significance in human cancers. Expert Rev Mol Diagn. 2019;19(11):979–95. https://doi.org/10.1080/14737159.2019.1673732.
Kerjaschki D, Bago-Horvath Z, Rudas M, Sexl V, Schneckenleithner C, Wolbank S, et al. Lipoxygenase mediates invasion of intrametastatic lymphatic vessels and propagates lymph node metastasis of human mammary carcinoma xenografts in mouse. J Clin Invest. 2011;121(5):2000–12. https://doi.org/10.1172/JCI44751.
Sohn YS, Tamir S, Song L, Michaeli D, Matouk I, Conlan AR, et al. NAF-1 and mitoNEET are central to human breast cancer proliferation by maintaining mitochondrial homeostasis and promoting tumor growth. Proc Natl Acad Sci U S A. 2013;110(36):14676–81. https://doi.org/10.1073/pnas.1313198110.
Salem AF, Whitaker-Menezes D, Howell A, Sotgia F, Lisanti MP. Mitochondrial biogenesis in epithelial cancer cells promotes breast cancer tumor growth and confers autophagy resistance. Cell Cycle. 2012;11(22):4174–80. https://doi.org/10.4161/cc.22376.
O'Neill S, Porter RK, McNamee N, Martinez VG, O'Driscoll L. 2-Deoxy-D-glucose inhibits aggressive triple-negative breast cancer cells by targeting glycolysis and the cancer stem cell phenotype. Sci Rep. 2019;9(1):3788. https://doi.org/10.1038/s41598-019-39789-9.
Song NY, Kim DH, Kim EH, Na HK, Kim NJ, Suh YG, et al. Multidrug resistance-associated protein 1 mediates 15-deoxy-Delta(12,14)-prostaglandin J2-induced expression of glutamate cysteine ligase expression via Nrf2 signaling in human breast cancer cells. Chem Res Toxicol. 2011;24(8):1231–41. https://doi.org/10.1021/tx200090n.
Seibt TM, Proneth B, Conrad M. Role of GPX4 in ferroptosis and its pharmacological implication. Free Radic Biol Med. 2019;133:144–52. https://doi.org/10.1016/j.freeradbiomed.2018.09.014.
Ge C, Cao B, Feng D, Zhou F, Zhang J, Yang N, et al. The down-regulation of SLC7A11 enhances ROS induced P-gp over-expression and drug resistance in MCF-7 breast cancer cells. Sci Rep. 2017;7(1):3791. https://doi.org/10.1038/s41598-017-03881-9.
Mele L, la Noce M, Paino F, Regad T, Wagner S, Liccardo D, et al. Glucose-6-phosphate dehydrogenase blockade potentiates tyrosine kinase inhibitor effect on breast cancer cells through autophagy perturbation. J Exp Clin Cancer Res. 2019;38(1):160.
Basu A, Ramamoorthi G, Jia Y, Faughn J, Wiener D, Awshah S, et al. Immunotherapy in breast cancer: current status and future directions. Adv Cancer Res. 2019;143:295–349. https://doi.org/10.1016/bs.acr.2019.03.006.
This work was supported by grants from National Natural Science Foundation of China (No. 81974049) and National key Research and development program (NO. 2017YFCO114210).
Ethics approval and consent to participate
Consent for publication
All authors state that there are no potential competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Wang, D., Wei, G., Ma, J. et al. Identification of the prognostic value of ferroptosis-related gene signature in breast cancer patients. BMC Cancer 21, 645 (2021). https://doi.org/10.1186/s12885-021-08341-2