One of the most prevalent hematological system cancers is acute myeloid leukemia (AML). Efferocytosis-related genes (ERGs) and N6-methyladenosine (m6A) have an important significance in the progression of cancer, and the metastasis of tumors.
The AML-related data were collected from The Cancer Genome Atlas (TCGA; TCGA-AML) database and Gene Expression Omnibus (GEO; GSE9476, GSE71014, and GSE13159) database. The “limma” R package and Venn diagram were adopted to identify differentially expressed ERGs (DE-ERGs). The m6A related-DE-ERGs were obtained by Spearman analysis. Subsequently, univariate Cox and Least Absolute Shrinkage and Selection Operator (LASSO) were used to construct an m6A related-ERGs risk signature for AML patients. The possibility of immunotherapy for AML was explored. The pRRophetic package was adopted to calculate the IC50 of drugs for the treatment of AML. Finally, the expression of characterized genes was validated by quantitative reverse transcription-PCR (qRT-PCR).
Based on m6A related-DE-ERGs, a prognostic model with four characteristic genes (UCP2, DOCK1, SLC14A1, and SLC25A1) was constructed. The risk score of model was significantly associated with the immune microenvironment of AML, with four immune cell types, 14 immune checkpoints, 20 HLA family genes and, immunophenoscore (IPS) all showing differences between the high- and low-risk groups. A total of 56 drugs were predicted to differ between the two groups, of which Erlotinib, Dasatinib, BI.2536, and bortezomib have been reported to be associated with AML treatment. The qRT-PCR results showed that the expression trends of DOCK1, SLC14A1 and SLC25A1 were consistent with the bioinformatics analysis.
In summary, 4 m6A related- ERGs were identified and the corresponding prognostic model was constructed for AML patients. This prognostic model effectively stratified the risk of AML patients.