A prognostic signature consisting of N6-methyladenosine modified mRNAs demonstrates clinical potential in prediction of biochemical recurrence and guidance on precision therapy in prostate cancer

•We presented the first transcriptome-wide profile of m6A modifications in prostate cancer.•What's more, we developed the MMM score, which utilizes m6A-modified differentially expressed genes to predict the biochemical recurrence of prostate cancer patients. This score was verified in nine external multi-center cohorts.•Additionally, our findings indicate that patients with a high MMM score have a poorer response to ARSI therapy and immunotherapy.•Two hub genes in the MMM score, TLE1 and PFKL, were found to be directly modified by m6A and promote invasion of prostate cancer cells. Novel biomarkers are urgently needed to improve the prediction of clinical outcomes and guide personalized treatment for prostate cancer (PCa) patients. However, the role of N6-methyladenosine (m6A) modifications in PCa initiation and progression remains largely elusive. In our study, we collected benign Prostate Hyperplasia (BPH), localized PCa, and metastatic PCa samples from patients and performed methylated RNA immunoprecipitation sequencing (MeRIP-Seq) to map m6A-methylated mRNAs. Furthermore, we developed a prognostic signature based on 239 differentially methylated RNAs and the TCGA-PRAD dataset, which can be used to calculate an m6A-modified mRNA (MMM) score for a PCa patient, validated by independent multi-center cohorts. Our findings revealed that differential m6A modifications were positively correlated with altered expressions of mapped m6A-modified mRNAs. Higher MMM scores were associated with shorter times to biochemical recurrence (BCR) in PCa patients, and the MMM scoring system outperformed three well-established signatures in nine independent validation cohorts, as demonstrated by Kaplan-Meier survival analysis, C-index and ROC. Patients who did not respond to androgen receptor signaling inhibitor (ARSI) therapy and immunotherapy were found to have high MMM scores. Two hub genes, TLE1 and PFKL, were confirmed to have m6A sites through MeRIP-qPCR, and their knockdown promoted PCa cell invasion. Bioinformatics analysis of single-cell databases identified cell types with high transcript abundance levels of these two genes. In summary, our study is the first to perform transcriptome-wide m6A mapping in prostate tissues. The translational potential of a prognostic signature, comprising m6A-methylated mRNAs, in predicting clinical outcomes and therapy responses for PCa patients, is demonstrated. [Display omitted]