Development and Validation of a Robust Ferroptosis-Related Gene Panel for Breast Cancer Disease-Specific Survival

Background: New biomarker combinations have been increasingly developed to improve the precision of current diagnostic and therapeutic modalities. Recently, researchers have found that tumor cells are more vulnerable to ferroptosis. Furthermore, ferroptosis-related genes (FRG) are promising therapeutic targets in breast cancer patients. Therefore, this study aimed to identify FRG that could predict disease-specific survival (DSS) in breast cancer patients. Methods: Gene expression matrix and clinical data were downloaded from public databases. We included 960, 1,900, and 234 patients from the TCGA, METABRIC, and GSE3494 cohorts, respectively. Data for FRG were downloaded from the FerrDb website. Differential expression of FRG was analyzed by comparing the tumors with adjacent normal tissues. Univariate Cox analysis of DSS was performed to identify prognostic FRG. The TCGA-BRCA cohort was used to generate a nine-gene panel with the LASSO cox regression. The METABRIC and GSE3494 cohorts were used to validate the panel. The panel's median cut-off value was used to divide the patients into high- or low-risk subgroups. Analyses of immune microenvironment, functional pathways, and clinical correlation were conducted via GO and KEGG analyses to determine the differences between the two subgroups. Results: The DSS of the low-risk subgroup was longer than that of the high-risk subgroup. The panel's predictive ability was confirmed by ROC curves (TCGA cohort AUC values were 0.806, 0.695, and 0.669 for 2, 3, and 5 years respectively, and the METABRIC cohort AUC values were 0.706, 0.734, and 0.7, respectively for the same periods). The panel was an independent DSS prognostic indicator in the Cox regression analyses. (TCGA cohort: HR = 3.51, 95% CI = 1.792-6.875, p < 0.001; METABRIC cohort: HR = 1.76, 95% CI = 1.283-2.413, p < 0.001). Immune-related pathways were enriched in the high-risk subgroup. The two subgroups that were stratified by the nine-gene panel were also associated with histology type, tumor grade, TNM stage, and Her2-positive and TNBC subtypes. The patients in the high-risk subgroup, whose CTLA4 and PD-1 statuses were both positive or negative, demonstrated a substantial clinical benefit from combination therapy with anti-CTLA4 and anti-PD-1. Conclusion: The new gene panel consisting of nine FRG may be used to assess the prognosis and immune status of patients with breast cancer. A precise therapeutic approach can also be possible with risk stratification.