Crosstalk of Eight Types of RNA Modification Regulators Defines Tumor Microenvironments, Cancer Hallmarks, and Prognosis of Lung Adenocarcinoma

Abstract

RNA modification has become an exciting underexplored field in recent years. In lung adenocarcinoma (LUAD), m6A was the best characterized and most studied RNA modification, while knowledge about other kinds of RNA modifications in LUAD is limited. In our study, we included a total of 100 RNA modification regulators of eight types of cancer-related RNA modifications (m6A, m1A, m5C, Nm, m7G, Ψ, A-to-I, and mcm5s2U) to systematically profile their specific roles in LUAD. By gene mutation and expression analysis, we identified extensive dysregulations and complicated interactions of 100 RNA modification regulators in LUAD. Based on unsupervised clustering analysis, gene set variation analysis (GSVA), and single-sample gene-set enrichment analysis (ssGSEA), two RNA modification patterns in LUAD were defined to show distinct biological characteristics. The favorable prognostic pattern was enriched with infiltrated immune cells, including activated B cells, CD8 T cells, eosinophil cells, dendritic cells, and natural killer cells, while the unfavorable prognostic pattern was enriched with cancer hallmarks, including hypoxia, epithelial-mesenchymal transition (EMT), angiogenesis, PI3K-AKT-mTOR pathway, MYC pathway, and glycolysis pathway. We also constructed an RNA modification score (RMScore) based on five critical genes (CYP17A1, NTSR1, PITX3, KRT6A, and ANLN) to evaluate the RNA modification status of individual LUAD patients. RMScore was revealed to be related to the infiltrated immune cells and cancer hallmarks and was an independent prognostic factor in the TCGA-LUAD cohort and two external GEO-LUAD cohorts. Our study was the first to comprehensively investigate the dysregulations, crosstalk, and potential prognostic value of eight types of RNA modifications in LUAD. Our results highlighted the significance of eight types of RNA modifications in tumor microenvironments and cancer hallmarks and provided novel prognostic biomarkers and potential therapeutic targets in the management of LUAD patients in the future.