- Open Access
M7G-related molecular subtypes can predict the prognosis and correlate with immunotherapy and chemotherapy responses in bladder cancer patients
European Journal of Medical Research volume 28, Article number: 55 (2023)
N7-methylguanosine (m7G) is closely associated with tumor prognosis and immune response in many cancer types. The correlation between m7G and bladder cancer (BC) needs further study. We aimed to orchestrate molecular subtypes and identify key genes for BC from the perspective of m7G.
RNA-seq and clinical data of BC patients were extracted from TCGA and GSE13507 datasets. The patients were subtyped by “ConsensusClusterPlus” and “limma.” The clusters were validated by the Kaplan‒Meier curves, univariable and multivariate Cox regression models, the concordance index, and calibration curves. The immunotherapy response was evaluated by immune checkpoints, immune infiltration, TIDE score, and IMvigor210 cohort. Genomics of Drug Sensitivity in Cancer was utilized to predict the chemotherapy response between the clusters.
The m7G-related cluster was ultimately established by EIF4G1, NUDT11, NUDT10, and CCNB1. The independent prognostic value of the m7G-related cluster was validated by the TCGA and GSE13507 datasets. The cluster was involved in immune-associated pathways, such as neutrophil degranulation, antigen processing cross-presentation, and signaling by interleukins pathways. Meanwhile, cluster 2 was positively correlated with many immune checkpoints, such as CD274, CTLA4, HAVCR2, LAG3, PDCD1, and PDCD1LG2. The cluster 2 was significantly correlated with a higher TIDE score than the cluster 1. Furthermore, in the IMvigor210 cohort, patients in the cluster 1 had a higher response rate than those in the cluster 2. Patients in the cluster 2 were sensitive to many chemotherapies.
We successfully determined molecular subtypes and identified key genes for BC from the perspective of m7G, thereby providing a roadmap for the evolution of immunotherapy and precision medicine.
Bladder cancer (BC) ranks as the 10th most commonly diagnosed carcinoma and caused 213,000 deaths globally in 2020 . According to the American Cancer Society estimation, BC is predicted to take fourth place among newly diagnosed cancer cases among American men . There are significant differences in morbidity between sexes. However, although the incidence of BC in females is low, the prognosis is worse than that in males . Radical cystectomy (RC) is the standard treatment for limited BC. However, the health expenditure of BC and considerable morbidity caused by RC have brought heavy economic, living, and mental burdens to BC patients. To improve these situations, many therapies have been performed clinically, such as surgery, chemotherapy, immunotherapy, and radiotherapy [4, 5]. Despite various therapies, BC is still prone to recurrence, progression and even death, which may be due to its high genetic heterogeneity . Under these circumstances, clinicians are striving to identify new powerful therapeutic targets. Novel immunotherapy is an outstanding representative of these efforts. Furthermore, to fully take advantage of the existing treatments, a powerful biomarker also requires much work from researchers, which can help determine the optimal treatment option . With a useful prognostic biomarker, doctors can decide whether aggressive therapies, such as RC, should be performed.
N7-methylguanosine (m7G), one of the prevalent posttranscriptional modifications of RNA, is positively charged and produced by the addition of a methyl group at position N7 of riboguanosine . In addition to participating in efficient pre-mRNA splicing and recruitment, the m7G cap is also required binding to the nuclear cap-binding complex . Moreover, m7G can also occur on the caps of mRNAs and internally in mRNAs. Compared with unm7G-modified mRNAs, internal m7G-modified mRNAs have higher translation efficiency . The core functional structure of m7G is the methyltransferase-like 1 (METTL1)/WD repeat domain (WDR4) complex, which can regulate RNA processing, metabolism, and function. Prior studies have noted the essential role of WDR4 in benign diseases, such as Galloway-Mowat syndrome , speech and language delay , and microcephalic primordial dwarfism . Recently, m7G has been reported to participate in tumorigenesis and can predict the prognosis and treatment outcome of tumors . Abnormal changes in m7G influence the occurrence and progression of malignant tumors. In glioblastoma, a study constructed a signature by four m7G-related genes (TMOD2, CACNG2, PLOD3, and TMSB10) and could predict the prognosis of glioblastoma patients . Li et al.  identified m7G-related clusters that could predict the prognosis of hepatocellular carcinoma patients. Furthermore, in the field of urological tumors, three studies demonstrated the prognostic value of m7G-related mRNAs, lncRNAs and miRNAs in clear cell renal cell carcinoma [16,17,18]. Meanwhile, M7G-related genes could influence the immune infiltration of prostate cancer [19, 20]. However, a few studies have discussed the correlation between m7G and BC through single m7G-related gene and review [13, 21]. The role of m7G in BC needs to be further explored. Therefore, this evidence inspired us to identify the relationship between m7G and BC.
To explore and prove this correlation, we collected BC data from online databases to orchestrate molecular subtypes and identify key genes for BC from the perspective of m7G to provide a roadmap for the evolution of precision medicine. Then, various methods were employed to validate the clusters.
Materials and methods
Microarray data and molecular subtypes
We extracted the RNA-seq and clinical data of BC patients from the Cancer Genome Atlas (www.gdc.cancer.gov, TCGA) to analyze the expression and function of m7G-related genes. In terms of prognosis analysis, patients with postoperative survival times shorter than 30 days were excluded. For external validation, GSE13507, including clinical information of BC patients, was used to evaluate the prognostic value of the m7G-related clusters. The m7G-related genes were collected from a combination of the literature , which selected 42 m7G-related genes by screening previous reports and the GSEA database (http://www.gsea-msigdb.org/gsea/index.jsp).
All 42 m7G-related genes were incorporated into the lasso regression model, which could apply penalties to these genes. The penalty parameter (λ) for the model was determined by tenfold cross-validation following the minimum criteria. Four genes (EIF4G1, NUDT11, NUDT10, and CCNB1) were selected. Then, the patients in the TCGA and GSE13507 datasets were subtyped by the R packages “ConsensusClusterPlus” and “limma” using the four genes. The consensus matrix k value denoted the number of m7G-related clusters. Then, based on a linear regression-based algorithm and a label propagation algorithm, GeneMANIA (www.genemania.org)  was employed to predict the interacting proteins of EIF4G1, NUDT11, NUDT10, and CCNB1.
Validation of the m7G-related cluster
Depending on the survival results of the TCGA and GSE13507 datasets, the prognostic value of the m7G-related cluster, including overall survival (OS), cancer-specific survival (CSS), and progression-free survival (PFS), was evaluated by Kaplan‒Meier curves. Due to lacking the data of PFS time and CSS time data, we did not plot these Kaplan‒Meier curves in GSE13507. Furthermore, the OS in different clinical subgroups, such as male, T2_4 stage, and age > 70 years subgroups, was also compared by the Kaplan‒Meier curves. The independent prognostic value of the cluster was estimated by univariable and multivariable Cox regression models based on the TCGA and GSE13507 datasets. In detail, according to the results of the univariable Cox regression model, a multivariable Cox regression model was constructed using factors with P < 0.1.
According to the results of above multivariable Cox regression models, two nomograms were carefully established based on the cluster and clinical parameters in the TCGA and GSE13507 datasets. For validation, the concordance index (C-index) and calibration curves were used to validate the reliability and accuracy of the nomograms.
Functional enrichment analysis
First, based on TCGA data, Gene Ontology (GO) enrichment analysis was used to explore the possible function of the cluster in molecular function (MF), biological process (BP), and cellular component (CC). Based on the GO results, GO terms with P < 0.05 and Q < 0.05 were selected and displayed in bubble plots. Similarly, the enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were screened with P < 0.05 and Q < 0.05. For further evaluation of the potential pathways, REACTOME pathways were carefully enriched by Gene Set Enrichment Analysis (GSEA) with P < 0.05 and FDR < 25%.
The correlation between the m7G cluster and immune checkpoints, immune infiltration, TCIA, immunotherapy response, and chemotherapy response
We analyzed and compared the expression of 20 checkpoints between the clusters. For estimation of immune cells in TCGA BC samples, the CIBERSORTx website (https://cibersortx.stanford.edu) was used to calculate the contents of 22 types of immune cells in BC samples. Then, we compared the contents of infiltrated immune cells between the clusters.
To further verify the above results, based on the TCGA BC cohort, the Tumor Immune Dysfunction and Exclusion (TIDE) algorithm was used to predict and compare the immunotherapy response between the cluster 1 and cluster 2. Meanwhile, we compared the RNA-seq and clinical data of BC patients receiving immunotherapy in the real world. The patients in the IMvigor210 cohort were subtyped by the R packages “ConsensusClusterPlus” and “limma” using the above four genes. Then, the IMvigor210 trial (EGAD00001003977), including patients with advanced or metastatic BC who were treated with an anti-PD-L1 agent (atezolizumab), was employed to assess the correlation between therapeutic outcomes and the clusters . Moreover, Genomics of Drug Sensitivity in Cancer (GDSC), a public pharmacogenomics database, was used to predict various chemotherapy responses for BC patients. The half-maximal inhibitory concentration (IC50) calculated by the “pRRophetic” package in R software was the main endpoint of the chemotherapy response and drug sensitivity, which was also compared between the clusters.
According to the normality and quality of variances of the data, one-way ANOVA or the Mann‒Whitney U test was used to perform statistical analysis of three or more continuous variables. Quantitative data in two groups were compared using Student’s t test. All analyzed data are displayed as the standard deviation (SD). A P < 0.05 was considered significant for all analyses, which were performed using R version 3.6.3 and relative packages. nsP ≥ 0.05; *P < 0.05; **P < 0.01; ***P < 0.001.
Gene selection and basal data
The workflow of our study is shown in Fig. 1. The 42 m7G-related genes displayed in Additional file 1: Table S1. As presented in Fig. 2A, the lasso regression model screened EIF4G1, NUDT11, NUDT10, and CCNB1. According to the results of the protein–protein interaction network, the four genes were mainly associated with CDK2, PPIP5K1, PPIP5K2, IP6K3, IP6K1, and CDK1 (Fig. 2B). Finally, these four m7G-related genes divided BC patients in the TCGA dataset into two clusters (Fig. 2C, consensus matrix k = 2). The expression of these four genes in the m7G-related clusters is shown in Additional file 2: Fig. S1. In the TCGA dataset, patients in the cluster 2 had worse OS (Fig. 2D, P = 0.002), CSS (Fig. 2E, P = 0.001), and PFS (Fig. 2F, P = 0.001) than those in the cluster 1. Similarly, patients in the GSE13507 dataset were classified into two clusters by the four m7G-related genes (Fig. 2G, consensus matrix k = 2). In the GSE13507 cohort, patients in the cluster 2 were significantly correlated with worse OS than patients in the cluster 1 (Fig. 2H, P = 0.014).
The basic information of the TCGA and GSE13507 datasets is displayed in Table 1 and Additional file 3: Table S2, respectively. In TCGA dataset, the cluster 2 was significantly associated with WHO high grade, poor OS, CSS, and PFS (Table 1).
Validation of the m7G-related cluster
The Kaplan‒Meier curves was utilized to estimate the prognostic value of the cluster in clinical subgroups. In the TCGA database, the cluster exhibited significant prognostic value in the male (Fig. 3A, P = 0.014), no lymphatic vessel invasion (Fig. 3C, P = 0.013), T2_4 stage (Fig. 3D, P = 0.021), no lymph node metastasis (Fig. 3E, P = 0.004), and no distant metastasis (Fig. 3F, P = 0.003) subgroups, except for the age > 70 years (Fig. 3B, P = 0.145) subgroup. In the GSE13507 dataset, the cluster exhibited significant prognostic value in the female (Fig. 3G, P = 0.043), Ta_1 stage (Fig. 3J, P = 0.044), no lymph node metastasis (Fig. 3K, P = 0.021), and no distant metastasis (Fig. 3L, P = 0.008) subgroups. No significant differences were found between the cluster 1 and cluster 2 in the age ≤ 70 years (Fig. 3H, P = 0.196) and WHO low-grade (Fig. 3I, P = 0.256) subgroups.
Construction and validation of a nomogram
Univariate Cox regression model was used to analyze the prognostic value of the cluster and clinical parameters in the TCGA (Fig. 4A) and GSE13507 (Fig. 4C). Then, the multivariate Cox regression models consisted of factors with P < 0.1 in these two datasets. The results identified that the cluster could independently predict the prognosis of BC patients in the TCGA (Fig. 4B, P = 0.007) and GSE13507 (Fig. 4D, P = 0.009) datasets. According to the results of the multivariate Cox regression models, we constructed two nomograms based on the factors with P < 0.05 in the TCGA (Fig. 4E) and GSE13507 (Fig. 4F). The TCGA nomogram had a C-index of 0.654 (0.632–0.677). In the GSE13507 nomogram, the C-index was 0.645 (0.614–0.677). In the TCGA nomogram, the calibration curves showed good agreement between the predicted value and the true value at 1, 3, and 5 years (Fig. 4G). Similarly, the calibration curves of GSE13507 also showed good agreement between the predicted value and the true value at 2, 3, and 5 years (Fig. 4H).
The results of GO, KEGG, and gene set enrichment analyses
The data from the TCGA dataset were used to perform GO, KEGG, and GSEA analyses. Figure 5A presents the GO results, showing that the cluster was associated with neutrophil chemotaxis, granulocyte chemotaxis, neutrophil migration in BP, cornified envelope, intermediate filament, intermediate filament cytoskeleton in CC, interleukin-1 receptor binding, cytokine activity, and chemokine receptor binding in MF. According to the KEGG results, the cluster was significantly involved in the IL-17 signaling pathway, chemical carcinogenesis, drug metabolism—cytochrome P450, cytokine‒cytokine receptor interaction, and metabolism of xenobiotics by cytochrome P450 pathways (Fig. 5B). In further GSEA, the cluster was significantly involved in immune-associated pathways, such as neutrophil degranulation, signaling by interleukins, immunoregulatory interactions between lymphoid and nonlymphoid cells, antigen processing cross-presentation, and downstream signaling events of b cell receptor BCR pathways (Fig. 5C). Moreover, the cluster also enriched metabolism-associated pathways, including synthesis of leukotrienes LT and eoxins ex, cytochrome p450 arranged by substrate type, fatty acid metabolism, keratan sulfate degradation, and fatty acids pathways (Fig. 5D).
The correlation between the risk score and immune checkpoints and immune cell infiltration
Based on TCGA data, we identified the presentation of the cluster in immune-related analysis. As shown in Fig. 5E, the cluster 2 was positively correlated with many immune checkpoints, such as CD274, CTLA4, HAVCR2, LAG3, PDCD1, and PDCD1LG2, while SIGLEC15 was highly expressed in the cluster 1. The results, as shown in Fig. 5F, indicated that samples in the cluster 2 had highly activated CD4+ memory T cells, NK resting cells, M0 macrophages, M1 macrophages, and resting mast cells infiltration. In the cluster 1, samples had more B memory cells, B plasma cells, resting memory CD4+ T cells, T regulatory cells, monocyte cells, and activated mast cells infiltration.
The correlation between the m7G-related risk score and immunotherapy and chemotherapy
As shown in Fig. 6A, the cluster 2 had a significantly higher TIDE score than the cluster 1 (P < 0.01). Meanwhile, the response rate of cluster 1 was 48.4%, which was higher than that of cluster 2 (39.8%) (Fig. 6B). In accordance with the TCGA and GSE13507 datasets, patients in the IMvigor210 cohort were also classified into two clusters (Fig. 6C, consensus matrix k = 2). According to the therapeutic outcome of IMvigor210, the response rate of cluster 1 was 29.1% (30/73), which was higher than the 18.4% (12/53) of cluster 2 (Fig. 6D). Both the TIDE and IMvigor210 results supported that patients in the cluster 1 were more likely to benefit from immunotherapy. In chemotherapy, patients in cluster2 were significantly sensitive to cisplatin, docetaxel, doxorubicin, gemcitabine, mitomycin C, and vinblastine based on the TCGA dataset (Fig. 6E, all P < 0.001).
As the key role of RNA modification, m7G is involved in the processing, metabolism, and function of mRNAs, microRNAs (miRNAs), tRNAs, and rRNAs . Recently, increasing evidence has found that m7G-related genes can promote or inhibit the processes of various cancers [14, 19, 21, 24,25,26]. Therefore, this study orchestrated molecular subtypes and identified key genes for BC from the perspective of m7G. In our results, the m7G-related cluster could independently predict the prognosis of BC patients. Furthermore, there was a significant correlation between the cluster and immunotherapy and chemotherapy. These results suggested that the m7G-related cluster mainly paved a roadmap for precision medicine.
As a part of a scaffold component of eukaryotic translation initiation factor 4F, eukaryotic translation initiation factor 4G 1 (EIF4G1) is closely associated with various tumors, especially in solid tumors . Nudix hydrolase 10 (NUDT10), located in Xp11.22, increases promoter methylation and is associated with OS in prostate cancer and gastric cancer [28, 29]. NUDT10 and Nudix hydrolase 11 (NUDT11) can promote prostate carcinogenesis and increase prostate cancer susceptibility by being involved in a variety of biological processes and mediating cellular stress responses . CyclinB1 (CCNB1) is mainly involved in cell cycle regulation and is associated with the epithelial–mesenchymal transition (EMT) process . In BC, CCNB1 is an oncogene and could predict the prognosis of nonmuscle-invasive patients .
In this study, we validated the prognostic value of the cluster in OS, CSS, and PFS for BC patients in the TCGA dataset. For further estimation, the GSE13507 cohort was used to identify the stability of this result, which also suggested that the cluster could independently predict the prognosis of BC patients. This finding was also reported by Ye et al. . They also identified that the m7G-related signature could predict the prognosis of prostate cancer. Similarly, Liu et al.  also identified the prognostic value of m7G-related lncRNAs in prostate cancer. Moreover, three studies demonstrated the prognostic value of m7G-related mRNAs, lncRNAs, and miRNAs in clear cell renal cell carcinoma [16,17,18]. Furthermore, we also assessed the prognostic value of the risk score in various subgroups. The subgroup results identified that the risk score could predict the prognosis of BC patients in TCGA and external validation cohorts. Thus, we speculate that the m7G-related cluster might be a potential biomarker for BC.
After assessing the prognostic value, we tried to explore the function of the m7G-related cluster. The GO, KEGG, and GSEA enrichment results indicated that the cluster was significantly associated with immune- and metabolism-associated pathways. In detail, the pathways involved in the regulation of immune factors, such as neutrophil degranulation and interleukins. Furthermore, these pathways could also influence the regulation of immune cells. Recently, a study reported that the cap structure of m7G can induce immune evasion by enhancing the identity of RNA as a ‘self-molecule’ for cellular recognition factors . This finding might explain why the cluster was mainly involved in immune-related pathways. These results inspired us to verify the correlation between the cluster and immunotherapy. In immune checkpoint analysis, CD274, CTLA4, HAVCR2, LAG3, PDCD1, PDCD1LG2, etc., while SIGLEC15 was highly expressed in cluster 1. No consensus has been reached on the predictive value of immune checkpoints in BC immunotherapy, such as PD-L1 and CTLA4 . In theory, an immunosuppressive microenvironment promotes the expression of immune checkpoints, as presented in the present study. To validate this hypothesis, we performed immune infiltration analysis. The results revealed that samples in cluster 2 had high NK resting cell, macrophage M0 cell, macrophage M1 cell, and mast resting cell infiltration. Of these, macrophages consist of 3 classes: M0, M1, and M2 macrophages. M0 macrophages can differentiate into M1 or M2 macrophages and were significantly associated with poor survival, which was consistent with our study . Based on the results of basic experiments, Bopp et al.  has reported that tumor-associated macrophages can promote immune evasion and tumor growth by expressing G-protein-coupled receptors. In addition to assessing the institution of the immune checkpoint and the tumor microenvironment, we also adopted the TIDE score to predict the response rate between cluster 1 and cluster 2. The immunotherapy response rate was negatively associated with the TIDE score . In the TIDE results, culster1 was significantly correlated with a lower TIDE score than cluster 2, which was consistent with the previous hypothesis that patients in cluster 1 were more likely to benefit from immunotherapy. For further validation of the results with real-world data, IMvigor210, including patients with advanced or metastatic BC who were treated with an anti-PD-L1 agent (atezolizumab), was employed to compare the therapeutic outcomes between cluster 1 and cluster 2. In accordance with the above results, a higher response rate was presented in cluster 1 patients. Taken together, we believe that patients in cluster 1 are more likely to benefit from immunotherapy, at least from anti-PD-L1 therapy. In addition to immunotherapy, we also explored the predictive function of the cluster in chemotherapy. In chemotherapy, patients in cluster 2 were associated with low IC50 values of cisplatin, docetaxel, doxorubicin, gemcitabine, mitomycin C, and vinblastine, which suggested that patients in cluster 2 might accept chemotherapy rather than immunotherapy. This evidence suggested that we successfully determined molecular subtypes and identified key genes for BC from the perspective of m7G to pave a roadmap for the evolution of immunotherapy and precision medicine.
A few limitations should be noted. Although we identified the expression and prognostic value of the signature in TCGA and external cohorts, basic experiments should also be performed, which is the next work of us. In this study, we also found some sensitive drugs to the four investigated genes, which need to be further studied in vivo and in vitro.
The results in TCGA and external validation cohorts identified that we successfully determined molecular subtypes and identified key genes for BC from the perspective of m7G, thereby paving a roadmap for the evolution of immunotherapy and precision medicine.
Availability of data and materials
All data from this study were downloaded from online databases. Therefore, everyone can get the data online. Further inquiries can be directed to the corresponding author.
WD Repeat Domain 4
Nudix Hydrolase 10
The Cancer Genome Atlas
Gene Set Enrichment Analysis
Kyoto Encyclopedia of Genes and Genomes
Genomics of Drug Sensitivity in Cancer
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global Cancer Statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71:209–49.
Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics. CA Cancer J Clin. 2022;72(2022):7–33.
Feng D, Li A, Hu X, Lin T, Tang Y, Han P. Comparative effectiveness of open, laparoscopic and robot-assisted radical cystectomy for bladder cancer: a systematic review and network meta-analysis. Ital J Urol Nephrol. 2020;72:251–64.
EAU Guidelines. Edn. Presented at the EAU Annual Congress Amsterdam 2022. ISBN 978-94-92671-16-5.
Li D-X, Wang X-M, Feng D-C, Han P. Neutrophil-to-lymphocyte ratio (NLR) during induction is a better predictor than preoperative NLR in non-muscle-invasive bladder cancer receiving Bacillus Calmette-GuÉRin. Asian J Surg. 2022. https://doi.org/10.1016/j.asjsur.2022.08.108.
Li DX, Wang XM, Tang Y, Yang YB, Feng DC, Li A, Zhang FC, Bai YJ, Han P. Prognostic value of preoperative neutrophil-to-lymphocyte ratio in histological variants of non-muscle-invasive bladder cancer. Investig Clin Urol. 2021;62:641–9.
Cowling VH. Regulation of mRNA cap methylation. Biochem J. 2009;425:295–302.
Fresco LD, Buratowski S. Conditional mutants of the yeast mRNA capping enzyme show that the cap enhances, but is not required for, mRNA splicing. RNA. 1996;2:584–96.
Malbec L, Zhang T, Chen YS, Zhang Y, Sun BF, Shi BY, Zhao YL, Yang Y, Yang YG. Dynamic methylome of internal mRNA N(7)-methylguanosine and its regulatory role in translation. Cell Res. 2019;29:927–41.
Braun DA, Shril S, Sinha A, Schneider R, Tan W, Ashraf S, Hermle T, Jobst-Schwan T, Widmeier E, Majmundar AJ, Daga A, Warejko JK, Nakayama M, Schapiro D, Chen J, Airik M, Rao J, Schmidt JM, Hoogstraten CA, Hugo H, Meena J, Lek M, Laricchia KM, Bagga A, Hildebrandt F. Mutations in WDR4 as a new cause of Galloway-Mowat syndrome. Am J Med Genet A. 2018;176:2460–5.
Chen X, Gao Y, Yang L, Wu B, Dong X, Liu B, Lu Y, Zhou W, Wang H. Speech and language delay in a patient with WDR4 mutations. Eur J Med Genet. 2018;61:468–72.
Shaheen R, Abdel-Salam GM, Guy MP, Alomar R, Abdel-Hamid MS, Afifi HH, Ismail SI, Emam BA, Phizicky EM, Alkuraya FS. Mutation in WDR4 impairs tRNA m(7)G46 methylation and causes a distinct form of microcephalic primordial dwarfism. Genome Biol. 2015;16:210.
Cheng W, Gao A, Lin H, Zhang W. Novel roles of METTL1/WDR4 in tumor via m(7)G methylation. Mol Ther Oncolytics. 2022;26:27–34.
Wu X, Li C, Wang Z, Zhang Y, Liu S, Chen S, Chen S, Liu W, Liu X. A bioinformatic analysis study of m(7)G regulator-mediated methylation modification patterns and tumor microenvironment infiltration in glioblastoma. BMC Cancer. 2022;22:729.
Chen J, Yao S, Sun Z, Wang Y, Yue J, Cui Y, Yu C, Xu H, Li L. The pattern of expression and prognostic value of key regulators for m(7)G RNA methylation in hepatocellular carcinoma. Front Genet. 2022;13: 894325.
Dong K, Gu D, Shi J, Bao Y, Fu Z, Fang Y, Qu L, Zhu W, Jiang A, Wang L. Identification and verification of m(7)G modification patterns and characterization of tumor microenvironment infiltration via multi-omics analysis in clear cell renal cell carcinoma. Front Immunol. 2022;13: 874792.
Ming J, Wang C. N7-methylguanosine-related lncRNAs: integrated analysis associated with prognosis and progression in clear cell renal cell carcinoma. Front Genet. 2022;13: 871899.
Hong P, Du H, Tong M, Cao Q, Hu D, Ma J, Jin Y, Li Z, Huang W, Tong G. A Novel M7G-related microRNAs risk signature predicts the prognosis and tumor microenvironment of kidney renal clear cell carcinoma. Front Genet. 2022;13: 922358.
Mei W, Jia X, Xin S, Liu X, Jin L, Sun X, Zhang JX, Zhang B, Yang G, Chen P, Ye L. A N(7)-methylguanine-related gene signature applicable for the prognosis and microenvironment of prostate cancer. J Oncol. 2022;2022:8604216.
Xin S, Deng Y, Mao J, Wang T, Liu J, Wang S, Song X, Song W, Liu X. Characterization of 7-methylguanosine identified biochemical recurrence and tumor immune microenvironment in prostate cancer. Front Oncol. 2022;12: 900203.
Ying X, Liu B, Yuan Z, Huang Y, Chen C, Jiang X, Zhang H, Qi D, Yang S, Lin S, Luo J, Ji W. METTL1-m(7) G-EGFR/EFEMP1 axis promotes the bladder cancer development. Clin Transl Med. 2021;11: e675.
Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, Franz M, Grouios C, Kazi F, Lopes CT, Maitland A, Mostafavi S, Montojo J, Shao Q, Wright G, Bader GD, Morris Q. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010;38:W214–20.
Mariathasan S, Turley SJ, Nickles D, Castiglioni A, Yuen K, Wang Y, Kadel EE III, Koeppen H, Astarita JL, Cubas R, Jhunjhunwala S, Banchereau R, Yang Y, Guan Y, Chalouni C, Ziai J, Şenbabaoğlu Y, Santoro S, Sheinson D, Hung J, Giltnane JM, Pierce AA, Mesh K, Lianoglou S, Riegler J, Carano RAD, Eriksson P, Höglund M, Somarriba L, Halligan DL, van der Heijden MS, Loriot Y, Rosenberg JE, Fong L, Mellman I, Chen DS, Green M, Derleth C, Fine GD, Hegde PS, Bourgon R, Powles T. TGFβ attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cells. Nature. 2018;554:544–8.
Zhang C, Zhou D, Wang Z, Ju Z, He J, Zhao G, Wang R. Risk model and immune signature of m7G-related lncRNA based on lung adenocarcinoma. Front Genet. 2022;13: 907754.
Zhou K, Yang J, Li X, Xiong W, Zhang P, Zhang X. N7-methylguanosine regulatory genes profoundly affect the prognosis, progression, and antitumor immune response of hepatocellular carcinoma. Front Surg. 2022;9: 893977.
Wu X, Tang J, Cheng B. Oral squamous cell carcinoma gene patterns connected with RNA methylation for prognostic prediction. Oral Dis. 2022. https://doi.org/10.1111/odi.14341.
Zhang S, Huang N, Pan X, Zang JL, Guan XX, Zhang JH, Liu LC, Lei XY. Role of eukaryotic translation initiation factor 4G in tumor. Acta Physiol Sin. 2016;68:207–14.
Li XY, Wang SL, Chen DH, Liu H, You JX, Su LX, Yang XT. Construction and validation of a m7G-related gene-based prognostic model for gastric cancer. Front Oncol. 2022;12: 861412.
Kamdar S, Isserlin R, Van der Kwast T, Zlotta AR, Bader GD, Fleshner NE, Bapat B. Exploring targets of TET2-mediated methylation reprogramming as potential discriminators of prostate cancer progression. Clin Epigenet. 2019;11:54.
Achilla C, Papavramidis T, Angelis L, Chatzikyriakidou A. The implication of X-linked genetic polymorphisms in susceptibility and sexual dimorphism of cancer. Anticancer Res. 2022;42:2261–76.
Kim SK, Roh YG, Park K, Kang TH, Kim WJ, Lee JS, Leem SH, Chu IS. Expression signature defined by FOXM1-CCNB1 activation predicts disease recurrence in non-muscle-invasive bladder cancer. Clin Cancer Res. 2014;20:3233–43.
Drazkowska K, Tomecki R, Warminski M, Baran N, Cysewski D, Depaix A, Kasprzyk R, Kowalska J, Jemielity J, Sikorski PJ. 2’-O-Methylation of the second transcribed nucleotide within the mRNA 5’ cap impacts the protein production level in a cell-specific manner and contributes to RNA immune evasion. Nucleic Acids Res. 2022;1862:272.
Labadie BW, Balar AV, Luke JJ. Immune checkpoint inhibitors for genitourinary cancers: treatment indications, investigational approaches and biomarkers. Cancers (Basel). 2021;13(21):5415.
Locati M, Curtale G, Mantovani A. Diversity, mechanisms, and significance of macrophage plasticity. Annu Rev Pathol. 2020;15:123–47.
Pan Y, Yu Y, Wang X, Zhang T. Tumor-associated macrophages in tumor immunity. Front Immunol. 2020;11: 583084.
Jiang P, Gu S, Pan D, Fu J, Sahu A, Hu X, Li Z, Traugh N, Bu X, Li B, Liu J, Freeman GJ, Brown MA, Wucherpfennig KW, Liu XS. Signatures of T cell dysfunction and exclusion predict cancer immunotherapy response. Nat Med. 2018;24:1550–8.
This study was supported by the Project of Health Commission of Sichuan Province (21PJ041).
Ethics approval and consent to participate
This study does not involve human participants and animals. Therefore, it does not require ethical review and approval.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1.
Additional file 2: Figure S1.
Comparison of four genes expression between the cluster1 and cluster2 in the TCGA dataset (A) and GSE13507 dataset (B).
Additional file 3: Table S2.
The clinicopathologic characteristics of the GSE13507 included patients.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Li, Dx., Feng, Dc., Wang, Xm. et al. M7G-related molecular subtypes can predict the prognosis and correlate with immunotherapy and chemotherapy responses in bladder cancer patients. Eur J Med Res 28, 55 (2023). https://doi.org/10.1186/s40001-023-01012-x
- Bladder cancer