Abstract
Aberrant regulation of m6A mRNA modification can lead to changes in gene expression, thus contributing to tumorigenesis in several types of solid tumors. In this study, by integrating analyses of m6A methylation and mRNA expression, we identified 84 m6A-regulated mRNAs in lung adenocarcinoma (LUAD). Although the m6A methylation levels of total RNA in LUAD patient tumor tissue were reduced, the majority (75.2%) of m6A-regulated mRNAs were hypermethylated. The m6A-hypermethylated mRNAs were mainly enriched in terms related to transcription factor activity. We established a 10-m6A-regulated-mRNA signature score system through least absolute shrinkage and selection operator Cox regression analysis, with its predictive value validated by Kaplan-Meier curve and time-dependent receiver operating characteristic curves. RFXAP and KHDRBS2 from the signature also exhibited an independent prognostic value. The co-expression and interaction network analyses demonstrated the strong correlation between m6A regulators and the genes in the signature, further supporting the results of the m6A methylation modification patterns. These findings highlight the potential utility of integrating multi-omics data (m6A methylation level and mRNA expression) to accurately obtain potential prognostic biomarkers, which may provide important insights into developing novel and effective therapies for LUAD.