An innovative model based on N7-methylguanosine-related lncRNAs for forecasting prognosis and tumor immune landscape in bladder cancer

Abstract

Abstract Background As a novel type of the prevalent post-transcriptional modifications, N7-methylguanosine (m7G) modification is essential in the tumorigenesis, progression, and invasion of many cancers, including bladder cancer (BCa). However, the integrated roles of m7G-related lncRNAs in BCa remain undiscovered. This study aims to develop a prognostic model based on the m7G-related lncRNAs and explore its predictive value of the prognosis and anti-cancer treatment sensitivity. Methods We obtained RNA-seq data and corresponding clinicopathological information from the TCGA database and collected m7G-related genes from previous studies and GSEA. Based on LASSO and Cox regression analysis, we developed a m7G prognostic model. The Kaplan–Meier (K-M) survival analysis and ROC curves were performed to evaluate the predictive power of the model. Gene set enrichment analysis (GSEA) was conducted to explore the molecular mechanisms behind apparent discrepancies between the low- and high-risk groups. We also investigated immune cell infiltration, TIDE score, TMB, the sensitivity of common chemotherapy drugs, and the response to immunotherapy between the two risk groups. Finally, we validated the expression levels of these ten m7G-related lncRNAs in BCa cell lines by qRT-PCR. Results We developed a m7G prognostic model (risk score) composed of 10 m7G-related lncRNAs that are significantly associated with the OS of BCa patients. The K-M survival curves revealed that the high-risk group patients had significantly worse OS than those in the low-risk group. The Cox regression analysis confirmed that the risk score was a significant independent prognostic factor for BCa patients. We found that the high-risk group had higher the immune scores and immune cell infiltration. Furthermore, the results of the sensitivity of common anti-BCa drugs showed that the high-risk group was more sensitive to neoadjuvant cisplatin-based chemotherapy and anti-PD1 immunotherapy. Finally, qRT-PCR revealed that AC006058.1, AC073133.2, LINC00677, and LINC01338 were significantly downregulated in BCa cell lines, while the expression levels of AC124312.2 and AL158209.1 were significantly upregulated in BCa cell lines compared with normal cell lines. Conclusion The m7G prognostic model can be applied to accurately predict the prognosis and provide robust directions for clinicians to develop better individual-based and precise treatment strategies for BCa patients.