Comprehensive Analysis of N6-Methylandenosine-Related Long Non-Coding RNAs Signature in Prognosis and Tumor Microenvironment of Bladder Cancer.

To investigate the role of N6-methyladenosine (m6A)- related long non-coding RNAs (lncRNAs) in bladder cancer (BC). 50 m6A-related lncRNAs were screened out and were correlated with prognosis from BC samples in The Cancer Genome Atlas (TCGA). The lncRNAs were subdivided into cluster 1 and cluster 2 with consensus cluster analysis, and it was found that lncRNAs in cluster 2 were associated with poor prognosis and increased PD-L1 expression. Gene set enrichment analysis (GSEA) revealed tumor-related pathways in cluster 2. Through least absolute shrinkage and selection operator (LASSO) Cox regression analysis, univariate and multivariate Cox regression, and ROC analyses, 14 prognostic lncRNAs were selected and used to construct the m6A-related lncRNA prognostic signature (m6A-LPS), furthermore, that m6A-LPS was as a valuable independent prognostic factor. Interestingly, the m6A-LPS risk score was positively correlated with the immune score, PD-L1 expression, and the infiltration of immune cell subtypes in BC. SNHG16, a member of the high-risk group based on m6A-LPS, was highly expressed in BC tissues and cell lines and interfered with siRNA resulted in suppressed proliferation, migration, and invasion in vitro. Our study illustrates the role of m6A-related lncRNAs in BC. The m6A-LPS may be an important regulatory target of the tumor microenvironment (TME) in BC.

Chen K ，Zhu S ，Yu W ，Xia Y ，Xing J ，Geng J ，Cheng F ... -  《Frontiers in Oncology》

Relationships of N6-Methyladenosine-Related Long Non-Coding RNAs With Tumor Immune Microenvironment and Clinical Prognosis in Lung Adenocarcinoma.

Background: Lung adenocarcinoma (LUAD) is the major subtype of lung cancer and is associated with very high mortality. Emerging studies have shown that N6-methyladenosine (m6A)-related long non-coding (lnc) RNAs play crucial roles in tumor prognosis and the tumor immune microenvironment (TME). We aimed to explore the expression patterns of different m6A-related lncRNAs concerning patient prognosis and construct an m6A-related lncRNA prognostic model for LUAD. Methods: The prognostic value of m6A-related lncRNAs was investigated in LUAD samples from The Cancer Genome Atlas (TCGA). Potential prognostic m6A-related lncRNAs were selected by Pearson's correlation and univariate Cox regression analysis. Patients were divided into clusters using principal component analysis and the m6A-related lncRNA prognostic signature was calculated using least absolute shrinkage and selection operator (LASSO) Cox regression analysis. Results: Based on 91 prognostic m6A-related lncRNAs, we identified two m6A-related-lncRNA pattern clusters with different overall survival (OS) and different TMEs. We subsequently verified our findings multidimensionally by constructing a 13 m6A-related lncRNA prognostic signature (m6A-LPS) to calculate the risk score, which was robust in different subgroups. The receiver operating characteristic (ROC) curves and concordance index demonstrated that m6A-LPS harbored a promising ability to predict OS in TCGA data set and independent GSE11969 cohort. The risk score was also related to OS, TME, and clinical stage, and the risk score calculated by our model was also identified as independent prognostic predictive factors for LUAD patients after adjustment for age, smoking, gender, and stage. Enrichment analysis indicated that malignancy and drug resistance-associated pathways were more common in cluster2 (LUAD-unfavorable m6A-LPS). Furthermore, the results indicated that the signaling pathway enriched by the target gene of 13 m6A-related lncRNAs may be associated with metastasis and progression of cancer according to current studies. Conclusion: The current results indicated that different m6A-related-lncRNA patterns could affect OS and TME in patients with LUAD, and the prognostic signature based on 13 m6A-related lncRNAs may help to predict the prognosis in LUAD patients.

Zhao J ，Lin X ，Zhuang J ，He F ... -  《Frontiers in Genetics》

An effective N6-methyladenosine-related long non-coding RNA prognostic signature for predicting the prognosis of patients with bladder cancer.

Bladder cancer (BLCA) typically has a poor prognosis due to high relapse and metastasis rates. A growing body of evidence indicates that N6-methyladenosine (m6A) and long non-coding RNAs (lncRNAs) play crucial roles in the progression of BLCA and the treatment response of patients with BLCA. Therefore, we conducted a comprehensive RNA-seq analysis of BLCA using data from The Cancer Genome Atlas (TCGA) to establish an m6A-related lncRNA prognostic signature (m6A-RLPS) for BLCA. Consensus clustering analysis was used to investigate clusters of BLCA patients with varying prognoses. The least absolute shrinkage and selection operator Cox regression were used to develop the m6A-RLPS. The ESTIMATE and CIBERSORT algorithms were used to evaluate the immune composition. A total of 745 m6A-related lncRNAs were identified using Pearson correlation analysis (|R| > 0.4, p < 0.001). Fifty-one prognostic m6A-related lncRNAs were screened using univariate Cox regression analysis. Through consensus clustering analysis, patients were divided into two clusters (clusters 1 and 2) with different overall survival rates and tumor stages based on the differential expression of the lncRNAs. Enrichment analysis demonstrated that terms related to tumor biological processes and immune-related activities were increased in patient cluster 2, which was more likely to exhibit low survival rates. Nine m6A-related prognostic lncRNAs were finally determined and subsequently used to construct the m6A-RLPS, which was verified to be an independent predictor of prognosis using univariate and multivariate Cox regression analyses. Further, a nomogram based on age, tumor stage, and the m6A-RLPS was generated and showed high accuracy and reliability with respect to predicting the survival outcomes of BLCA patients. The prognostic signature was found to be strongly correlated to tumor-infiltrating immune cells and immune checkpoint expression. We established a novel m6A-RLPS with a favorable prognostic value for patients with BLCA. We believe that this prognostic signature can provide new insights into the tumorigenesis of BLCA and predict the treatment response in patients with BLCA.

Ma T ，Wang X ，Meng L ，Liu X ，Wang J ，Zhang W ，Tian Z ，Zhang Y ... -  《BMC CANCER》

Dissecting the Role of N6-Methylandenosine-Related Long Non-coding RNAs Signature in Prognosis and Immune Microenvironment of Breast Cancer.

Breast cancer (BC) represents a molecularly and clinically heterogeneous disease. Recent progress in immunotherapy has provided a glimmer of hope for several BC subtypes. The relationship between N6-methyladenosine (m6A) modification and long non-coding RNAs (LncRNAs) is still largely unexplored in BC. Here, with the intention to dissect the landscape of m6A-related lncRNAs and explore the immunotherapeutic value of the m6A-related lncRNA signature, we identified m6A-related lncRNAs by co-expression analysis from The Cancer Genome Atlas (TCGA) and stratified BC patients into different subgroups. Furthermore, we generated an m6A-related lncRNA prognostic signature. Four molecular subtypes were identified by consensus clustering. Cluster 3 preferentially had favorable prognosis, upregulated immune checkpoint expression, and high level of immune cell infiltration. Twenty-one m6A-related lncRNAs were applied to construct the m6A-related lncRNA model (m6A-LncRM). Survival analysis and receiver operating characteristic (ROC) curves further confirmed the prognostic value and prediction performance of m6A-LncRM. Finally, high- and low-risk BC subgroups displayed significantly different clinical features and immune cell infiltration status. Overall, our study systematically explored the prognostic value of the m6A-related LncRNAs and identified a high immunogenicity BC subtype. The proposed m6A-related LncRNA model might serve as a robust prognostic signature and attractive immunotherapeutic targets for BC treatment.

Zhang J ，Shan B ，Lin L ，Dong J ，Sun Q ，Zhou Q ，Chen J ，Han X ... -  《Frontiers in Cell and Developmental Biology》

Comprehensive analysis of tumor immune microenvironment and prognosis of m6A-related lncRNAs in gastric cancer.

N6-methyladenosine (m6A) modification and long non-coding RNAs (lncRNAs) play pivotal roles in gastric cancer (GC) progression. The emergence of immunotherapy in GC has created a paradigm shift in the approaches of treatment, whereas there is significant heterogeneity with regard to degree of treatment responses, which results from the variability of tumor immune microenvironment (TIME). How the interplay between m6A and lncRNAs enrolling in the shaping of TIME remains unclear. The RNA sequencing and clinical data of GC patients were collected from TCGA database. Pearson correlation test and univariate Cox analysis were used to screen out m6A-related lncRNAs. Consensus clustering method was implemented to classify GC patients into two clusters. Survival analysis, the infiltration level of immune cells, Gene set enrichment analysis (GSEA) and the mutation profiles were analyzed and compared between two clusters. A competing endogenous RNA (ceRNA) network and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were applied for the identification of pathways in which m6A-related lncRNAs enriched. Then least absolute shrinkage and selection operator (LASSO) COX regression was implemented to select pivotal lncRNAs, and risk model was constructed accordingly. The prognosis value of the risk model was explored. In addition, the response to immune checkpoint inhibitors (ICIs) therapy were compared between different risk groups. Finally, we performed qRT-PCR to detect expression patterns of the selected lncRNAs in the 35 tumor tissues and their paired adjacent normal tissues, and validated the prognostic value of risk model in our cohort (N = 35). The expression profiles of 15 lncRNAs were included to cluster patients into 2 subtypes. Cluster1 with worse prognosis harbored higher immune score, stromal score, ESTIMATE score and lower mutation rates of the genes. Different immune cell infiltration patterns were also displayed between the two clusters. GSEA showed that cluster1 preferentially enriched in tumor hallmarks and tumor-related biological pathways. KEGG pathway analysis found that the target mRNAs which m6A-related lncRNAs regulated by sponging miRNAs mainly enriched in vascular smooth muscle contraction, cAMP signaling pathway and cGMP-PKG signaling pathway. Next, eight lncRNAs were selected by LASSO regression algorithm to construct risk model. Patients in the high-risk group had poor prognoses, which were consistent in our cohort. As for predicting responses to ICIs therapy, patients from high-risk group were found to have lower tumor mutation burden (TMB) scores and account for large proportion in the Microsatellite Instability-Low (MSI-L) subtype. Moreover, patients had distinct immunophenoscores in different risk groups. Our study revealed that the interplay between m6A modification and lncRNAs might have critical role in predicting GC prognosis, sculpting TIME landscape and predicting the responses to ICIs therapy.

Wang Y ，Zhu GQ ，Tian D ，Zhou CW ，Li N ，Feng Y ，Zeng MS ... -  《BMC CANCER》