Comprehensive Analysis of m5C Methylation Regulatory Genes and Tumor Microenvironment in Prostate Cancer

BACKGROUND: 5-Methylcytidine (m5C) methylation is an emerging epigenetic modification in recent years, which is associated with the development and progression of various cancers. However, the prognostic value of m5C regulatory genes and the correlation between m5C methylation and the tumor microenvironment (TME) in prostate cancer remain unknown. METHODS: In the current study, the genetic and transcriptional alterations and prognostic value of m5C regulatory genes were investigated in The Cancer Genome Atlas and Gene Expression Omnibus datasets. Then, an m5C prognostic model was established by LASSO Cox regression analysis. Gene set variation analyses (GSVA), gene set enrichment analysis (GSEA), clinical relevance, and TME analyses were conducted to explain the biological functions and quantify the TME scores between high-risk and low-risk subgroups. m5C regulatory gene clusters and m5C immune subtypes were identified using consensus unsupervised clustering analysis. The Cell-type Identification By Estimating Relative Subsets of RNA Transcripts algorithm was used to calculate the contents of immune cells. RESULTS: TET3 was upregulated at transcriptional levels in PCa compared with normal tissues, and a high TET3 expression was associated with poor prognosis. An m5C prognostic model consisting of 3 genes (NSUN2, TET3, and YBX1) was developed and a nomogram was constructed for improving the clinical applicability of the model. Functional analysis revealed the enrichment of pathways and the biological processes associated with RNA regulation and immune function. Significant differences were also found in the expression levels of m5C regulatory genes, TME scores, and immune cell infiltration levels between different risk subgroups. We identified two distinct m5C gene clusters and found their correlation with patient prognosis and immune cell infiltration characteristics. Naive B cells, CD8+ T cells, M1 macrophages and M2 macrophages were obtained and 2 m5C immune subtypes were identified. CTLA4, NSUN6, TET1, and TET3 were differentially expressed between immune subtypes. The expression of CTLA4 was found to be correlated with the degree of immune cell infiltration. CONCLUSIONS: Our comprehensive analysis of m5C regulatory genes in PCa demonstrated their potential roles in the prognosis, clinical features, and TME. These findings may improve our understanding of m5C regulatory genes in the tumor biology of PCa.