Skip to main content

Landscape of internal N7-methylguanosine of long non-coding RNA modifications in resistant acute myeloid leukemia



Growing evidence indicates that RNA methylation plays a fundamental role in epigenetic regulation, which is associated with the tumorigenesis and drug resistance. Among them, acute myeloid leukemia (AML), as the top acute leukemia for adults, is a deadly disease threatening human health. Although N7-methylguanosine (m7G) has been identified as an important regulatory modification, its distribution has still remained elusive.


The present study aimed to explore the long non-coding RNA (lncRNA) functional profile of m7G in AML and drug-resistant AML cells. The transcriptome-wide m7G methylation of lncRNA was analyzed in AML and drug-resistant AML cells. RNA MeRIP-seq was performed to identify m7G peaks on lncRNA and differences in m7G distribution between AML and drug-resistant AML cells. The Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were conducted to predict the possible roles and m7G-associated pathway.


Using m7G peak sequencing, it was found that a sequence motif was necessary for m7G methylation in drug-resistant AML lncRNA. Unsupervised hierarchical cluster analysis confirmed that lncRNA m7G methylation occurred more frequently in drug-resistant AML cells than in AML cells. RNA sequencing demonstrated that more genes were upregulated by methylation in drug-resistant AML cells, while methylation downregulated more genes in AML cells. The GO and KEGG pathway enrichment analyses revealed that genes having a significant correlation with m7G sites in lncRNA were involved in drug-resistant AML signaling pathways.


Significant differences in the levels and patterns of m7G methylation between drug-resistant AML cells and AML cells were revealed. Furthermore, the cellular functions potentially influenced by m7G in drug-resistant AML cells were predicted, providing evidence implicating m7G-mediated lncRNA epigenetic regulation in the progression of drug resistance in AML. These findings highlight the involvement of m7G in the development of drug resistance in AML.

Peer Review reports


Acute myeloid leukemia (AML) is a type of blood cancer that originates in the bone marrow and is characterized by the uncontrolled growth of myeloid cells. Drug resistance in AML remains a significant challenge to effective treatment, with several patients relapsing after initial therapy. Understanding the mechanisms underlying drug resistance is therefore an essential step towards developing novel therapies [1]. Emerging evidence suggests that RNA methylation, particularly m7G methylation, of long non-coding RNAs (lncRNAs) plays a critical role in the development of drug resistance in AML. Alterations in RNA methylation patterns, including aberrant m7G methylation of lncRNAs, have been implicated in the acquisition and maintenance of drug resistance in AML [2].

Dysregulated m7G methylation of specific lncRNAs has been associated with altered gene expression patterns, dysfunctions in cellular signaling pathways, and the development of drug-resistant phenotypes in AML cells. Aberrant m7G methylation of lncRNAs has been found to be associated with the modulation of drug efflux pumps, drug metabolism pathways, DNA repair mechanisms, and anti-apoptotic processes, contributing to the resistance of AML cells to therapeutic interventions [3]. Targeting dysregulated m7G methylation of drug-resistant AML-associated lncRNAs presents a potential strategy for overcoming treatment resistance and improving patient outcomes. Comprehensive profiling of RNA methylation patterns, including specific analysis of m7G methylation in lncRNAs, may enable the identification of novel biomarkers and therapeutic targets for combating drug resistance in AML [4]. Integration of RNA methylation data, including the assessment of m7G methylation in lncRNAs, with clinical and genomic information can enhance our understanding of the mechanisms driving drug resistance in AML and guide the development of personalized treatment strategies.

As a post-transcriptional modification, RNA methylation plays a fundamental role in epigenetic regulation. To date, over 100 RNA methylations have been reported [5]. Examples are m7G (7-methyl guanosine), m5C (5-methylcytosine), m1A (N1-methyladenosine), m6Am (2-O-dimethyladenosine), and m6A (N6-methyladenosine). During epigenetic regulation, RNA methylation influences various bio-processes, including RNA splicing, translation, nuclear export, stability, DNA damage repair, initiation of miRNA biogenesis and immunogenicity, thereby affecting cellular differentiation, embryonic development, spermatogenesis, sex determination, learning and memory, immune response and occurrence and development of cancer [6]. For instance, SRY-box transcription factor 2 modification by m6A in glioblastoma is positively correlated with radiation resistance and maintenance of cells similar to glioma stem via METTL3, and m1A in tRNA drives liver tumorigenesis by regulating cholesterol metabolism[7, 8].

Notably, m7G is a widely recognized RNA modification on mRNA caps and inside rRNAs, tRNAs and mRNAs[9]. Despite its abundance, m7G is highly difficult to be investigated by using standard sequencing-based technologies owing to high stability in neutral aqueous solutions and neutrality to Watson-Crick base pair of m7G in RNA, without affecting reverse transcription. In virtue of advances in techniques such as high-throughput sequencing, m7G modifications in RNA species with low contents of m7G, including lncRNAs and mRNA, can be identified and quantified. Previous studies have reported functions and distribution of m7G in different classes of RNA. Nevertheless, understanding of distribution over the transcriptome and prevalence of m7G in lncRNA remains at its infancy [10,11,12].

LncRNAs are RNA molecules with over 200 nucleotides but no or limited protein-coding potential. Indeed, LncRNAs serve as key gene expression regulators in multiple pathological and physiological processes, including tumorigenesis, cell differentiation, and inflammation [13]. Nowadays, lncRNAs have attracted great attention in post-transcriptional regulation, as well as epigenetic modification. Specifically, lncRNAs interact with RNAs, proteins and chromatins in cytoplasm/nucleus, which is associated with multiple diseases (e.g., diverse types of cancer) [14,15,16]. For instance, lncRNA MEG3 contributes to drug resistance in AML by positively regulating ALG9 through sponging miR-155[17].

According to the American Cancer Society, around 20,000 cases of AML were diagnosed in 2020in the United States. Despite some breakthroughs in AML diagnosis and treatment over the recent decades, relapsed and refractory nature of AML remains significant owing to its heterogeneity and malignancy [18]. Nevertheless, recent studies revealed a key role of methylation modifications in AML progression, demonstrating it as a novel effective AML therapy [19,20,21,22].

In the present study, the influence of m7G methylation on lncRNAs in AML-resistant cells was investigated. Specifically, m7G was globally mapped via RNA MeRIP-seq in AML-resistant and non-resistant cells, while distribution and prevalence in the two cell types and cellular compartments were determined and compared with each other. The results demonstrated that the m7G modification degree was significantly higher in AML-resistant cells than that in non-resistant cells. This could be reflected by intra-cell consistency, as well as inter-cells differences, involving all chromosomes. This study could facilitate the understanding of epigenetic regulation of lncRNAs in AML by m7G, as well as new therapies for AML.

Materials and methods

Cell lines and cell culture

HL60/MX2 and HL60 cells were acquired from the American Type Culture Collection and cultured in RPMI-1640 containing 10% thermally deactivated fetal calf serum, 1% (v/v) penicillin streptomycin in a humidified atmosphere containing 5% CO2 at 37 °C.

Development of cDNA libraries

In the present study, m7G peaks were acquired from each sample by using the MeRIP-seq technique. RIP was executed by a commercial service provider. High-quality reads were aligned to the UCSC HG19 by utilizing hisat2 (ver. 2.0.4). Then, fragments per kilobase of exon per million fragments mapped values from transcript-level LncRNAs were obtained by utilizing cuffdiff (ver. 2.2.1) under guidance by the gtf files. Additionally, the P-value and fold-change of the two groups were determined to screen for lncRNAs that were expressed differentially.

Analysis of sequencing data

Sequencing and visualization of the libraries developed were achieved by using a sequencer (Illumina NovaSeq 6000). The mapping quality is reflected by ‘Q30’: a read above 30 was considered to be reliable. Meanwhile, library mapping was achieved by STAR (ver. 2.5.1b). Sequencing datasets were aligned to UCSC HG19 by utilizing Hisat2 (ver. 2.0.4). Methylation peaks of the two cells were identified by using MACS, while their differences were determined by using DiffReps.

The GO and KEGG pathway enrichment analyses

The GO and KEGG pathway enrichment analyses were carried out using the DAVID online tool to clarify the properties of the above-mentioned genes. The GO analysis aimed to analyze MFs (molecular functions), CCs (cellular components), and BPs (biological processes) of these genes based on differentially methylated genes. P < 0.05 was set to be the critical level of significant GO terms. Additionally, the KEGG pathway analysis [23,24,25] was conducted to investigate functions and biological pathways. Similarly, P < 0.05 indicated a significant difference.


Identification of LncRNAs related to transcriptome-wide m7G methylation in human AML-resistant cells

m7G methylation over transcriptome in the two cells was investigated by RNA MeRIP-seq of lncRNA. As observed, a total of 7,049 and 6,170 m7G peaks were observed in HL-60 cells and HL60/MX2 cells, respectively (Fig. 1a). Up to 6491 and 5,757 annotated genes were mapped in HL-60 and HL60/MX2 cells, respectively (Fig. 1b). Among them, 462 m7G peaks appeared in both HL-60 cells and HL60/MX2 cells. A fold-change cut-off of 2.0 was employed to identify differentially methylated peaks. Additionally, 99% and 84% of sites identified in this study were indeed specific to HL60/MX2 and HL-60 cells, respectively, suggesting that the possibility of m7G methylation is higher in AML-resistant cells than in AML cells.

Fig. 1
figure 1

m7G methylation (transcriptome-wide) and role of lncRNAs in the two cells. (A) Venn diagram of m7G methylation sites recognized in lncRNAs from the two cells. (B) Venn diagram of m7G genes in the two cells. (C) The sequence motif of m7G sites in HL-60 cells. (D) The sequence motif of m7G sites in MX2/ HL60 cells. (E) Percentages of lncRNAs containing different quantities of m7G peaks in the two cells; most lncRNAs contain one m7G peak only. (F) A schematic showing the pathways and biological functions affected and how they differ in drug-resistant versus sensitive AML [26]

The sequences of m7G-methylated peaks were investigated by scanning to verify that an m7G motif is present. URUUUY (R = A/G, Y = U/C) was determined to be the most conservative motif in HL-60 cells and its E-value was 1.1e-009 (Fig. 1C). Meanwhile, the most conservative MOTIF for m7G peaks in HL60/MX2 cells was YMCAG (Y = C/U, M = C/A) and its E-value was 9.7e-020 (Fig. 1D), indicating that the MOTIF sequences of m7G methylation exhibited significant difference in the two groups.

Additionally, the density of m7G peaks in lncRNA was determined on the basis of the cell type. The results showed that 50% of lncRNAs comprised of a m7G peak in the CDS region in HL-60 cells, compared to 45.9% of single peak lncRNA in MX2/ HL60 cells. Percentages for the 3’UTR region (20.1% in HL-60 cells vs. 17.7% in HL60/MX2 cells), 5’UTR region (10.2% in HL-60 cells vs. 13.3% in HL60/MX2 cells), stopC (11.5% in HL-60 cells vs. 11.1% in HL60/MX2 cells) and startC (8.2% in HL-60 cells vs. 11.9% in HL60/MX2 cells) were also determined (Fig. 1E). A schematic showing the pathways and biological functions affected and how they differ in drug-resistant versus sensitive AML is presented in Fig. 1F [26].

Analysis of sources of lncRNAs methylation in the two cells

The source data of methylated lncRNAs were summarized and a pie chart was plotted (Fig. 2A and B). The results showed that most of the lncRNAs were derived from exon sense overlapping regions, followed by intergenic regions. Compared with HL60 cells, HL60/MX2 cells had fewer methylated lncRNAs derived from intergenic regions, and more methylated lncRNAs from exon sense overlapping regions. Meanwhile, more lncRNAs were derived from exons in HL60/MX2 cells, which may lead to an increase of the number of lncRNAs that can perform silencing and affect cell function. Additionally, the distributions of lncRNA methylation sites on the chromosomes were analyzed by using Circos (Fig. 2C). The results showed that the lncRNA methylation site distribution on each chromosome was different between the two groups. Indeed, the X chromosome and chromosome 2 were the more evident ones. Also, there are various types of leukemia and tumors associated with specific regions on chromosome 2. Compared with autosomes, sex chromosomes in both groups were less methylated.

Fig. 2
figure 2

 A and B. Pie chart of the source of methylated lncRNA in the two cells. C. Visualized chromosome-level distribution of m7G in lncRNAs in the two cells

Analysis of differential expression of lncRNAs on the basis of m7G methylation

The expressions of lncRNAs in the two cells were investigated by RNA sequencing according to m7G methylation. It was revealed that 301 lncRNAs were upregulated and exhibited increased expression levels in the HL60/MX2 cells. On the other hand, 355 lncRNAs were downregulated in HL60/MX2 cells compared with HL-60 cells (Fig. 3A). The top 10 downregulated and upregulated lncRNAs are listed in Table 1.

Fig. 3
figure 3

The differentiated expressions of lncRNAs categorized based on m7G methylation. (A) Differential expression of lncRNAs in the two cells. (B) Cumulative distribution of lncRNA expressions in the two cells for genes with upregulated m7G (red) and downregulated m5C (green); blue denotes others

Table 1 Top 10 upregulated and downregulated lncRNAs in HL60 and HL60/MX2 cells screened by microarray

The top 10 lncRNAs with statistically different transcription levels of hypermethylated genes are listed in Table 2. Functional abbreviations and names of such hypermethylated genes are summarized in Table 3. These genes may be correlated with AML resistance pathogenesis and shall be further investigated. Additionally, the number of downregulated m7G modified lncRNAs was higher than that of upregulated ones in the HL60/MX2 cells (Fig. 3B). Consequently, m7G may be associated with gene expression in various transcripts of AML resistant cells.

Table 2 Top 10 differentially expressed lncRNAs in HL60/MX2 and HL60 cells
Table 3 Functional roles of top 10 differentially expressed lncRNAs

Analysis of m7G genes in lncRNAs of MX2/ HL60 cells by KEGG and GO enrichment

The biological classifications of m7G genes in lncRNAs of MX2/ HL60 cells were investigated by functional and pathway enrichment analyses.

The GO analysis revealed that BPs of m7G genes that were up-methylated in MX2/ HL60 cells were significantly enriched in cellular metabolic process, nitrogen compound metabolic process, organic substance metabolic process, and primary metabolic process, while the genes with m7G that were down-methylated were mainly related to nitrogen compound metabolic process and primary metabolic process, as well as their regulation. Enrichment of MF of the up-methylated m7G genes in MX2/ HL60 cells was basically in protein binding, ion binding, and heterocyclic compound binding, while the genes with m7G that were down-methylated were mainly related to metal ion binding, cation binding, enzyme binding, and protein binding. Enrichment of CC of the up-methylated m7G genes in MX2/ HL60 cells was basically in cytoplasm, intracellular, intracellular membrane-bounded organelles, intracellular organelles, membrane-bounded organelles, and organelles, which is consistent with methylation sites upregulated in HL60 cells (Fig. 4).

Fig. 4
figure 4

GO analysis of m7G genes in lncRNAs of HL60/MX2 cells. (A-C) The leading ten significantly enriched GO terms for (A) BPs, (B) CCs, and (C) MFs in up-methylated m7G genes of HL60/MX2 cells. (D–E) The leading ten significantly enriched GO terms of (D) BPs, (E) CCs, and (F) MFs in down-methylated m7G genes of MX2/ HL60 cells

The KEGG pathway analysis [23,24,25] revealed significant enrichment of lncRNAs with up-methylated m7G in HL60-MX2 cells primarily in the fanconi anemia pathway, other types of O-glycan biosynthesis, and non-homologous end-joining. Conversely, lncRNAs with down-methylated m7G showed significant enrichment in the Fc gamma R-mediated phagocytosis, B cell receptor signaling pathway, and osteoclast differentiation (Fig. 5).

Fig. 5
figure 5

KEGG analysis of m7G genes in lncRNAs of MX2/ HL60 cells. (A) Top 10 significantly enriched pathways for up-methylated m7g genes in MX2/ HL60 cells. (B) Top 10 significantly enriched pathways for down-methylated m7G genes in MX2/ HL60 cells


Risk stratification in AML patients was conducted based on the results of cytogenetic and molecular analyses [27]. Specifically, the AML cases can be categorized into three sub-groups: favorable, intermediate, and poor outcomes. In the favorable sub-group, inv(16), t(8;21) and t(15;17) are most common cytogenetic abnormalities, while those with adverse prognosis exhibited abnormal 3q and 11q23 abnormalities besides t(9;11), as well as complex karyotypes. In the intermediate sub-group, most patients exhibit normal karyotypes and this case is regarded as CN-AML (cytogenetically normal AML). CN-AML accounts for 50% of adult AML patients, and it refers to patients with significantly broad diversity of clinical outcomes and molecular characterization. Meanwhile, CN-AML patients exhibited a wide range of molecular features and a wide diversity of clinical outcomes. Despite the FDA approval of several medicine for the treatment of AML, primary and secondary drug resistances remain a major issue in recent years. In one specific sub-group, some patients exhibited complete remission, while others showed drug resistance. Although the mechanisms of AML resistance to multiple drugs have been thorough investigated (genetic biomarkers), it is increasingly difficult for a single target to achieve AML treatment. The current AML classification faces huge challenges in prognosis and effective treatment due to the heterogeneity in each subgroup [28,29,30].

Epigenetic alterations, especially methylation, have emerged as effective molecular markers for evaluation of prognosis or drug therapeutic potential of cancers. Genes methylation has a high incidence (60-90%) in gastric cancer, Hodgkin, non-Hodgkin’s lymphoma, and breast cancer [31,32,33,34]. Located in cytoplasm and nuclei, eIF4E plays a key role in export and translation of specific m(7)G-capped mRNAs, respectively [35]. Nuclear accumulation of eIF4E in AML patients is related to increased export of transcripts, which are responsible for encoding oncoproteins, dependent on eIF4E. Additionally, its subcellular localization is directly related to responses of the patient. However, the expression level of m7G, as well as its role in AML resistance, remain poorly clarified.

As a potential critical layer of biological regulation, lncRNAs are aberrantly expressed in AML by transcription factors, DNA methylation or histone modifications, as well as by post-transcriptional modifications. Dysregulated lncRNAs may alter activity and expression of mRNAs, proteins and miRNAs, leading to epigenetic alterations in AML. For instance, AC026150.8, which is upregulated in AML, is correlated with poor prognosis, and its overexpression is positively related to drug resistance of AML cells, and confirms its scaffolding effect combined with splicing factors [36]. WT1 was upregulated in HOX transcript antisense RNA by targeting miR-20a-5p, and the regulation of resistance of Adriamycin by curcumin in AML cells can be reversed by inhibiting miR-20a-5p [37]. The transcriptional expression of lncRNAs in AML has been extensively investigated. Nevertheless, the mechanism by which m7G methylation in lncRNAs promotes the progression of AML drug resistance remains unclear, and further exploration is needed. In our study, the extent of cell methylation was estimated by examining the ratio of m7G in AML cells and AML-resistant cells. RNA MeRIP-seq was employed to globally locate m7G in human AML resistant and non-resistant cells; the distribution and prevalence of m7G in both cells and corresponding cellular compartments were compared with each other to clarify the m7G methylation levels of lncRNA in AML resistant cells. A substantial increase of m7G modification in AML resistant cells compared to non-resistant cells was observed. Meanwhile, 764 lncRNAs were upregulated, with high expression levels in HL60/MX2 cells. On the other hand, 870 lncRNAs were downregulated in HL 60/MX2 cells. The results suggested a correlation of molecular heterogeneous genomic events in AML with drug resistance. Specifically, the m7G site and the lncRNA set of genes localized by the m7G site in the AML-resistant cells far exceeded those in the AML cells. This was observed in all chromosomes. Cluster analysis showed that the methylation level is an effective index to distinguish AML cells from drug-resistant cells, confirming the correlation of m7G with AML resistance.

The present study utilized RNA MeRIP-seq to investigate m7G methylation over the transcriptome in human AML-resistant cells. The analysis revealed a total of 7,049 m7G peaks in HL-60 cells and 6,170 m7G peaks in HL60/MX2 cells. These peaks were associated with annotated genes, with 6,491 genes mapped in HL-60 cells and 5,757 genes mapped in HL60/MX2 cells. The study identified specific m7G motifs in the two cell types. In HL-60 cells, the most conservative motif was determined to be URUUUY (R = A/G, Y = U/C), while in HL60/MX2 cells, the most conservative motif was YMCAG (Y = C/U, M = C/A). The presence of these motifs indicated significant differences in the sequences of m7G methylation between the two cell groups. The density of m7G peaks in lncRNAs was also analyzed based on the cell type. The results showed that approximately 50% of lncRNAs in HL-60 cells and 45.9% of lncRNAs in HL60/MX2 cells contained m7G peaks in the CDS region. The percentages varied in other regions, such as 3’UTR, 5’UTR, stop codon (stopC), and start codon (startC), suggesting differential distribution patterns of m7G methylation in lncRNAs between the two cell types. Furthermore, the study analyzed the sources of methylated lncRNAs and found that most of the methylated lncRNAs were derived from exon sense overlapping regions, followed by intergenic regions. The distributions of lncRNA methylation sites across chromosomes were also examined, revealing differences between the two groups. The X chromosome and chromosome 2 exhibited notable differences in the distribution of lncRNA methylation sites. The differential expression of lncRNAs based on m7G methylation was assessed, and it was observed that 301 lncRNAs were upregulated and 355 lncRNAs were downregulated in HL60/MX2 cells compared to HL-60 cells. The study also provided a list of the top 10 upregulated and downregulated lncRNAs. Functional and pathway enrichment analyses were performed to investigate the biological classifications of m7G genes in lncRNAs of HL60/MX2 cells. GO analysis revealed enriched BPs and MFs that were associated with the up-methylated and down-methylated m7G genes. The up-methylated genes were enriched in cellular metabolic processes, protein binding, and ion binding, while the down-methylated genes were related to nitrogen compound metabolic processes and metal ion binding, among others. The KEGG pathway analysis demonstrated that the lncRNAs with up-methylated m7G in HL60/MX2 cells were enriched in pathways, such as the fanconi anemia pathway, other types of O-glycan biosynthesis, and non-homologous end-joining. On the other hand, the lncRNAs with down-methylated m7G were enriched in pathways including Fc gamma R-mediated phagocytosis, B cell receptor signaling pathway, and osteoclast differentiation.

Uptake of amino acid, catabolism and homeostate levels are elevated in LSCs (leukemia stem cells). LSCs that are obtained from naive AML patients are dependent on the metabolism of amino acids for OP (oxidative phosphorylation) and survival [35]. LSCs isolated from relapsed AML patients are independent of metabolism of amino acids owing to compensation by increased metabolism of fatty acids (FAs). It is known that cancer cells primarily rely on glycolysis for energy production rather than oxidative phosphorylation (OP). However, recent studies have demonstrated that cancer stem cells (CSCs) in different types of tumors depend on oxidative phosphorylation and display lower glycolytic reserves compared to non-stem cancer cells. The metabolic shift known as the Warburg effect, characterized by aerobic glycolysis, plays a pivotal role in cancer cells, leading to increased lactate production and heightened glucose consumption [35]. Metabolic plasticity and reprogramming are correlated with drug resistance, metastasis, and tumorigenesis [38]. For instance, mutations of Flt3-ITD elicit AKT-mediated PHOS of mitochondrial HK2 (hexokinase 2) and facilitate aerobic glycolysis in AML cells. 3-Brop (3-Bromopyruvate) leads to increased cell necrosis by sensitizing tumor cells to therapeutic performances of doxorubicin and sorafenib. mTORC1 is constitutively activated in AML cells and plays a significant role in the sensitivity of glycolytic inhibitors. The biosynthesis of FAs can be generated by reprogramming of glycolysis and glutaminolysis, thus facilitating survival and growth of cancer cells. In summary, the metabolic reprogramming signatures of AML were observed, and the inhibition of related pathways may contribute to the sensitivity of chemotherapeutic medicine. Interestingly, our Subsequent biological function predictions revealed that lncRNAs upregulated by m7G modification methylation were mainly enriched in complement activation and glycolysis/gluconeogenesis-related metabolic pathways, which may suggest that in AML cells, regulation of lncRNAs by m7G methylation could positively regulate glycolysis/gluconeogenesis-related metabolic pathways leading to cell drug resistance.

Additionally, it was noted that certain genes, which exhibited low expression levels in AML cells, demonstrated high expression in AML-resistant cells accompanied by the elevated levels of methylation. However, the precise role of this methylation in AML resistance necessitates further investigation. The present study provided a new perspective for the future study on the role of m7G methylation in drug resistance of AML cells to discover new therapeutic targets.


In conclusion, the present study explored the regulatory mechanisms of m7G methylation in AML resistance. The findings shed light on the identification of differentially methylated lncRNAs, specific m7G motifs, and their distribution patterns in AML-resistant cells compared with AML cells. The results suggest a potential role for m7G methylation in modulating gene expression in AML-resistant cells. However, it is important to note that further investigation is warranted to fully understand the epigenetic alterations associated with AML resistance. Future research in this area may provide valuable insights into the underlying mechanisms and potential therapeutic targets for overcoming AML resistance.

Data Availability

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the following link:


  1. Kirtonia A, Ashrafizadeh M, Zarrabi A, Hushmandi K, Zabolian A, Bejandi AK, Rani R, Pandey AK, Baligar P, Kumar V, Das BC, Garg M. Long noncoding RNAs: a novel insight in the leukemogenesis and drug resistance in acute myeloid leukemia. J Cell Physiol. 2022;237(1):450–65.

    Article  CAS  PubMed  Google Scholar 

  2. Izadirad M, Jafari L, James AR, Unfried JP, Wu ZX, Chen ZS. Long noncoding RNAs have pivotal roles in chemoresistance of acute myeloid leukemia. Drug Discovery Today. 2021;26(7):1735–43.

    Article  CAS  PubMed  Google Scholar 

  3. Zhang J, Gao J, Hu M, Xu S, Cheng C, Zheng W, Zhang J. Integrated investigation of the clinical implications and targeted landscape for RNA methylation modifications in hepatocellular carcinoma. Eur J Med Res. 2023;28(1):46.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Wang Y, Fang Z, Hong M, Yang D, Xie W. Long-noncoding RNAs (lncRNAs) in drug metabolism and disposition, implications in cancer chemo-resistance. Acta Pharm Sinica B. 2020;10(1):105–12.

    Article  CAS  Google Scholar 

  5. Yang B, Wang J-Q, Tan Y, Yuan R, Chen Z-S, Zou C. RNA methylation and cancer treatment. Pharmacol Res. 2021;174:105937.

    Article  CAS  PubMed  Google Scholar 

  6. Song P, Tayier S, Cai Z, Jia G. RNA methylation in mammalian development and cancer. Cell Biol Toxicol. 2021;37(6):811–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Maimaiti A, Tuersunniyazi A, Meng X, Pei Y, Ji W, Feng Z, Jiang L, Wang Z, Kasimu M, Wang Y, et al. N6-methyladenosine RNA methylation regulator-related alternative splicing gene signature as prognostic predictor and in immune microenvironment characterization of patients with low-grade glioma. Front Genet. 2022;13:872186.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Qi L, Wang Y, Hu H, Li P, Hu H, Li Y, Wang K, Zhao Y, Feng M, Lyu H, et al. mA methyltransferase METTL3 participated in sympathetic neural remodeling post-MI via the TRAF6/NF-κB pathway and ROS production. J Mol Cell Cardiol. 2022;170:87–99.

    Article  CAS  PubMed  Google Scholar 

  9. Cheng W, Gao A, Lin H, Zhang W. Novel roles of METTL1/WDR4 in tumor via mG methylation. Mol Ther Oncolytics. 2022;26:27–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Chen Y, Lin H, Miao L, He J. Role of N7-methylguanosine (mG) in cancer. Trends Cell Biol 2022.

  11. Wu X, Li C, Wang Z, Zhang Y, Liu S, Chen S, Chen S, Liu W, Liu X. A bioinformatic analysis study of mG regulator-mediated methylation modification patterns and tumor microenvironment infiltration in glioblastoma. BMC Cancer. 2022;22(1):729.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Han H, Zheng S, Lin S. N-methylguanosine (mG) tRNA modification: a novel autophagy modulator in cancer. Autophagy 2022:1–3.

  13. Herman AB, Tsitsipatis D, Gorospe M. Integrated lncRNA function upon genomic and epigenomic regulation. Mol Cell. 2022;82(12):2252–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Xu Q, Guo T. Somatic mutation-associated risk index based on lncRNA expression for predicting prognosis in acute myeloid leukemia. Hematol (Amsterdam Netherlands). 2022;27(1):659–71.

    CAS  Google Scholar 

  15. Kong W, Yin G, Zheng S, Liu X, Zhu A, Yu P, Zhang J, Shan Y, Ying R, Jin H. Long noncoding RNA (lncRNA) HOTAIR: pathogenic roles and therapeutic opportunities in gastric cancer. Genes Dis. 2022;9(5):1269–80.

    Article  CAS  PubMed  Google Scholar 

  16. Chen Z, Zhang M, Lu Y, Ding T, Liu Z, Liu Y, Zhou Z, Wang L. Overexpressed lncRNA FTX promotes the cell viability, proliferation, migration and invasion of renal cell carcinoma via FTX/miR–4429/UBE2C axis. Oncol Rep 2022, 48(3).

  17. Yu Y, Kou D, Liu B, Huang Y, Li S, Qi Y, Guo Y, Huang T, Qi X, Jia L. LncRNA MEG3 contributes to drug resistance in acute myeloid leukemia by positively regulating ALG9 through sponging miR-155. Int J Lab Hematol. 2020;42(4):464–72.

    Article  PubMed  Google Scholar 

  18. Weigert N, Rowe JM, Lazarus HM, Salman MY. Consolidation in AML: abundant opinion and much unknown. Blood Rev. 2022;51:100873.

    Article  PubMed  Google Scholar 

  19. Lejman M, Dziatkiewicz I, Jurek M. Straight to the point-the novel strategies to Cure Pediatric AML. Int J Mol Sci 2022, 23(4).

  20. Jones LM, Tarlock K, Cooper T. Targeted therapy in Pediatric AML: an Evolving Landscape. Paediatr Drugs. 2021;23(5):485–97.

    Article  PubMed  Google Scholar 

  21. Li M, Zhang D. DNA methyltransferase-1 in acute myeloid leukaemia: beyond the maintenance of DNA methylation. Ann Med. 2022;54(1):2011–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Zhao X, Liu H-Q, Wang L-N, Yang L, Liu X-L. Current and emerging molecular and epigenetic disease entities in acute myeloid leukemia and a critical assessment of their therapeutic modalities. Sem Cancer Biol. 2022;83:121–35.

    Article  CAS  Google Scholar 

  23. Kanehisa M, Goto S. KEGG: Kyoto Encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28:1947–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51:D587–92.

    Article  CAS  PubMed  Google Scholar 

  26. Premnath N, Madanat YF. Novel investigational agents and pathways that May influence the Future Management of Acute myeloid leukemia. Cancers (Basel). 2023;15(11):2958.

    Article  CAS  PubMed  Google Scholar 

  27. Wang M, Yang C, Zhang L, Schaar DG. Molecular mutations and their Cooccurrences in Cytogenetically Normal Acute myeloid leukemia. Stem Cells Int. 2017;2017:6962379.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Yuan X-Q, Peng L, Zeng W-J, Jiang B-Y, Li G-C, Chen X-P. DNMT3A R882 mutations predict a poor prognosis in AML: a Meta-analysis from 4474 patients. Medicine. 2016;95(18):e3519.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Papaioannou D, Ozer HG, Nicolet D, Urs AP, Herold T, Mrózek K, Batcha AMN, Metzeler KH, Yilmaz AS, Volinia S, et al. Clinical and molecular relevance of genetic variants in the non-coding transcriptome of patients with cytogenetically normal acute myeloid leukemia. Haematologica. 2022;107(5):1034–44.

    Article  CAS  PubMed  Google Scholar 

  30. Yang L, Zhang H, Yang X, Lu T, Ma S, Cheng H, Yen K, Cheng T. Prognostic prediction of cytogenetically normal Acute Myeloid Leukemia based on a gene expression model. Front Oncol. 2021;11:659201.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Hong YH, Aziz N, Park JG, Lee D, Kim JK, Kim SA, Choi W, Lee CY, Lee HP, Huyen Trang HT, et al. The EEF1AKMT3/MAP2K7/TP53 axis suppresses tumor invasiveness and metastasis in gastric cancer. Cancer Lett. 2022;544:215803.

    Article  CAS  PubMed  Google Scholar 

  32. Liu H, Gu J, Huang Z, Han Z, Xin J, Yuan L, Du M, Chu H, Wang M, Zhang Z. Fine particulate matter induces METTL3-mediated mA modification of BIRC5 mRNA in bladder cancer. J Hazard Mater. 2022;437:129310.

    Article  CAS  PubMed  Google Scholar 

  33. Liu R, Jia Y, Kong G, He A. Novel insights into roles of N6-methyladenosine reader YTHDF2 in cancer progression. J Cancer Res Clin Oncol. 2022;148(9):2215–30.

    Article  CAS  PubMed  Google Scholar 

  34. Gonzalez ME, Naimo GD, Anwar T, Paolì A, Tekula SR, Kim S, Medhora N, Leflein SA, Itkin J, Trievel R, et al. EZH2 T367 phosphorylation activates p38 signaling through lysine methylation to promote breast cancer progression. iScience. 2022;25(8):104827.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Takahashi S, Shibutani S, Iwata H. Nuclear-targeted 4E-BP1 is dephosphorylated, induces nuclear translocation of eIF4E, and alters mRNA translation. Exp Cell Res. 2022;418(1):113246.

    Article  CAS  PubMed  Google Scholar 

  36. Zhang H, Zhao Y, Liu X, Liu Y, Wang X, Fu Y, Fu S, Zhang J. A novel upregulated LncRNA-AC026150.8 promotes chemo-resistance and predicts poor prognosis in acute myeloid leukemia. Cancer Med. 2021;10(23):8614–29.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Zhao W, Geng D, Li S, Chen Z, Sun M. LncRNA HOTAIR influences cell growth, migration, invasion, and apoptosis via the miR-20a-5p/HMGA2 axis in breast cancer. Cancer Med. 2018;7(3):842–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Pouysségur J, Marchiq I, Parks SK, Durivault J, Ždralević M, Vucetic M. ‘Warburg effect’ controls tumor growth, bacterial, viral infections and immunity - Genetic deconstruction and therapeutic perspectives. Sem Cancer Biol 2022.

Download references


Not applicable.


This study was supported by the National Natural Science Foundation of China (Grant No. 81473484).

Author information

Authors and Affiliations



JYH conceptualized and designed the study, and drafted the initial manuscript. QQL, YZ, and DL analyzed the data and revised the manuscript. DL and RW designed the data collection instruments, coordinated and supervised data collection, and critically reviewed the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Ran Wang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Han, J., Liu, Q., Zhou, Y. et al. Landscape of internal N7-methylguanosine of long non-coding RNA modifications in resistant acute myeloid leukemia. BMC Genomics 24, 425 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: