Comprehensive analyses indicated the association between m6A related long non-coding RNAs and various pathways in glioma.

Glioma is one of the most malignant brain tumors and diseases. N6-methyladenosine modification (m6A) is the most abundant and prevalent internal chemical modification of mRNA and long non-coding RNAs (lncRNAs) in eukaryotes. Nevertheless, the correlated pathways and clinical utilization of m6A-related lncRNAs have not been fully evaluated in glioma. Public RNA-sequencing and clinical annotation data were retrieved from TCGA, CGGA and GEO database. Differential expression analysis and univariate Cox regression analysis were performed to identify the m6A-related and differentially expressed lncRNAs with prognostic function (m6A-DELPF). The consensus clustering was performed to identify the expression pattern of m6A-DELPF. LASSO Cox regression analysis was performed to construct the lncRNA-based signature. The CIBERSORT and ESTIMATE algorithms were performed to analyze immune infiltration and tumor microenvironment, respectively. Immunotherapy sensitivity analysis was performed using data from TCIA. The small molecule drugs prediction analysis was performed using The Connectivity Map (CMap) database and STITCH database. A competing endogenous RNAs (ceRNA) network was constructed based on miRcode, miRDB, miRTarBase, TargetScan database. Two clusters (cluster1 and cluster2) were identified after unsupervised cluster analysis based on m6A-DELPF. Additionally, a 15-gene prognostic signature namely m6A-DELPFS was constructed. Analyses of epithelial-mesenchymal-transition score, tumor microenvironment, immune infiltration, clinical characterization analysis, and putative drug prediction were performed to confirm the clinical utility and efficacy of m6A-DELPFS. The potential mechanisms including tumor immune microenvironment of m6A-DELPF influence the initiation and progression of glioma. A clinically accessible nomogram was also constructed based on the m6A-DELPF and other survival-relevant clinical parameters. Two miRNAs and 114 mRNAs were identified as the downstream of seven m6A-related lncRNAs in a ceRNA network. Our present research confirmed the clinical value of m6A related lncRNAs and their high correlation with tumor immunity, tumor microenvironment, tumor mutation burden and drug sensitivity in glioma.

Chen Z ，Zhang W ，Yan Z ，Zhang M ... -  《Cancer Medicine》

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》

N6-methyladenosine Modification-Related Long Non-Coding RNAs are Potential Biomarkers for Predicting the Prognosis of Patients With Osteosarcoma.

Yang K ，Wang F ，Li K ，Peng G ，Yang H ，Xu H ，Xiang Y ，Sun H ... -  《-》

The prognostic value and immune landscaps of m6A/m5C-related lncRNAs signature in the low grade glioma.

N6-methyladenosine (m6A) and 5-methylcytosine (m5C) are the main RNA methylation modifications involved in the oncogenesis of cancer. However, it remains obscure whether m6A/m5C-related long non-coding RNAs (lncRNAs) affect the development and progression of low grade gliomas (LGG). We summarized 926 LGG tumor samples with RNA-seq data and clinical information from The Cancer Genome Atlas and Chinese Glioma Genome Atlas. 105 normal brain samples with RNA-seq data from the Genotype Tissue Expression project were collected for control. We obtained a molecular classification cluster from the expression pattern of sreened lncRNAs. The least absolute shrinkage and selection operator Cox regression was employed to construct a m6A/m5C-related lncRNAs prognostic signature of LGG. In vitro experiments were employed to validate the biological functions of lncRNAs in our risk model. The expression pattern of 14 sreened highly correlated lncRNAs could cluster samples into two groups, in which various clinicopathological features and the tumor immune microenvironment were significantly distinct. The survival time of cluster 1 was significantly reduced compared with cluster 2. This prognostic signature is based on 8 m6A/m5C-related lncRNAs (GDNF-AS1, HOXA-AS3, LINC00346, LINC00664, LINC00665, MIR155HG, NEAT1, RHPN1-AS1). Patients in the high-risk group harbored shorter survival times. Immunity microenvironment analysis showed B cells, CD4 + T cells, macrophages, and myeloid-derived DC cells were significantly increased in the high-risk group. Patients in high-risk group had the worse overall survival time regardless of followed TMZ therapy or radiotherapy. All observed results from the TCGA-LGG cohort could be validated in CGGA cohort. Afterwards, LINC00664 was found to promote cell viability, invasion and migration ability of glioma cells in vitro. Our study elucidated a prognostic prediction model of LGG by 8 m6A/m5C methylated lncRNAs and a critical lncRNA regulation function involved in LGG progression. High-risk patients have shorter survival times and a pro-tumor immune microenvironment.

Li R ，Chen H ，Li C ，Qi Y ，Zhao K ，Wang J ，You C ，Huang H ... -  《BMC BIOINFORMATICS》

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》