Identification of miR-342-5p/MDM4/p53 network in acute myeloid leukemia

Abstract

Acute myeloid leukemia (AML) is one of the most prevalent hematological malignancies. miRNAs play roles in cancer initiation and progression in various cancer types by post-transcriptional regulation of gene expression. The aim of this study is to investigate the mechanisms in the development and progression of acute myeloid leukemia and to identify potential target genes and miRNAs by bioinformatic analysis. miRNA expression profiles were obtained from the GSE51908 dataset on the Gene Expression Omnibus (GEO). GEO2R was used to identify differentially expressed miRNAs. The diagnostic and overall survival effects of the identified miRNA were determined using ROC analysis and Kaplan-Meier curve, respectively. Putative miRNA targets were determined based on miRWalk and miRDB tools. The expression change and overall survival analysis of the identified target gene were analyzed by Gene Expression Profiling Interactive Analysis (GEPIA). Protein-protein interaction (PPI) networks of the target gene were determined using STRING and GeneMANIA. Functional enrichment analysis was performed using the DAVID program. 24 DE-miRNAs were identified, including 16 upregulated and 8 downregulated genes. miR-342-5p expression had significantly shorter survival than those in higher expression control group (p = 0.0001), and its AUC value to discriminate AML from control groups was 0.795. High expression of MDM4 predicts an unfavorable prognosis in AML patients. The MDM4 gene was determined to be associated with decreased survival rates. According to KEGG results, microRNAs, p53 signaling pathway, and cell cycle are associated with AML development. The current study based on the GEO database, miR-342-5p/MDM4/p53 axis AML may provide new therapeutic targets.

1 Introduction

Acute myeloid leukemia (AML) is characterized by the accumulation of immature blast cells due to the inhibition of normal hematopoiesis, and aberrant proliferation and differentiation of immature myeloid hematopoietic cells in the bone marrow (BM) (Xin et al., 2022; Cheng et al., 2020). It is the most prevalent variety of acute leukemia in adults and has a high morbidity and mortality rate (Cheng et al., 2020; Juliusson et al., 2012). Tumor formation in AML is highly complex and is associated with various gene mutations (Cheng et al., 2020). In recent years, there have been advances in genomic, transcriptomic and epigenomic studies of AML (Kirtonia et al., 2022). Although these studies have improved our understanding of AML, its pathogenesis has not yet been completely elucidated (Liu et al., 2019). AML affects numerous organs, including the lymph nodes, liver, spleen, and central nervous system, and has a poor prognosis (Cheng et al., 2020). AML can be treated with chemotherapy, hematopoietic stem cell transplantation, cell therapy, biological immunotherapy, and gene-targeted therapy, although the success rate of treatment is not that high (Xin et al., 2022). Despite the evolving treatment options, the ineffectiveness of chemotherapy and the lack of information on the molecular pathogenesis of AML of the disease still remain a problem (Kirtonia et al., 2022). The majority of patients experience a relapse after obtaining complete remission, and some even dies (Cheng et al., 2020). Accurate diagnosis, prognosis, and efficient treatment of AML depend on the identification and development of therapeutic targets (Liu et al., 2019). It is predicted that studies on the pathogenesis and therapeutic targets of AML may have important clinical application value (Cheng et al., 2020). For this purpose, revealing the potential roles of non-coding RNAs, especially in AML, will be crucial in terms of understanding the molecular pathways related to AML development, chemotherapeutic response, and relapse. Previous studies have shown the role of non-coding RNAs (ncRNAs) in the pathogenesis of cancers, including AML (Kirtonia et al., 2022).

MicroRNAs (miRNAs) are highly conserved small non-coding RNA molecules, usually 18–24 nucleotides long, acting as negative post transcriptional regulation of gene expression (He and Hannon, 2004). The abnormal expression of miRNAs is associated with a variety of diseases, including cancer, and offers crucial information on the molecular pathophysiology of the diseases (O’Brien et al., 2018). miRNAs are dysregulated levels in AML and are crucial for the development or repression of the disease (Fletcher et al., 2022). Through alterations in proximity to the oncogenic genomic area, epigenetic modifications, inappropriate targeting of miRNA promoter regions by altered transcription factors or oncoproteins, and processing of deregulated miRNAs, these molecules could be involved in the pathogenesis of AML (Liu et al., 2019). Depending on the molecular pathways of the mRNA they target, miRNAs exhibit either oncogenic or tumor-suppressive properties (Otmani and Lewalle, 2021). There are many miRNAs involved in the development of AML, one of which is miR-342 (Fletcher et al., 2022; Journal et al., 2019; Expression et al., 2022). Abnormal expression of miR-342 has been observed in various types of cancer (hepatocellular carcinoma, cervical cancer, lung cancer, etc.) and has been associated with tumor development and proliferation (Gao et al., 2017; Li X. et al., 2014; Zhao and Zhang, 2015). miR-342 functions as a tumor suppressor in AML (Journal et al., 2019). However, the molecular mechanisms involved in miR-342 underlying the pathogenesis of AML remain unclear. Identifying the potential target of miR-342 and the associated biological pathways may suggest the miR-342 is a new potential therapeutic target in AML patients.

Murine double minute (MDM) proteins are overexpressed in a various of cancer types, and they primarily exhibit their oncogenic properties by inhibiting the p53 tumor suppressor (Li and Lozano, 2013). It also modulates and responds to many signaling networks (Wade et al., 2013). Both MDM2 (murine double minute 2) and MDM4 (murine double minute 4, also know as MDMX) form a heterodimer that closely regulates p53 function (Eskandari et al., 2021). MDM4 binds p53 to induce transcriptional inactivation and thus inhibits p53 function. Experimental evidence showed that p53 is kept inactive in the cytoplasm at high levels by interactions with MDM4. Further investigation of p53 in AML cell lines OCI/AML-2 (AML2) is mostly due to cytoplasmic and MDM4 binding (Tan B. X. et al., 2014). It is crucial to disrupt the MDM-p53 relationship in order to activate p53 while developing effective therapeutic approaches for tumor therapy. Therefore, scientists have focused on the roles of MDM4 as a negative regulators of the p53 tumor suppressor (Li and Lozano, 2013).

Microarray analysis can be used to identify miRNA changes in AML (Candia et al., 2015). Based on these data, new and functional miRNAs can be identified by bioinformatic analysis. In recent years, bioinformatics analysis has been used in oncology research to identify genetic alterations and new potential cancer-related biomarkers (Sheng et al., 2021). Few studies have been conducted to analyze the prognostic lncRNA-miRNA-mRNA (Wang et al., 2019) and circRNA-miRNA-mRNA (Lv et al., 2018) ceRNA network in AML (Cheng et al., 2020).

The aim of this study is to investigate the mechanisms underlying the development of AML, potential target genes, and new potential biomarkers for AML prognosis through miRNA-associated mRNA network in silico analysis. We focus on the mir342-5p/MDM4/p53 network given these findings because it may present exciting opportunities for comprehending the underlying molecular process and for cancer therapy.

2 Materials and methods

2.1 Microarray data

The dataset we used for this study was downloaded from the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) (http://www.ncbi.nlm.nih.gov/geo) (Edgar et al., 2002). In the present study miRNAs expression profile data of GSE51908 based on the platform of GPL8786 (Affymetrix Multispecies miRNA-1 Array) were obtained from the GEO database, which was deposited by CivinCandia et al. (2015), Tan Y. S. et al. (2014). GSE51908 dataset includes 190 samples (AML cell lines and patient samples, B ALL cell lines and Patient samples, T ALL cell lines and patient samples, normal B cells, granulocytes, normal monocytes, T cells and CD34⁺ cells). In this study, we selected 42 sample for AML (AML cell lines and patient sample) and 50 control (normal B cells, normal granulocytes, normal monocytes, normal T cells and normal CD34⁺ cells from human peripheral blood of normal healthy individuals) from GSE51908 datasets.

2.2 Data processing

GEO2R (www.ncbi.nlm.nih.gov/geo/geo2r) was used to identify differentially expressed miRNAs (DE-miRNAs) in the microarray datasets between the AML and control groups (Bao and Jiang, 2019). A large number of experimental datasets are included in GEO2R, and false-positive rates are modified using an adjusted P-value (adj. P). For the purpose of choosing DE-miRNAs, the adjusted p-value cut-off was set at p < 0.05 and llogFCl>2.

2.3 Prediction of miRNA targets

The miRWalk database (http://mirwalk.umm.uni-heidelberg.de/search_mirnas Version 3.0) was used to estimate the targets of the overlapping genes among miRNA datasets. MiRWalk is an intuitive interface that generates predicted and verified miRNA binding sites. To improve the precision of miRNA target prediction, the miRDB and miRTarBase methods in the miRWalk database were used (Sticht et al., 2018). Only target genes that were recognized by both databases were used to choose putative mRNAs. Additionally, AML cell lines (HL-60 and KG-1 cells) used the online database miRDB to estimate the targets of miRNA (Chen and Wang, 2020).

2.4 Differential expression of genes in AML and healthy tissue

Based on information from the Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx), the interactive web program Gene Expression Profiling Interactive Analysis (GEPIA) provides a variety of visualization and analysis capabilities for gene expression. Using the gene expression profiling database GEPIA (http://gepia.cancer-pku.cn/), we identified changes in the expression of identified genes between AML and healthy tissue (Tang et al., 2017).

2.5 Survival analysis

Using TCGA data from AML, GEPIA’s online tool was used to carry out a survival analysis of MDM4 (Tang et al., 2017). In order to assess and compare the survival rates of AML to death, Kaplan-Meier analysis and the log-rank test were employed to generate survival curves. To calculate and compare the overall survival (OS) of AML with control groups, the Cox proportional hazards regression model was employed. The hazard ratio (HR) with 95% confidence intervals (CIs) and log rank P-value were evaluated.

2.6 Construction of protein‐protein interactions (PPI) network

The STRING database is a protein-protein association network that provides all known and predicted information about the direct (physical) and indirect (functional) relationships that occur between different proteins (Szklarczyk et al., 2011). Another online resource, GeneMANIA is a database consisting of an intuitive interface for gene function predictions and interactions of genes with each other (Warde-farley et al., 2010). The PPI network of MDM4-associated target genes was established by the STRING (http://www.string-db.org/) (Szklarczyk et al., 2011), and GeneMANIA (http://genemania.org/) database (Warde-farley et al., 2010).

2.7 Functional enrichment analysis

The Gene Ontology (GO) resource is a free public resource that explains how genes move in biological systems and provides information about cellular location, molecular functions, and biological processes (Gene et al., 2011). The Kyoto Encyclopedia of Genes and Genomes (KEGG) contains a wide variety of databases, including analysis of genome sequences, biological functions, diseases, drugs, and chemicals (Kanehisa and Goto, 2000). The Database for Annotation Visualization and Integrated Discovery (DAVID) (https://david.ncifcrf.gov/ Version 6.8) software was used to analyze GO and KEGG pathway enrichment analysis of the potential activities of chosen genes (MDM4, MDMD2, and p53) (Huang et al., 2007). The cut-off value was set as P < 0.05.

2.8 Statistical analysis

IBM SPSS Statistics for Windows (version 23.0; IBM Corp. Armonk, NY, United States) was used to conduct the statistical analyses. The Kolmogorov-Smirnov test was used to evaluate the distribution of the variables. The examination of two independent variables with non-normal distributions was determined using the Mann-Whitney U test. On Medcalc software version 19.1 (Medcalc software BVBA, Belgium), receiver operating characteristic (ROC) curves were examined.

In this study, the values for the area under the curve (AUC) was used to evaluate the diagnostic accuracy of predictors. In comparing the prognostic accuracy for AML response with miR-342-5p expression levels, AUC value was calculated from the ROC curve using the Medcalc program. P values were two-sided, and statistical significance was defined as ≤0.05.

3 Results

To identify DE-miRNAs, a comparison between AML and normal control groups were performed by GEO2R, and a total of 250 miRNAs were identified between AML and normal samples. 24 De-miRNAs (16 upregulated and 8 downregulated) (Table 1) were identified from the GSE51908, using adjusted p < 0.05 and llogFCl>2 for DE-miRNAs selection. Volcano plot, visually demonstrating of DE-miRNAs is shown in Figure 1A. Box data normalization (Figure 1C) and scatter plots of miR-342-5p (Figure 1B) are shown in Figure 1.

FIGURE 1

The levels of miR-342-5p were significantly different between control and AML groups (p = 0.0001). To determine the power of miR-342-5p expression levels to distinguish patients with AML from control, the ROC analysis was performed. Values for AUC, sensitivity, specificity and cut-off points for miR-342-5p shown in Figure 2A. miR-342-5p showed significantly higher AUC values (AUC = 0.795, p = 0.0001). The cut-off of miR-342-5p was 6.85.

FIGURE 2

To further evaluate the clinical relevance of AML, the overall survival of AML patients was assessed by the Kaplan-Meier test using data from the article of Bullinger et al. (2010). Interestingly, as shown in Figure 2B, AML patients with low miR-342-5p expression had significantly shorter survival than those in higher expression control group (p = 0.0001), and its AUC value to discriminate AML from control groups was 0.795 (Figure 2A).

miR-342-5p in AML, potential target genes were predicted by the mirWALK database (miRDB and miRTarBase algorithms) (Table 2). Furthermore, the target genes of hsa-mir-342-5p at HL-60 and KG-1 cells were determined using miRDB (Figures 3A, B).

FIGURE 3

Bullinger et al. (2010).

Common target genes identified in HL-60 and KG-1 cells, and the mirWALK database were selected. Target gene expression values in AML were determined using GEPIA (LAML; Acute Myeloid Leukemia) (Table 3). And also, Multi-Gene comparison was performed according to AML and TCGA normal and GTEx data (Figure 4A). MDM4 expression level is increased in tumor tissue compared to healthy tissue (Figures 4B, C).

FIGURE 4

High expression of MDM4 predicts an unfavorable prognosis in AML patients. This highlights the important prognostic value of MDM4 in AML. The overall survival analysis of MDM4 is shown in Figure 5.

FIGURE 5

The PPI network of MDM4-related target genes was obtained from the STRING and GeneMANIA database (Figures 6A, B). STRING database demonstrated that the MDM4 gene primarily interacts with TP53, MDMD2, CNSK1A1, UBE2S2, UBB, UBC, USP7, CDKN2A, ATM, and CHEK2. GeneMANIA demonstrated that mdm4 interacts with MDM2, TP53, TP73,USP2,CBLC, HGF, ABL1, CD44, CDKN2A, USP7, SMARCD2, SMARCD3, SMARCD1, BECN1, MAP2K5, PELI1, GLIS2, CDKN2B, CDK4, and FAM193A. Among the genes, especially MDMD2 and p53, which are related to cell cycle and apoptosis, were selected.

FIGURE 6

The Gene Ontology and KEGG pathway analyses were performed to better understand the pathway and process affected by the identified genes. GO terms enriched by MDM4, MDM2, and p53 genes are shown in Table 4. Three terms biological process (BP), cellular component (CC), and molecular function (MF) are included in the GO ontology. The biological processes of these differentially expressed genes were primarily involved in the cellular response to hypoxia, macromolecular complex assembly, regulation of cell cycle, negative regulation of the apoptotic process, negative regulation of transcription from RNA polymerase II promoter, negative regulation of transcription, DNA-templated, positive regulation of transcription from RNA polymerase II promoter, and regulation of transcription from RNA polymerase II promoter. Meanwhile, the genes related to cellular components were mostly involved in transcriptional repressor complex, nucleoplasm, and nucleus. In terms of molecular function, these differential genes were mostly enriched in p53 binding, enzyme binding, zinc ion binding, metal ion binding, and protein binding.

In addition, the top 20 results of KEGG pathway enrichment analysis of MDM4, MDM2, and p53 genes are shown in Table 5. Among these pathways, the PI3K-Akt signaling pathway, microRNAs in cancer, proteoglycans in cancer, p53 signaling pathway, cellular senescence, and cell cycle are closely correlated with the development of AML.

4 Discussion

Acute myeloid leukemia is one of the most prevalent hematopoietic malignancies in adults (Short et al., 2018). Although significant advancements in the clinical treatment of AML, including stem cell transplantation and chemotherapies (Deng et al., 2017), the treatment is not at the desired level (Wang et al., 2022). The 5-year survival rate for AML patients is less than 15% (Zhao et al., 2015). Early diagnosis and timely intervention in patients with AML can improve their chances of survival. Cytogenetic and molecular analysis are very important in predicting the remission and survival rates of AML patients (Article et al., 2017). Understanding the molecular mechanisms underlying the pathogenesis of AML, identifying new biomarkers for diagnosis and prognosis, exploring new therapeutic strategies and improving the clinical outcome of AML is important (Journal et al., 2019; Article et al., 2017). Therefore, an improved understanding of the pathogenesis of AML will provide new insights into the diagnosis and treatment of AML. Recent studies have revealed that miRNAs are associated with the formation and development of AML, and play a role in the pathogenesis of AML by regulating their target genes (Fletcher et al., 2022; Wang et al., 2022). In order to comprehend the molecular mechanism of AML and develop new therapeutic targets for AML, it is necessary to identify unregulated miRNAs or mRNAs.

Microarray and High-throughput RNA sequencing (RNA-Seq) are widely used in cancer research to understand gene functions, find target molecules, and identify biomarkers for disease classification and diagnosis (Tang et al., 2017). Data obtained through experimental analysis are used in bioinformatics approaches. In particular, it is seen that bioinformatics data mining is a new trend in cancer research (Çakmak, 2022). Information on gene expression variations between normal and malignant tissues, chemical interactions, and the discovery of novel potential biomarkers for diagnosis and treatment can be obtained using bioinformatics tools (Çakmak, 2022).

The goal of the study is to identify potential target genes and important miRNAs for the diagnosis and therapy of AML by bioinformatics analysis. To the best of our knowledge, this is the first study to investigate the relationship between AML and miR-342-5p and target genes by bioinformatics analysis in the literature. Microarray data of GSE51908 were obtained from the GEO database to identify DE-miRNAs between AML and normal control (healthy). It was determined that several miRNAs were significantly differentially expressed, especially miR-342-5p′s expression was decreased in AML (Table 1; Figure 1). miRNA expression profiles are associated with prognosis, highlighting the potential significance of miRNAs in this disease (Garzon et al., 2008).

The Kaplan-Meier general database showed that MDM4 was significantly associated with the prognosis of AML patients (Figure 5). In addition, Kaplan-Meier analysis was performed to determine the survival times of the AML and control groups over miR-342-5p expression levels. Low miR-342-5p strongly predicted overall survival in AML patients. Moreover, Cox regression analysis showed that miR-342-5p expression independently provided prognostic information. Similarly, studies have shown that some miRs have a lower mortality rate and a huge effect on survival. Among the target genes of miR-10a-3p, SLC14A1, ARHGAP5, PIK3CA, it has been shown to show significant negative association with it, leading to longer survival in AML (Bullinger et al., 2010; Thanh et al., 2021). Similar to these studies, it has been shown that miR-342-5p is associated with the diagnosis and prognosis of patients with AML and has a positive effect on survival.

The findings of the ROC and overall survival studies show that miR-342-5p had diagnostic and prognostic significance in AML (Figure 2). The potential use of miRNAs as disease biomarkers has been the most common clinical use to far (Wallace and Connell, 2017). ROC analysis was performed to confirm whether the miR-342-5p can be used as one of the markers for AML. The miR-342-5p levels exhibited a significantly higher AUC value (AUC = 0.795, p = 0.0001), (Figure 2A). Based on this finding, it is concluded that miR-342-5p levels can be used as a reliable marker to evaluate AML. Previous studies revealed that plasma miR-92a, miR-143, and miR-342 could be promising biomarkers for AML (Article et al., 2017). In AML, miRNAs act as tumor suppressors or oncomiRs to affect a wide variety of leukemic processes, such as self-renewal, differentiation, proliferation, and epigenetic regulation. These molecules influence leukemic development and progression by targeting directly at the mRNA level or promoting malignancy (Wallace and Connell, 2017).

In addition to being a promising biomarker of miR-342-5p, it provides information about the molecular pathways involved in AML pathology. For this purpose, the target genes of miR-342-5p were determined using the miRwalk database. In addition, target genes identified in HL-60 and KG-1 cells and genes that are match as a result of mirwalk data were selected. Using the GEPIA database, expression analyses of these genes in AML were performed. In the GEPIA analysis, TCGA normal and database based on GTEx data in 173 AML patients and 70 healthy subjects, MDM4 expression level was higher in AML patients (Table 3; Figure 4). In addition, high expression levels of MDM4 were associated with poor overall survival (Figure 5).

Mouse Double Minute 4, which is among the target genes of miR-342-5p, has been reported to increase its expression in glioma, soft tissue sarcoma, melanoma, retinoblastoma, breast cancers and hematological malignancies (Li L. et al., 2014; Touqan et al., 2013; Furgason et al., 2015; Haupt et al., 2015; Carvajal et al., 2018). It has been stated that patients with AML have very high MDM4 expression and exhibit oncogenic activity compared to many other tumor types (Carvajal et al., 2018).

Functional enrichment analysis and PPI network were performed to identify the proteins that MDM4 is associated with in the development of AML and to reveal how these proteins function. According to the match data of the STRING, and GeneMANIA database, p53 and MDMD2 proteins were determined, especially in relation to apoptosis and cell cyle (Figures 6A, B). According to the literature, MDM4 inhibit p53 transcriptional activity by interacting with the p53 transactivation domain through its N-terminal domain (Tan B. X. et al., 2014). The tumor suppressor p53 protein, which guards the genome, triggers cell DNA repair, cell cycle arrest, senescence, and/or apoptosis pathways (Vousden and Lane, 2007). Downregulate p53, MDM2 and MDM4 are mutually dependent on one another (Badciong and Haas, 2002). MDM4 forms heterodimers with MDM2 via its RING domains and enhances the degradation of p53 by stimulating MDM2 E3 ubiquitin ligase activity (Wang et al., 2011). Loss of p53 function in stem and myeloid progenitor cells promotes the onset of leukemia (Zhao et al., 2010). In AML, TP53 mutations occur in less than 10% of patients (Kadia et al., 2016). Despite the TP53 mutations found in some subtypes of AML, MDM4 or MDM2 overexpression is the primary cause of p53 inactivation. Given the critical roles of MDM4 or MDM2, disruption of MDM4/MDM2/p53 interactions may offer a new way to restore p53 activity. Despite recent efforts to develop MDM4/MDM2 inhibitors, the desired success has not been achieved as a result of lack of affinity and dose-limiting toxicity (Carvajal et al., 2018). Therefore, targeting MDM4 with miRNA replacement therapy may be an attractive therapeutic strategy in AML.

Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis revealed that MDM4, MDM2, and p53 genes are associated with PI3K-Akt signaling pathway in cancer, MicroRNAs in cancer, Proteoglycans in cancer, p53 signaling pathway, cellular senescence and cell cycle. Previous studies have reported the involvement of PI3K-Akt signaling pathway (Bertacchini et al., 2015), p53 signaling (Chen and Sun, 2022), cellular senescence (Mao et al., 2022), cell cycle (Zhou et al., 2021), microRNAs in cancer (Trino et al., 2018), and proteoglycans in AML (Wu et al., 2020).

When the studies with miR-342 in the literature are examined; It has been suggested that dysregulation of the miR-342 Naa10p axis may be involved in the aggressive progression of pediatric AML. Naa10p-siRNA was reported to significantly suppress cell proliferation and increase cell apoptosis (Journal et al., 2019). Wu et al. showed that miR-342-5p is downregulated and has a significant inhibitory effect on cell proliferation in CML. miR-342-5p was shown to target the 3′-UTR domain of CCND1 and reduce its expression, and overexpression of miR-342-5p increased imatinib-induced DNA double-strand breaks and apoptosis (Wu et al., 2021). In another study, it was determined that miR-342-3p expression was decreased in the plasma of AML patients. Further studies have shown that miR-342-3p targets SOX12 and thus its mimics suppress AML cell growth, accelerate apoptosis and induce G0/G1 cell cycle arrest (Wang et al., 2022).

Given all these data, miR-342-5p is thought to act as a tumor suppressor. Given that the expression of miR-342-5p is reduced in AML, miRNA replacement therapy may provide an opportunity to restore the previously lost contribution of miR-342-5p. MiR-342-5p mimics can suppress cell proliferation in AML. miR-342-5p negatively regulates MDM4 expression, modulates DNA replication signaling pathways by contributing to the hsa-mir342-5p/MDM4/p53 network and may promote apoptosis. miR-342-5p and MDM4 are thought to be targets for the treatment of AML.

This study has a few limitations; only one microarray profile was analyzed as the database of miRNA profiles is rare. Experimental validation was not performed for the biological functions and targets of the identified miRNA. Further studies are needed for the mechanism of miR-342-5p and its potential targets in the development of AML. We believe that it can be a potential biomarker for the early diagnosis of AML patients and that cohort studies should be conducted on this topic.

miRNA and candidate genes associated with the pathogenesis and therapeutic targets of AML were identified using integrated bioinformatics approaches. A new miRNA-mRNA regulatory was axis identified that could improve understanding of the molecular mechanisms underlying the development of AML. These target molecules can serve as prognostic biomarkers and therapeutic targets. However, further studies are needed to confirm these predictive results.

References

• 1

  ArticleO.ElhamamsyA. R.SuleimanM.SharkawyE.ZanatyA. F.AhmedM.et al (2017). Circulating miR-92a, miR-143 and miR-342 in plasma are novel potential biomarkers for acute myeloid leukemia. Int. J. Mol. Cell Med.6 (2), 77–86. 10.22088/acadpub.BUMS.6.2.2

  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  BadciongJ. C.HaasA. L. (2002). MdmX is a RING finger ubiquitin ligase capable of synergistically enhancing Mdm2 ubiquitination. J. Biol. Chem.277 (51), 49668–49675. 10.1074/jbc.M208593200

  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  BaoM.JiangG. (2019). Differential expression and functional analysis of lung cancer gene expression datasets: a systems biology perspective. Oncol. Lett.18, 776–782. 10.3892/ol.2019.10362

  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  BertacchiniJ.HeidariN.MedianiL.CapitaniS.ShahjahaniM.AhmadzadehA.et al (2015). Targeting PI3K/AKT/mTOR network for treatment of leukemia. Cell Mol. Life Sci.72, 2337–2347. 10.1007/s00018-015-1867-5

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  BullingerL.EhrichM.DoK.SchlenkR. F.DoH.NelsonM. R.et al (2010). Quantitative DNA methylation predicts survival in adult acute myeloid leukemia. Blood115 (3), 636–642. 10.1182/blood-2009-03-211003

  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  ÇakmakE. (2022). A bioinformatics approach to identify potential biomarkers in non-small cell lung cancer. Cell Lung Cancer43 (1), 6–13. 10.17776/csj.976510

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  CandiaJ.CherukuriS.GuoY.DoshiK. A.CivinC. I.LosertW.et al (2015). Uncovering low-dimensional, miR-based signatures of acute myeloid and lymphoblastic leukemias with a machine-learning-driven network approach. Converg. Sci. Phys. Oncol.1 (2), 025002. 10.1088/2057-1739/1/2/025002

  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  CarvajalL. A.NeriahD. B.SenecalA.BenardL.ThiruthuvanathanV.YatsenkoT.et al (2018). Dual inhibition of MDMX and MDM2 as a therapeutic strategy in leukemia. Sci. Transl. Med.10 (436), eaao3003. 10.1126/scitranslmed.aao3003

  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  ChenY.SunL. (2022). Inhibition of NEDD8 NEDDylation induced apoptosis in acute myeloid leukemia cells via p53 signaling pathway. Biosci. Rep.42 (8), BSR20220994. 10.1042/BSR20220994

  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  ChenY.WangX. (2020). miRDB: an online database for prediction of functional microRNA targets. Nucleic Acids Res.48, D127–D131. 10.1093/nar/gkz757

  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  ChengY.SuY.WangS.LiuY.JinL.WanQ.et al (2020). Identification of circRNA-lncRNA-miRNA-mRNA competitive endogenous RNA network as novel prognostic markers for acute myeloid leukemia. Genes (Basel)11, 868. 10.3390/genes11080868

  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  DengL.JiangL.LinX.TsengK.LiuY.ZhangX.et al (2017). The PI3K/mTOR dual inhibitor BEZ235 suppresses proliferation and migration and reverses multidrug resistance in acute myeloid leukemia. Nat. Publ. Gr38, 382–391. 10.1038/aps.2016.121

  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  EdgarR.DomrachevM.LashA. E. (2002). Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res.30 (1), 207–210. 10.1093/nar/30.1.207

  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  EskandariM.ShiY.LiuJ.AlbaneseJ.GoelS.WangY.et al (2021). The expression of MDM2, MDM4, p53 and p21 in myeloid neoplasms and the effect of MDM2/MDM4 dual inhibitor. Leuk. Lymphoma62 (1), 167–175. 10.1080/10428194.2020.1817441

  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  ExpressionG.MyeloidA.SayedA.AhmedS. (2022). Clinical significance of MicroRNA-29a and MicroRNA-100 gene expression in pediatric acute myeloid leukemia. J. Pediatr. Hematol. Oncol.44 (2), 391–395. 10.1097/MPH.0000000000002168

  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  FletcherD.GuinnB.BrownE.JavadalaJ.Uysal-onganerP. (2022). microRNA expression in acute myeloid leukaemia: new targets for therapy. EJHaem3 (3), 596–608. 10.1002/jha2.441

  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  FurgasonJ. M.KoncarR. F.MichelhaughS. K.SarkarF. H.MittalS.SloanA. E.et al (2015). Whole genome sequence analysis links chromothripsis to EGFR, MDM2, MDM4, and CDK4 amplification in glioblastoma. Oncoscience2 (7), 618–628. 10.18632/oncoscience.178

  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  GaoY.ZhangS.WangZ.LiaoJ. (2017). Down-regulation of miR-342-3p in hepatocellular carcinoma tissues and its prognostic significance. Eur. Rev. Med. Pharmacol. Sci.21 (9), 2098–2102.

  • Google Scholar
  • View reference in article

• 19

  GarzonR.VoliniaS.LiuC.Fernandez-cymeringC.PalumboT.PichiorriF.et al (2008). MicroRNA signatures associated with cytogenetics and prognosis in acute myeloid leukemia. Blood111 (6), 3183–3189. 10.1182/blood-2007-07-098749

  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  GeneT.ConsortiumO.AshburnerM.BallC. A.BlakeJ. A.BotsteinD.et al (2011). Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat. Genet.25 (1), 25–29. 10.1038/75556

  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  HauptS.BuckleyD.PangJ.PanimayaJ.PaulP. J.GamellC.et al (2015). Targeting Mdmx to treat breast cancers with wild-type p53. Cell Death Dis.6 (7), e1821. 10.1038/cddis.2015.173

  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  HeL.HannonG. J. (2004). MicroRNAs: small RNAs with a big role in gene regulation. Nat. Rev. Genet.5 (7), 522–531. 10.1038/nrg1379

  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  HuangD. W.ShermanB. T.TanQ.KirJ.LiuD.BryantD.et al (2007). DAVID Bioinformatics Resources: expanded annotation database and novel algorithms to better extract biology from large gene lists. Nucleic Acids Res.35, 169–175. 10.1093/nar/gkm415

  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  JournalA. I.WangH.HeH.YangC. (2019). miR-342 suppresses the proliferation and invasion of acute myeloid leukemia by targeting Naa10p. Nanomedicine, Biotechnol.47 (1), 3671–3676. 10.1080/21691401.2019.1596930

  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  JuliussonG.LazarevicV.HoA.HagbergO.HoM.Swedish Acute Leukemia Registry Group (2012). Acute myeloid leukemia in the real world: why population-based registries are needed. Blood119 (17), 3890–3899. 10.1182/blood-2011-12-379008

  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  KadiaT. M.JainP.RavandiF.Garcia-maneroG.AndreefM.TakahashiK.et al (2016). TP53 mutations in newly diagnosed acute myeloid leukemia: clinicomolecular characteristics, response to therapy, and outcomes. Cancer122 (22), 3484–3491. 10.1002/cncr.30203

  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  KanehisaM.GotoS. (2000). KEGG: kyoto Encyclopedia of genes and genomes. Nucleic Acids Res.28 (1), 27–30. 10.1093/nar/28.1.27

  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  KirtoniaA.AshrafizadehM.ZarrabiA.HushmandiK.ZabolianA.BejandiA. K.et al (2022). Long noncoding RNAs: a novel insight in the leukemogenesis and drug resistance in acute myeloid leukemia. J. Cell Physiol.237 (1), 450–465. 10.1002/jcp.30590

  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  LiL.TanY.ChenX.XuZ.YangS.RenF.et al (2014). MDM4 overexpressed in acute myeloid leukemia patients with complex karyotype and wild-type TP53. PLoS One9 (11), e113088. 10.1371/journal.pone.0113088

  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  LiQ.LozanoG. (2013). Molecular pathways: targeting Mdm2 and Mdm4 in cancer therapy. Clin. Cancer Res.19 (1), 34–41. 10.1158/1078-0432.CCR-12-0053

  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  LiX.ChuH.LvT.WangL.KongS.DaiS. (2014). miR-342-3p suppresses proliferation, migration and invasion by targeting FOXM1 in human cervical cancer. FEBS Lett.588 (17), 3298–3307. 10.1016/j.febslet.2014.07.020

  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  LiuY.ChengZ.PangY.CuiL.QianT.QuanL.et al (2019). Role of microRNAs, circRNAs and long noncoding RNAs in acute myeloid leukemia. J. Hematol. Oncol.12 (1), 51. 10.1186/s13045-019-0734-5

  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  LvC.SunL.GuoZ.LiH.KongD.XuB.et al (2018). Circular RNA regulatory network reveals cell – cell crosstalk in acute myeloid leukemia extramedullary infiltration. J. Transl. Med.16, 361. 10.1186/s12967-018-1726-x

  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  MaoY.XuJ.XuX.QiuJ.HuZ.JiangF.et al (2022). Comprehensive analysis for cellular senescence-related immunogenic characteristics and immunotherapy prediction of acute myeloid leukemia. Front. Pharmacol.13, 987398. 10.3389/fphar.2022.987398

  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  O’BrienJ.HayderH.ZayedY.PengC. (2018). Overview of microRNA biogenesis, mechanisms of actions, and circulation. Front. Endocrinol. (Lausanne)9 (AUG), 402–412. 10.3389/fendo.2018.00402

  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  OtmaniK.LewalleP. (2021). Tumor suppressor miRNA in cancer cells and the tumor microenvironment: mechanism of deregulation and clinical implications. Front. Oncol.15 (11), 708765. 10.3389/fonc.2021.708765

  • CrossRef
  • Google Scholar
  • View reference in article

• 37

  ShengH.PanH.YaoM.XuL.LuJ.LiuB.et al (2021). Integrated analysis of circular RNA-associated ceRNA network reveals potential circRNA biomarkers in human breast cancer. Comput. Math. Methods Med.2021, 1732176. 10.1155/2021/1732176

  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  ShortN. J.RyttingM. E.CortesJ. E. (2018). Acute myeloid leukaemia. Lancet.392 (10147), 593–606. 10.1016/S0140-6736(18)31041-9

  • CrossRef
  • Google Scholar
  • View reference in article

• 39

  StichtC.De La TorreC.ParveenA.GretzN. (2018). miRWalk: an online resource for prediction of microRNA binding sites. PLoS One13, e0206239. 10.1371/journal.pone.0206239

  • CrossRef
  • Google Scholar
  • View reference in article

• 40

  SzklarczykD.FranceschiniA.KuhnM.SimonovicM.RothA.MinguezP.et al (2011). The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored. Nucleic Acids Res.39, 561–568. 10.1093/nar/gkq973

  • CrossRef
  • Google Scholar
  • View reference in article

• 41

  TanB. X.KhooK. H.LimT. M.LaneD. P. (2014). High Mdm4 levels suppress p53 activity and enhance its half-life in acute myeloid leukaemia. Oncotarget5 (4), 933–943. 10.18632/oncotarget.1559

  • CrossRef
  • Google Scholar
  • View reference in article

• 42

  TanY. S.KimM.KingsburyT. J.CivinC. I.ChengW. (2014). Regulation of RAB5C is important for the growth inhibitory effects of MiR-509 in human precursor-B acute lymphoblastic leukemia. PLoS One.9 (11), e111777. 10.1371/journal.pone.0111777

  • CrossRef
  • Google Scholar
  • View reference in article

• 43

  TangZ.LiC.KangB.GaoG.LiC.ZhangZ. (2017). GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res.45, W98–W102. 10.1093/nar/gkx247

  • CrossRef
  • Google Scholar
  • View reference in article

• 44

  ThanhT.FriedrichV.KristyS.ChristophW. W.MelindaR.MolloyT. J.et al (2021). miR-10a as a therapeutic target and predictive biomarker for MDM2 inhibition in acute myeloid leukemia. Leukemia35, 1933–1948. 10.1038/s41375-020-01095-z

  • CrossRef
  • Google Scholar
  • View reference in article

• 45

  TouqanN.DiggleC. P.VergheseE. T.PerryS.HorganK.MerchantW.et al (2013). An observational study on the expression levels of MDM2 and MDMX proteins, and associated effects on P53 in a series of human liposarcomas. BMC Clin. Pathol.13 (1), 32. 10.1186/1472-6890-13-32

  • CrossRef
  • Google Scholar
  • View reference in article

• 46

  TrinoS.LamorteD.CaivanoA.LaurenzanaI.TagliaferriD.FalcoG.et al (2018). MicroRNAs as new biomarkers for diagnosis and prognosis, and as potential therapeutic targets in acute myeloid leukemia. Int. J. Mol. Sci.19 (2), 460. 10.3390/ijms19020460

  • CrossRef
  • Google Scholar
  • View reference in article

• 47

  VousdenK. H.LaneD. P. (2007). p53 in health and disease. Nat. Rev. Mol. Cell Biol.8, 275–283. 10.1038/nrm2147

  • CrossRef
  • Google Scholar
  • View reference in article

• 48

  WadeM.LiY.-C.WahlG. M. (2013). MDM2, MDMX and p53 in oncogenesis and cancer therapy. Nat. Rev. Cancer13 (2), 83–96. 10.1038/nrc3430

  • CrossRef
  • Google Scholar
  • View reference in article

• 49

  WallaceJ. A.ConnellR. M. O. (2017). MicroRNAs and acute myeloid leukemia: therapeutic implications and emerging concepts. Blood130 (11), 1290–1301. 10.1182/blood-2016-10-697698

  • CrossRef
  • Google Scholar
  • View reference in article

• 50

  WangJ.ZhouH.TuX.HeY.LiuQ.LiuQ.et al (2019). Prediction of competing endogenous RNA coexpression network as prognostic markers in AML. Aging (Albany NY)11 (10), 3333–3347. 10.18632/aging.101985

  • CrossRef
  • Google Scholar
  • View reference in article

• 51

  WangX.WangJ.JiangX. (2011). MdmX protein is essential for Mdm2 protein-mediated p53 polyubiquitination. J. Biol. Chem.286 (27), 23725–23734. 10.1074/jbc.M110.213868

  • CrossRef
  • Google Scholar
  • View reference in article

• 52

  WangY.GuoX.WangL.XingL.ZhangX.RenJ. (2022). miR-342-3p inhibits acute myeloid leukemia progression by targeting SOX12. Oxid. Med. Cell Longev.2022, 1275141. 10.1155/2022/1275141

  • CrossRef
  • Google Scholar
  • View reference in article

• 53

  Warde-farleyD.DonaldsonS. L.ComesO.ZuberiK.BadrawiR.ChaoP.et al (2010). The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res.38, 214–220. 10.1093/nar/gkq537

  • CrossRef
  • Google Scholar
  • View reference in article

• 54

  WuJ.XiaoY.SunJ.SunH.ChenH.ZhuY.et al (2020). A single-cell survey of cellular hierarchy in acute myeloid leukemia. J. Hematol. Oncol.13 (1), 128. 10.1186/s13045-020-00941-y

  • CrossRef
  • Google Scholar
  • View reference in article

• 55

  WuY.LaiH.HuangT.ChenY.YeR.ChangP.et al (2021). Aberrantly reduced expression of miR-342-5p contributes to CCND1-associated chronic myeloid leukemia progression and imatinib resistance. Cell Death Dis.12 (10), 908. 10.1038/s41419-021-04209-2

  • CrossRef
  • Google Scholar
  • View reference in article

• 56

  XinX.XuZ.WeiJ.ZhangY. (2022). MiR-376a-3p increases cell apoptosis in acute myeloid leukemia by targeting MT1X. Cancer Biol. Ther.23 (1), 234–242. 10.1080/15384047.2022.2054243

  • CrossRef
  • Google Scholar
  • View reference in article

• 57

  ZhaoL.ZhangY. (2015). miR-342-3p affects hepatocellular carcinoma cell proliferation via regulating nF-κB pathway. Biochem. Biophys. Res. Commun.457 (3), 370–377. 10.1016/j.bbrc.2014.12.119

  • CrossRef
  • Google Scholar
  • View reference in article

• 58

  ZhaoY.ZhangX.ZhaoY.KongD.QinF.SunJ.et al (2015). Identification of potential therapeutic target genes, key miRNAs and mechanisms in acute myeloid leukemia based on bioinformatics analysis. Med. Oncol.32, 152. 10.1007/s12032-015-0572-4

  • CrossRef
  • Google Scholar
  • View reference in article

• 59

  ZhaoZ.ZuberJ.Diaz-floresE.LintaultL.KoganS. C.ShannonK.et al (2010). p53 loss promotes acute myeloid leukemia by enabling aberrant self-renewal. Genes Dev.24, 1389–1402. 10.1101/gad.1940710

  • CrossRef
  • Google Scholar
  • View reference in article

• 60

  ZhouC.DuJ.ZhaoL.LiuW.ZhaoT.LiangH.et al (2021). GLI1 reduces drug sensitivity by regulating cell cycle through PI3K/AKT/GSK3/CDK pathway in acute myeloid leukemia. Cell Death Dis.12 (3), 231. 10.1038/s41419-021-03504-2

  • CrossRef
  • Google Scholar
  • View reference in article