Exploration of prognostic genes and risk signature in breast cancer patients based on RNA binding proteins associated with ferroptosis

Background: Breast cancer (BRCA) is a life-threatening malignancy in women with an unsatisfactory prognosis. The purpose of this study was to explore the prognostic biomarkers and a risk signature based on ferroptosis-related RNA-binding proteins (FR-RBPs). Methods: FR-RBPs were identified using Spearman correlation analysis. Differentially expressed genes (DEGs) were identified by the “limma” R package. The univariate Cox and multivariate Cox analyses were executed to determine the prognostic genes. The risk signature was constructed and verified with the training set, testing set, and validation set. Mutation analysis, immune checkpoint expression analysis in high- and low-risk groups, and correlation between risk signature and chemotherapeutic agents were conducted using the “maftools” package, “ggplot2” package, and the CellMiner database respectively. The Human Protein Atlas (HPA) database was employed to confirm protein expression trends of prognostic genes in BRCA and normal tissues. The expression of prognostic genes in cell lines was verified by Real-time quantitative polymerase chain reaction (RT-qPCR). Kaplan-meier (KM) plotter database analysis was applied to predict the correlation between the expression levels of signature genes and survival statuses. Results: Five prognostic genes (GSPT2, RNASE1, TIPARP, TSEN54, and SAMD4A) to construct an FR-RBPs-related risk signature were identified and the risk signature was validated by the International Cancer Genome Consortium (ICGC) cohort. Univariate and multivariate Cox regression analysis demonstrated the risk score was a robust independent prognostic factor in overall survival prediction. The Tumor Mutational Burden (TMB) analysis implied that the high- and low-risk groups responded differently to immunotherapy. Drug sensitivity analysis suggested that the risk signature may serve as a chemosensitivity predictor. The results of GSEA suggested that five prognostic genes might be related to DNA replication and the immune-related pathways. RT-qPCR results demonstrated that the expression trends of prognostic genes in cell lines were consistent with the results from public databases. KM plotter database analysis suggested that high expression levels of GSPT2, RNASE1, and SAMD4A contributed to poor prognoses. Conclusion: In conclusion, this study identified the FR-RBPs-related prognostic genes and developed an FR-RBPs-related risk signature for the prognosis of BRCA, which will be of great significance in developing new therapeutic targets and prognostic molecular biomarkers for BRCA.