AUTHOR=Shao Dongqi , Li Yu , Wu Junyong , Zhang Binbin , Xie Shan , Zheng Xialin , Jiang Zhiquan 

TITLE=An m6A/m5C/m1A/m7G-Related Long Non-coding RNA Signature to Predict Prognosis and Immune Features of Glioma

JOURNAL=Frontiers in Genetics

VOLUME=13

YEAR=2022

URL=https://www.frontiersin.org/journals/genetics/articles/10.3389/fgene.2022.903117

DOI=10.3389/fgene.2022.903117

ISSN=1664-8021

ABSTRACT=<p><bold>Background:</bold> Gliomas are the most common and fatal malignant type of tumor of the central nervous system. RNA post-transcriptional modifications, as a frontier and hotspot in the field of epigenetics, have attracted increased attention in recent years. Among such modifications, methylation is most abundant, and encompasses N6-methyladenosine (m6A), 5-methylcytosine (m5C), N1 methyladenosine (m1A), and 7-methylguanosine (m7G) methylation.</p><p><bold>Methods:</bold> RNA-sequencing data from healthy tissue and low-grade glioma samples were downloaded from of The Cancer Genome Atlas database along with clinical information and mutation data from glioblastoma tumor samples. Forty-nine m6A/m5C/m1A/m7G-related genes were identified and an m6A/m5C/m1A/m7G-lncRNA signature of co-expressed long non-coding RNAs selected. Least absolute shrinkage and selection operator Cox regression analysis was used to identify 12 m6A/m5C/m1A/m7G-related lncRNAs associated with the prognostic characteristics of glioma and their correlation with immune function and drug sensitivity analyzed. Furthermore, the Chinese Glioma Genome Atlas dataset was used for model validation.</p><p><bold>Results:</bold> A total of 12 m6A/m5C/m1A/m7G-related genes (AL080276.2, AC092111.1, SOX21-AS1, DNAJC9-AS1, AC025171.1, AL356019.2, AC017104.1, AC099850.3, UNC5B-AS1, AC006064.2, AC010319.4, and AC016822.1) were used to construct a survival and prognosis model, which had good independent prediction ability for patients with glioma. Patients were divided into low and high m6A/m5C/m1A/m7G-LS groups, the latter of which had poor prognosis. In addition, the m6A/m5C/m1A/m7G-LS enabled improved interpretation of the results of enrichment analysis, as well as informing immunotherapy response and drug sensitivity of patients with glioma in different subgroups.</p><p><bold>Conclusion:</bold> In this study we constructed an m6A/m5C/m1A/m7G-LS and established a nomogram model, which can accurately predict the prognosis of patients with glioma and provides direction toward promising immunotherapy strategies for the future.</p>