Single-cell and spatial transcriptomics reveal 5-methylcytosine RNA methylation regulators immunologically reprograms tumor microenvironment characterizations, immunotherapy response and precision treatment of clear cell renal cell carcinoma

•We first revealed the relationship between the estimation of stromal, immune cells score and 5-methylcytosine (m5C) regulators at the pan-cancer level (33 types of cancer), and then depicted the association between m5C regulators and tumor-infiltrating immune cells in a meta-cohort of more than 1000 patients with RCC.•Next, we established a novel 5-methylcytosine RNA modification related gene index (M5CRMRGI).•Additionally, we studied the effect of m5C on tumor microenvironment (TME) distribution, using single-cell sequencing data for in-depth analysis, and finally verified it using spatial sequencing data.•Notably, we predicted the efficacy of m5C on adjuvant therapy, which was verified in real-world immunological and targeted therapy cohorts. Clear cell Renal Cell Carcinoma (ccRCC) is a highly heterogeneous disease, making it challenging to predict prognosis and therapy efficacy. In this study, we aimed to explore the role of 5-methylcytosine (m5C) RNA modification in ccRCC and its potential as a predictor for therapy response and overall survival (OS). We established a novel 5-methylcytosine RNA modification-related gene index (M5CRMRGI) and studied its effect on the tumor microenvironment (TME) using single-cell sequencing data for in-depth analysis, and verified it using spatial sequencing data. Our results showed that M5CRMRGI is an independent predictor of OS in multiple datasets and exhibited outstanding performance in predicting the OS of ccRCC. Distinct mutation profiles, hallmark pathways, and infiltration of immune cells in TME were observed between high- and low-M5CRMRGI groups. Single-cell/spatial transcriptomics revealed that M5CRMRGI could reprogram the distribution of tumor-infiltrating immune cells. Moreover, significant differences in tumor immunogenicity and tumor immune dysfunction and exclusion (TIDE) were observed between the two risk groups, suggesting a better response to immune checkpoint blockade therapy of the high-risk group. We also predicted six potential drugs binding to the core target of the M5CRMRGI signature via molecular docking. Real-world treatment cohort data proved once again that high-risk patients were appropriate for immune checkpoint blockade therapy, while low-risk patients were appropriate for Everolimus. Our study shows that the m5C modification landscape plays a role in TME distribution. The proposed M5CRMRGI-guided strategy for predicting survival and immunotherapy efficacy, we reported here, might also be