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

Objective: To identify ferroptosis-related molecular clusters, and to develop and validate a ferroptosis-based molecular signature for predicting biochemical recurrence-free survival (BCRFS) and tumor immune microenvironment of prostate cancer (PCa). Materials and Methods: The clinical data and transcriptome data of PCa were downloaded from TCGA and GEO database. Ferroptosis-related genes (FRGs) were obtained from FerrDb database. We performed consensus clustering analysis to identify ferroptosis-related molecular subtypes for PCa. Univariate and multivariate Cox regression analysis were used to establish a ferroptosis-based signature for predicting BCRFS. Internal verification, external verification and subgroup survival analysis were then successfully performed. Results: There was a total of 40 differentially expressed FRGs in PCa. We then identified three ferroptosis-related molecular clusters of PCa, which have significantly different immune infiltrating cells, tumor immune microenvironment and PD-L1 expression level. More importantly, a novel ferroptosis-based signature for predicting BCRFS of PCa based on four FRGs (including ASNS, GPT2, NFE2L2, RRM2) was developed. Internal and external verifications were then successfully performed. Patients with high-risk score were associated with significant poor BCRFS compared with those with low-risk score in training cohort, testing cohort and validating cohort, respectively. The area under time-dependent Receiver Operating Characteristic (ROC) curve were 0.755, 0.705 and 0.726 in training cohort, testing cohort and validating cohort, respectively, indicating the great performance of this signature. Independent prognostic analysis indicated that this signature was an independent predictor for BCRFS of PCa. Subgroup analysis revealed that this signature was particularly suitable for younger or stage T III-IV or stage N0 or cluster 1-2 PCa patients. Patients with high-risk score have significantly different tumor immune microenvironment in comparison with those with low-risk score. The results of qRT-PCR successfully verified the mRNA expression levels of ASNS, GPT2, RRM2 and NFE2L2 in DU-145 and RWPE-1 cells while the results of IHC staining exactly verified the relative protein expression levels of ASNS, GPT2, RRM2 and NFE2L2 between PCa and BPH tissues. Conclusions: This study successfully identified three ferroptosis-related molecular clusters. Besides, we developed and validated a novel ferroptosis-based molecular signature, which performed well in predicting BCRFS and tumor immune microenvironment of PCa.