Comprehensive analysis of m6A/m5C/m1A-related gene expression, immune infiltration, and sensitivity of antineoplastic drugs in glioma

Abstract

This research aims to develop a prognostic glioma marker based on m⁶A/m⁵C/m¹A genes and investigate the potential role in the tumor immune microenvironment. Data for patients with glioma were downloaded from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA). The expression of genes related to m⁶A/m⁵C/m¹A was compared for normal and glioma groups. Gene Ontology and Kyoto Encyclopedia of Genes and Gene enrichment analysis of differentially expressed genes were conducted. Consistent clustering analysis was performed to obtain glioma subtypes and complete the survival analysis and immune analysis. Based on TCGA, Lasso regression analysis was used to obtain a prognostic model, and the CGGA database was used to validate the model. The model-based risk scores and the hub genes with the immune microenvironment, clinical features, and antitumor drug susceptibility were investigated. The clinical glioma tissues were collected to verify the expression of hub genes via immunohistochemistry. Twenty genes were differentially expressed, Consensus cluster analysis identified two molecular clusters. Overall survival was significantly higher in cluster 2 than in cluster 1. Immunological analysis revealed statistically significant differences in 26 immune cells and 17 immune functions between the two clusters. Enrichment analysis detected multiple meaningful pathways. We constructed a prognostic model that consists of WTAP, TRMT6, DNMT1, and DNMT3B. The high-risk and low-risk groups affected the survival prognosis and immune infiltration, which were related to grade, gender, age, and survival status. The prognostic value of the model was validated using another independent cohort CGGA. Clinical correlation and immune analysis revealed that four hub genes were associated with tumor grade, immune cells, and antitumor drug sensitivity, and WTAP was significantly associated with microsatellite instability(MSI). Immunohistochemistry confirmed the high expression of WTAP, DNMT1, and DNMT3B in tumor tissue, but the low expression of TRMT6. This study established a strong prognostic marker based on m⁶A/m⁵C/m¹A methylation regulators, which can accurately predict the prognosis of patients with gliomas. m⁶A/m⁵C/m¹A modification mode plays an important role in the tumor microenvironment, can provide valuable information for anti-tumor immunotherapy, and have a profound impact on the clinical characteristics.

Keywords: N1-adenylate methylation (m1A); N6-adenylate methylation (m6A); RNA methylation modification; cytosine hydroxylation (m5C); glioma; tumor immune microenvironment.