key: cord-0031812-96nxtmie
authors: Xu, Yuhao; Du, Yaoqiang; Zheng, Qinghui; Zhou, Tao; Ye, Buyun; Wu, Yihao; Xu, Qiuran; Meng, Xuli
title: Identification of Ferroptosis-Related Prognostic Signature and Subtypes Related to the Immune Microenvironment for Breast Cancer Patients Receiving Neoadjuvant Chemotherapy
date: 2022-05-04
journal: Front Immunol
DOI: 10.3389/fimmu.2022.895110
sha: 38668eff76ebd3854f19991e0882d9d3b065a2b6
doc_id: 31812
cord_uid: 96nxtmie

PURPOSE: To identify molecular clusters associated with ferroptosis and to develop a ferroptosis-related signature for providing novel potential targets for the recurrence-free survival and treatment of breast cancer. METHODS: Ferroptosis-related gene (FRG) signature was constructed by univariate and multivariate Cox regression and least absolute shrinkage and selection operator (LASSO). Receiver operating characteristic curves, Kaplan–Meier survival analysis, principal component analysis, and univariate and multivariate Cox regression analyses in the training and test cohorts were used to evaluate the application of this signature. Quantitative reverse transcriptase–PCR (qRT-PCR) was employed to detect the expression of FRGs in the model. Furthermore, the correlations between the signature and immune microenvironment, somatic mutation, and chemotherapeutic drugs sensitivity were explored. RESULTS: Internal and external validations affirmed that relapse-free survival differed significantly between the high-risk and low-risk groups. Univariate and multivariate Cox regression analyses indicated that the riskScore was an independent prognostic factor for BRCA. The areas under the curve (AUCs) for predicting 1-, 2-, and 3-year survival in the training and test cohorts were satisfactory. Significant differences were also found in the immune microenvironment and IC50 of chemotherapeutic drugs between different risk groups. Furthermore, we divided patients into three clusters based on 18 FRGs to ameliorate the situation of immunotherapy failure in BRCA. CONCLUSIONS: The FRG signature functions as a robust prognostic predictor of the immune microenvironment and therapeutic response, with great potential to guide individualized treatment strategies in the future.

Breast cancer has surpassed lung cancer as being the most commonly diagnosed cancer with approximately 2.3 million new cases in 2020, accounting for 11.7% of all new cancer cases (1) . Another scary truth is the drop in average onset age (2) . Because breast cancer is a highly heterogeneous systemic disease, advancements in therapy are particularly crucial (3) .

Neoadjuvant chemotherapy (NAC) is seen as the standard and first-line treatment for locally advanced breast cancer (4, 5) , which not only is beneficial to breast-conserving surgery but also can detect tumor sensitivity to anticancer therapy for locally advanced breast cancer (6) , and it could also be employed as a bridge to other therapies (7, 8) . Anthracyclines and taxanes serve as the backbone of NAC regimens and are widely used clinically (9) .

Ferroptosis is an emerging form of programmed cell death featured by the iron-dependent accumulation of lipid reactive oxygen species (ROS) of metabolic dysfunctions, iron accumulation, and antioxidant vulnerability (10) (11) (12) . Accumulating evidence showed that the role of ferroptosis in carcinogenesis, progression, and chemoresistance had made progress. Fascin regulates SLC7A11 stability to induce ferroptosis (13) . Renovation of SLC7A11 rescues miR-5096mediated ferroptosis and antitumor effects of breast cancer (14) . Ferroptosis-related gene (FRG) GPX4 promotes chemoresistance in nasopharyngeal carcinoma (15) . Bufotalin induces ferroptosis by facilitating the ubiquitination and degradation of GPX4 in non-small cell lung cancer cells (16) . Via ferroptosis, ETS1/miR-23a-3p/ACSL4 axis stimulates sorafenib resistance in HCC (17) .

Ferroptosis has the characteristics of inhibiting chemoresistance and enhancing antitumor immunity (18) , which may be a potential strategy to overcome the drug resistance mechanism of traditional cancer treatments (12) . Previous studies prove the feasibility of ferroptosis-related prognostic markers to predict overall survival and immune characteristics. FRG signatures were constructed to predict overall survival in lung adenocarcinoma (19) , colorectal cancer (20) , and pancreatic adenocarcinoma (21) . However, as far as we are aware, studies focusing on the correlation of ferroptosis with biochemical recurrence and antitumor immunology of BRCA were rather limited. Thus, it is an urgent need to discover a robust biomarker to predict relapse-free survival (RFS) in BRCA.

In this study, we constructed an FRG prognostic signature and identified three ferrClusters in predicting the RFS internally and externally, exploring the status of immune infiltrates and drug sensitivity of BRCA patients receiving NAC for guiding clinical practice. This signature may also serve as a novel and robust prediction tool for evaluating whether BRCA patients can benefit from immunotherapy.

Open expression matrix of mRNA (FPKM values) and clinical files of BRCA samples were downloaded from The Cancer Genome Atlas (TCGA) database. Datasets GSE25055 in the Gene Expression Omnibus (GEO) database were used to acquire RNA-sequencing (RNA-Seq) and clinical data of BRCA patients receiving NAC as a training cohort and GSE16446 and GSE25065 as test cohorts. Gene expression file of GSE25055 and GSE25065 was collected using platform GPL96 [HG-U133A] Affymetrix Human Genome U133A Array, and GSE16446 using platform GPL570 [HG-U133_Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array. Batch effects and other unwanted variations in high-throughput experiments were eliminated using the "combat" function in the "sva" package (22) in R 4.1.1. Copy number variation (CNV) data were collected from the University of California, Santa Cruz (UCSC) website.

FRGs including 150 drivers, 109 suppressors, and 123 markers were collected from FerrDb (19, 23, 24) . A univariate Cox proportional hazards regression analysis was conducted to filtrate prognostic FRGs in the GSE25055 cohort with p < 0.05 considered to be statistically significant using the "coxph" function. Subsequently, with the help of the "cv.glmnet" function, the least absolute shrinkage and selection operator (LASSO) was performed for the dimension reduction and Kfold cross-validation, which was multiplied by ten, and the optimal parameter was the l value that corresponded to the lowest deviation. The optimal penalty parameter was defined as the value within one SD of the minimum cross-validated partial likelihood deviance to obtain the best model. The proteins with non-zero regression coefficients were chosen for subsequent 

First, in the GSE25055 dataset, the Kaplan-Meier (K-M) survival analysis using the "Surv" function in the "survival" package and univariate and multivariate Cox regression analyses between gene expression and clinical characters using the "coxph" function in the "survival" package were performed to certify that riskScore served as an independent predictor in predicting recurrence-free survival (RFS). Principal component analysis (PCA) using the "prcomp" function was used to visualize sample distribution. Receiver operating characteristic (ROC) using the "timeROC" package was done, and area under the curve (AUC) plots were generated for the 1-year, 2-year, and 3year survival rates to assess the sensitivity and specificity of the prognostic model. Then, the prognostic signature was validated in the GSE25065 and GSE16446 datasets via the above methods. Quantitative Reverse Transcriptase-PCR Total RNAs were isolated from cells using the TRIzol reagent (Invitrogen, Carlsbad, CA, USA). PrimeScript ™ RT reagent Kit (Takara, Maebashi, Japan) was employed to reverse transcribe into cDNA following the manufacturer's protocol. Then SYBR Green PCRMaster Mix (Applied TaKaRa, Otsu, Japan) was used to conduct Real-time PCR on Applied Biosystems 7500 Fast Real-Time RCR System (Applied Biosystems, Foster City, CA, USA). The primers of FRGs for qRT-PCR utilized in this research were as follows:

Primer sequence (5′ to 3′)

With the use of the "pRRophetic" package, the halfmaximal inhibitory concentration (IC50) of BRCA patients was calculated on Genomics of Drug Sensitivity in Cancer (GDSC) (25) (https://www.cancerrxgene.org/) based on the given gene expression profiles in these datasets to evaluate the drug sensitivities (26) (27) (28) .

The CIBERSORT algorithm was used to explore the proportion of different types of immune cells in BRCA patients using CIBERSORT R script v1.04 (29) (30) (31) . Based on the expression level of immune cell-related genes, the ESTIMATE algorithm was conducted to calculate the stromal score (SS), estimate score (ES), and immune score (IS) (the SS represents the level of stroma content in a tumor; the IS reflects the infiltration of immune cells in a tumor; the estimated score infers tumor purity) among the high-and low-risk groups using the "estimate" package (32). Single-sample gene set enrichment analysis (ssGSEA) was performed to calculate scores for antitumor immunity and protumor suppression for each sample (33) using "GSEABase" and "GSVA" packages.

Consensus clustering based on Euclidean distance and Ward's linkage was performed for hierarchical clustering to identify different subtypes using the "ConsensusClusterPlus" package and repeated the procedures 1,000 times to guarantee the stability of the classification (34) . In consideration of a high consistency of clusters, a low coefficient of variation, and no significant increase in the CDF curve, the optimum cluster number could be determined (35) .

Correlation coefficients were calculated by Spearman's and distance correlation analyses. For comparison of more than two groups, the Kruskal-Wallis and one-way ANOVAs were chosen as nonparametric and parametric methods, while Wilcoxon's t-test was used for two groups. Student's t-test was used to explore the statistical significance of quantitative data. The K-M and log-rank tests were employed to confirm the significance of prognostic differences (22) . R 4.1.1 software was the main tool to conduct the statistical analysis. For all statistical results, a p-value of <0.05 was considered to be statistically significant.

GSE25055 dataset was used as a training cohort; meanwhile, GSE25065 and GSE16446 datasets were used as test cohorts. Batch effects were removed for further study ( Figures 1A, B) . First, we performed a univariate Cox regression analysis in GSE25055. Among 382 FRGs retrieved from the FerrDB database, 76 FRGs were identified to be associated with RFS, with the standard of p < 0.01 ( Figure 1C ). Pearson's correlation analysis revealed a correlation among these genes ( Figure 1D ). Then, LASSO regression analysis was used to establish the FRGs prognostic signature ( Figures 1E, F) :

Genomic mutations were common in these genes with 40 (4.07%) of 983 patients having experienced genetic changes, and a mutation frequency of 1% was observed in ENPP2 and NEDD4L ( Figure 1G ). We also found that CNV is prevalent among the 18 FRGs. ENPP2, CHMP6, MAFG, VEGFA, VDAC2, LPCAT3, CYP4F8, SLC1A4, XBP1, PIR, and ACO1 showed copy number amplification, while deletion happened in the other FRGs ( Figure 1H ). The location of the 18 FRGs in human chromosomes could be seen in Figure 1I . The result of differential analysis in normal breast tissue and tumor tissue showed that ACO1, CHMP6, ENPP2, MAFG, NEDD4L, PIR, SLC1A4, SLC7A5, and VEGFA had significant differential expression in breast cancer with p-value <0.001; ACADSB and NOX3 with p-value <0.01; and PEX12 with p-value <0.05 ( Figure 1J ). SLC7A5 was seen as a significantly upregulated gene, while ENPP2 and ACO1 were seen as significantly downregulated genes with |log FC| > 1 ( Figure 1K ). The result of RT-PCR provided strong support for our conclusion ( Figures 1L-N) . As described above, FRGs had significant heterogeneity of genetic variation and transcriptomic alteration landscape in BRCA patients, which played an important part in regulating the happening, aggravation, and prognosis of BRCA.

After the riskScore of each patient based on the risk model was calculated, we divided patients into the high-and low-risk groups with the standard median score in GSE25055 (Figure 2A) . With the use of the median of GSE25055, patients in GSE25065 ( Figure 2H ) and GSE16446 (Figures 2H, O) were separated into the high-risk and low-risk groups in the same manner. The result of PCA showed significant heterogeneity between high-risk and low-risk patients in GSE25055 ( Figures 2B, C) , GSE25065 ( Figures 2I, J) , and GSE16446 ( Figures 2P, Q) , which certified the superior discrimination of the FRG model. For the purpose of exploring whether the signature could represent its prognostic value independently of other clinical factors, we conducted univariate and multivariate Cox regression analyses in the training and test cohorts. In univariate analyses, this risk score was able to independently predict survival outcomes in GEO cohorts (GSE25055, hazard ratio (HR) = 4.690, p < 0.001; GSE25065, HR = 6.350, p < 0.001; GSE16446, HR = 7.648, p < 0.001) ( Figures 2D,  K, R) . The same conclusion could be drawn in multivariate analyses ( Figures 2E, L, S) . The results revealed that riskScore and pathologic response served as independent factors affecting receiving NAC BRCA patients' prognosis. The AUCs of the time-dependent ROC curves at 1, 2, and 3 years were 0.818, 0.824, and 0.783 in GSE25055 ( Figure 2F) ; 0.812, 0.824, and 0.783 in GSE25065 ( Figure 2M) ; and 0.715, 0.725, and 0.723 in GSE16446 ( Figure 2T) . The AUCs in different years and cohorts were relatively high compared with those of other published literature, which suggested high sensitivity and specificity of the signature for predicting RFS. The K-M survival curve showed that patients in the high-risk group had a higher recurrence rate than those in the low-risk group using log-rank tests with p < 0.001 ( Figure 2G) , p < 0.001( Figure 2N ), p = 0.029 ( Figure 2U ). Ferroptosis is a recently recognized form of regulated cell death that is characterized by lipid peroxidation, which mediates cell death in breast cancer. Among genes in our signature, ferroptosis driver genes such as NOX3 and PEBP1 had negative coefficients, while ferroptosis suppressor genes such as PIR and VDAC2 had positive coefficients. Therefore, high riskScore indicated that ferroptosis was suppressed in breast cancer, which might imply a worse prognosis.

We further anatomized the association between riskScore and clinical parameters of BRCA patients. The detailed results depicted that the riskScore had a positive correlation with T stage, N stage, American Joint Committee on Cancer (AJCC) stage, and grade ( Figure 3 ). BRCA patients with higher T, N, AJCC stage, and grade, combined with lower age, and negative status of progesterone receptor and estrogen receptor seemed to have higher riskScore, indicating a higher incidence rate of relapse, which was consistent with the conclusions of current accumulated literature. In other words, the results implied that the riskScore had a correlation with clinicopathological parameters.

In order to improve the therapeutic benefit of BRCA patients from neoadjuvant therapy, we further explored whether FRG signature could predict the sensitivity to several chemotherapy drugs widely used in BRCA between two groups. According to the results calculated based on the GDSC database, IC50 values of chemotherapy drugs covering axitinib, bicalutamide, bleomycin, bortezomib, dasatinib, doxorubicin, gefitinib, lapatinib, and paclitaxel were evaluated. Compared with the low-risk group, IC50 values of paclitaxel, gefitinib, doxorubicin, bleomycin, and bortezomib were lower in the high-risk groups, which indicated that high-risk patients were more sensitive to these drugs ( Figures 4A-I) . The above results demonstrated that the riskScore had potential predictive value for chemotherapy and targeted therapy in breast cancer.

We calculated the constitution of tumor-infiltrating immune cells in BRCA through the CIBERSORT algorithm ( Figure 5A) . Compared with the low-risk groups, the proportion of resting mast cells was lower in the high-risk groups ( Figures 5B, C) . Then, the IS, SS, and ES of patients were evaluated using the ESTIMATE algorithm. Based on the optimum cutoff value of ISs or SSs respectively, BRCA patients were divided into the high and low IS/SS/ES groups. The K-M curves showed that patients with high IS/SS/ES exhibited significantly worse RFS as compared to the ones with low IS/SS/ES (Figures 5D-F) . We further explored the relationships between the IS/SS/ES and riskScore. The result of Wilcoxon's rank-sum test displayed that there is no significant difference between the high-risk and low-risk groups in SS (p = 0.53, Figure 5G ) but significant in IS (p = 5.7e−09, Figure 5H ) and ES (p = 1.9e−05, Figure 5I ). Pearson's correlation analysis showed that riskScore was positively associated with IS (R = 0.26, p = 1.6e−10, Figure 5K ) and ES (R = 0.2, p = 5.1e−07, Figure 5L ). However, the riskScore was not significantly correlated to the SS (P = 0.14, Figure 5J ) Furthermore, we dissected the role of riskScore in immune checkpoint blockade (ICB) treatment. We noticed that the expression levels of all immune checkpoints were significantly higher in the high riskScore group (Figures 5M, N) . Taken together, the prognostic signature could predict the potential response to immunotherapy in BRCA patients, which provided guidance on whether or what to use for immunotherapy in clinical practice.

On the basis of the expression of 18 FRGs in the signature, we employed the "Partition Around Medoids" algorithm, along with Pearson's distance to estimate similarity among patients to identify three clusters. We noticed that K = 3 seemed to be an optimal selection by clustering variable (k) increasing from 2 to 9, in which the intergroup correlations were the lowest and the intragroup correlations were the greatest ( Figure 6B) , indicating the optimal clustering stability of the three molecular phenotypes. The consensus cumulative distribution function (CDF) diagram showed that when k = 3, distribution reached an approximate maximum ( Figure 6C) , implying robust clustering for all samples ( Figure 6A ). The delta area plot depicts the relative change compared to k − 1 showing that the delta area was optimum when k = 3 ( Figure 6D ). Prognostic analysis of the three clusters revealed that patients in ferrCluster A were the least likely to relapse, while in ferrCluster B, they were the most likely to relapse ( Figure 6E) .

We then performed an ssGSEA to quantify the scores of various immune cell subpopulations to further compare the differences in the number of immune cells among the three types of ferrClusters. The results indicated that the contents of monocyte cells were not significantly different. The proportion of immune cells was significantly different among the three clusters. Contents of nearly all types of immune cells in ferrCluster A seemed to be the poorest. The levels of activated CD4+ T cells, CD8+ T cells, dendritic cells, CD56 bright and dim NK cells, gd-T cells, Tregs, and T helper cells were relatively the highest in ferrCluster B. Hence, we could draw the conclusion that ferrCluster A was a type of immune failure, ferrCluster B was a type of immune-activated characterized by T-cell subset enrichment, and ferrCluster C was a type of immune-activated characterized by B-cell subset enrichment ( Figure 6F) . These results indicated that the FRGs play key roles in immune cell infiltration and characteristic tumor immune microenvironment (TME) formation and affect the prognosis of BRCA patients.

With a view to the individual heterogeneity and complexity of BRCA patients, we calculated ferrScore to assess the ferroptosis pattern of each patient based on the PCA on the 18 FRGs in the model. The scoring framework was defined as ferrScore = PC1 + PC2 to quantify individual ferroptosis patterns of BRCA patients (36) , further facilitating precise treatment. As indicated from the K-M curve, patients with lower ferrScore had a lower probability of relapse ( Figure 7A ). The ferrScore was closely related to immune cells ( Figure 7B ). We also observed that ferrScores of patients in ferrCluster A were significantly lower than those in ferrCluster B and C, while there was no significant difference between ferrClusters B and C ( Figure 7D ). The Sankey diagram shows the attribute changes in riskScore, ferrCluster, ferrScore, and recurrence status, indicating that the higher the riskScore and ferrScore, the higher the risk of relapse after receiving NAC ( Figure 7C ). The above results enriched treatment strategies for BRCA patients not only in targeted therapy and chemotherapy but also in immunotherapy. At last, the expression of CTLA4 was examined to elucidate a potential response to immunotherapy, and the high ferrScore group showed relatively high levels of expression (p = 8.6e−11, Figure 7B ).

Ferroptosis is a newfound programmed cell death pattern distinguished from traditional cell death such as apoptosis, necrosis, and autophagy (37) . Accumulating evidence demonstrated that dysregulated expression and genetic variations of FRGs were closely related to cell death, tumor carcinogenesis, and progression (22, 38) . TME is a cradle for tumorigenesis and cancer progression, in which immune infiltrating cells affect therapeutic outcomes (39) . The relationships between TME infiltration immune cells and ferroptosis modifications have become a hotspot in the mechanism of tumorigenesis and development (40, 41) . MIF secreted by nasopharyngeal carcinoma could suppress ferroptosis of macrophages and then increase the rate of metastasis (42) . BEBT-908 induces immunogenic ferroptosis to potentiate cancer immune checkpoint therapy (43) . SCD1 and FABP4 could drive ferroptosis, thereby leading to tumor resistance (44) . Ferroptotic cells could also release chemotaxis to interact with immune cells, such as CD + T cells, and then modulate the anticancer immunity (45) .

High-throughput genomic studies provided cutting-edge sights into the molecular mechanisms and identified new potential targets of breast cancer. Our research developed and verified a stepwise multivariate Cox regression model including 18 FRGs using LASSO and multivariate Cox regression for removing redundant factors to forecast the RFS of individual patients in GSE25055. The expression of FRGs in the signature was higher in BRCA tissues than in adjacent normal tissues, which was verified in several breast cancer cell lines using realtime PCR. Meanwhile, CNVs and mutation frequencies of FRGs were prevalent. Internal and external validations exhibited an excellent ability to predict the prognosis of BRCA patients. Specifically, a higher riskScore indicated a higher rate of recurrence. Moreover, riskScore was associated closely with clinicopathological features. With a view to the significance of the immune system in antiviral and antitumor responses, we calculated the proportion of different types of tumor-infiltrating immune cells in BRCA using CIBERSORT and used ESTIMATE to explore IS, SS, and tumor purity. Higher SSs and ISs were observed in high-risk patients, leading to an unfavorable prognosis, which was consistent with a line of evidence from previous research (46, 47) .

Extensive interest in cancer immunotherapy is reported according to the clinical importance of CTLA-4 and PD-1/PD-L1 in immune checkpoint therapies (48) . The main immune checkpoints for breast cancer include CTLA-4, PD-1/PD-L1, lymphocyte activation gene 3 (LAG-3), T-cell immunoglobulin domain and mucin 3 (TIM-3), and other molecules (49) . Clinical trials like SOLTI-1503 PROMETEO TRIAL (50), KEYNOTE-086 (51), NIMBUS (52), KEYNOTE-173 (53) , and KEYNOTE-522 (54) showed that immunological checkpoint inhibitors have made significant progress in breast cancer immunotherapy, which is expected to become a new treatment for breast cancer.

Furthermore, for the purpose of exploring the response to chemotherapy sensitivity of patients, we calculated the IC50 value. The sensitivities of chemotherapeutic drugs widely used in BRCA showed a significant difference between the two groups.

In accordance with the expression matrix of the 18 FRGs in the signature, we identified three ferroptosis-related molecular clusters via consensus clustering analysis. The rate of relapse was significantly different among the three clusters. ssGSEA identified that the three ferrClusters as three immune types of immune failure, immune-activated characterized by T-cell subset enrichment, and immune-activated characterized by B-cell subset enrichment.

Inevitably, numerous limitations of our study should be included in the consideration. First, although our conclusion came through internal and external validation in TCGA, GSE25055, GSE25065, and GSE16446 cohorts, when it comes to its clinical application, caution is advised. Multicenter large-scale prospective clinical studies were needed rather than only retrospective data from public open databases to verify the signature. Second, the expression matrix of patients in GSE25055 and GSE25065 was extracted via platform GPL96 [HG-U133A] Affymetrix Human Genome U133A Array in 2010, which only included 12,549 genes, while GPL570 [HG-U133_Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array for GSE16446 contained 21,655 genes. Due to the relatively small number of detectable genes, bias may be amplified. Finally, detailed molecular mechanisms in the BRCA of the FRGs in the signature had not been fully revealed. Further in-depth studies were required to confirm relationships between FRGs and tumor microenvironment, and between ferroptosis and chemoresistance.

In brief, we constructed a novel FRG signature and identified three molecular subtypes for predicting the RFS of BRCA patients, which could predict the immune status of the tumor microenvironment and RFS of patients. It is worth noting that our conclusions provided more clues for the rational choices of chemotherapeutic drugs for patients with BRCA, provided a new immunological perspective and a new basis for immunotherapy of BRCA in the clinic, and had the potential possibility to coach and guide individualized healthcare decisions.

The expression and clinical datasets presented in this study are available in databases of TCGA and Gene Expression Omnibus (GSE25055, GSE25065 and GSE16446), and the copy number variation data are presented in the University of California, Santa Cruz (UCSC) website https://xenabrowser.net/datapages/? dataset=TCGA-BRCA.gistic.tsv&host=https%3A%2F%2Fgdc. xenahubs.net&removeHub=https%3A%2F%2Fxena.treehouse. gi.ucsc.edu%3A443.

XM and QX substantially contributed to the conception of the work. YX, QZ, and TZ contributed to the data collection. YX, YD, and QX wrote the manuscript. QZ and YD helped to perform the enrichment and network analysis. YX, YD, and XM drafted and revised the manuscript. All authors contributed to the article and approved the submitted version. 

Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries

Bilateral Intra-Areolar Polythelia

Current Landscape of Immunotherapy in Breast Cancer: A Review

Early Breast Cancer Trialists' Collaborative G. Long-Term Outcomes for Neoadjuvant Versus Adjuvant Chemotherapy in Early Breast Cancer: Meta-Analysis of Individual Patient Data From Ten Randomised Trials

Pathological Complete Response and Long-Term Clinical Benefit in Breast Cancer: The CTNeoBC Pooled Analysis

Stromal Lymphocyte Infiltration After Neoadjuvant Chemotherapy Is Associated With Aggressive Residual Disease and Lower Disease-Free Survival in HER2-Positive Breast Cancer

Impact of 70-Gene Signature Use on Adjuvant Chemotherapy Decisions in Patients With Estrogen Receptor-Positive Early Breast Cancer: Results of a Prospective Cohort Study

Breast Cancer

Early Stage Triple Negative Breast Cancer: Management and Future Directions

Process and Function

Ferroptosis: A Regulated Cell Death Nexus Linking Metabolism

Ferroptosis: A Promising Target for Cancer Immunotherapy

Fascin Enhances the Vulnerability of Breast Cancer to Erastin-Induced Ferroptosis

SLC7A11/ xCT is a Target of miR-5096 and its Restoration Partially Rescues miR-5096-Mediated Ferroptosis and Anti-Tumor Effects in Human Breast Cancer Cells

EBV Infection-Induced GPX4 Promotes Chemoresistance and Tumor Progression in Nasopharyngeal Carcinoma

Bufotalin Induces Ferroptosis in non-Small Cell Lung Cancer Cells by Facilitating the Ubiquitination and Degradation of GPX4

Epigenetic Regulation of Ferroptosis via ETS1/miR-23a-3p/ACSL4 Axis Mediates Sorafenib Resistance in Human Hepatocellular Carcinoma

The Evolving Role of Ferroptosis in Breast Cancer: Translational Implications Present and Future

Establishment and Validation of a Ferroptosis-Related Gene Signature to Predict Overall Survival in Lung Adenocarcinoma

Prognostic and Therapeutic Significance of a Novel Ferroptosis Related Signature in Colorectal Cancer Patients

Development and Validation of a Novel 3-Gene Prognostic Model for Pancreatic Adenocarcinoma Based on Ferroptosis-Related Genes

Ferroptosis Regulator Modification Patterns and Tumor Microenvironment Immune Infiltration Characterization in Hepatocellular Carcinoma

A Novel Ferroptosis-Based Molecular Signature Associated With Biochemical Recurrence-Free Survival and Tumor Immune Microenvironment of Prostate Cancer

A Novel Prognostic Signature Based on Ferroptosis-Related Genes Predicts the Prognosis of Patients With Advanced Bladder Urothelial Carcinoma

A Pyroptosis-Based Prognostic Model for Immune Microenvironment Estimation of Hepatocellular Carcinoma

Using Immune-Related lncRNAs to Construct Novel Biomarkers and Investigate the Immune Landscape of Breast Cancer

Necroptosis-Related lncRNAs: Predicting Prognosis and the Distinction Between the Cold and Hot Tumors in Gastric Cancer

Identification of a Four-Gene-Based SERM Signature for Prognostic and Drug Sensitivity Prediction in Gastric Cancer

Comprehensive Analysis of Immune Cell Infiltration and Significant Genes in Head and Neck Squamous Cell Carcinoma

Identification of Driver Genes and Interaction Networks Related to Brain Metastasis in Breast Cancer Patients

Profiling Tumor Infiltrating Immune Cells With CIBERSORT

Inferring Tumour Purity and Stromal and Immune Cell Admixture From Expression Data

Identification and Verification of Tumor Immune Microenvironment-Related Prognostic Genes in Kidney Renal Clear Cell Carcinoma

ConsensusClusterPlus: A Class Discovery Tool With Confidence Assessments and Item Tracking

Construction of a Risk Prediction Model Using M6a RNA Methylation Regulators in Prostate Cancer: Comprehensive Bioinformatic Analysis and Histological Validation

Pyroptosis Regulators and Tumor Microenvironment Infiltration Characterization in Clear Cell Renal Cell Carcinoma

Ferroptosis: An Iron-Dependent Form of Nonapoptotic Cell Death

Activation of the P62-Keap1-NRF2 Pathway Protects Against Ferroptosis in Hepatocellular Carcinoma Cells

Vanden Berghe T. Targeting Ferroptosis to Iron Out Cancer

Ferroptosis in the Tumor Microenvironment: Perspectives for Immunotherapy

Ferroptosis in Carcinoma: Regulatory Mechanisms and New Method for Cancer Therapy

Exosomal MIF Derived From Nasopharyngeal Carcinoma Promotes Metastasis by Repressing Ferroptosis of Macrophages

A Dual PI3K/HDAC Inhibitor Induces Immunogenic Ferroptosis to Potentiate Cancer Immune Checkpoint Therapy

Tumor Resistance to Ferroptosis Driven by Stearoyl-CoA Desaturase-1 (SCD1) in Cancer Cells and Fatty Acid Biding Protein-4 (FABP4) in Tumor Microenvironment Promote Tumor Recurrence

CD8(+) T Cells Regulate Tumour Ferroptosis During Cancer Immunotherapy

Immunological Hallmarks of Stromal Cells in the Tumour Microenvironment

A Novel Prognostic Signature for Survival Prediction and Immune Implication Based on SARS-CoV-2-Related Genes in Kidney Renal Clear Cell Carcinoma

Immune Evasion Mechanism and AXL

Research Progresses in Immunological Checkpoint Inhibitors for Breast Cancer Immunotherapy

SOLTI-1503 PROMETEO TRIAL: Combination of Talimogene Laherparepvec With Atezolizumab in Early Breast Cancer

Pembrolizumab Monotherapy for Previously Untreated, PD-L1-Positive, Metastatic Triple-Negative Breast Cancer: Cohort B of the Phase II KEYNOTE-086 Study

Effect of Costimulatory Blockade With Abatacept After Ustekinumab Withdrawal in Patients With Moderate to Severe Plaque Psoriasis: The PAUSE Randomized Clinical Trial

Pembrolizumab Plus Chemotherapy as Neoadjuvant Treatment of High-Risk, Early-Stage Triple-Negative Breast Cancer: Results From the Phase 1b Open-Label, Multicohort KEYNOTE-173 Study

Pembrolizumab for Early Triple-Negative Breast Cancer

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.Copyright © 2022 Xu, Du, Zheng, Zhou, Ye, Wu, Xu and Meng. This is an openaccess 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.