Background
N6-methyladenosine (m6A) has been confirmed to function critically in acute myeloid leukemia (AML) progression. Hitherto, the subtyping and prognostic predictive significance of m6A-correlated genes in AML is unclear.
Method
From The Cancer Genome Atlas (TCGA-LAML), Therapeutically Applicable Research to Generate Effective Treatments (TARGET-AML) and Gene Expression Omnibus (GEO, GSE71014) databases, we collected the sequencing data of AML patients. The batch effect was removed via limma package for TCGA-LAML and TARGET-AML, and the aggregated samples were AML cohorts. Samples in the AML cohort identified m6A models in AML by consensus clustering based on 23-m6A-related modulators. M6A-related differentially expressed genes (m6ARDEGs) influencing the overall survival (OS) of AML were determined by performing differential expression analysis and univariate COX analysis, and consensus-based clustering was utilized to access AML molecular subtypes. LASSO and multivariate COX analyses were performed to obtain the optimized m6ARDEGs to construct the m6A Prognostic Risk Score (m6APR_Score). Whether the model was robust was evaluated according to Kaplan–Meier (K-M) and receiver operator characteristic (ROC) curves. Further, the abundance of immune cell infiltration was explored in different m6A modification patterns and molecular subtypes and m6APR_Score groupings. Finally, nomogram was constructed to predict OS in AML. Quantitative real-time polymerase chain reaction (RT-qPCR) and cell counting kit-8 (CCK-8) assay were used to validate the genes in m6APR_Score in AML cells.
Results
The m6A models (m6AM1, m6AM2, m6AM3) and molecular subtypes (C1, C2, C3) were identified in the AML cohort, exhibiting different prognosis and immunoreactivity. We recognized novel prognostic biomarkers of AML such as CD83, NRIP1, ACSL1, METTL7B, OGT, and C4orf48. AML patients were grouped into high-m6APR_Score and low-m6APR_Score groups, with the later group showing a better prognosis than former one. Both the AML cohort and the validation cohort GSE71014 demonstrated excellent prediction. Finally, the nomogram accurately predicted the survival of patients suffering from AML. Further, the decision curves showed that both nomogram and m6APR_Score showed excellent prediction. It was confirmed in vitro experiments that mRNA expressions of NRIP1, ACSL1, METTL7B and OGT were elevated, while CD83 and C4orf48 mRNA expressions downregulated in AML cells. A significant increase in the viability of U937 and THP-1 cell lines after inhibition of CD83, while siMETTL7B had contrast results.
Conclusion
Our study demonstrated that m6APR_Score and CD83, NRIP1, ACSL1, METTL7B, OGT, and C4orf48 potentially provided novel and promising prognostic support for AML patients.