Transcriptome-wide analysis of RNA N6-methyladenosine in Adriamycin-resistant acute myeloid leukemia cells

N6-Methyladenosine (m6A) is the most abundant post-transcriptional modification in eukaryotes, the imbalance of which is reported to be associated with various pathological processes, including drug resistance. In this study, we analyzed the methylated RNA immunoprecipitation combined with next-generation sequencing (MeRIP-seq) data of AML cell line HL60 and its adriamycin-resistant cell line HL60/ADR. We found a total of 40550 peaks, representing 15640 genes in HL60, and a total of 38834 peaks, representing 15285 genes in HL60/ADR. A total of 4437 differentially methylated m6A peaks within 3461 genes have been found between HL60 and HL60/ADR. Among them, 3587 differentially m6A peaks within 2790 genes were hyper-methylated, and 850 m6A peaks within 671 genes were hypo-methylated. KEGG pathway analysis showed that pathways were enriched in tumor and drug-resistant related signaling pathway. Results of MeRIP-seq showed that fold enrichment of global m6A peaks was higher in HL60/ADR compared to HL60. This study provides a framework for the application of comprehensive mRNA m6A profiling towards acute myeloid leukemia cell line (HL60) and its adriamycin-resistant acute myeloid leukemia cell line (HL60/ADR). mRNA N6-methyladenosine profiles of acute myeloid leukemia cell line (HL60) and its adriamycin-resistant acute myeloid leukemia cell line (HL60/ADR)