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

BACKGROUND

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.

METHODS

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.

RESULTS

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.

CONCLUSION

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.