N6-methyladenosine-related lncRNAs identified as potential biomarkers for predicting the overall survival of Asian gastric cancer patients.

Gastric cancer (GC) is one of the most prevalent malignant tumors in Asian countries. Studies have proposed that lncRNAs can be used as diagnostic and prognostic indicators of GC due to the high specificity of lncRNAs expression involvement in GC. Recently, N6-methyladenosine (m6A) has also emerged as an important modulator of the expression of lncRNAs in GC. This study aimed at establishing a novel m6A-related lncRNAs prognostic signature that can be used to construct accurate models for predicting the prognosis of GC in the Asian population. First, the levels of m6A modification and m6A methyltransferases expression in GC samples were determined using dot blot and western blot analyses. Next, we evaluated the lncRNAs expression profiles and the corresponding clinical data of 88 Asian GC patients retrieved from The Cancer Genome Atlas (TCGA) database. Differential expression of m6A-related lncRNAs between GC and normal tissues was investigated. The relationship between these target lncRNAs and potential immunotherapeutic signatures was also analyzed. Gene set enrichment analysis (GSEA) was performed to identify the malignancy-associated pathways. Univariate Cox regression, LASSO regression, and multivariate Cox regression analyses were performed to establish a novel prognostic m6A-related lncRNAs prognostic signature. Moreover, we constructed a predictive nomogram and determined the expression levels of nine m6A-related lncRNAs in 12 pairs of clinical samples. We found that m6A methylation levels were significantly increased in GC tumor samples compared to adjacent normal tissues, and the increase was positively correlated with tumor stage. Patients were then divided into two clusters (cluster 1 and cluster 2) based on the differential expression of the m6A-related lncRNAs. Results showed that there was a significant difference in survival probability between the two clusters (p = 0.018). Notably, the low survival rate in cluster 2 may be associated with high expression of immune cells (resting memory CD4+ T cells, p = 0.027; regulatory T cells, p = 0.0018; monocytes, p = 0.00095; and resting dendritic cells, p = 0.015), and low expression of immune cells (resting NK cells, p = 0.033; and macrophages M1, p = 0.045). Enrichment analysis indicated that malignancy-associated biological processes were more common in the cluster 2 subgroup. Finally, the risk model comprising of six m6A-related lncRNAs was identified as an independent predictor of prognoses, which could divide patients into high- or low-risk groups. Time-dependent ROC analysis suggested that the risk score could accurately predict the prognosis of GC patients. Patients in the high-risk group had worse outcomes compared to patients in the low-risk group, and the risk score showed a positive correlation with immune cells (resting memory CD4+ T cells, R = 0.31, P = 0.038; regulatory T cells, R = 0.42, P = 0.0042; monocytes, R = 0.42, P = 0.0043). However, M1 macrophages (R = -0.37, P = 0.012) and resting NK cells (R = -0.31, P = 0.043) had a negative correlation with risk scores. Furthermore, analysis of clinical samples validated the weak positive correlation between the risk score and tumor stage. The risk model described here, based on the six m6A-related lncRNAs signature, and may predict the clinical prognoses and immunotherapeutic response in Asian GC patients.

Xu S ，Chen W ，Wang Y ，Zhang Y ，Xia R ，Shen J ，Gong X ，Liang Y ，Xu J ，Tang H ，Zhao T ，Zhang Y ，Chen T ，Wang C ... -  《BMC CANCER》

Identification of a novel m6A-related lncRNA pair signature for predicting the prognosis of gastric cancer patients.

Accumulating studies have demonstrated that lncRNAs play vital roles in the prognosis of gastric cancer (GC); however, the prognostic value of N6-methyladenosine-related lncRNAs has not been fully reported in GC. This study aimed to construct and validate an m6A-related lncRNA pair signature (m6A-LPS) for predicting the prognosis of GC patients. GC cohort primary data were downloaded from The Cancer Genome Atlas. We analysed the coexpression of m6A regulators and lncRNAs to identify m6A-related lncRNAs. Based on cyclical single pairing along with a 0-or-1 matrix and least absolute shrinkage and selection operator-penalized regression analyses, we constructed a novel prognostic signature of m6A-related lncRNA pairs with no dependence upon specific lncRNA expression levels. All patients were divided into high-risk and low-risk group based on the median risk score. The predictive reliability was evaluated in the testing dataset and whole dataset with receiver operating characteristic (ROC) curve analysis. Gene set enrichment analysis was used to identify potential pathways. Fourteen m6A-related lncRNA pairs consisting of 25 unique lncRNAs were used to construct the m6A-LPS. Kaplan-Meier analysis showed that the high-risk group had poor prognosis. The area under the curve for 5-year overall survival was 0.906, 0.827, and 0.882 in the training dataset, testing dataset, and whole dataset, respectively, meaning that the m6A-LPS was highly accurate in predicting GC patient prognosis. The m6A-LPS served as an independent prognostic factor for GC patients after adjusting for other clinical factors (p < 0.05). The m6A-LPS had more accuracy and a higher ROC value than other prognostic models for GC. Functional analysis revealed that high-risk group samples mainly showed enrichment of extracellular matrix receptor interactions and focal adhesion. Moreover, N-cadherin and vimentin, known biomarkers of epithelial-mesenchymal transition, were highly expressed in high-risk group samples. The immune infiltration analysis showed that resting dendritic cells, monocytes, and resting memory CD4 T cells were significantly positively related to the risk score. Thus, m6A-LPS reflected the infiltration of several types of immune cells. The signature established by pairing m6A-related lncRNAs regardless of expression levels showed high and independent clinical prediction value in GC patients.

Wang JM ，Li X ，Yang P ，Geng WB ，Wang XY ... -  《BMC GASTROENTEROLOGY》

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 ... -  《-》

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》

Analysis and prognostic significance of tumour immune infiltrates and immune microenvironment of m6A-related lncRNAs in patients with gastric cancer.

Studies have shown that long noncoding RNAs and N6-methyladenosine play important roles in gastric cancer. The purpose of this study was to determine the correlation and prognostic value of m6A-related lncRNAs and immune infiltration in gastric cancer. We downloaded the clinically related information and RNA-Seq transcriptome data of gastric cancer patients from the TCGA database. Univariate Cox regression analysis and Pearson analysis were used to screen out m6A-related lncRNAs. Consensus cluster analysis was used to divide the sample into two clusters, and LASSO analysis and Cox regression analysis were used to construct a risk scoring model. A total of 25 lncRNA expression profiles were screened, and gastric cancer patients were divided into different subtypes. Cluster 2 had a better prognosis, but its stromal score, ESTIMATE score and immune score were low. Cluster 1 was rich in resting memory CD4 T cells, regulatory T cells, monocytes, and resting mast cells, and Cluster 2 was rich in activated memory CD4 T cells and follicular helper T cells. Thirteen lncRNAs were selected to construct a risk model, and the prognosis of gastric cancer patients in the high-risk group was poor. The expression of PD-L1 in tumours is significantly higher than that in normal tissues. Univariate and multivariate Cox regression analysis results showed that the overall survival rate was significantly related to stage and the risk score, which can be used as an independent prognostic factor. The results of the heatmap and scatter plot showed that clusters (P = 0.0045) and grade (G1-2, G3, P = 0.0037) were significantly related to prognosis. The relationship between the risk score and immune cell infiltration showed that memory B cells, resting dendritic cells, M0 macrophages, and M2 macrophages were positively correlated with the risk score, while resting mast cells, monocytes, activated NK cells, and follicular helper T cells were negatively correlated with the risk score. The results of this study indicate that m6A-related lncRNAs may play an important role in the prognosis of gastric cancer patients and the tumour immune microenvironment and may provide help for the treatment of gastric cancer patients.

Huang J ，Song J ，Li X ，Liu S ，Huang W ，Shen Z ，Cheng Y ，Kou S ，Gao Z ，Tian Y ，Hu J ... -  《BMC Medical Genomics》