Roles of m5C RNA Modification Patterns in Biochemical Recurrence and Tumor Microenvironment Characterization of Prostate Adenocarcinoma.

ABSTRACT

Background: Prostate cancer is the second most common cancer with a high risk of biochemical recurrence (BCR) among men. Recently, 5-methylcytosine (m5C) modification has attracted more attention as a new layer of RNA post-transcriptional regulation. Hence, we aimed at investigating the potential roles of m5C modification regulators in the BCR of prostate adenocarcinoma (PRAD). Methods: CNV data, mutation annotation data, mRNA expression profiles, and clinical data were downloaded from TCGA and GEO databases. Kaplan-Meier curves analysis, log-rank test, univariate and multivariate Cox regression, and time-dependent ROC curves analysis were performed to evaluate the prognostic factors. Principal components analysis (PCA) was applied to validate the distinction between subgroups. Gene set variation analysis (GSVA) was used to investigate the underlying pathways associated with m5C modification patterns. Single sample gene set enrichment analysis (ssGSEA) was utilized to assess the infiltration of distinct immune cells. Tumor Immune Dysfunction and Exclusion (TIDE) prediction was carried out to assess the potential response to immune checkpoint blockade (ICB) therapy. The m5C modification signature was constructed via LASSO Cox's proportional hazards regression method. Results: After comprehensively analyzing various types of data from TCGA dataset, and exploring the differential expression and prognostic value of each m5C regulator, we identified m5C modification patterns based on 17 m5C regulators. Two patterns presented a significant difference in the risk of BCR, the tumor microenvironment (TME), and immunotherapy response in PRAD. We found that TET2, which was highly expressed in adjacent normal tissues compared to tumor tissues, was closely associated with many infiltrating immune cells. The m5C modification signature was constructed for the clinical application. Risk score calculated by m5C signature was associated with T stage, N stage, Gleason score, and the possibility of BCR (HR, 4.197; 95% CI, 3.016-5.842; p < 0.001). A higher risk score also represented the possibility of immunotherapy response. Finally, the potential roles of m5C modification signature were validated in the testing dataset. Conclusions: Our study revealed the potential roles of m5C modification in the PRAD BCR and TME diversity, which may provide new insight into the field of prostate cancer in future research.