Abstract
RNA modifications implicate pathological and prognosis significance in cancer development and progression, of which, m6A and m5C are representative regulators. These RNA modifications could produce effects on the function of other RNA by regulating gene expression. Thus, in this study, we aimed to explore the correlation between m6A/m5C regulators and early-stage lung adenocarcinoma (LUAD). Only the early-stage LUAD samples were included in this investigation, and the RNA-seq dataset of The Cancer Genome Atlas (TCGA) cohort was utilized to evaluate the expression of 37 m6A/m5C regulated genes. Based on the expression level of these 37 genes, early-stage LUAD patients were divided into 2 clusters, which were performed by consensus clustering, and the m6A/m5C subtypes had significantly different prognostic outcomes (p < 0.001). Cluster1, which has a better prognosis, was characterized by the C3 (inflammatory) immune subtype, low immune infiltration, chemokine expression, major histocompatibility complex (MHC) expression, and immune checkpoint molecule expression. Furthermore, compared with cluster1, cluster2 showed a T cell exhaustion state, characterized by a high expression of immune checkpoint genes, and immune cells, such as T cells, CD8+ T cells, cytotoxic lymphocytes, NK cells, and so on. In addition, patients in cluster2 were with high tumor mutational burden (TMB) and numerous significant mutated oncogene and tumor suppressor genes, such as WNT10B, ERBB4, SMARCA4, TP53, and CDKN2A (p < 0.001). A total of 19 genes were mostly related to the prognosis of LUAD and were upregulated in cluster2 (p < 0.05), showing a positive correlation with the mRNA expression of 37 m6A/m5C regulated genes. The predictive risk model was constructed using Cox and LASSO (least absolute shrinkage and selection operator) regression analysis. Finally, a seven-gene m6A/m5C risk model, comprising of METTL3, NPLOC4, RBM15, YTHDF1, IGF2BP1, NSUN3, and NSUN7, was constructed to stratify the prognosis of early-stage LUAD (p = 0.0049, AUC = 0.791). The high-risk score was associated with a poorer prognosis. This model was also validated using two additional GEO datasets: GSE72094 (p = 0.011, AUC = 0.736) and GSE50081 (p = 0.012, AUC = 0.628). In summary, it was established that the m6A/m5C-regulated genes performed a crosstalk function in the mRNA expression of early-stage LUAD. By interacting with other mRNA genes, m6A/m5C modification disturbs DNA replication and the tumor immune microenvironment (TIME). The seven-gene risk model may be a critical tool for the prognostic assessment of early-stage LUAD.