- Primary research
- Open Access
High expression of miR-25 predicts favorable chemotherapy outcome in patients with acute myeloid leukemia
Cancer Cell Internationalvolume 19, Article number: 122 (2019)
Acute myeloid leukemia (AML) pertains to a hematologic malignancy with heterogeneous therapeutic responses. Improvements in risk stratification in AML patients are warranted. MicroRNAs have been associated with the pathogenesis of AML.
To examine the prognostic value of miR-25, 162 cases with de novo AML were classified into two groups according to different treatment regimens.
In the chemotherapy group, cases with upregulated miR-25 expression showed relatively longer overall survival (OS; P = 0.0086) and event-free survival (EFS; P = 0.019). Multivariable analyses revealed that miR-25 upregulation is an independent predictor for extended OS (HR = 0.556, P = 0.015) and EFS (HR = 0.598, P = 0.03). In addition, allogeneic hematopoietic stem cell transplantation (allo-HSCT) circumvented the poor prognosis that was related to miR-25 downregulation with chemotherapy. The expression level pattern of miR-25 coincided with AML differentiation and proliferation, which included HOXA and HOXB cluster members, as well as the HOX cofactor MEIS1. The MYH9 gene was identified as a direct target of miR-25.
The miR-25 levels are correlated with prognosis in AML independently of other powerful molecular markers. The expression of miR-25 may contribute to the selection of the optimal treatment regimen between chemotherapy and allo-HCST for AML patients.
Acute myeloid leukemia (AML) is a group of clonal malignant diseases that derive from the hematopoietic stem cells. AML is characterized by a large group of germinal cells, which leads to a loss of normal hematopoietic function . The clinical prognosis of patients with AML is various. The differences in outcomes among AML patients depend on multiple intrinsic factors [2, 3]. With the development of methodologies of massive sequencing, it has been demonstrated that somatic mutations in NPM1, FLT3, CEBPA, IDH1 and IDH2 are connected to prognosis in AML . To be specific, patients with mutated FLT3 have a dismal outcome, while mutations in NPM1 and CEBPA are related with favorable prognosis. The advent of chemotherapy and allogeneic hematopoietic stem cell transplantation (allo-HSCT) has significantly improved AML treatment outcomes . Relapse and refractory of leukemia remain the most disturbing problems in AML patients . Thus, it is urge to explore more reliable and effective prognostic biomarkers to enhance the capacity of prediction and thus improve the outcome of AML by choosing optimal therapeutic approach.
MicroRNAs are short non-coding RNAs, which are implicated in a diverse group of critical cellular mechanisms, such as apoptosis, differentiation, cell cycle progression, and immune responses . Recently, more and more attentions have been focused on the prognostic role of microRNAs in AML. A recent study has shown that the upregulation of miR-181a facilitates better survival of AML patients who are cytogenetically normal . However, AML patients who are cytogenetically normal and upregulated miR-212 and miR-3151 have shorter overall and disease-free survival [9, 10]. However, most microRNA analyses did not differentiate the AML patients treated with chemotherapy and allo-HSCT. Thus, microRNAs may have varied prognostic roles in chemotherapy and allo-HSCT treatment group, respectively.
MiR-25, a member of miR-106b-25 cluster, is located on human chromosome 7q22.1 . Previous studies revealed that miR-25 was involved in many kinds of cancers . It has been identified that miR-25 is a potential biomarker for pediatric AML based on Pipeline of Outlier MicroRNA Analysis (POMA) model . More important, Garzon et al.  reported that miR-25 is significantly down-regulated in 122 newly diagnosed AML samples compared with CD34+ normal cells. However, clinical and prognostic role of miR-25 in AML are still unclear. A total of 162 recently diagnosed de novo AML patients were enrolled in this evaluation. The cases were placed into two groups based on the treatment that they received. The present study suggested that miR-25 is a solitary AML prognostic biomarker. Furthermore, our study revealed that allo-HSCT would be more beneficial to patients showing downregulated miR-25.
Patients and methods
Approximately 162 patients with a diagnosis of de novo AML were included in this study. The data sets used in this investigation were acquired from The Cancer Genome Atlas (TCGA). A single-institution tissue banking strategy endorsed by the human studies committee of Washington University was used in this study. All of the patients provided their written informed consent. AML diagnosis and classification were made according to the French–American–British (FAB) and the World Health Organization (WHO) criteria. The cases were placed into two groups according to the clinical treatment received. Ninety patients accepted chemotherapy, and the rest accepted allo-HSCT.
The samples from 155 patients both had been obtained mRNA and microRNA expression data. These data were applied to identify the mRNA-expression signature associated with miR-25 expression. The sequencing read count for each miRNA was normalized to Reads per million reads (RPM). The mRNA expression values were logged (base 2) prior to analysis . Spearman correlation was used to correlate the mRNA-expression profile with miR-25 expression. Hierarchical clustering analysis was used to reorder the gene rows. To screen for target genes of miR-25, Targetscan, miRNApath and miRDB website tools were implemented. Gene Ontology enrichment assessment of genes in miR-340 related signature was performed with the Database for Annotation, Visualization, and Integrated Discovery (DAVID).
The clinical endpoints of this investigation on treatment outcomes included overall survival (OS) and event-free survival (EFS). OS pertains to the time interval from diagnosis to death or last follow-up of the patient. EFS is described as the time interval from diagnosis to disease progression, relapse, or death attributed to any cause. The patients were assigned to the high or low expression groups based on the median miR-25 expression. Descriptive statistics (median and/or range) were used to summarize patients’ clinical and molecular characteristics. To elucidate the role of miR-25 expression in AML clinical and molecular features, the Pearson Chi-square and Fisher’s exact tests were used to screen for significant differences between two categorical variables. In addition, the Mann–Whitney’s U test was used for continuous variables. For univariable and multivariable analysis, a Cox proportional hazards model was employed to determine the effect of various risk factors on patient OS and EFS. The limited backward elimination procedure was applied to assess hazard ratios (HRs) and P values. Kaplan–Meier analysis was performed to determine the impact of miR-25 expression on OS and EFS. Statistical analysis was conducted with SPSS and GraphPad Prism. Differences among variables were determined to be statistically significant when the P value was < 0.05.
Correlation analysis of miR-25 expression and clinical characteristics
To establish the correlation among miR-25 expression and various clinical profiles, we assigned the patients who underwent chemotherapy and allo-HSCT to one of two groups according to median miR-25 expression levels, respectively. The associations of the clinical features with miR-25 expression levels are summarized in Table 1. In the chemotherapy group, subjects who exhibited upregulated miR-25 had a higher percentage of RUNX1-RUNX1T1 compared to those with downregulated expression (P = 0.026). In addition, high miR-25 expresser involved in more good risk cases of AML (P = 0.002). However, no significant differences were observed in gender, age, WBC count, BM blast, PB blast, FAB subtypes, FLT3-ITD, NPM1, DNMT3A, RUNX1, MLL-PTD, TP53, IDH1 and IDH2 among the high and low miR-25 expression group. In the allo-HSCT group, study participants with upregulated miR-25 exhibited a lower frequency for FLT3-ITD mutations (P = 0.045) compared to those with downregulated miR-25. No significant differences in as far as gender, age, WBC count, BM blast, PB blast, and mutations in the NPM1, RUNX1, DNMT3A, MLL-PTD, IDH1, IDH2, and TP53 genes were observed among the upregulated and downregulated miR-25 groups.
Prognostic value of miR-25 profiles in AML patients
We performed Kaplan–Meier analysis and log-rank test to assess the prognostic value of miR-25 profiles in AML patients. The chemotherapy group showed that AML patients with upregulated miR-25 were connected to better EFS (P = 0.019) and OS (P = 0.0086) relative to those with downregulated miR-25 (Fig. 1a, b). However, AML patients who received allo-HSCT did not exhibit any connection among prognosis and miR-25 expression (Fig. 1c, d). These findings revealed that miR-25 may be utilized as a chemotherapy-specific prognostic marker for AML.
High level of miR-25 is independently associated with favorable prognosis
To determine whether miR-25 expression could be used as an independent predictor for AML patient survival, we conducted univariate and multivariate Cox analyses. For the chemotherapy group, univariate analysis revealed that the upregulation of miR-25 was connected with longer EFS (HR = 0.598, 95% CI 0.376–0.951, P = 0.030) and OS (HR = 0.556, 95% CI 0.347–0.890, P = 0.015). Furthermore, multivariate cox analysis indicated that miR-25 upregulation was independently connected with longer EFS (HR = 0.561, 95% CI 0.333–0.943, P = 0.029) and OS (HR = 0.502, 95% CI 0.296–0.851, P = 0.011) after adjustment of mutation status for the FLT3-ITD, NPM1, DNMT3A, RUNX1, IDH1, and IDH2 genes and WBC count (Table 2).
Univariate analysis of the allo-HSCT group suggested that AML cases harboring FLT3-ITD mutations had shorter EFS (HR = 1.873, 95% CI 1.020–3.437, P = 0.043) and OS (HR = 1.998, 95% CI 1.053–3.788, P = 0.034). Patients with mutations only in the RUNX1 gene exhibited shorter OS (HR = 2.253, 95% CI 1.046–4.849, P = 0.038). Multivariate analysis indicated that FLT3-ITD and RUNX1 remained independent outcome predictors after adjusting for all other prognostic factors (Table 3). However, allo-HSCT patients did not show any significant differences between upregulated and downregulated miR-25 expression.
Allo-HSCT may circumvent poor patient outcomes that are related to downregulated miR-25 expression
To determine whether allo-HSCT therapy could circumvent the severe prognosis that was associated with downregulated miR-25, the whole cohort of 162 cases was split into two groups according to the median miR-25 expression levels. In the downregulated miR-25 group, the AML cases who received allo-HSCT showed significantly longer EFS (HR = 0.515, 95% CI 0.327–0.831, P = 0.0069) and OS (HR = 0.405, 95% CI 0.250–0.639, P = 0.0002) relative to cases who underwent standard chemotherapy alone (Fig. 2a, b). For the upregulated miR-25 group, no obvious differences in EFS (P = 0.969) and OS (P =0.364) were observed among the allo-HSCT and chemotherapy regimens. Thus, the AML patients showing downregulated miR-25 may benefit from treatment with allo-HSCT.
Biological insights into miR-25 profiles in AML
To generate insights into the molecular mechanism of miR-25, we analyzed a gene expression signature that was connected with miR-25 expression among AML cases. An association between the expression of 205 genes and miR-25 was observed. Among these genes, 145 were negatively correlated and 60 were positively correlated with the expression of miR-25 (Fig. 3). MiR-25 expression was inversely correlated with the expression of HOXA and HOXB, as well as the HOX cofactor MEIS1. Notably, these genes are crucial for the leukemogenesis and self-renewal capacities of AML [8, 16, 17]. Furthermore, we discovered that the expression of miR-25 was negatively connected with the levels of the PRDM16, Which involved in AML translocation ; CD97, an EGF-TM7 receptor ; IRAK1, which activates NF-κB pathways by the interaction with TRAF6 ; NFKB2, a pro-inflammatory response gene ; MYH9, which predicts unfavorable outcome of AML ; HDAC11, a epigenetic regulator. Notably, MYH9 was a in silico predicted target of miR-25.
Gene Ontology showed that genes involved in cellular metabolic process, system development, immune system process, transcription, hematopoietic or lymphoid organ development, hemopoiesis and myeloid cell differentiation were markedly overrepresented among differentially expressed genes associated with miR-25 expression (Table 4).
AML has been considered to occur as the result of genetic abnormalities, including chromosomal rearrangements, gene deregulations and mutations . The deregulated expression of microRNAs in AML can influence cell proliferation, survival and hematopoietic differentiation . The association of microRNAs with prognosis in heterogeneous patients with AML is still largely unclear. In this evaluation, the upregulated of miR-25 was determined to be an independently favorable prognosticator of AML cases who were administered chemotherapy. Furthermore, allo-HSCT may overcome the poor prognosis of AML cases with low miR-25 expression.
A correlation between aberrant miRNA expression and AML prognosis has been established [25, 26]. However, most of previous microRNA markers is restricted to AML without cytogenetic abnormalities. In our study, univariate and multivariate analyses demonstrated that miR-25 is an independently biomarker for cases administered chemotherapy. High miR-25 expression can predict favorable outcome. The prognostic role of miR-25 is different with previously established prognostic factors in a heterogeneous population of AML. MiR-25, as an independent outcome predictor, may improve the current clinical risk-based classification of patients with AML.
To further understand the biological insight into the molecular mechanism underlying miR-25, we identified genes significantly correlated with miR-25 expression. We discovered that the expression of miR-25 negatively connected with the levels of PRDM16, HOXAs, HOXBs, MEIS1, CD97, IRAK1, NFKB2 and MYH9. HOXA and HOXB gene clusters are the common characters of AML [27, 28]. Of these genes, HOXB4 is positively involved in the renewal of hematopoietic stem cell [29, 30]. A previous study has shown that HOXA9 contributes to the proliferation, apoptosis, and differentiation processes of leukemia . In addition, HOXA9 has been correlated with poor AML prognosis . Prior evaluations have revealed that IRAK1 may be utilized as a therapeutic target for AML, and TRAF6 may be used to activate pathways such as NFKB, MAPK, and AKT [20, 33]. PRDM16, also known as MEL1, is highly homologous to MDS1/EVI1. High expression of PRDM16 can predict the adverse outcome of AML . Moreover, MYH9 has also been predicted as a direct target of miR-25. High expression of MYH9 can induce resistant to chemotherapy and predict poor clinical outcome in AML . Taken together, the miR-25-associated gene-expression profiling analyses provide insights into the leukemogenic role of genes that are either direct or indirect targets of miR-25. Therefore, the miR-25-associated gene-expression signature analysis give novel insights into the oncogenic role of these genes. These miR-25-related genes could contribute to the chemotherapeutic responses of AML patients.
The FMS-like tyrosine kinase 3 (FLT3) gene is pivotal to hematopoietic stem cell proliferation and differentiation . FLT3 mutations take a great account of most frequent genetic aberrations in AML . FLT3-ITD mutation is one of FLT3 mutations, which can keep the tyrosine kinase persistently active, and result in the abnormal proliferation of leukemic cells. Mutations in the FLT-ITD gene have been associated with higher risk for relapse and poor OS and EFS . Consistent with the conclusion, our data suggested that FLT3-ITD mutation is a poor outcome marker in patients undergoing allo-HSCT. These analysis results indicate that allo-HSCT cannot overcome all adverse prognosis of molecular markers. The findings of this study have revealed that allo-HSCT circumvents the poor chemotherapy outcomes that are related to downregulated miR-25 expression. Thus, low miR-25 expression may be employed as a predictor of adverse prognoses among patients who received chemotherapy, as well as identify patients who require strategies in selecting the best treatment regimen, i.e., chemotherapy and/or allo-HCST.
In conclusion, high expression of miR-25 was identified to independently predict favorable survival in a highly heterogeneous population of patients with AML. Our findings may offer more information for the therapeutic strategies and the prediction of patients with AML, which may improve the survival and reduce the relapse of them. More importantly, allo-HSCT circumvents poor chemotherapeutic outcomes in cases with downregulated miR-25. The expression levels of miR-25 may thus be utilized in determining whether chemotherapy or allo-HSCT is the optimal treatment regimen for a specific AML patient.
Availability of data and materials
The datasets supporting the conclusions of this article are included within this article. Clinical data for all patients, including the treatment approach and outcomes data, are publicly accessible from the TCGA website (https://cancergenome.nih.gov).
acute myeloid leukemia
The Cancer Genome Atlas
World Health Organization
allogeneic hematopoietic stem cell transplantation
Ding L, Ley TJ, Larson DE, Miller CA, Koboldt DC, Welch JS, et al. Clonal evolution in relapsed acute myeloid leukaemia revealed by whole-genome sequencing. Nature. 2012;481(7382):506–10.
Papaemmanuil E, Gerstung M, Bullinger L, Gaidzik VI, Paschka P, Roberts ND, et al. Genomic classification and prognosis in acute myeloid leukemia. N Engl J Med. 2016;374(23):2209–21.
Niu M, Shen Y, Qi J, Liu X, Sang W, Wu Q, et al. Effects of realgar (As4S4) on degradation of PML-RARA harboring acquired arsenic-resistance mutations. Ann Hematol. 2017;96(11):1945–8.
Bullinger L, Dohner K, Dohner H. Genomics of acute myeloid leukemia diagnosis and pathways. J Clin Oncol. 2017;35(9):934–46.
Chen X, Xie H, Wood BL, Walter RB, Pagel JM, Becker PS, et al. Relation of clinical response and minimal residual disease and their prognostic impact on outcome in acute myeloid leukemia. J Clin Oncol. 2015;33(11):1258–64.
Khwaja A, Bjorkholm M, Gale RE, Levine RL, Jordan CT, Ehninger G, et al. Acute myeloid leukaemia. Nat Rev Dis Primers. 2016;2:16010.
Nowek K, Sun SM, Dijkstra MK, Bullinger L, Dohner H, Erkeland SJ, et al. Expression of a passenger miR-9* predicts favorable outcome in adults with acute myeloid leukemia less than 60 years of age. Leukemia. 2016;30(2):303–9.
Schwind S, Maharry K, Radmacher MD, Mrozek K, Holland KB, Margeson D, et al. Prognostic significance of expression of a single microRNA, miR-181a, in cytogenetically normal acute myeloid leukemia: a Cancer and Leukemia Group B study. J Clin Oncol. 2010;28(36):5257–64.
Sun SM, Rockova V, Bullinger L, Dijkstra MK, Dohner H, Lowenberg B, et al. The prognostic relevance of miR-212 expression with survival in cytogenetically and molecularly heterogeneous AML. Leukemia. 2013;27(1):100–6.
Eisfeld AK, Marcucci G, Maharry K, Schwind S, Radmacher MD, Nicolet D, et al. miR-3151 interplays with its host gene BAALC and independently affects outcome of patients with cytogenetically normal acute myeloid leukemia. Blood. 2012;120(2):249–58.
Qu J, Li M, Zhong W, Hu C. Prognostic role of microRNA-25 in cancers: evidence from a meta-analysis. Int J Clin Exp Med. 2015;8(8):12921–7.
Wang M, Yang YO, Jin Q, Shang L, Zhang J. Function of miR-25 in the invasion and metastasis of esophageal squamous carcinoma cells and bioinformatical analysis of the miR-106b-25 cluster. Exp Ther Med. 2018;15(1):440–6.
Yan W, Xu L, Sun Z, Lin Y, Zhang W, Chen J, et al. MicroRNA biomarker identification for pediatric acute myeloid leukemia based on a novel bioinformatics model. Oncotarget. 2015;6(28):26424–36.
Garzon R, Volinia S, Liu CG, Fernandez-Cymering C, Palumbo T, Pichiorri F, et al. MicroRNA signatures associated with cytogenetics and prognosis in acute myeloid leukemia. Blood. 2008;111(6):3183–9.
Zhang H, Zhang N, Wang R, Shao T, Feng Y, Yao Y, et al. High expression of miR-363 predicts poor prognosis and guides treatment selection in acute myeloid leukemia. J Transl Med. 2019;17(1):106.
Argiropoulos B, Humphries RK. Hox genes in hematopoiesis and leukemogenesis. Oncogene. 2007;26(47):6766–76.
Eklund EA. The role of HOX genes in malignant myeloid disease. Curr Opin Hematol. 2007;14(2):85–9.
Yamato G, Yamaguchi H, Handa H, Shiba N, Kawamura M, Wakita S, et al. Clinical features and prognostic impact of PRDM16 expression in adult acute myeloid leukemia. Genes Chromosomes Cancer. 2017;56(11):800–9.
Wobus M, Bornhauser M, Jacobi A, Krater M, Otto O, Ortlepp C, et al. Association of the EGF-TM7 receptor CD97 expression with FLT3-ITD in acute myeloid leukemia. Oncotarget. 2015;6(36):38804–15.
Beverly LJ, Starczynowski DT. IRAK1: oncotarget in MDS and AML. Oncotarget. 2014;5(7):1699–700.
Morikawa M, Nakano S, Mitsui N, Murasawa H, Masuki S, Nose H. Effects of dried tofu supplementation during interval walking training on the methylation of the NFKB2 gene in the whole blood of older women. J Physiol Sci. 2017;68:749–57.
Yu M, Wang J, Zhu Z, Hu C, Ma Q, Li X, et al. Prognostic impact of MYH9 expression on patients with acute myeloid leukemia. Oncotarget. 2017;8(1):156–63.
Grimwade D, Ivey A, Huntly BJ. Molecular landscape of acute myeloid leukemia in younger adults and its clinical relevance. Blood. 2016;127(1):29–41.
Dell’Aversana C, Giorgio C, D’Amato L, Lania G, Matarese F, Saeed S, et al. miR-194-5p/BCLAF1 deregulation in AML tumorigenesis. Leukemia. 2018;32(2):573.
Shivarov V, Dolnik A, Lang KM, Kronke J, Kuchenbauer F, Paschka P, et al. MicroRNA expression-based outcome prediction in acute myeloid leukemia: novel insights through cross-platform integrative analyses. Haematologica. 2016;101(11):e454–6.
Yang C, Shao T, Zhang H, Zhang N, Shi X, Liu X, et al. MiR-425 expression profiling in acute myeloid leukemia might guide the treatment choice between allogeneic transplantation and chemotherapy. J Transl Med. 2018;16(1):267.
Spencer DH, Young MA, Lamprecht TL, Helton NM, Fulton R, O’Laughlin M, et al. Epigenomic analysis of the HOX gene loci reveals mechanisms that may control canonical expression patterns in AML and normal hematopoietic cells. Leukemia. 2015;29(6):1279–89.
Kontro M, Kumar A, Majumder MM, Eldfors S, Parsons A, Pemovska T, et al. HOX gene expression predicts response to BCL-2 inhibition in acute myeloid leukemia. Leukemia. 2017;31(2):301–9.
Antonchuk J, Sauvageau G, Humphries RK. HOXB4-induced expansion of adult hematopoietic stem cells ex vivo. Cell. 2002;109(1):39–45.
Miyake N, Brun AC, Magnusson M, Miyake K, Scadden DT, Karlsson S. HOXB4-induced self-renewal of hematopoietic stem cells is significantly enhanced by p21 deficiency. Stem Cells. 2006;24(3):653–61.
Collins C, Wang J, Miao H, Bronstein J, Nawer H, Xu T, et al. C/EBPalpha is an essential collaborator in Hoxa9/Meis1-mediated leukemogenesis. Proc Natl Acad Sci USA. 2014;111(27):9899–904.
Collins CT, Hess JL. Role of HOXA9 in leukemia: dysregulation, cofactors and essential targets. Oncogene. 2016;35(9):1090–8.
Rhyasen GW, Bolanos L, Fang J, Jerez A, Wunderlich M, Rigolino C, et al. Targeting IRAK1 as a therapeutic approach for myelodysplastic syndrome. Cancer Cell. 2013;24(1):90–104.
Shlush LI, Zandi S, Mitchell A, Chen WC, Brandwein JM, Gupta V, et al. Identification of pre-leukaemic haematopoietic stem cells in acute leukaemia. Nature. 2014;506(7488):328–33.
Stone RM, Mandrekar SJ, Sanford BL, Laumann K, Geyer S, Bloomfield CD, et al. Midostaurin plus chemotherapy for acute myeloid leukemia with a FLT3 mutation. N Engl J Med. 2017;377(5):454–64.
Brunet S, Labopin M, Esteve J, Cornelissen J, Socie G, Iori AP, et al. Impact of FLT3 internal tandem duplication on the outcome of related and unrelated hematopoietic transplantation for adult acute myeloid leukemia in first remission: a retrospective analysis. J Clin Oncol. 2012;30(7):735–41.
The research was supported by National Natural Science Foundation of China (81670142, 81772658, 81870163); Jiangsu Provincial Key Research and Development Program (BE2017638, BE2017636); Natural Science Foundation of Jiangsu Province (BK20180104); The Foundation of Jiangsu Province Six Talents Peak (2017-WSN-120); Jiangsu Qing Lan Project for Mingshan Niu; Postgraduate Research & Practice Innovation Program of Jiangsu Province (SJCX18_0704, SJKY19_2110).
Ethics approval and consent to participate
All data sets were derived from The Cancer Genome Atlas (TCGA) data portal and are publicly accessible from the TCGA website. The patients were enrolled in a single-institution tissue banking protocol approved by the human studies committee at Washington University. Written informed consent was obtained from all patients.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.