Comprehensive Analysis of the PD-L1 and Immune Infiltrates of N6-Methyladenosine Related Long Non-Coding RNAs in Bladder Cancer

Comprehensive Analysis of the PD-L1 and Immune Infiltrates of N6-Methyladenosine Related Long Non-Coding RNAs in Bladder Cancer | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Comprehensive Analysis of the PD-L1 and Immune Infiltrates of N6-Methyladenosine Related Long Non-Coding RNAs in Bladder Cancer MQ Xue, YL Wang, JC Wang, XD Wang, XJ Wang, YQ Zhang This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-1341357/v1 This work is licensed under a CC BY 4.0 License Status: Under Review Version 1 posted 9 You are reading this latest preprint version Abstract Background & Aim: Bladder cancer (BLCA) is one of the most frequent genitourinary cancers, with a high rate of morbidity and mortality. The connection of m6A-related lncRNAs with PD-L1 and tumor immune microenvironment (TIME) in BLCA prognosis was extensively investigated in this study, which could suggest novel therapeutic targets for further investigation. Results: 30 m6A-associated lncRNAs with predictive values from the TCGA data set were identified with co-expression analysis. Cluster2 was correlated with a poor prognosis,upregulated PD-L1 expression, and higher immune ratings. Cluster2 had larger amounts of resting CD4 memory-activated T cells, M2 macrophages, neutrophils, and NK cells infiltration. “CHEMOKINE SIGNALING PATHWAY” was the most significantly enriched signaling pathway according to GSEA，which may play an important role in the diffrent immune cell infiltrates between cluster1/2. The risk model for m6A-related lncRNAs could be employed in a prognostic model to predict BLCA prognosis, regardless of other clinical features. Collectively, m6A-related lncRNAs were linked to PD-L1 and TIME, which would dynamically affect the number of tumor-infiltrating immune cells. Conclusion: m6A-related lncRNAs may be key mediators of PD-L1 expression and immune cells infiltration and may strongly affect the TIME of BLCA. Bioinformatics analysis N6-methyladenosine lncRNAs Immune cells infiltration Bladder cancer PD-L1 Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 Figure 9 Figure 10 Figure 11 Figure 12 Figure 13 Figure 14 Figure 15 Figure 16 Introduction Bladder cancer is the ninth most frequent malignancy in the world. It is mainly represented by bladder urothelial cancer (BUC), accounting for more than 90% of bladder cancer, and smoking is considered the leading risk factor [ 1 , 2 ]. Furthermore, BLCA has a significant recurrence rate, and approximately half of patients relapsing following major surgery and developing metastases. [ 3 , 4 ]. Patients with metastatic urothelial carcinoma have reaped extraordinary gains from platinum-based chemotherapy and Immune Checkpoint Inhibitors (ICIs), but the diverse features of BLCA prompt to variable clinical outcomes for BLCA patients even though they received standard therapy. [ 5 ]. Bladder cancer develops through a complex process including aberrant genetic changes and epigenetic aberrations. Epigenetic anomalies can be seen on a range of levels, including DNA. [ 6 ], RNA[ 7 ]and histone modifications [ 8 – 10 ]. Researchers must find novel biomarkers to better predict the prognosis and therapeutic response to BLCA in order to enhance survival and minimize the burden of BLCA patients. Long noncoding RNAs (lncRNAs) are a class of noncoding RNA with a length of more than 200 nucleotides that have been implicated in post-transcriptional regulatory elements that regulate mRNA splicing, stability, and translation [ 11 , 12 ]. LncRNAs may also play a key role in regulating genes that code for cancer-fighting proteins, according to growing evidence[ 13 ]. NKILA lncRNA, for example, aids tumor immune evasion by making T cells more susceptible to activation-induced cell death[ 14 ]. By generating an RNA-protein complex with enhancer of zeste homolog2(EZH2) and boosting the binding of EZH2 and H3K27me3 on the E-cadherin promoter region, lncRNA MRPL23-AS1 facilitated adenoid cystic carcinoma lung metastasis [ 15 ]. However, the potential role of immune-related lncRNA signatures as an effective therapeutic strategy in BLCA is unknown. N6-methyladenosine (m6A) alteration has been shown to regulate lncRNAs in several research [ 16 – 19 ]. M6A, the methylation modification at the sixth N atom of adenine, is the most frequent post-transcriptional alteration on mRNA, mediating > 60% RNA methylation [ 20 , 21 ]. The m6A enzyme's abnormal expression affects tumor cell function and the tumor microenvironment(TIME)[ 22 – 24 ]. Tumor stem cell self-renewal is accelerated by abrupt m6A alteration, which plays a vital role in tumorigenesis [ 25 ]. m6A enzymes and lncRNAs are both excellent diagnostic and prognostic indicators. Evidence is overwhelming indicating m6A-related mRNAs and lncRNAs could be served as promising targets for estimating prognosis in a multitude of malignancies [ 26 – 29 ]. Degradation of m6A reader YTHDF2-modified lncRNA FENDER, for example, enhanced endometrioid endometrial tumorigenesis dramatically [ 30 ]. Wang [ 31 ]et al discovered that m6A-induced lncRNA RP11 expression triggers the malignant and immunosuppressive gastric cancer cells via upregulation of YAP1 expression. Therefore, Identification of novel and reliable prognostic molecular characteristics from various dimensions is critical for identifying the suitable therapy methods and improving the dismal prognosis in BLCA patients. As tumors develop, the body's immune system is triggered to resist tumor development [ 32 ]. Several immunological checkpoints, such as programmed death ligand-1 (PD-L1), have been identified as one of the immune escape mechanisms [ 33 ]. The TIME may influence a patient's response to immune checkpoint inhibitors, hence PD-L1 expression in BLCA should be taken into account while analyzing tumor immunity [ 34 ]. This study focus on systematically assessing the relationships of m6A-related lncRNAs with prognosis, programmed death ligand 1 (PD-L1), and TIME in BLCA. Results Identification Of M6a Related Lncrnas To distinguish mRNAs and lncRNAs, data on m6A-related gene expression were retrieved from transcriptome data. Next, The connection between m6A-related gene expression and lncRNAs was analyzed via co-expression analysis (table S1). A network plot was generated to visualize this correlation (Figure 1A). A forest plot was used to visualize the univariate Cox regression analysis (Figure 1B, table S2). heatmaps and box plots were constructed using the differential expression of prognosis-related m6A-related lncRNAs between tumor and normal tissues, as shown in Figure 2A and 2B. 30 m6A prognosis-related lncRNAs were identified that were differently expressed between tumor and surrounding normal tissues. LINC02604、AC116914.2、 Z84485.1、 ZNF32−AS2、 AC004076.2、 AL138921.1、 TMEM147−AS1、 SNHG20、 PTOV1−AS2 and AC004148.1 were dramatically higher in BLCA tissues than in normal adjacent tissues (p < 0. 001). The expression level of ATP1B3−AS1、AC025280.1、AC012568.1、BDNF−AS was markedly lower in BLCA tissues than in normal tissues (p < 0. 01). Figure 2 . Different expression of m6A prognostic related lncRNAs. A: Boxplot of the difference expression of m6A prognostic related lncRNAs among tumor and normal tissue. *:P < 0.05; **:P < 0.01; ***:P < 0.001. B: Heatmap of the difference expression of m6A prognostic related lncRNAs among tumor and normal tissue. *:P < 0.05; **:P < 0.01; ***:P < 0.001. Red represents high expression, while blue represents low expression. The abscissa represents the sample, while the ordinate represents prognostic related lncRNA. The Role of m6A-related lncRNAs Based on the similarity demonstrated by the expression levels of m6A-related lncRNAs and the proportion of ambiguous clustering measure, the k = 2 was identified as having optimal clustering stability from k = 2 to 9. A total of 433 patients with BLCA were clustered into two subtypes, namely, cluster1 (n = 132) and cluster2 (n = 271), based on the expression levels of the m6A-related lncRNAs (Fig. 3 , table S3). To assess the potential role of m6A-related lncRNAs, a survival study based on lncRNA subtypes was conducted, and the overall survival (OS) of cluster1 was higher than that of cluster2 (p = 0. 007), as shown in Fig. 4 A. The expression of most m6A-related lncRNAs was lower in the cluster2 than in the cluster1, especially the expression levels of LINC0260、AC104532.2、AL022328.2 and EHMT2 − AS1 (Fig. 4 B). While the expression of individual m6A-related lncRNAs was higher in the cluster2 than in the cluster1, such as AC097359.2 and AC025280.1 (Fig. 4 B). Additionally, The clinicopathological features between the two subtypes were then compared (Fig. 4 B). The cluster2 was preferentially associated with Stage III − IV (p < 0. 001). Association Of Pd-l1 With M6a-related Lncrnas To explore the correlation of PD-L1 with m6A-related lncRNAs, the differential expression of PD-L1 in two subtypes were assessed. The expression level of PD-L1 in the cluster2 was distinctly higher than that in the cluster1 (p<0.01; Fig.5A). However, there was no distinct difference in the expression of prognosis-related lncRNAs in the different tissues (Fig.5B), probably due to small sample size. In BLCA, a gene correlation study was performed to see if there was a link between the target gene and the prognostic m6A-related lncRNAs (Fig.6), and we found that the expression of PD-L1 had a significantly negative association with m6A-lncRNA LINC0260、AC104532. 2、AL022328.2 and EHMT2−AS1. Figure 5 . Difference expression of target gene in related subtype. The expression level of PD-L1 was upregulated in the cluster2 compared with cluster1(p 0. 05). The Role of Immune cells Infiltration and the Tumor Microenvironment Differential analyses of immune cell infiltration and immunological scores in distinct clusters were carried out to evaluate the effect of m6A-related lncRNAs on the TIME of BLCA (table S4). The infiltration abundance of 22 types of immune cells in each cluster was shown in Figure 7A. Additionally, cluster2 exhibits high levels of CD4 memory-activated T cells, macrophages M2, neutrophils and NK cells resting (P < 0. 05; Figure 7B-7E), whereas cluster1 was more correlated with plasma cells (P < 0. 05; Figure7F) and T cells regulatory (Tregs) (P < 0. 001; Figure 7A and 7G). As shown in Figure 8A-8C, All of the scores were higher in cluster2, indicating a higher level of immune-related cells in the tumor microenvironment (P < 0. 05) (table S5). We applied gene set enrichment analysis (GSEA) to explore the probable regulatory mechanisms underlying the variations in TIME between the two groupings. The hallmark pathways and functions involved in cluster2 included Chemokine Signaling Pathway、 Cytokine-Cytokine Receptor Interaction、Cell Adhesion Molecules (CAMs)、Natural Killer Cell Mediated Cytotoxicity、Antigen Processing and Presentation and ECM-Receptor Interaction(Figure 9A-9F). According to the results, the most enriched signaling pathway was “Chemokine signaling pathway”. These functions and pathways were up-regulated in cluster2, Both the FDR q-value and the FWER p-values were < 0. 05. These functions and pathways might be implicated in the distinct TIME of cluster1/2. Figure 9 . Part of the GSEA results are listed above. high risk of m6A-related lncRNAs were enriched in multiple cancer-related functions and pathways, which were upregulated in class C2. Both FDR q-value and FWER p-Value < 0.01. Construction Of Prognostic Signatures Based On M6a-related Lncrnas A train group (70%) and a test group (30%) were utilized for the Kaplan-Meier survival analysis and associated plots to accurately predict the clinical outcome of m6A-related lncRNAs in BLCA patients, as shown in Fig. 10 A - 10 B and Table S6 and S7. The results showed that the OS of patients in the high-risk group was significantly lower than that in the low-risk group, as shown in Fig. 10 C- 10 D. The corresponding ROC curve was acquired using the time-ROC package to test the accuracy of our model for predicting the survival of patients with the disease (Fig. 10 E-F). ROC curve analysis showed that the AUC value of test group was 0. 666, while the AUC value of train group was 0. 680, indicating the considerable accuracy of our model for predicting the survival of patients with BLCA. Univariate and multivariate Cox analyses were performed to evaluate whether risk score was an independent prognostic factor for BLCA patients. The distribution of the risk scores, OS, OS status, and expression profiles of the eleven m6A-related lncRNAs in TCGA training and validation cohorts was displayed in Fig. 11A-F. which revealed that the higher the immune scores, the higher the risk. Both in the test group and train group, after univariate analysis obtained factors related to OS in BLCA patients, multivariate Cox regression analysis showed that risk score, stage and age were identified as independent prognostic factors (P 1), in Fig. 12 A-D. The link between risk score and clinical features, as well as cluster subgroups, was studied further. In the high-risk and low-risk groups, the heatmap showed the expression levels of eleven m6A-related lncRNAs(Fig. 13 ). AL136295.2 、BDNF − AS 、AC073575. 4 、AC074117. 1 、AC073534. 2 、AL022322. 1 、AC004076. 2 、AC104564. 3 、TMEM147 − AS1 、AC104532. 2 、AC116914. 2 were highly expressed in cluster1 and low risk group, while ATP1B3 − AS1 and AC097359. 2 were highly expressed in cluster2 and high-risk group. Furthermore, the different risk scores among subtype(P < 0.001), immune scores(P < 0.001), stage(P < 0.001) and N stage(P < 0.05) were examined in this heatmap. The results show that the risk score of cluster2、higher immune scores 、stage III − IV and N1–N3 is significantly higher. These findings indicate that the risk score of BLCA patients may profoundly influence clinical outcomes. To evaluate and verify if our model could be applied to diverse clinical parameters, we utilized model validation for clinical groups, as shown in Fig. 14 . The result indicated that our model could be applied to the following different clinical parameters: age, gender, lymph node metastasis, stage and T stage (P < 0.05). In our BLCA model, genetic differential analysis was used to examine the expression differences of target genes in the different risk groups, as shown in Fig. 15 . The expression level of PD-L1 was higher in high-risk group(P < 0. 001). The influence of eleven m6A-related lncRNAs on the BLCA immune microenvironment was estimated using the relationship between the risk score and immune infiltration levels(Fig. 16). The risk score was positively correlated with the infiltration levels of Macrophages M0, Macrophages M1, Macrophages M2, Mast cells activated, Neutrophils and NK resting cells with R > 0 and P < 0. 05. A significantly negative correlation was observed between the risk score and infiltration levels of Plasma cells, T cells gamma delta and T cells regulatory (Tregs) with R < 0 and P < 0.05. Risk signatures based on m6A-related lncRNAs are likely to have a role in BLCA's immune microenvironment regulation. Figure 11. Distribution of risk score, OS, and OS status of the 11 prognostic m6A-related lncRNAs in the TCGA test cohort (A+C+E) and TCGA train cohort (B+D+F) Figure 16 . Relationships between the Risk Score and Infiltration Abundances of ten Immune Cell Types. (A) Macrophages M0, (B) Macrophages M1, (C) Macrophages M2, (D) Mast cells activated (E) Neutrophils and (F) NK resting cells are positive related with risk score, R > 0 and P < 0.05. (G) Plasma cells, (H) T cells gamma delta and (I) T regulatory cells, R < 0 and P < 0.05. Discussion If not treated effectively, BLCA is a complicated disease with high morbidity and fatality rates [ 35 ]. M6A modification regulates mRNA splicing, stability, nuclear export, and translation[ 36 ]. Furthermore, The alteration of m6A has an impact on lncRNA processing, as well as the processes that control cell proliferation and maturation[ 37 ]. LncRNA-based signatures have been shown to accurately predict BLCA survival and recurrence. [ 38 , 39 ]. For instance, NKILA lncRNA boots tumor immune evasion by sensitizing T cells to activation-induced cell death[ 40 ]. Several studies have found that using m6A-related lncRNAs to build prognostic models was effective in predicting tumor prognosis [ 43 – 45 ]. In the current study, co-expression analysis was carried out to identify the correlation between m6A-related gene expression and lncRNAs. This was followed by univariate Cox regression analysis for the identification of m6A-related prognostic lncRNAs. As a result, 30 m6A-associated lncRNAs were found to be strongly related to overall survival outcomes. Some m6A-related lncRNAs are highly expressed in tumor tissues, while others are highly expressed in normal tissues(P < 0.05). Therefore, the expression, prognostic value, and immune significance of m6A-related lncRNAs in BLCA were researched which could guide future research. Next, two molecular subtypes (clusters 1/2) were identified by consensus clustering of 30 m6A-related lncRNAs. The prognosis and different clinicopathological aspects of BLCA were determined by the cluster1/2 subtype, which was directly associated to PD-L1, immunological scores, and immune cell infiltration levels. Patients in cluster2 showed a worse survival rate than those in cluster1, which could be due to the greater importance of immunological scores and PD-L1 expression levels in cluster2. These findings backed up the findings of a previous study, which found that bladder cancer patients with high immunological and stromal scores had a lower overall survival rate[ 41 ]. We also discovered that PD-L1 expression was considerably higher in cluster2 than in cluster1. In addition, it has been observed that PD-L1 is significantly negatively correlated with LINC0260、AC104532. 2、AL022328.2 and EHMT2 − AS1. Further research is required to see whether the regulators can estimate the effectiveness of immunotherapy in BLCA patients. Increasing evidence also suggests that targeting lncRNA UCA1 with CRISPR-Cas9, which blocks PD-1 function, can improve antitumor activity in BLCA patients[ 42 ] . However, research on m6A alteration of lncRNAs is currently small in number. Thus, there is an urgent need for further research on lncRNA m6A modification and recognition corroborate our findings. Moreover, There was a substantial variation in TIME between cluster1/2 subtypes. The results revealed that CD4 memory-activated T cells, macrophages M2, neutrophils, and NK cells resting in cluster2 were greater than those resting in cluster1, but plasma cells and T cells regulatory (Tregs) were strongly clustered in cluster1. Furthermore, in agreement with our findings, tumor-associated macrophages are critical immunosuppressive cells driving tumorigenesis and metastasis[ 43 , 44 ]. Macrophages that infiltrate the tumor microenvironment are usually referred to as tumor-associated macrophages (TAMs) [ 45 ]. TAMs have been shown to promote tumor cell proliferation, invasion, and metastasis, strengthen angiogenesis, and suppress the anti-tumor response[ 46 ]. It has been observed that BLCA patients with larger numbers of macrophages, particularly macrophages M2, had poorer prognoses [ 47 , 48 ]. Next, we carried out GSEA. “CHEMOKINE SIGNALING PATHWAY” was the most significantly enriched signaling pathway. CCL2, CCL3, CCL4, and CCL5 are CC chemokines that are well-known chemotactic elements for macrophage populations in tumors [ 49 , 50 ]. CCL2 stimulation shifts human peripheral blood CD11b + cells toward a CD206 + M2-polarized phenotype according to Roca et al [ 51 ]. Tripathi et al. also discovered that oncostatin M and the chemokine CCL11/eotaxin, both generated by hypoxic cancer cells, shifted macrophages toward an M2 phenotype [ 52 , 53 ]. Additionaly ,CCL17 and CCL22 can promote Treg migration by interacting with the CCR4 receptor[ 54 ]. Thus, attraction of immunosuppressive immune cells through chemokine production is one of the pro-tumoral characteristics of TAMs. Taken together, our findings demonstrate that m6A methylation of lncRNA plays a crucial role in shaping the TIME for immune evasion, providing novel insights for effective cancer immunotherapy. Lasso regression was used to establish a statistical model for m6A-related lncRNAs. The OS of the patients in high-risk group was shorter than that of the patients in low-risk group in the TCGA test and train cohorts. Cluster2 had considerably higher risk ratings, immunological scores, and PD-L1 levels, according to further research. This observation was consistent with previous results. previous research have found that PD-L1 can bind to PD-1 on T cells, B cells, and macrophages activated on the surface of tumor cells, thus showing immunosuppressive effects[ 55 ]. A high level of PD-L1 expression has been linked to malignancy and a poor prognosis in BLCA patients, and such individuals have a greater rate of recurrence following surgery[ 56 – 58 ]. The results further confirming that the PD-L1 gene may be an oncogene of bladder cancer. Furthermore, this conducted model could be applied to different clinical groups. The risk score was found to be an independent prognostic factor for BLCA patients in both univariate and multivariate Cox regression models. Further experimental studies is necessary to determine the regulation of these lncRNAs by m6A modification and to clarify the associated processes in bladder cancer progression and immune evasion. Studies with larger sample sizes are also needed to confirm the prognostic value of m6A-related lncRNA risk score. Conclusion In conclusion, this study thoroughly identified the expression of m6A-related lncRNAs in BLCA, their correlation with PD-L1, effects on the TIME and potential regulatory mechanisms. Two BLCA subtypes (cluster1/2) were identified, each of which classified the prognosis of BLCA patients while also presenting considerably distinct TIME. By modifying TIME and PD-L1 expression, m6A-related lncRNAs may promote the responsiveness of BLCA patients to immunotherapy. GSEA further identified that the m6A-related lncRNAs might be involved in the regulation of BLCA immune microenvironment in synergy with the “CHEMOKINE SIGNALING PATHWAY”. More importantly, we developed and validated an prognostic model based on m6A-related lncRNA risk score with robust prognostic value and the ability to predict response to immunotherapy in patients with bladder cancer. The risk score is highly correlated with the malignant clinicopathological features of BLCA, which enhanced the understanding of m6A-related lncRNAs in TIME cells infiltration and immune evasion, providing novel insights for guiding more effective immunotherapy strategies in BLCA. Methods Ethics statement Genes were investigated exclusively using sequences available in public databases. This research did not include human or animal subjects. The all study were performed in accordance with the relevant guidelines and regulations and adhered to the Declaration of Helsinki. Sample Data Acquisition and Collation Transcriptome RNA sequencing data of 433 colorectal cancer samples were gained from the TCGA website (https://portal.gdc.cancer. gov/)[59], which contains 414 tumor samples and 19 adjacent normal samples. Then, utilising PERL software (https://www.perl.org/), we created the mRNA matrix and used the associated script to organize transcriptomic data and transform gene IDs. Finally, the clinical statistical analyses were performed using the same software and a specific script. Identification of m6A Related lncRNAs To distinguish mRNAs and lncRNAs by employing the collated transcriptome data., wee stablished a gene expression matrix and human conFigure uration file including the expression levels of m6A-related lncRNAs gene profiles. Then, The gene IDs were then converted to gene names using data from the Ensembl database (http://asia.ensembl.org/info/data/index. html). Next, co-expression analysis was carried out to illustrate the connection between m6A-related gene expression and lncRNAs and a network plot was depicted via the igraph package to visualize the correlation. Furthermore, prognostic-related lncRNAs were obtained, and the survival package estimated the confidence interval and hazard ratio. Forest plots were used to visualize the results of the univariate Cox regression analysis. The differential expression data of prognosis-related m6A-related lncRNAs between tumor and normal tissues were acquired through the limma package, pheatmap package, reshape2 package and ggpubr package in R software. Differences with P < 0. 05 were deemed statistically meaningful. To visualize the differences of expression, heatmaps and boxplots were drafted. The Role of m6A-related lncRNAs We used the "ConsensusClusterPlus" tool to examine the biological features of m6A regulators in BLCA (http://www. bioconductor. org/, 1000 iterations and resampling rate of 80%) to categorize BLCA patients from the TCGA database into different groups. Moreover, Survminer and survival programs were used to do a survival analysis based on lncRNA subtypes. A heatmap was depicted to demonstrate the correlation between lncRNAs and clinicopathological indicators. the standard names of the target genes was identified via NCBI (https://www.ncbi.nlm.nih. gov/). Differences with P < 0. 05 were deemed statistically meaningful. The Role of Immune cells Infiltration and the Tumor Microenvironment CIBERSORT (http://cibersort.stanford.edu/), a deconvolution algorithm based on gene expression was computed the composition of immune cells[60]. Based on the expression profile data of BLCA in the TCGA database, the CIBERSORT software was used to calculate the infiltration level of 22 immune cells. Subsequently, was used the ESTIMATE algorithm to calculate the immune scores of each patient, and the different immune scores between the two cluster subgroups was evaluated[61]. The biological processes involved in distinct subgroups were studied using gene set enrichment analysis (GSEA). GSEA hallmarks were used to identify predetermined gene sets, and p-values were determined via 5000 permutations according to the gene set. As referenced in the Results section, a pathway or function with a p-value < 0. 05 and a false discovery rate (FDR) < 0.05 was considered relevant. m6A-related lncRNAs Prognostic Model The prognostic signature of m6A-related lncRNAs were identified by Lasso regression analysis. The following computational equation was utilized to calculate the coefficients for each bladder cancer case: risk score = sum of coefficients × the lncRNA expression. In both the training and test cohorts, the risk score of all BlCA patients was calculated. The patients were then divided into high-and low-risk groups, with the median value of the risk score serving as the cutoff point. Moreover, the correlation between the risk score and the abundance of immune cells was also calculated and “ggplot2” was used to visualize the result. Declarations Author contributions statement X.M.Q designed the research, conducted statistical analysis and drafted the manuscript. W.Y.L. designed the research and reviewed the article. W.J.C. and W.X.D reviewed the results. Z.Y.Q and W.X.J reviewed the article and data analysis. All authors reviewed the manuscript. Availability of Data and Materials Publicly available datasets were analyzed in this study. This data can be found here: TCGA database (http://www.cancer.gov/tcga). The raw data and relative software/code can be accessed in a supplementary file. The data generated during the current study are available from the corresponding author on reasonable request. Competing interests The author(s) declare no competing interests. References Antoni, S., et al., Bladder Cancer Incidence and Mortality: A Global Overview and Recent Trends. Eur Urol, 2017. 71 (1): p. 96-108. Ferlay, J., et al., Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer, 2015. 136 (5): p. E359-86. Witjes, J.A., et al., Updated 2016 EAU Guidelines on Muscle-invasive and Metastatic Bladder Cancer. European Urology, 2016: p. S0302283816302901. Cambier, S., et al., EORTC Nomograms and Risk Groups for Predicting Recurrence, Progression, and Disease-specific and Overall Survival in Non–Muscle-invasive Stage Ta–T1 Urothelial Bladder Cancer Patients Treated with 1–3 Years of Maintenance Bacillus Calmette-Guérin. European Urology, 2016: p. 60-69. Szabados, B., et al., Response Rate to Chemotherapy After Immune Checkpoint Inhibition in Metastatic Urothelial Cancer. European Urology, 2017: p. 149. Akhavan-Niaki, H. and A.A. Samadani, DNA methylation and cancer development: molecular mechanism. Cell Biochem Biophys, 2013. 67 (2): p. 501-13. Roignant, J.Y. and M. Soller, m(6)A in mRNA: An Ancient Mechanism for Fine-Tuning Gene Expression. Trends Genet, 2017. 33 (6): p. 380-390. Morgan, M.A.J. and A. Shilatifard, Reevaluating the roles of histone-modifying enzymes and their associated chromatin modifications in transcriptional regulation. Nat Genet, 2020. 52 (12): p. 1271-1281. Litt, M.D., et al., Correlation between histone lysine methylation and developmental changes at the chicken beta-globin locus. Science, 2001. 293 (5539): p. 2453-5. Noma, K., C.D. Allis, and S.I. Grewal, Transitions in distinct histone H3 methylation patterns at the heterochromatin domain boundaries. Science, 2001. 293 (5532): p. 1150-5. Yao, R.W., Y. Wang, and L.L. Chen, Cellular functions of long noncoding RNAs. Nature Cell Biology, 2019. 21 (5): p. 542. Di Martino, M.T., et al., miRNAs and lncRNAs as Novel Therapeutic Targets to Improve Cancer Immunotherapy. Cancers (Basel), 2021. 13 (7). Denaro, N., M.C. Merlano, and C.L. Nigro, Long noncoding RNAs as regulators of cancer immunity. Molecular Oncology, 2019. 13 (1). Huang, D., et al., NKILA lncRNA promotes tumor immune evasion by sensitizing T cells to activation-induced cell death. Nature Immunology, 2018. Shuai, et al., Synergistic Antitumor Effect on Bladder Cancer by Rational Combination of Programmed Cell Death 1 Blockade and CRISPR-Cas9-Mediated Long Non-Coding RNA Urothelial Carcinoma Associated 1 Knockout. Human Gene Therapy, 2018. Fazi, F. and A. Fatica, Interplay Between N6-Methyladenosine (m6A) and Non-coding RNAs in Cell Development and Cancer. Frontiers in Cell and Developmental Biology, 2019. 7 : p. 116. Dai, F., et al., Crosstalk between RNA m6A Modification and Non-coding RNA Contributes to Cancer Growth and Progression - ScienceDirect. 2020. He, R.Z., J. Jiang, and D.X. Luo, The functions of N6-methyladenosine modification in lncRNAs. Genes & Diseases, 2020. Yi, Y.C., et al., Novel insights into the interplay between m6A modification and noncoding RNAs in cancer. Molecular Cancer, 2020. 19 (1). Du, K., et al., m(6)A RNA Methylation Controls Neural Development and Is Involved in Human Diseases. Mol Neurobiol, 2019. 56 (3): p. 1596-1606. Wang, S., et al., Roles of RNA methylation by means of N(6)-methyladenosine (m(6)A) in human cancers. Cancer Lett, 2017. 408 : p. 112-120. Han, X., et al., RNA methylations in human cancers. Seminars in Cancer Biology, 2020(70–81). Cancer stem cells in progression of colorectal cancer. Oncotarget, 2017. 9 (70). Zhang, L., et al., Immune Landscape of Colorectal Cancer Tumor Microenvironment from Different Primary Tumor Location. Frontiers in Immunology, 2018. 9 : p. 1578-. Jin-Zhao, et al., METTL14 suppresses the metastatic potential of hepatocellular carcinoma by modulating N6 -methyladenosine-dependent primary MicroRNA processing. Hepatology (Baltimore, Md.), 2017. Barros-Silva, D., et al., VIRMA-dependent N6-methyladenosine modifications regulate the expression of long non-coding RNAs CCAT1 and CCAT2 in prostate cancer. Cancers, 2020. 12 (4): p. 771. Li, Z., et al., Identification of three m6A‐related mRNAs signature and risk score for the prognostication of hepatocellular carcinoma. Cancer Medicine, 2020. 9 . Meng, Z., et al., The m6A-Related mRNA Signature Predicts the Prognosis of Pancreatic Cancer Patients. Molecular Therapy — Oncolytics, 2020. 17 . Jin, Y., et al., Analysis of m6A-Related Signatures in the Tumor Immune Microenvironment and Identification of Clinical Prognostic Regulators in Adrenocortical Carcinoma. Frontiers in Immunology, 2021. 12 : p. 637933. Shen, J., et al., N-methyladenosine reader YTHDF2-mediated long noncoding RNA FENDRR degradation promotes cell proliferation in endometrioid endometrial carcinoma. Lab Invest, 2021. 101 (6): p. 775-784. Wang, J., et al., RP11-323N12.5 promotes the malignancy and immunosuppression of human gastric cancer by increasing YAP1 transcription. Gastric Cancer, 2021. 24 (1): p. 85-102. Sanmamed, M.F. and L. Chen, A Paradigm Shift in Cancer Immunotherapy: From Enhancement to Normalization. Cell, 2018. 175 (2): p. 313-326. Kim, J.W., et al., Emerging immunotherapies for bladder cancer. Curr Opin Oncol, 2015. 27 (3): p. 191-200. Yi, L., et al., Comprehensive Analysis of the PD-L1 and Immune Infiltrates of m(6)A RNA Methylation Regulators in Head and Neck Squamous Cell Carcinoma. Mol Ther Nucleic Acids, 2020. 21 : p. 299-314. Kamat, A.M., et al., Bladder cancer. Lancet, 2016. 388 (10061): p. 2796-2810. Lan, Q., et al., The Critical Role of RNA m(6)A Methylation in Cancer. Cancer Res, 2019. 79 (7): p. 1285-1292. Dai, D., et al., N6-methyladenosine links RNA metabolism to cancer progression. Cell Death & Disease, 2018. 9 (2): p. 124. He, A., et al., Prognostic value of long non-coding RNA signatures in bladder cancer. Aging (Albany NY), 2019. 11 (16): p. 6237-6251. Lian, P., et al., An eight-long non-coding RNA signature as a candidate prognostic biomarker for bladder cancer. Aging (Albany NY), 2019. 11 (17): p. 6930-6940. Huang, D., et al., NKILA lncRNA promotes tumor immune evasion by sensitizing T cells to activation-induced cell death. Nat Immunol, 2018. 19 (10): p. 1112-1125. Li, F., et al., Profiles of tumor-infiltrating immune cells and prognostic genes associated with the microenvironment of bladder cancer. Int Immunopharmacol, 2020. 85 : p. 106641. Zhen, S., et al., Synergistic Antitumor Effect on Bladder Cancer by Rational Combination of Programmed Cell Death 1 Blockade and CRISPR-Cas9-Mediated Long Non-Coding RNA Urothelial Carcinoma Associated 1 Knockout. Hum Gene Ther, 2018. 29 (12): p. 1352-1363. Poh, A.R. and M. Ernst, Targeting Macrophages in Cancer: From Bench to Bedside. Front Oncol, 2018. 8 : p. 49. Ostuni, R., et al., Macrophages and cancer: from mechanisms to therapeutic implications. Trends Immunol, 2015. 36 (4): p. 229-39. Noy, R. and J.W. Pollard, Tumor-associated macrophages: from mechanisms to therapy. Immunity, 2014. 41 (1): p. 49-61. Grivennikov, S.I., F.R. Greten, and M. Karin, Immunity, inflammation, and cancer. Cell, 2010. 140 (6): p. 883-99. Chung, I.H., et al., The long non-coding RNA LINC01013 enhances invasion of human anaplastic large-cell lymphoma. Sci Rep, 2017. 7 (1): p. 295. Xue, Y., et al., Tumor ‑infiltrating M2 macrophages driven by specific genomic alterations are associated with prognosis in bladder cancer. Oncol Rep, 2019. 42 (2): p. 581-594. De, I., et al., CSF1 Overexpression Promotes High-Grade Glioma Formation without Impacting the Polarization Status of Glioma-Associated Microglia and Macrophages. Cancer Res, 2016. 76 (9): p. 2552-60. Wu, Y., et al., CCL3-CCR5 axis regulates intratumoral accumulation of leukocytes and fibroblasts and promotes angiogenesis in murine lung metastasis process. J Immunol, 2008. 181 (9): p. 6384-93. Roca, H., et al., CCL2 and interleukin-6 promote survival of human CD11b+ peripheral blood mononuclear cells and induce M2-type macrophage polarization. J Biol Chem, 2009. 284 (49): p. 34342-54. Jose, P.J., et al., Eotaxin: a potent eosinophil chemoattractant cytokine detected in a guinea pig model of allergic airways inflammation. J Exp Med, 1994. 179 (3): p. 881-7. Tripathi, C., et al., Macrophages are recruited to hypoxic tumor areas and acquire a pro-angiogenic M2-polarized phenotype via hypoxic cancer cell derived cytokines Oncostatin M and Eotaxin. Oncotarget, 2014. 5 (14): p. 5350-68. Ishida, T. and R. Ueda, CCR4 as a novel molecular target for immunotherapy of cancer. Cancer Sci, 2006. 97 (11): p. 1139-46. Bellmunt, J., T. Powles, and N.J. Vogelzang, A review on the evolution of PD-1/PD-L1 immunotherapy for bladder cancer: The future is now. Cancer Treat Rev, 2017. 54 : p. 58-67. Nakanishi, J., et al., Overexpression of B7-H1 (PD-L1) significantly associates with tumor grade and postoperative prognosis in human urothelial cancers. Cancer Immunol Immunother, 2007. 56 (8): p. 1173-82. Inman, B.A., et al., PD-L1 (B7-H1) expression by urothelial carcinoma of the bladder and BCG-induced granulomata: associations with localized stage progression. Cancer, 2007. 109 (8): p. 1499-505. Wang, Y., et al., Immune-related signature predicts the prognosis and immunotherapy benefit in bladder cancer. Cancer Med, 2020. 9 (20): p. 7729-7741. Wang, Z., M.A. Jensen, and J.C. Zenklusen, A Practical Guide to The Cancer Genome Atlas (TCGA). Methods Mol Biol, 2016. 1418 : p. 111-41. Newman, A.M., et al., Robust enumeration of cell subsets from tissue expression profiles. Nat Methods, 2015. 12 (5): p. 453-7. Yoshihara, K., et al., Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun, 2013. 4 : p. 2612. Additional Declarations No competing interests reported. Supplementary Files rawdataandrelativesoftwareandcode.zip supplementarymaterials.xlsx Cite Share Download PDF Status: Under Review Version 1 posted Editorial decision: Major revision 06 Apr, 2022 Reviews received at journal 26 Mar, 2022 Reviewers agreed at journal 17 Mar, 2022 Reviewers agreed at journal 14 Mar, 2022 Reviewers invited by journal 14 Mar, 2022 Editor assigned by journal 14 Mar, 2022 Editor invited by journal 08 Mar, 2022 Submission checks completed at journal 08 Mar, 2022 First submitted to journal 09 Feb, 2022 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-1341357","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":89110955,"identity":"1b0f0dd9-8ca2-46e8-9bfb-ef63df451eff","order_by":0,"name":"MQ Xue","email":"","orcid":"","institution":"The First Affiliated Hospital of Zhengzhou University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"MQ","middleName":"","lastName":"Xue","suffix":""},{"id":89110956,"identity":"8ba58de9-899c-49e8-be41-bf97f22a1f47","order_by":1,"name":"YL Wang","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAAp0lEQVRIiWNgGAWjYBACA2YQWWDDw8/eQJIWgzQZyZ4DxGqBkIdtDG44EKnFnJ35mDSPwXkehhsMjB8+5hChxbKZLQ2o5TYP4+wGZsmZ24hx2GEes9sgLcwyB9iYeUnQco6HTSKBNC0HeHiI1gL0S/rPOQbJPBI8B5uJ84s5/+HDBm8q7Oztjzcf/PCRGC0gwMQDphgbiFQPUvuDeLWjYBSMglEwEgEATCwu6Cioo+QAAAAASUVORK5CYII=","orcid":"","institution":"The First Affiliated Hospital of Zhengzhou University","correspondingAuthor":true,"submittingAuthor":false,"prefix":"","firstName":"YL","middleName":"","lastName":"Wang","suffix":""},{"id":89110957,"identity":"9bfbd76a-2014-497a-abf6-8f6d97ad49e0","order_by":2,"name":"JC Wang","email":"","orcid":"","institution":"Henan Bioengineering Technology Research Center","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"JC","middleName":"","lastName":"Wang","suffix":""},{"id":89110958,"identity":"207b6609-e70b-4214-94c2-e2771edb5c03","order_by":3,"name":"XD Wang","email":"","orcid":"","institution":"Henan General Hospital","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"XD","middleName":"","lastName":"Wang","suffix":""},{"id":89110959,"identity":"00cc3d16-4dd9-4353-8c8e-ed04105a2c85","order_by":4,"name":"XJ Wang","email":"","orcid":"","institution":"Henan General Hospital","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"XJ","middleName":"","lastName":"Wang","suffix":""},{"id":89110960,"identity":"e57a2519-6eff-479c-9b23-cfa502cc4de7","order_by":5,"name":"YQ Zhang","email":"","orcid":"","institution":"Zhengzhou Technical College","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"YQ","middleName":"","lastName":"Zhang","suffix":""}],"badges":[],"createdAt":"2022-02-09 06:29:08","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-1341357/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-1341357/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":19099595,"identity":"1e76741f-bd2a-4304-96cb-fd923c65a531","added_by":"auto","created_at":"2022-03-10 21:08:08","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":364763,"visible":true,"origin":"","legend":"\u003cp\u003eThe expression of m6A-long noncoding RNAs (lncRNAs) and their function in the prognosis of BLCA patients. A: Network plot of correlation between m6A related gene expression and lncRNAs. B: forest plot of univariate cox regression analysis. The confidence interval and hazard ratio were determined with data from prognostic associated lncRNAs. Red represents high risk, while green represents low risk.\u003c/p\u003e","description":"","filename":"fig1.png","url":"https://assets-eu.researchsquare.com/files/rs-1341357/v1/4afa8fa7f537b1585ed81d07.png"},{"id":19099848,"identity":"1c4f68dd-3d3d-4ea0-ada1-5c628bbda904","added_by":"auto","created_at":"2022-03-10 21:11:08","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":580041,"visible":true,"origin":"","legend":"\u003cp\u003eDifferent expression of m6A prognostic related lncRNAs. A: Boxplot of the difference expression of m6A prognostic related lncRNAs among tumor and normal tissue. *:P \u0026lt; 0.05; **:P \u0026lt; 0.01; ***:P \u0026lt; 0.001. B: Heatmap of the difference expression of m6A prognostic related lncRNAs among tumor and normal tissue. *:P \u0026lt; 0.05; **:P \u0026lt; 0.01; ***:P \u0026lt; 0.001. Red represents high expression, while blue represents low expression. The abscissa represents the sample, while the ordinate represents prognostic related lncRNA.\u003c/p\u003e","description":"","filename":"fig2.png","url":"https://assets-eu.researchsquare.com/files/rs-1341357/v1/03937f23fad279d1e4f7c24e.png"},{"id":19100027,"identity":"b8ead546-1cc9-43ff-930e-09cb3e72a7c3","added_by":"auto","created_at":"2022-03-10 21:14:08","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":70670,"visible":true,"origin":"","legend":"\u003cp\u003eThe type of prognostic m6A-related lncRNAs. We divided lncRNAs into two categories according to their expression: cluster1 and cluster2. When K = 2, there was the least cross-mixing part between the two types and the CDF value was lowest, thus we classified them into two types: cluster1 and cluster2.\u003c/p\u003e","description":"","filename":"fig3.png","url":"https://assets-eu.researchsquare.com/files/rs-1341357/v1/7ed96fff9e8a91366ed488fd.png"},{"id":19100024,"identity":"b32c4f64-3ecf-41f3-a83d-a8606c7530ca","added_by":"auto","created_at":"2022-03-10 21:14:08","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":404669,"visible":true,"origin":"","legend":"\u003cp\u003ePrognostic associated m6A lncRNAs survival analyses and correlation with clinicopathological factors A: Survival analysis according to subtypes of lncRNAs, the survival rate of cluster1 was higher, P = 0.007. B: Heatmap of difference expression of prognostic related lncRNAs and relationship with clinicopathological parameters in different cluster. *: P \u0026lt; 0.05. Red represents high expression, while blue represents low expression. The abscissa represents the sample, while the ordinate represents prognostic related lncRNA.\u003c/p\u003e","description":"","filename":"fig4.png","url":"https://assets-eu.researchsquare.com/files/rs-1341357/v1/149ab3dc6e163c0d4db21b4b.png"},{"id":19099597,"identity":"3e2f6b9e-d2ff-4c1a-8711-e17877d7cc9d","added_by":"auto","created_at":"2022-03-10 21:08:08","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":112233,"visible":true,"origin":"","legend":"\u003cp\u003eDifference expression of target gene in related subtype. The expression level of PD-L1 was upregulated in the cluster2 compared with cluster1(p \u0026lt; 0. 001). There was not significant difference between BLCA tissues and para-carcinoma tissues (p \u0026gt; 0. 05).\u003c/p\u003e","description":"","filename":"fig5.png","url":"https://assets-eu.researchsquare.com/files/rs-1341357/v1/16a55ee697e6ae9ed7750093.png"},{"id":19099598,"identity":"8988a60d-319e-47e3-91fe-25d62ab9aaa2","added_by":"auto","created_at":"2022-03-10 21:08:08","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":389617,"visible":true,"origin":"","legend":"\u003cp\u003eAssociation of PD-L1 with m6A-related lncRNAs. Red means positive correlation, while blue means negative correlation, * means the difference is statistically significant (P \u0026lt; 0.05).\u003c/p\u003e","description":"","filename":"fig6.png","url":"https://assets-eu.researchsquare.com/files/rs-1341357/v1/d41c2fd941a574b191934eef.png"},{"id":19099605,"identity":"597a689d-6a17-4d8b-97ae-3dc7f91165c4","added_by":"auto","created_at":"2022-03-10 21:08:08","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":244071,"visible":true,"origin":"","legend":"\u003cp\u003eDifferent analysis of immune cell infiltration in different cluster. Vioplot: (A); Boxplot: (B) T cells CD4 memory activated, (C) macrophages M2, (D) Neutrophils, (E) NK cells resting, (F) plasma cells and (G) T cells regulatory (Tregs).\u003c/p\u003e","description":"","filename":"fig7.png","url":"https://assets-eu.researchsquare.com/files/rs-1341357/v1/8b2aeed176f2e254e0b3c52a.png"},{"id":19100151,"identity":"d025ee9d-4229-4573-b62a-06737270ae32","added_by":"auto","created_at":"2022-03-10 21:17:08","extension":"png","order_by":8,"title":"Figure 8","display":"","copyAsset":false,"role":"figure","size":62157,"visible":true,"origin":"","legend":"\u003cp\u003eImmune scores in the cluster1/2 subtypes. (A) ImmuneScore, (B) ESTIMATEScore, (C) StromalScore. All of the scores are higher in cluster1 (p \u0026lt; 0. 001).\u003c/p\u003e","description":"","filename":"fig8.png","url":"https://assets-eu.researchsquare.com/files/rs-1341357/v1/fe4780f5e567618f3c1ea9fa.png"},{"id":19100025,"identity":"6fc79e08-4f14-4280-8497-8e94b7277e9f","added_by":"auto","created_at":"2022-03-10 21:14:08","extension":"png","order_by":9,"title":"Figure 9","display":"","copyAsset":false,"role":"figure","size":201208,"visible":true,"origin":"","legend":"\u003cp\u003ePart of the GSEA results are listed above. high risk of m6A-related lncRNAs were enriched in multiple cancer-related functions and pathways, which were upregulated in class C2. Both FDR q-value and FWER p-Value \u0026lt; 0.01.\u003c/p\u003e","description":"","filename":"fig9.png","url":"https://assets-eu.researchsquare.com/files/rs-1341357/v1/433c9c53698173e05006d0aa.png"},{"id":19099856,"identity":"a63e5268-e5bd-45ec-86b8-f89bdd9b70a6","added_by":"auto","created_at":"2022-03-10 21:11:09","extension":"png","order_by":10,"title":"Figure 10","display":"","copyAsset":false,"role":"figure","size":285699,"visible":true,"origin":"","legend":"\u003cp\u003eConstruction of a novel prognostic risk signature for BLCA. (A + B): Lasso regression of 11 m6A-related prognostic lncRNAs was carried out to construct a prognostic model. (C + D): Overall survival analysis for patients in high/low risk (C represents test group, P \u0026lt; 0.01, while D represents train group, P \u0026lt; 0.001). (E + F): ROC curve to evaluate the accuracy of the predictive model (E represents test group, while F represents train group, AUC \u0026gt; 0.5).\u003c/p\u003e","description":"","filename":"fig10.png","url":"https://assets-eu.researchsquare.com/files/rs-1341357/v1/6f13bde44d503b761f831dc8.png"},{"id":19099857,"identity":"395f9195-ce56-422b-903d-e150d013035d","added_by":"auto","created_at":"2022-03-10 21:11:09","extension":"png","order_by":11,"title":"Figure 11","display":"","copyAsset":false,"role":"figure","size":349179,"visible":true,"origin":"","legend":"\u003cp\u003eDistribution of risk score, OS, and OS status of the 11 prognostic m6A-related lncRNAs in the TCGA test cohort (A+C+E) and TCGA train cohort (B+D+F)\u003c/p\u003e","description":"","filename":"fig11.png","url":"https://assets-eu.researchsquare.com/files/rs-1341357/v1/2900317163cd3e21467524e3.png"},{"id":19099602,"identity":"bd947f2b-5c56-4679-8fa2-52ceb68b8e10","added_by":"auto","created_at":"2022-03-10 21:08:08","extension":"png","order_by":12,"title":"Figure 12","display":"","copyAsset":false,"role":"figure","size":79228,"visible":true,"origin":"","legend":"\u003cp\u003eMultivariate and univariate analysis of independent prognostic analysis. (A + B represents test; C + D represents train, Stage, age and risk Score were risk factors for the prognosis of BLCA; A + C: multivariate, B + D: univariate analysis), P \u0026lt; 0.05.\u003c/p\u003e","description":"","filename":"fig12.png","url":"https://assets-eu.researchsquare.com/files/rs-1341357/v1/66b28db7ff0193fc7b9fd41f.png"},{"id":19100152,"identity":"d7ddb035-8e56-4b65-8fba-354824bce0da","added_by":"auto","created_at":"2022-03-10 21:17:09","extension":"png","order_by":13,"title":"Figure 13","display":"","copyAsset":false,"role":"figure","size":418077,"visible":true,"origin":"","legend":"\u003cp\u003eHeatmap and clinicopathologic features of high- and low-risk groups. ∗p \u0026lt; 0.05, ∗∗p \u0026lt; 0.01, and ∗∗∗p \u0026lt; 0.001.\u003c/p\u003e","description":"","filename":"fig13.png","url":"https://assets-eu.researchsquare.com/files/rs-1341357/v1/23b23b0b3eae79d175cb8f3c.png"},{"id":19099611,"identity":"74fcf29e-ab14-4a6e-af48-24ec1874c06b","added_by":"auto","created_at":"2022-03-10 21:08:09","extension":"png","order_by":14,"title":"Figure 14","display":"","copyAsset":false,"role":"figure","size":603728,"visible":true,"origin":"","legend":"\u003cp\u003eSurvival curves for model validation. The model could be applied to different clinical groups: age (A+B), gender (C+D), M stage (E+F), N stage (G+H), T stage (I+J), and stage (K+L), (P \u0026lt; 0.05).\u003c/p\u003e","description":"","filename":"fig14.png","url":"https://assets-eu.researchsquare.com/files/rs-1341357/v1/5960ce51d1fdacb050c9248a.png"},{"id":19099851,"identity":"61d37d2f-1f24-4197-83ba-11398ae42f11","added_by":"auto","created_at":"2022-03-10 21:11:08","extension":"png","order_by":15,"title":"Figure 15","display":"","copyAsset":false,"role":"figure","size":109340,"visible":true,"origin":"","legend":"\u003cp\u003eThe PD-L1 expression level by risk score group. PD-L1 expression was higher in high-risk group (P \u0026lt;0. 001).\u003c/p\u003e","description":"","filename":"fig15.png","url":"https://assets-eu.researchsquare.com/files/rs-1341357/v1/c301362dfd670b5a5bdafdfa.png"},{"id":19099853,"identity":"23f6d0f9-620a-4c90-b401-a27ebd94e3f7","added_by":"auto","created_at":"2022-03-10 21:11:08","extension":"png","order_by":16,"title":"Figure 16","display":"","copyAsset":false,"role":"figure","size":140287,"visible":true,"origin":"","legend":"\u003cp\u003eRelationships between the Risk Score and Infiltration Abundances of ten Immune Cell Types. (A) Macrophages M0, (B) Macrophages M1, (C) Macrophages M2, (D) Mast cells activated (E) Neutrophils and (F) NK resting cells are positive related with risk score, R \u0026gt; 0 and P \u0026lt; 0.05. (G) Plasma cells, (H) T cells gamma delta and (I) T regulatory cells, R \u0026lt; 0 and P \u0026lt; 0.05.\u003c/p\u003e","description":"","filename":"fig16.png","url":"https://assets-eu.researchsquare.com/files/rs-1341357/v1/91b03bb9e7d7a9259bcb5f0d.png"},{"id":19100153,"identity":"d278240f-60f7-4b51-a288-7442176f9df1","added_by":"auto","created_at":"2022-03-10 21:17:11","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":476143,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-1341357/v1/164f6730-14c7-4a88-a631-3f9519b4857f.pdf"},{"id":19099847,"identity":"35fab1f3-e0e2-4f4c-8659-4eaee4db2b56","added_by":"auto","created_at":"2022-03-10 21:11:08","extension":"zip","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":25,"visible":true,"origin":"","legend":"","description":"","filename":"rawdataandrelativesoftwareandcode.zip","url":"https://assets-eu.researchsquare.com/files/rs-1341357/v1/92c339527d1c68e46defdbe7.zip"},{"id":19099609,"identity":"96529e77-7e36-4f8e-a185-c882b2a011da","added_by":"auto","created_at":"2022-03-10 21:08:09","extension":"xlsx","order_by":3,"title":"","display":"","copyAsset":false,"role":"supplement","size":226405,"visible":true,"origin":"","legend":"","description":"","filename":"supplementarymaterials.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-1341357/v1/fe4155098bac9f3d0f5a882b.xlsx"}],"financialInterests":"No competing interests reported.","formattedTitle":"Comprehensive Analysis of the PD-L1 and Immune Infiltrates of N6-Methyladenosine Related Long Non-Coding RNAs in Bladder Cancer","fulltext":[{"header":"Introduction","content":"\u003cp\u003eBladder cancer is the ninth most frequent malignancy in the world. It is mainly represented by bladder urothelial cancer (BUC), accounting for more than 90% of bladder cancer, and smoking is considered the leading risk factor [\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e, \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e]. Furthermore, BLCA has a significant recurrence rate, and approximately half of patients relapsing following major surgery and developing metastases. [\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e, \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e]. Patients with metastatic urothelial carcinoma have reaped extraordinary gains from platinum-based chemotherapy and Immune Checkpoint Inhibitors (ICIs), but the diverse features of BLCA prompt to variable clinical outcomes for BLCA patients even though they received standard therapy. [\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e]. Bladder cancer develops through a complex process including aberrant genetic changes and epigenetic aberrations. Epigenetic anomalies can be seen on a range of levels, including DNA. [\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e], RNA[\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e]and histone modifications [\u003cspan additionalcitationids=\"CR9\" citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e]. Researchers must find novel biomarkers to better predict the prognosis and therapeutic response to BLCA in order to enhance survival and minimize the burden of BLCA patients.\u003c/p\u003e \u003cp\u003eLong noncoding RNAs (lncRNAs) are a class of noncoding RNA with a length of more than 200 nucleotides that have been implicated in post-transcriptional regulatory elements that regulate mRNA splicing, stability, and translation [\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e, \u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e]. LncRNAs may also play a key role in regulating genes that code for cancer-fighting proteins, according to growing evidence[\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e]. NKILA lncRNA, for example, aids tumor immune evasion by making T cells more susceptible to activation-induced cell death[\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e]. By generating an RNA-protein complex with enhancer of zeste homolog2(EZH2) and boosting the binding of EZH2 and H3K27me3 on the E-cadherin promoter region, lncRNA MRPL23-AS1 facilitated adenoid cystic carcinoma lung metastasis [\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e]. However, the potential role of immune-related lncRNA signatures as an effective therapeutic strategy in BLCA is unknown.\u003c/p\u003e \u003cp\u003eN6-methyladenosine (m6A) alteration has been shown to regulate lncRNAs in several research [\u003cspan additionalcitationids=\"CR17 CR18\" citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e]. M6A, the methylation modification at the sixth N atom of adenine, is the most frequent post-transcriptional alteration on mRNA, mediating\u0026thinsp;\u0026gt;\u0026thinsp;60% RNA methylation [\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e, \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e]. The m6A enzyme's abnormal expression affects tumor cell function and the tumor microenvironment(TIME)[\u003cspan additionalcitationids=\"CR23\" citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e]. Tumor stem cell self-renewal is accelerated by abrupt m6A alteration, which plays a vital role in tumorigenesis [\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e]. m6A enzymes and lncRNAs are both excellent diagnostic and prognostic indicators. Evidence is overwhelming indicating m6A-related mRNAs and lncRNAs could be served as promising targets for estimating prognosis in a multitude of malignancies [\u003cspan additionalcitationids=\"CR27 CR28\" citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e]. Degradation of m6A reader YTHDF2-modified lncRNA FENDER, for example, enhanced endometrioid endometrial tumorigenesis dramatically [\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e]. Wang [\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e]et al discovered that m6A-induced lncRNA RP11 expression triggers the malignant and immunosuppressive gastric cancer cells via upregulation of YAP1 expression. Therefore, Identification of novel and reliable prognostic molecular characteristics from various dimensions is critical for identifying the suitable therapy methods and improving the dismal prognosis in BLCA patients.\u003c/p\u003e \u003cp\u003eAs tumors develop, the body's immune system is triggered to resist tumor development [\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e]. Several immunological checkpoints, such as programmed death ligand-1 (PD-L1), have been identified as one of the immune escape mechanisms [\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e]. The TIME may influence a patient's response to immune checkpoint inhibitors, hence PD-L1 expression in BLCA should be taken into account while analyzing tumor immunity [\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e]. This study focus on systematically assessing the relationships of m6A-related lncRNAs with prognosis, programmed death ligand 1 (PD-L1), and TIME in BLCA.\u003c/p\u003e"},{"header":"Results","content":"\u003ch2\u003eIdentification Of M6a Related Lncrnas\u003c/h2\u003e\n\u003cp\u003e To distinguish mRNAs and lncRNAs, data on m6A-related gene expression were retrieved from transcriptome data. Next, The connection between m6A-related gene expression and lncRNAs was analyzed via co-expression analysis (table S1). A network plot was generated to visualize this correlation (Figure 1A). A forest plot was used to visualize the univariate Cox regression analysis (Figure 1B, table S2). heatmaps and box plots were constructed using the differential expression of prognosis-related m6A-related lncRNAs between tumor and normal tissues, as shown in Figure 2A and 2B. 30 m6A prognosis-related lncRNAs were identified that were differently expressed between tumor and surrounding normal tissues. LINC02604、AC116914.2、 Z84485.1、 ZNF32\u0026minus;AS2、 AC004076.2、 AL138921.1、 TMEM147\u0026minus;AS1、 SNHG20、 PTOV1\u0026minus;AS2 and AC004148.1 were dramatically higher in BLCA tissues than in normal adjacent tissues (p \u0026lt; 0. 001). The expression level of ATP1B3\u0026minus;AS1、AC025280.1、AC012568.1、BDNF\u0026minus;AS was markedly lower in BLCA tissues than in normal tissues (p \u0026lt; 0. 01).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cb\u003eFigure 2\u003c/b\u003e. Different expression of m6A prognostic related lncRNAs. A: Boxplot of the difference expression of m6A prognostic related lncRNAs among tumor and normal tissue. *:P\u0026thinsp;\u0026lt;\u0026thinsp;0.05; **:P\u0026thinsp;\u0026lt;\u0026thinsp;0.01; ***:P\u0026thinsp;\u0026lt;\u0026thinsp;0.001. B: Heatmap of the difference expression of m6A prognostic related lncRNAs among tumor and normal tissue. *:P\u0026thinsp;\u0026lt;\u0026thinsp;0.05; **:P\u0026thinsp;\u0026lt;\u0026thinsp;0.01; ***:P\u0026thinsp;\u0026lt;\u0026thinsp;0.001. Red represents high expression, while blue represents low expression. The abscissa represents the sample, while the ordinate represents prognostic related lncRNA.\u003c/p\u003e \u003cdiv id=\"Sec8\" class=\"Section2\"\u003e \u003ch2\u003eThe Role of m6A-related lncRNAs\u003c/h2\u003e \u003cp\u003eBased on the similarity demonstrated by the expression levels of m6A-related lncRNAs and the proportion of ambiguous clustering measure, the k\u0026thinsp;=\u0026thinsp;2 was identified as having optimal clustering stability from k\u0026thinsp;=\u0026thinsp;2 to 9. A total of 433 patients with BLCA were clustered into two subtypes, namely, cluster1 (n\u0026thinsp;=\u0026thinsp;132) and cluster2 (n\u0026thinsp;=\u0026thinsp;271), based on the expression levels of the m6A-related lncRNAs (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e3\u003c/span\u003e, table S3). To assess the potential role of m6A-related lncRNAs, a survival study based on lncRNA subtypes was conducted, and the overall survival (OS) of cluster1 was higher than that of cluster2 (p\u0026thinsp;=\u0026thinsp;0. 007), as shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e4\u003c/span\u003eA. The expression of most m6A-related lncRNAs was lower in the cluster2 than in the cluster1, especially the expression levels of LINC0260、AC104532.2、AL022328.2 and EHMT2\u0026thinsp;\u0026minus;\u0026thinsp;AS1 (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e4\u003c/span\u003eB). While the expression of individual m6A-related lncRNAs was higher in the cluster2 than in the cluster1, such as AC097359.2 and AC025280.1 (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e4\u003c/span\u003eB). Additionally, The clinicopathological features between the two subtypes were then compared (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e4\u003c/span\u003eB). The cluster2 was preferentially associated with Stage III\u0026thinsp;\u0026minus;\u0026thinsp;IV (p\u0026thinsp;\u0026lt;\u0026thinsp;0. 001).\u003c/p\u003e\u003c/div\u003e\n\u003ch2\u003eAssociation Of Pd-l1 With M6a-related Lncrnas\u003c/h2\u003e\n\u003cp\u003e To explore the correlation of PD-L1 with m6A-related lncRNAs, the differential expression of PD-L1 in two subtypes were assessed. The expression level of PD-L1 in the cluster2 was distinctly higher than that in the cluster1 (p\u0026lt;0.01; Fig.5A). However, there was no distinct difference in the expression of prognosis-related lncRNAs in the different tissues (Fig.5B), probably due to small sample size. In BLCA, a gene correlation study was performed to see if there was a link between the target gene and the prognostic m6A-related lncRNAs (Fig.6), and we found that the expression of PD-L1 had a significantly negative association with m6A-lncRNA LINC0260、AC104532. 2、AL022328.2 and EHMT2\u0026minus;AS1.\u003c/p\u003e \u003cp\u003e \u003cb\u003eFigure 5\u003c/b\u003e. Difference expression of target gene in related subtype. The expression level of PD-L1 was upregulated in the cluster2 compared with cluster1(p\u0026thinsp;\u0026lt;\u0026thinsp;0. 001). There was not significant difference between BLCA tissues and para-carcinoma tissues (p\u0026thinsp;\u0026gt;\u0026thinsp;0. 05).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cdiv id=\"Sec10\" class=\"Section2\"\u003e \u003ch2\u003eThe Role of Immune cells Infiltration and the Tumor Microenvironment\u003c/h2\u003e \u003cp\u003eDifferential analyses of immune cell infiltration and immunological scores in distinct clusters were carried out to evaluate the effect of m6A-related lncRNAs on the TIME of BLCA (table S4). The infiltration abundance of 22 types of immune cells in each cluster was shown in Figure 7A. Additionally, cluster2 exhibits high levels of CD4 memory-activated T cells, macrophages M2, neutrophils and NK cells resting (P \u0026lt; 0. 05; Figure 7B-7E), whereas cluster1 was more correlated with plasma cells (P \u0026lt; 0. 05; Figure7F) and T cells regulatory (Tregs) (P \u0026lt; 0. 001; Figure 7A and 7G). As shown in Figure 8A-8C, All of the scores were higher in cluster2, indicating a higher level of immune-related cells in the tumor microenvironment (P \u0026lt; 0. 05) (table S5). We applied gene set enrichment analysis (GSEA) to explore the probable regulatory mechanisms underlying the variations in TIME between the two groupings. The hallmark pathways and functions involved in cluster2 included Chemokine Signaling Pathway、 Cytokine-Cytokine Receptor Interaction、Cell Adhesion Molecules (CAMs)、Natural Killer Cell Mediated Cytotoxicity、Antigen Processing and Presentation and ECM-Receptor Interaction(Figure 9A-9F). According to the results, the most enriched signaling pathway was \u0026ldquo;Chemokine signaling pathway\u0026rdquo;. These functions and pathways were up-regulated in cluster2, Both the FDR q-value and the FWER p-values were \u0026lt; 0. 05. These functions and pathways might be implicated in the distinct TIME of cluster1/2.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cb\u003eFigure 9\u003c/b\u003e. Part of the GSEA results are listed above. high risk of m6A-related lncRNAs were enriched in multiple cancer-related functions and pathways, which were upregulated in class C2. Both FDR q-value and FWER p-Value\u0026thinsp;\u0026lt;\u0026thinsp;0.01.\u003c/p\u003e \u003c/div\u003e\n\u003ch2\u003eConstruction Of Prognostic Signatures Based On M6a-related Lncrnas\u003c/h2\u003e\n\u003cp\u003eA train group (70%) and a test group (30%) were utilized for the Kaplan-Meier survival analysis and associated plots to accurately predict the clinical outcome of m6A-related lncRNAs in BLCA patients, as shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e10\u003c/span\u003eA -\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e10\u003c/span\u003eB and Table S6 and S7. The results showed that the OS of patients in the high-risk group was significantly lower than that in the low-risk group, as shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e10\u003c/span\u003eC-\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e10\u003c/span\u003eD. The corresponding ROC curve was acquired using the time-ROC package to test the accuracy of our model for predicting the survival of patients with the disease (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e10\u003c/span\u003eE-F). ROC curve analysis showed that the AUC value of test group was 0. 666, while the AUC value of train group was 0. 680, indicating the considerable accuracy of our model for predicting the survival of patients with BLCA. Univariate and multivariate Cox analyses were performed to evaluate whether risk score was an independent prognostic factor for BLCA patients. The distribution of the risk scores, OS, OS status, and expression profiles of the eleven m6A-related lncRNAs in TCGA training and validation cohorts was displayed in Fig.\u0026nbsp;11A-F. which revealed that the higher the immune scores, the higher the risk. Both in the test group and train group, after univariate analysis obtained factors related to OS in BLCA patients, multivariate Cox regression analysis showed that risk score, stage and age were identified as independent prognostic factors (P\u0026thinsp;\u0026lt;\u0026thinsp;0. 01, HR\u0026thinsp;\u0026gt;\u0026thinsp;1), in Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e12\u003c/span\u003eA-D. The link between risk score and clinical features, as well as cluster subgroups, was studied further. In the high-risk and low-risk groups, the heatmap showed the expression levels of eleven m6A-related lncRNAs(Fig.\u0026nbsp;\u003cspan refid=\"Fig9\" class=\"InternalRef\"\u003e13\u003c/span\u003e). AL136295.2 、BDNF\u0026thinsp;\u0026minus;\u0026thinsp;AS 、AC073575. 4 、AC074117. 1 、AC073534. 2 、AL022322. 1 、AC004076. 2 、AC104564. 3 、TMEM147\u0026thinsp;\u0026minus;\u0026thinsp;AS1 、AC104532. 2 、AC116914. 2 were highly expressed in cluster1 and low risk group, while ATP1B3\u0026thinsp;\u0026minus;\u0026thinsp;AS1 and AC097359. 2 were highly expressed in cluster2 and high-risk group. Furthermore, the different risk scores among subtype(P\u0026thinsp;\u0026lt;\u0026thinsp;0.001), immune scores(P\u0026thinsp;\u0026lt;\u0026thinsp;0.001), stage(P\u0026thinsp;\u0026lt;\u0026thinsp;0.001) and N stage(P\u0026thinsp;\u0026lt;\u0026thinsp;0.05) were examined in this heatmap. The results show that the risk score of cluster2、higher immune scores 、stage III\u0026thinsp;\u0026minus;\u0026thinsp;IV and N1\u0026ndash;N3 is significantly higher. These findings indicate that the risk score of BLCA patients may profoundly influence clinical outcomes. To evaluate and verify if our model could be applied to diverse clinical parameters, we utilized model validation for clinical groups, as shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig10\" class=\"InternalRef\"\u003e14\u003c/span\u003e. The result indicated that our model could be applied to the following different clinical parameters: age, gender, lymph node metastasis, stage and T stage (P\u0026thinsp;\u0026lt;\u0026thinsp;0.05). In our BLCA model, genetic differential analysis was used to examine the expression differences of target genes in the different risk groups, as shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig11\" class=\"InternalRef\"\u003e15\u003c/span\u003e. The expression level of PD-L1 was higher in high-risk group(P\u0026thinsp;\u0026lt;\u0026thinsp;0. 001). The influence of eleven m6A-related lncRNAs on the BLCA immune microenvironment was estimated using the relationship between the risk score and immune infiltration levels(Fig.\u0026nbsp;16). The risk score was positively correlated with the infiltration levels of Macrophages M0, Macrophages M1, Macrophages M2, Mast cells activated, Neutrophils and NK resting cells with R\u0026thinsp;\u0026gt;\u0026thinsp;0 and P\u0026thinsp;\u0026lt;\u0026thinsp;0. 05. A significantly negative correlation was observed between the risk score and infiltration levels of Plasma cells, T cells gamma delta and T cells regulatory (Tregs) with R\u0026thinsp;\u0026lt;\u0026thinsp;0 and P\u0026thinsp;\u0026lt;\u0026thinsp;0.05. Risk signatures based on m6A-related lncRNAs are likely to have a role in BLCA's immune microenvironment regulation.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cb\u003eFigure\u003c/b\u003e 11. Distribution of risk score, OS, and OS status of the 11 prognostic m6A-related lncRNAs in the TCGA test cohort (A+C+E) and TCGA train cohort (B+D+F)\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cb\u003eFigure 16\u003c/b\u003e. Relationships between the Risk Score and Infiltration Abundances of ten Immune Cell Types. (A) Macrophages M0, (B) Macrophages M1, (C) Macrophages M2, (D) Mast cells activated (E) Neutrophils and (F) NK resting cells are positive related with risk score, R\u0026thinsp;\u0026gt;\u0026thinsp;0 and P\u0026thinsp;\u0026lt;\u0026thinsp;0.05. (G) Plasma cells, (H) T cells gamma delta and (I) T regulatory cells, R\u0026thinsp;\u0026lt;\u0026thinsp;0 and P\u0026thinsp;\u0026lt;\u0026thinsp;0.05.\u003c/p\u003e"},{"header":"Discussion","content":"\u003cp\u003eIf not treated effectively, BLCA is a complicated disease with high morbidity and fatality rates [\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e]. M6A modification regulates mRNA splicing, stability, nuclear export, and translation[\u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e]. Furthermore, The alteration of m6A has an impact on lncRNA processing, as well as the processes that control cell proliferation and maturation[\u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e]. LncRNA-based signatures have been shown to accurately predict BLCA survival and recurrence. [\u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e38\u003c/span\u003e, \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e39\u003c/span\u003e]. For instance, NKILA lncRNA boots tumor immune evasion by sensitizing T cells to activation-induced cell death[\u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e]. Several studies have found that using m6A-related lncRNAs to build prognostic models was effective in predicting tumor prognosis [\u003cspan additionalcitationids=\"CR44\" citationid=\"CR43\" class=\"CitationRef\"\u003e43\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e]. In the current study, co-expression analysis was carried out to identify the correlation between m6A-related gene expression and lncRNAs. This was followed by univariate Cox regression analysis for the identification of m6A-related prognostic lncRNAs. As a result, 30 m6A-associated lncRNAs were found to be strongly related to overall survival outcomes. Some m6A-related lncRNAs are highly expressed in tumor tissues, while others are highly expressed in normal tissues(P\u0026thinsp;\u0026lt;\u0026thinsp;0.05). Therefore, the expression, prognostic value, and immune significance of m6A-related lncRNAs in BLCA were researched which could guide future research.\u003c/p\u003e \u003cp\u003eNext, two molecular subtypes (clusters 1/2) were identified by consensus clustering of 30 m6A-related lncRNAs. The prognosis and different clinicopathological aspects of BLCA were determined by the cluster1/2 subtype, which was directly associated to PD-L1, immunological scores, and immune cell infiltration levels. Patients in cluster2 showed a worse survival rate than those in cluster1, which could be due to the greater importance of immunological scores and PD-L1 expression levels in cluster2. These findings backed up the findings of a previous study, which found that bladder cancer patients with high immunological and stromal scores had a lower overall survival rate[\u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e]. We also discovered that PD-L1 expression was considerably higher in cluster2 than in cluster1. In addition, it has been observed that PD-L1 is significantly negatively correlated with LINC0260、AC104532. 2、AL022328.2 and EHMT2\u0026thinsp;\u0026minus;\u0026thinsp;AS1. Further research is required to see whether the regulators can estimate the effectiveness of immunotherapy in BLCA patients. Increasing evidence also suggests that targeting lncRNA UCA1 with CRISPR-Cas9, which blocks PD-1 function, can improve antitumor activity in BLCA patients[\u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e] \u003cb\u003e.\u003c/b\u003e However, research on m6A alteration of lncRNAs is currently small in number. Thus, there is an urgent need for further research on lncRNA m6A modification and recognition corroborate our findings.\u003c/p\u003e \u003cp\u003eMoreover, There was a substantial variation in TIME between cluster1/2 subtypes. The results revealed that CD4 memory-activated T cells, macrophages M2, neutrophils, and NK cells resting in cluster2 were greater than those resting in cluster1, but plasma cells and T cells regulatory (Tregs) were strongly clustered in cluster1. Furthermore, in agreement with our findings, tumor-associated macrophages are critical immunosuppressive cells driving tumorigenesis and metastasis[\u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e43\u003c/span\u003e, \u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e]. Macrophages that infiltrate the tumor microenvironment are usually referred to as tumor-associated macrophages (TAMs) [\u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e]. TAMs have been shown to promote tumor cell proliferation, invasion, and metastasis, strengthen angiogenesis, and suppress the anti-tumor response[\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e]. It has been observed that BLCA patients with larger numbers of macrophages, particularly macrophages M2, had poorer prognoses [\u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e47\u003c/span\u003e, \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e]. Next, we carried out GSEA. \u0026ldquo;CHEMOKINE SIGNALING PATHWAY\u0026rdquo; was the most significantly enriched signaling pathway. CCL2, CCL3, CCL4, and CCL5 are CC chemokines that are well-known chemotactic elements for macrophage populations in tumors [\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e, \u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e50\u003c/span\u003e]. CCL2 stimulation shifts human peripheral blood CD11b\u0026thinsp;+\u0026thinsp;cells toward a CD206\u0026thinsp;+\u0026thinsp;M2-polarized phenotype according to Roca et al [\u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e51\u003c/span\u003e]. Tripathi et al. also discovered that oncostatin M and the chemokine CCL11/eotaxin, both generated by hypoxic cancer cells, shifted macrophages toward an M2 phenotype [\u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e52\u003c/span\u003e, \u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e53\u003c/span\u003e]. Additionaly ,CCL17 and CCL22 can promote Treg migration by interacting with the CCR4 receptor[\u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e]. Thus, attraction of immunosuppressive immune cells through chemokine production is one of the pro-tumoral characteristics of TAMs. Taken together, our findings demonstrate that m6A methylation of lncRNA plays a crucial role in shaping the TIME for immune evasion, providing novel insights for effective cancer immunotherapy.\u003c/p\u003e \u003cp\u003eLasso regression was used to establish a statistical model for m6A-related lncRNAs. The OS of the patients in high-risk group was shorter than that of the patients in low-risk group in the TCGA test and train cohorts. Cluster2 had considerably higher risk ratings, immunological scores, and PD-L1 levels, according to further research. This observation was consistent with previous results. previous research have found that PD-L1 can bind to PD-1 on T cells, B cells, and macrophages activated on the surface of tumor cells, thus showing immunosuppressive effects[\u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e55\u003c/span\u003e]. A high level of PD-L1 expression has been linked to malignancy and a poor prognosis in BLCA patients, and such individuals have a greater rate of recurrence following surgery[\u003cspan additionalcitationids=\"CR57\" citationid=\"CR56\" class=\"CitationRef\"\u003e56\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR58\" class=\"CitationRef\"\u003e58\u003c/span\u003e]. The results further confirming that the PD-L1 gene may be an oncogene of bladder cancer. Furthermore, this conducted model could be applied to different clinical groups. The risk score was found to be an independent prognostic factor for BLCA patients in both univariate and multivariate Cox regression models. Further experimental studies is necessary to determine the regulation of these lncRNAs by m6A modification and to clarify the associated processes in bladder cancer progression and immune evasion. Studies with larger sample sizes are also needed to confirm the prognostic value of m6A-related lncRNA risk score.\u003c/p\u003e"},{"header":"Conclusion","content":"\u003cp\u003eIn conclusion, this study thoroughly identified the expression of m6A-related lncRNAs in BLCA, their correlation with PD-L1, effects on the TIME and potential regulatory mechanisms. Two BLCA subtypes (cluster1/2) were identified, each of which classified the prognosis of BLCA patients while also presenting considerably distinct TIME. By modifying TIME and PD-L1 expression, m6A-related lncRNAs may promote the responsiveness of BLCA patients to immunotherapy. GSEA further identified that the m6A-related lncRNAs might be involved in the regulation of BLCA immune microenvironment in synergy with the \u0026ldquo;CHEMOKINE SIGNALING PATHWAY\u0026rdquo;. More importantly, we developed and validated an prognostic model based on m6A-related lncRNA risk score with robust prognostic value and the ability to predict response to immunotherapy in patients with bladder cancer. The risk score is highly correlated with the malignant clinicopathological features of BLCA, which enhanced the understanding of m6A-related lncRNAs in TIME cells infiltration and immune evasion, providing novel insights for guiding more effective immunotherapy strategies in BLCA.\u003c/p\u003e"},{"header":"Methods","content":"\u003cp\u003e\u003cstrong\u003eEthics statement\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eGenes were investigated exclusively using sequences available in public databases. This research did not include human or animal subjects. The all study were performed in accordance with the relevant guidelines and regulations and adhered to the Declaration of Helsinki.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eSample Data Acquisition and Collation\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTranscriptome RNA sequencing data of 433 colorectal cancer samples were gained from the TCGA website (https://portal.gdc.cancer. gov/)[59], which contains 414 tumor samples and 19 adjacent normal samples. Then, utilising PERL software (https://www.perl.org/), we created the mRNA matrix and used the associated script to organize transcriptomic data and transform gene IDs. Finally, the clinical statistical analyses were performed using the same software and a specific script.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eIdentification of m6A Related lncRNAs\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTo distinguish mRNAs and lncRNAs by employing the collated transcriptome data., wee stablished \u0026nbsp;a gene expression matrix and human conFigure uration file including the expression levels of m6A-related lncRNAs gene profiles. Then, The gene IDs were then converted to gene names using data from the Ensembl database (http://asia.ensembl.org/info/data/index. html). Next, co-expression analysis was carried out to illustrate the connection between m6A-related gene expression and lncRNAs and a network plot was depicted via the igraph package to visualize the correlation. Furthermore, prognostic-related lncRNAs were obtained, and the survival package estimated the confidence interval and hazard ratio. Forest plots were used to visualize the results of the univariate Cox regression analysis. The differential expression data of prognosis-related m6A-related lncRNAs between tumor and normal tissues were acquired through the limma package, pheatmap package, reshape2 package and ggpubr package in R software. Differences with P\u0026thinsp;\u0026lt;\u0026thinsp;0. 05 were deemed statistically meaningful. To visualize the differences of expression, heatmaps and boxplots were drafted.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eThe Role of m6A-related lncRNAs\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWe used the \u0026quot;ConsensusClusterPlus\u0026quot; tool to examine the biological features of m6A regulators in BLCA (http://www. bioconductor. org/, 1000 iterations and resampling rate of 80%) to categorize BLCA patients from the TCGA database into different groups. Moreover, Survminer and survival programs were used to do a survival analysis based on lncRNA subtypes. A heatmap was depicted to demonstrate the correlation between lncRNAs and clinicopathological indicators. the standard names of the target genes was identified via NCBI (https://www.ncbi.nlm.nih. gov/). Differences with P\u0026thinsp;\u0026lt;\u0026thinsp;0. 05 were deemed statistically meaningful. \u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eThe Role of Immune cells Infiltration and the Tumor Microenvironment\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eCIBERSORT (http://cibersort.stanford.edu/), a deconvolution algorithm based on gene expression was computed the composition of immune cells[60]. Based on the expression profile data of BLCA in the TCGA database, the CIBERSORT software was used to calculate the infiltration level of 22 immune cells. Subsequently, was used the ESTIMATE algorithm to calculate the immune scores of each patient, and the different immune scores between the two cluster subgroups was evaluated[61]. The biological processes involved in distinct subgroups were studied using gene set enrichment analysis (GSEA). GSEA hallmarks were used to identify predetermined gene sets, and p-values were determined via\u0026nbsp;5000 permutations according to the gene set. As referenced in the Results section, a pathway\u0026nbsp;or function with a p-value \u0026lt;\u0026nbsp;0. 05 and a false discovery rate (FDR) \u0026lt;\u0026nbsp;0.05 was considered relevant.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003em6A-related lncRNAs Prognostic Model\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe prognostic signature of m6A-related lncRNAs were identified by Lasso regression analysis. The following computational equation was utilized to calculate the coefficients for each bladder cancer case: risk score\u0026thinsp;=\u0026thinsp;sum of coefficients\u0026thinsp;\u0026times;\u0026thinsp;the lncRNA expression. In both the training and test cohorts, the risk score of all BlCA patients was calculated. The patients were then divided into high-and low-risk groups, with the median value of the risk score serving as the cutoff point. Moreover, the correlation between the risk score and the abundance of immune cells was also calculated and \u0026ldquo;ggplot2\u0026rdquo; was used to visualize the result.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eAuthor contributions statement\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e\u0026nbsp;\u003c/strong\u003eX.M.Q designed the research, conducted statistical analysis and drafted the manuscript. W.Y.L. designed the research and reviewed the article. W.J.C. and W.X.D reviewed the results. Z.Y.Q and W.X.J reviewed the article and data analysis. All authors reviewed the manuscript.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAvailability of Data and Materials\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003ePublicly available datasets were analyzed in this study. This data can be found here: TCGA database (http://www.cancer.gov/tcga). The raw data and relative software/code can be accessed in a supplementary file. The data generated during the current study are available from the corresponding author on reasonable request.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe author(s) declare no competing interests.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eAntoni, S., et al., \u003cem\u003eBladder Cancer Incidence and Mortality: A Global Overview and Recent Trends.\u003c/em\u003e Eur Urol, 2017. \u003cstrong\u003e71\u003c/strong\u003e(1): p. 96-108.\u003c/li\u003e\n\u003cli\u003eFerlay, J., et al., \u003cem\u003eCancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012.\u003c/em\u003e Int J Cancer, 2015. \u003cstrong\u003e136\u003c/strong\u003e(5): p. E359-86.\u003c/li\u003e\n\u003cli\u003eWitjes, J.A., et al., \u003cem\u003eUpdated 2016 EAU Guidelines on Muscle-invasive and Metastatic Bladder Cancer.\u003c/em\u003e European Urology, 2016: p. S0302283816302901.\u003c/li\u003e\n\u003cli\u003eCambier, S., et al., \u003cem\u003eEORTC Nomograms and Risk Groups for Predicting Recurrence, Progression, and Disease-specific and Overall Survival in Non\u0026ndash;Muscle-invasive Stage Ta\u0026ndash;T1 Urothelial Bladder Cancer Patients Treated with 1\u0026ndash;3 Years of Maintenance Bacillus Calmette-Gu\u0026eacute;rin.\u003c/em\u003e European Urology, 2016: p. 60-69.\u003c/li\u003e\n\u003cli\u003eSzabados, B., et al., \u003cem\u003eResponse Rate to Chemotherapy After Immune Checkpoint Inhibition in Metastatic Urothelial Cancer.\u003c/em\u003e European Urology, 2017: p. 149.\u003c/li\u003e\n\u003cli\u003eAkhavan-Niaki, H. and A.A. Samadani, \u003cem\u003eDNA methylation and cancer development: molecular mechanism.\u003c/em\u003e Cell Biochem Biophys, 2013. \u003cstrong\u003e67\u003c/strong\u003e(2): p. 501-13.\u003c/li\u003e\n\u003cli\u003eRoignant, J.Y. and M. Soller, \u003cem\u003em(6)A in mRNA: An Ancient Mechanism for Fine-Tuning Gene Expression.\u003c/em\u003e Trends Genet, 2017. \u003cstrong\u003e33\u003c/strong\u003e(6): p. 380-390.\u003c/li\u003e\n\u003cli\u003eMorgan, M.A.J. and A. Shilatifard, \u003cem\u003eReevaluating the roles of histone-modifying enzymes and their associated chromatin modifications in transcriptional regulation.\u003c/em\u003e Nat Genet, 2020. \u003cstrong\u003e52\u003c/strong\u003e(12): p. 1271-1281.\u003c/li\u003e\n\u003cli\u003eLitt, M.D., et al., \u003cem\u003eCorrelation between histone lysine methylation and developmental changes at the chicken beta-globin locus.\u003c/em\u003e Science, 2001. \u003cstrong\u003e293\u003c/strong\u003e(5539): p. 2453-5.\u003c/li\u003e\n\u003cli\u003eNoma, K., C.D. Allis, and S.I. Grewal, \u003cem\u003eTransitions in distinct histone H3 methylation patterns at the heterochromatin domain boundaries.\u003c/em\u003e Science, 2001. \u003cstrong\u003e293\u003c/strong\u003e(5532): p. 1150-5.\u003c/li\u003e\n\u003cli\u003eYao, R.W., Y. Wang, and L.L. Chen, \u003cem\u003eCellular functions of long noncoding RNAs.\u003c/em\u003e Nature Cell Biology, 2019. \u003cstrong\u003e21\u003c/strong\u003e(5): p. 542.\u003c/li\u003e\n\u003cli\u003eDi Martino, M.T., et al., \u003cem\u003emiRNAs and lncRNAs as Novel Therapeutic Targets to Improve Cancer Immunotherapy.\u003c/em\u003e Cancers (Basel), 2021. \u003cstrong\u003e13\u003c/strong\u003e(7).\u003c/li\u003e\n\u003cli\u003eDenaro, N., M.C. Merlano, and C.L. Nigro, \u003cem\u003eLong noncoding RNAs as regulators of cancer immunity.\u003c/em\u003e Molecular Oncology, 2019. \u003cstrong\u003e13\u003c/strong\u003e(1).\u003c/li\u003e\n\u003cli\u003eHuang, D., et al., \u003cem\u003eNKILA lncRNA promotes tumor immune evasion by sensitizing T cells to activation-induced cell death.\u003c/em\u003e Nature Immunology, 2018.\u003c/li\u003e\n\u003cli\u003eShuai, et al., \u003cem\u003eSynergistic Antitumor Effect on Bladder Cancer by Rational Combination of Programmed Cell Death 1 Blockade and CRISPR-Cas9-Mediated Long Non-Coding RNA Urothelial Carcinoma Associated 1 Knockout.\u003c/em\u003e Human Gene Therapy, 2018.\u003c/li\u003e\n\u003cli\u003eFazi, F. and A. Fatica, \u003cem\u003eInterplay Between N6-Methyladenosine (m6A) and Non-coding RNAs in Cell Development and Cancer.\u003c/em\u003e Frontiers in Cell and Developmental Biology, 2019. \u003cstrong\u003e7\u003c/strong\u003e: p. 116.\u003c/li\u003e\n\u003cli\u003eDai, F., et al., \u003cem\u003eCrosstalk between RNA m6A Modification and Non-coding RNA Contributes to Cancer Growth and Progression - ScienceDirect.\u003c/em\u003e 2020.\u003c/li\u003e\n\u003cli\u003eHe, R.Z., J. Jiang, and D.X. Luo, \u003cem\u003eThe functions of N6-methyladenosine modification in lncRNAs.\u003c/em\u003e Genes \u0026amp; Diseases, 2020.\u003c/li\u003e\n\u003cli\u003eYi, Y.C., et al., \u003cem\u003eNovel insights into the interplay between m6A modification and noncoding RNAs in cancer.\u003c/em\u003e Molecular Cancer, 2020. \u003cstrong\u003e19\u003c/strong\u003e(1).\u003c/li\u003e\n\u003cli\u003eDu, K., et al., \u003cem\u003em(6)A RNA Methylation Controls Neural Development and Is Involved in Human Diseases.\u003c/em\u003e Mol Neurobiol, 2019. \u003cstrong\u003e56\u003c/strong\u003e(3): p. 1596-1606.\u003c/li\u003e\n\u003cli\u003eWang, S., et al., \u003cem\u003eRoles of RNA methylation by means of N(6)-methyladenosine (m(6)A) in human cancers.\u003c/em\u003e Cancer Lett, 2017. \u003cstrong\u003e408\u003c/strong\u003e: p. 112-120.\u003c/li\u003e\n\u003cli\u003eHan, X., et al., \u003cem\u003eRNA methylations in human cancers.\u003c/em\u003e Seminars in Cancer Biology, 2020(70\u0026ndash;81).\u003c/li\u003e\n\u003cli\u003e\u003cem\u003eCancer stem cells in progression of colorectal cancer.\u003c/em\u003e Oncotarget, 2017. \u003cstrong\u003e9\u003c/strong\u003e(70).\u003c/li\u003e\n\u003cli\u003eZhang, L., et al., \u003cem\u003eImmune Landscape of Colorectal Cancer Tumor Microenvironment from Different Primary Tumor Location.\u003c/em\u003e Frontiers in Immunology, 2018. \u003cstrong\u003e9\u003c/strong\u003e: p. 1578-.\u003c/li\u003e\n\u003cli\u003eJin-Zhao, et al., \u003cem\u003eMETTL14 suppresses the metastatic potential of hepatocellular carcinoma by modulating N6 -methyladenosine-dependent primary MicroRNA processing.\u003c/em\u003e Hepatology (Baltimore, Md.), 2017.\u003c/li\u003e\n\u003cli\u003eBarros-Silva, D., et al., \u003cem\u003eVIRMA-dependent N6-methyladenosine modifications regulate the expression of long non-coding RNAs CCAT1 and CCAT2 in prostate cancer.\u003c/em\u003e Cancers, 2020. \u003cstrong\u003e12\u003c/strong\u003e(4): p. 771.\u003c/li\u003e\n\u003cli\u003eLi, Z., et al., \u003cem\u003eIdentification of three m6A‐related mRNAs signature and risk score for the prognostication of hepatocellular carcinoma.\u003c/em\u003e Cancer Medicine, 2020. \u003cstrong\u003e9\u003c/strong\u003e.\u003c/li\u003e\n\u003cli\u003eMeng, Z., et al., \u003cem\u003eThe m6A-Related mRNA Signature Predicts the Prognosis of Pancreatic Cancer Patients.\u003c/em\u003e Molecular Therapy \u0026mdash; Oncolytics, 2020. \u003cstrong\u003e17\u003c/strong\u003e.\u003c/li\u003e\n\u003cli\u003eJin, Y., et al., \u003cem\u003eAnalysis of m6A-Related Signatures in the Tumor Immune Microenvironment and Identification of Clinical Prognostic Regulators in Adrenocortical Carcinoma.\u003c/em\u003e Frontiers in Immunology, 2021. \u003cstrong\u003e12\u003c/strong\u003e: p. 637933.\u003c/li\u003e\n\u003cli\u003eShen, J., et al., \u003cem\u003eN-methyladenosine reader YTHDF2-mediated long noncoding RNA FENDRR degradation promotes cell proliferation in endometrioid endometrial carcinoma.\u003c/em\u003e Lab Invest, 2021. \u003cstrong\u003e101\u003c/strong\u003e(6): p. 775-784.\u003c/li\u003e\n\u003cli\u003eWang, J., et al., \u003cem\u003eRP11-323N12.5 promotes the malignancy and immunosuppression of human gastric cancer by increasing YAP1 transcription.\u003c/em\u003e Gastric Cancer, 2021. \u003cstrong\u003e24\u003c/strong\u003e(1): p. 85-102.\u003c/li\u003e\n\u003cli\u003eSanmamed, M.F. and L. Chen, \u003cem\u003eA Paradigm Shift in Cancer Immunotherapy: From Enhancement to Normalization.\u003c/em\u003e Cell, 2018. \u003cstrong\u003e175\u003c/strong\u003e(2): p. 313-326.\u003c/li\u003e\n\u003cli\u003eKim, J.W., et al., \u003cem\u003eEmerging immunotherapies for bladder cancer.\u003c/em\u003e Curr Opin Oncol, 2015. \u003cstrong\u003e27\u003c/strong\u003e(3): p. 191-200.\u003c/li\u003e\n\u003cli\u003eYi, L., et al., \u003cem\u003eComprehensive Analysis of the PD-L1 and Immune Infiltrates of m(6)A RNA Methylation Regulators in Head and Neck Squamous Cell Carcinoma.\u003c/em\u003e Mol Ther Nucleic Acids, 2020. \u003cstrong\u003e21\u003c/strong\u003e: p. 299-314.\u003c/li\u003e\n\u003cli\u003eKamat, A.M., et al., \u003cem\u003eBladder cancer.\u003c/em\u003e Lancet, 2016. \u003cstrong\u003e388\u003c/strong\u003e(10061): p. 2796-2810.\u003c/li\u003e\n\u003cli\u003eLan, Q., et al., \u003cem\u003eThe Critical Role of RNA m(6)A Methylation in Cancer.\u003c/em\u003e Cancer Res, 2019. \u003cstrong\u003e79\u003c/strong\u003e(7): p. 1285-1292.\u003c/li\u003e\n\u003cli\u003eDai, D., et al., \u003cem\u003eN6-methyladenosine links RNA metabolism to cancer progression.\u003c/em\u003e Cell Death \u0026amp; Disease, 2018. \u003cstrong\u003e9\u003c/strong\u003e(2): p. 124.\u003c/li\u003e\n\u003cli\u003eHe, A., et al., \u003cem\u003ePrognostic value of long non-coding RNA signatures in bladder cancer.\u003c/em\u003e Aging (Albany NY), 2019. \u003cstrong\u003e11\u003c/strong\u003e(16): p. 6237-6251.\u003c/li\u003e\n\u003cli\u003eLian, P., et al., \u003cem\u003eAn eight-long non-coding RNA signature as a candidate prognostic biomarker for bladder cancer.\u003c/em\u003e Aging (Albany NY), 2019. \u003cstrong\u003e11\u003c/strong\u003e(17): p. 6930-6940.\u003c/li\u003e\n\u003cli\u003eHuang, D., et al., \u003cem\u003eNKILA lncRNA promotes tumor immune evasion by sensitizing T cells to activation-induced cell death.\u003c/em\u003e Nat Immunol, 2018. \u003cstrong\u003e19\u003c/strong\u003e(10): p. 1112-1125.\u003c/li\u003e\n\u003cli\u003eLi, F., et al., \u003cem\u003eProfiles of tumor-infiltrating immune cells and prognostic genes associated with the microenvironment of bladder cancer.\u003c/em\u003e Int Immunopharmacol, 2020. \u003cstrong\u003e85\u003c/strong\u003e: p. 106641.\u003c/li\u003e\n\u003cli\u003eZhen, S., et al., \u003cem\u003eSynergistic Antitumor Effect on Bladder Cancer by Rational Combination of Programmed Cell Death 1 Blockade and CRISPR-Cas9-Mediated Long Non-Coding RNA Urothelial Carcinoma Associated 1 Knockout.\u003c/em\u003e Hum Gene Ther, 2018. \u003cstrong\u003e29\u003c/strong\u003e(12): p. 1352-1363.\u003c/li\u003e\n\u003cli\u003ePoh, A.R. and M. Ernst, \u003cem\u003eTargeting Macrophages in Cancer: From Bench to Bedside.\u003c/em\u003e Front Oncol, 2018. \u003cstrong\u003e8\u003c/strong\u003e: p. 49.\u003c/li\u003e\n\u003cli\u003eOstuni, R., et al., \u003cem\u003eMacrophages and cancer: from mechanisms to therapeutic implications.\u003c/em\u003e Trends Immunol, 2015. \u003cstrong\u003e36\u003c/strong\u003e(4): p. 229-39.\u003c/li\u003e\n\u003cli\u003eNoy, R. and J.W. Pollard, \u003cem\u003eTumor-associated macrophages: from mechanisms to therapy.\u003c/em\u003e Immunity, 2014. \u003cstrong\u003e41\u003c/strong\u003e(1): p. 49-61.\u003c/li\u003e\n\u003cli\u003eGrivennikov, S.I., F.R. Greten, and M. Karin, \u003cem\u003eImmunity, inflammation, and cancer.\u003c/em\u003e Cell, 2010. \u003cstrong\u003e140\u003c/strong\u003e(6): p. 883-99.\u003c/li\u003e\n\u003cli\u003eChung, I.H., et al., \u003cem\u003eThe long non-coding RNA LINC01013 enhances invasion of human anaplastic large-cell lymphoma.\u003c/em\u003e Sci Rep, 2017. \u003cstrong\u003e7\u003c/strong\u003e(1): p. 295.\u003c/li\u003e\n\u003cli\u003eXue, Y., et al., \u003cem\u003eTumor\u003c/em\u003e\u003cem\u003e‑infiltrating M2\u0026nbsp;macrophages driven by specific genomic alterations are associated with prognosis in bladder cancer.\u003c/em\u003e Oncol Rep, 2019. \u003cstrong\u003e42\u003c/strong\u003e(2): p. 581-594.\u003c/li\u003e\n\u003cli\u003eDe, I., et al., \u003cem\u003eCSF1 Overexpression Promotes High-Grade Glioma Formation without Impacting the Polarization Status of Glioma-Associated Microglia and Macrophages.\u003c/em\u003e Cancer Res, 2016. \u003cstrong\u003e76\u003c/strong\u003e(9): p. 2552-60.\u003c/li\u003e\n\u003cli\u003eWu, Y., et al., \u003cem\u003eCCL3-CCR5 axis regulates intratumoral accumulation of leukocytes and fibroblasts and promotes angiogenesis in murine lung metastasis process.\u003c/em\u003e J Immunol, 2008. \u003cstrong\u003e181\u003c/strong\u003e(9): p. 6384-93.\u003c/li\u003e\n\u003cli\u003eRoca, H., et al., \u003cem\u003eCCL2 and interleukin-6 promote survival of human CD11b+ peripheral blood mononuclear cells and induce M2-type macrophage polarization.\u003c/em\u003e J Biol Chem, 2009. \u003cstrong\u003e284\u003c/strong\u003e(49): p. 34342-54.\u003c/li\u003e\n\u003cli\u003eJose, P.J., et al., \u003cem\u003eEotaxin: a potent eosinophil chemoattractant cytokine detected in a guinea pig model of allergic airways inflammation.\u003c/em\u003e J Exp Med, 1994. \u003cstrong\u003e179\u003c/strong\u003e(3): p. 881-7.\u003c/li\u003e\n\u003cli\u003eTripathi, C., et al., \u003cem\u003eMacrophages are recruited to hypoxic tumor areas and acquire a pro-angiogenic M2-polarized phenotype via hypoxic cancer cell derived cytokines Oncostatin M and Eotaxin.\u003c/em\u003e Oncotarget, 2014. \u003cstrong\u003e5\u003c/strong\u003e(14): p. 5350-68.\u003c/li\u003e\n\u003cli\u003eIshida, T. and R. Ueda, \u003cem\u003eCCR4 as a novel molecular target for immunotherapy of cancer.\u003c/em\u003e Cancer Sci, 2006. \u003cstrong\u003e97\u003c/strong\u003e(11): p. 1139-46.\u003c/li\u003e\n\u003cli\u003eBellmunt, J., T. Powles, and N.J. Vogelzang, \u003cem\u003eA review on the evolution of PD-1/PD-L1 immunotherapy for bladder cancer: The future is now.\u003c/em\u003e Cancer Treat Rev, 2017. \u003cstrong\u003e54\u003c/strong\u003e: p. 58-67.\u003c/li\u003e\n\u003cli\u003eNakanishi, J., et al., \u003cem\u003eOverexpression of B7-H1 (PD-L1) significantly associates with tumor grade and postoperative prognosis in human urothelial cancers.\u003c/em\u003e Cancer Immunol Immunother, 2007. \u003cstrong\u003e56\u003c/strong\u003e(8): p. 1173-82.\u003c/li\u003e\n\u003cli\u003eInman, B.A., et al., \u003cem\u003ePD-L1 (B7-H1) expression by urothelial carcinoma of the bladder and BCG-induced granulomata: associations with localized stage progression.\u003c/em\u003e Cancer, 2007. \u003cstrong\u003e109\u003c/strong\u003e(8): p. 1499-505.\u003c/li\u003e\n\u003cli\u003eWang, Y., et al., \u003cem\u003eImmune-related signature predicts the prognosis and immunotherapy benefit in bladder cancer.\u003c/em\u003e Cancer Med, 2020. \u003cstrong\u003e9\u003c/strong\u003e(20): p. 7729-7741.\u003c/li\u003e\n\u003cli\u003eWang, Z., M.A. Jensen, and J.C. Zenklusen, \u003cem\u003eA Practical Guide to The Cancer Genome Atlas (TCGA).\u003c/em\u003e Methods Mol Biol, 2016. \u003cstrong\u003e1418\u003c/strong\u003e: p. 111-41.\u003c/li\u003e\n\u003cli\u003eNewman, A.M., et al., \u003cem\u003eRobust enumeration of cell subsets from tissue expression profiles.\u003c/em\u003e Nat Methods, 2015. \u003cstrong\u003e12\u003c/strong\u003e(5): p. 453-7.\u003c/li\u003e\n\u003cli\u003eYoshihara, K., et al., \u003cem\u003eInferring tumour purity and stromal and immune cell admixture from expression data.\u003c/em\u003e Nat Commun, 2013. \u003cstrong\u003e4\u003c/strong\u003e: p. 2612.\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"scientific-reports","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"scirep","sideBox":"Learn more about [Scientific Reports](http://www.nature.com/srep/)","snPcode":"","submissionUrl":"","title":"Scientific Reports","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Scientific Reports","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"Bioinformatics analysis, N6-methyladenosine, lncRNAs, Immune cells infiltration, Bladder cancer, PD-L1","lastPublishedDoi":"10.21203/rs.3.rs-1341357/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-1341357/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003e\u003cstrong\u003eBackground \u0026amp; Aim:\u003c/strong\u003e\u003c/p\u003e\u003cp\u003eBladder cancer (BLCA) is one of the most frequent genitourinary cancers, with a high rate of morbidity and mortality.\u0026nbsp;The connection of m6A-related lncRNAs with PD-L1 and tumor immune microenvironment (TIME) in BLCA prognosis was extensively investigated in this study, which could suggest novel therapeutic targets for further investigation. \u003c/p\u003e\u003cp\u003e\u003cstrong\u003eResults:\u003c/strong\u003e\u003c/p\u003e\u003cp\u003e30 m6A-associated lncRNAs with predictive values from the TCGA data set were identified with co-expression analysis. Cluster2 was correlated with a poor prognosis,upregulated\u0026nbsp;\u0026nbsp;PD-L1 expression, and higher immune ratings. Cluster2 had larger amounts of resting CD4 memory-activated T cells, M2 macrophages, neutrophils, and NK cells infiltration. \u0026nbsp;“CHEMOKINE SIGNALING PATHWAY” was the most significantly enriched signaling pathway according to GSEA，which may play an important role in the diffrent immune cell infiltrates between cluster1/2. The risk model for m6A-related lncRNAs could be employed in a prognostic model to predict BLCA prognosis, regardless of other clinical features. Collectively, m6A-related lncRNAs were linked to PD-L1 and TIME, which would dynamically affect the number of tumor-infiltrating immune cells. \u003c/p\u003e\u003cp\u003e\u003cstrong\u003eConclusion:\u003c/strong\u003e\u003c/p\u003e\u003cp\u003em6A-related lncRNAs may be key mediators of PD-L1 expression and immune cells infiltration and may strongly affect the TIME of BLCA.\u003c/p\u003e","manuscriptTitle":"Comprehensive Analysis of the PD-L1 and Immune Infiltrates of N6-Methyladenosine Related Long Non-Coding RNAs in Bladder Cancer","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2022-03-10 21:08:06","doi":"10.21203/rs.3.rs-1341357/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Major revision","date":"2022-04-06T15:39:10+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2022-03-27T03:15:23+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"c481083e-2b58-4d9b-a743-ad8827886cf4","date":"2022-03-17T07:12:57+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"8f0be8d9-40ac-4b44-ba06-0613acb19374","date":"2022-03-14T07:38:21+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2022-03-14T06:49:21+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2022-03-14T06:41:51+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"","date":"2022-03-08T10:07:50+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2022-03-08T09:38:23+00:00","index":"","fulltext":""},{"type":"submitted","content":"Scientific Reports","date":"2022-02-09T06:20:43+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"scientific-reports","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"scirep","sideBox":"Learn more about [Scientific Reports](http://www.nature.com/srep/)","snPcode":"","submissionUrl":"","title":"Scientific Reports","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Scientific Reports","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"aa58c472-0738-4f45-a4d8-171ebb01217e","owner":[],"postedDate":"March 10th, 2022","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"under-review","subjectAreas":[],"tags":[],"updatedAt":"2022-06-01T05:29:21+00:00","versionOfRecord":[],"versionCreatedAt":"2022-03-10 21:08:06","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-1341357","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-1341357","identity":"rs-1341357","version":["v1"]},"buildId":"-HB7Z8yhvgn0wM9Nzuekk","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}