Abstract
BACKGROUND: N6-methyladenosine (m6A) is the most frequent RNA modification in mammals, and its role in bladder cancer (BC) remains rarely revealed. OBJECTIVE: To predict the value of m6A-related genes in prognosis and immunity in BC. METHODS: We performed multiple omics analysis of 618 TCGA and GEO patients and used principal component analysis (PCA) to calculate the m6A score for BC patients. RESULTS: We described the multiple omics status of 23 m6A methylation-related genes (MRGs), and four m6A clusters were identified, which showed significant differences in immune infiltration and biological pathways. Next, we intersected the differential genes among m6A clusters, and 11 survival-related genes were identified, which were used to calculate the m6A score for the patients. We found that the high-score (HS) group showed lower tumor mutation burden (TMB) and TP53 mutations and better prognosis than the low-score (LS) group. Lower immune infiltration, higher expression of PD-L1, PD-1, and CTLA4, and higher immune dysfunction and immune exclusion scores were identified in the LS group, suggesting a higher possibility of immune escape. Finally, the experimental verification shows that the m6A related genes, such as IGFBP1, plays an important role in the growth and metastasis of bladder cancer. CONCLUSIONS: These findings revealed the important roles of m6A MRGs in predicting prognosis, TMB status, TP53 mutation, immune functions and immunotherapeutic response in BC.