- Open Access
High expression of miR-363 predicts poor prognosis and guides treatment selection in acute myeloid leukemia
Journal of Translational Medicinevolume 17, Article number: 106 (2019)
Acute myeloid leukemia (AML) is a highly heterogeneous malignancy with various outcomes, and therefore needs better risk stratification tools to help select optimal therapeutic options.
In this study, we identify miRNAs that could predict clinical outcome in a heterogeneous AML population using TCGA dataset.
We found that MiR-363 is a novel prognostic factor in AML patients undergoing chemotherapy. In multivariable analyses, high miR-363 remained predictive for shorter OS (HR = 2.349, P = 0.012) and EFS (HR = 2.082, P = 0.001) independent of other well-known prognostic factors. More importantly, allogeneic hematopoietic stem cell transplantation (allo-HSCT) overcame the adverse outcomes related to high miR-363 expression. In gene expression profiling, high miR-363 expression was positively correlated with the amounts of leukemogenic transcription factors, including Myb, RUNX3, GATA3, IKZF3, ETS1 and MLLT3. Notably, we found that the in silico predicted target genes (EZH2, KLF6 and PTEN) of miR-363 were downregulated in association with high miR-363 expression.
In summary, miR-363 expression may help identify patients in need of strategies to select the optimal therapy between chemotherapeutic and allo-HCST regimens. AML patients with high miR-363 expression may be highly recommended for early allo-HSCT regimen.
Acute myeloid leukemia (AML), the most frequent form of acute leukemia in adults, is caused by a rapid clonal proliferation of neoplastic myeloblasts . Patients with AML manifest complex and heterogeneous outcomes after receiving different treatments. Conventional cytotoxic treatment with chemotherapeutics is the first-line therapy in AML . High-risk patients could receive effectively antileukemic action and potential cure after accepting allogeneic hematopoietic stem cell transplantation (allo-HSCT). The variation in AML patient prognosis is related to various inherent factors, including cytogenetics and genetic alterations. Somatic mutations in NPM1, CEBPA, FLT3, IDH1, IDH2 and TET2 are associated with outcomes of patients and served as prognostic markers in AML . Despite the molecular mechanisms of leukemogenesis are well known, most AML patients are not cured. Notably, the currently available risk stratification systems are not completely accurate. Therefore, novel prognostic markers are needed to improve AML risk classification and select optimal therapeutic schedule.
MicroRNAs (miRNAs) represent short noncoding RNAs, which hybridize to target mRNAs with high specificity and decrease protein levels through translation inhibition . Dysregulation of miRNAs expression in AML can affect cell proliferation, survival and hematopoietic differentiation . More importantly, abnormal expression of miRNAs is related to clinical outcomes of AML patients. For instance, high miR-181a level has been confirmed to predict favorable survival in cytogenetically normal AML cases . Patients with high miR-212 level tend to have better outcome independently of cytogenetic subtype . Moreover, high miR-3151 expression is associated with worse overall survival and disease-free survival in patients cytogenetically normal AML . However, the majority of previous studies did not distinguish the various effects of chemotherapy and allo-HSCT treatment on the therapeutic outcome. As is well known, the prognostic impact of a marker is treatment-dependent in AML. Consequently, miRNAs may have different prognostic effects in chemotherapy and allo-HSCT treatment, respectively.
In this study, we identify miR-363 that could predict clinical outcome in a heterogeneous AML population using genome-wide screening. The prognostic role of miR-363 is independent of known potent clinical and molecular predictors. The miR-363 expression contributed to risk classification in AML patients undergoing chemotherapy. We also evaluated whether allo-HCST could overcome the poor prognostic effects of high miR-363 level in the same cohort. In order to evaluate biological insights of miR-363, we performed genome-wide gene and miRNA expression analyses.
Materials and methods
We studied a total of 162 patients with newly diagnosed AML according to the WHO classification. The RNA-Seq expression data of these AML patients were provided by The Cancer Genome Atlas (TCGA) . This study has been approved by Human Studies Committee of the Washington University. Patients with AML were included in a single center’s tissue protocol and followed NCCN guidelines to receive treatment. Patients with unfavorable risk underwent allo-HSCT if they were medically fit for the risks of transplantation, and if a suitably matched donor was available. In this cohort, 90 patients were only treated by chemotherapy and another 72 patients accepted both chemotherapy and allo-HSCT. All clinical data are available on the TCGA website.
Of the 162 patients, only 155 had both microRNA and mRNA expression data. For mRNA-seq data, genes expressed at or below a noise threshold of RPKM (Reads per kilobase per million mapped reads) ≤ 0.2 in at least 75% of samples were removed. For miRNA-seq data, read counts were normalized to RPM (Reads per million reads). The expression data were log2 transformed before analysis. The gene/microRNA expression signatures were derived by Spearman correlation analysis (Benjamini–Hochberg adjusted P value < 0.01). Finally, gene rows were reordered using hierarchical clustering analysis. The miRBase Targets Version 7 and Targetscan Release 7.1 were employed to predict the targets of miR-363. Gene Ontology enrichment assessment of genes in miR-363 related signature was conducted with the Database for Annotation, Visualization, and Integrated Discovery (DAVID).
A comparison of baseline characteristics between patients with high and low miRNA expression patients was conducted. The median miR-363 level was used to identify patients with low and high miRNA expression, respectively. Mann–Whitney U test was performed to test relations between two continuous variables. Fisher’s exact and Chi square tests were determined for categorical variables. Overall survival (OS) was the time from patient diagnosis to death at the final follow-up. Event-free survival (EFS) was the time from patient diagnosis to adverse events, including relapse and death. Kaplan–Meier method was performed to evaluate OS and EFS distributions and the log-rank test was employed to compare survival curves.
Univariable Cox proportional hazards models were constructed for assessing correlations of miR-363 expressions with OS and EFS, respectively. We establish multivariable Cox proportional hazards models to identify factors affecting OS and EFS. The factors included in the evaluation model contained miR-363 expression levels, FLT3–ITD, NPM1, DNMT3A, RUNX1, TP53, TET2, MLL–PTD, IDH1/IDH2 and NRAS/KRAS mutation statuses, and WBC involvement. Factors showing significance with α = 0.20 in univariable analysis were entered into limited backward selection to generate multivariable models. Variables remaining in the final models were significant at α = 0.05. The R software 3.1.5, GraphPad Prism and SPSS were used for statistical analysis, with P < 0.05 indicating statistical significance.
Association of miR-363 level with clinico-molecular properties
The patients were divided into two groups, chemotherapy and allo-HSCT groups. Subsequently, each group was subdivided into two groups in accordance with the median of miR-363. The relationship between clinical-genetic characteristics and miR-363 expression is shown in Table 1. In patients who underwent chemotherapy, cases with high miR-363 levels showed higher relapse rate (P = 0.001), and lower WBC count (P = 0.001) and circulating blast amounts (P = 0.007) at initial diagnosis in comparison with those expressing low miR-363 amounts. Patients with elevated miR-363 expression comprised less cases with favorable risk (P = 0.002), but more with poor risk of AML (P = 0.018). Furthermore, patients with high miR-363 expression included 92% of all cases with complex karyotypes and all cases with TP53 gene mutation. Meanwhile, Low miR-363 expressers encompassed all cases with the CBFβ-MYH11 fusion gene and 69% of all cases with NPM1 gene mutation.
Prognostic value of miR-363 expression in patients treated with chemotherapy or allo-HSCT
We performed genome-wide screening of miRNAs in AML cases in order to acquire prognostic marker to improve the classification of AML. MiR-363 was identified as a new prognostic marker for chemotherapy in AML patients. In order to evaluate survival of patients, we employed the Kaplan–Meier method and log-rank test. The expression level distribution of miR-363 was shown in Fig. 1a. In the chemotherapy group, cases highly expressing miR-363 showed reduced OS (HR = 2.28, P = 0.0004) and EFS (HR = 2.14, P = 0.0012) in comparison with low expressers (Fig. 1b). We further performed a survival analysis in the good/intermediate group, patients with high miR-363 expression had significantly shorter OS (P = 0.0009) and EFS (P = 0.0019) compared with patients with low miR-363 expression (Fig. 1c). However, miR-363 expression level was not associated with outcome in AML patients treated with allo-HCST (Fig. 1d). These data suggested that high expression of miR-363 was a poor prognostic factor in AML patients treated with chemotherapy.
MiR-363 is associated with clinical outcome in AML
Univariate and multivariate cox analyses were performed to assess whether miR-363 level is an independent predictor of survival in AML. Univariate analysis (Table 2) showed that high miR-363 had an adverse prognostic value for predicting OS (HR = 2.389, P < 0.001) and EFS (HR = 2.224, P = 0.001) in cases administered chemotherapy. In multivariable analysis, miR-363 and multiple demonstrated prognostic factors were assessed (Table 2). High miR-363 remained an independent predictor of shorter OS (HR = 2.349, 95% CI 1.305–4.229, P = 0.012) and EFS (HR = 2.082, 95% CI 1.172–3.699, P = 0.001).
In patients receiving allo-HSCT, univariate analysis showed that adverse OS in patients with TP53-mutant. However, miR-363 expression status was not associated with OS and EFS in the allo-HSCT group (Table 3). Multivariable analysis revealed that TP53 and FLT3–ITD mutations independently predict adverse OS (P = 0.002 and P = 0.049, respectively). The miR-363 expression status did not persist as OS and EFS predictors in multivariable analysis.
Allo-HSCT overcomes the adverse prognostic role of miR-363 expression
Next, we investigated whether allo-HSCT could overcome the adverse outcomes of miR-363 expression. The 162 patients were divided into 2 groups according to the median level of miR-363. In the high miR-363 group, cases administered allo-HSCT showed markedly improved OS (HR = 0.361, 95% CI 0.225–0.588, P < 0.0001) and EFS (HR = 0.447, 95% CI 0.287–0.751, P = 0.002) in comparison with cases administered chemotherapy (Fig. 2a). In patients with lower expression of miR-363, no marked differences in OS (P = 0.127) and EFS (P = 0.226) were found between the chemotherapy and allo-HSCT groups (Fig. 2b). These results suggested that miR-363 may be considered as a prognostic marker for the detection of patients requiring optimal therapeutic schedules.
Biological insight of miR-363 expression in AML
To further investigate the biological function of miR-363, gene expression signature associated with miRNA expression was determined in AML cases. We observed that the levels of 178 genes were strongly associated with miR-363 expression, including 130 and 48 with positive and negative correlations, respectively (Fig. 3). Differentially upregulated genes in patients with high miR-363 expression included leukemogenic transcription factors (Myb, RUNX3, GATA3, IKZF3, HMGA2 and ETS1) [10, 11]. Notably, MLLT3/AF9 was up-regulated in the high miR-363 group, as a frequent fusion partner of the MLL gene in translocations t(9;11)(p22;q23) related to AML . Among downregulated genes, we found that miR-363 expression showed negative correlations with the levels of tumor suppressor genes (EZH2, KLF6 and PTEN). Interestingly, these three genes were predicted miR-363 targets according to in silico analysis. Gene Ontology showed that genes associated with cell migration, T cell activation, system development, cell differentiation, response to chemicals and immune response were significantly correlated with miR-363 expression (Table 4). Thus, the miR-363 associated gene-expression profiling signature supported clinical finding in AML obtained by miRNA analysis.
As current molecular stratification schemes do not fully grasp the heterogeneity of prognosis in patients with AML, the identification of novel prognostic markers is urgent . In heterogeneous cohorts of AML patients, the correlation of miRNAs as predictive molecular markers remains largely unknown. In this study, miR-363 was determined as an independent prognostic factor in AML cases undergoing chemotherapy. Meanwhile, the miR-363 expression provides a powerful tool for risk stratification of AML patients. More importantly, allo-HSCT can overcome miR-363 expression-associated adverse outcomes.
We showed that miR-363 expression levels constitute independent prognostic marker of AML in a heterogeneous cohort administered chemotherapy. High miR-363 levels could still predict adverse outcome after consideration of other molecular prognostic factors in multivariable analysis. Thus, miR-363 could increase the prognostic value of previously defined molecular factors in a highly heterogeneous cohort of AML cases. Strikingly, patients with high miR-363 expression levels showed markedly poor OS and EFS. These findings suggest that miR-363 independently influences treatment outcomes and may synergistically drive leukemogenesis. More importantly, miR-363 expression levels could be useful to the identification of patients with adverse outcome in AML patients administrated chemotherapy.
Conventional chemotherapy and allo-HCST constitute the standard post-remission treatment strategies for AML . However, there is a lack of efficient prognostic markers for guiding rational treatment options. We found that high miR-363 expressers administered allo-HSCT showed markedly improved OS and EFS in comparison with cases administered chemotherapy. In cases lowly expressing miR-363, there was no advantage for those administered allo-HSCT in comparison with the chemotherapy group. These findings indicate that patients with low miR-363 expression may not benefit from allo-HSCT as first-line therapy. Therefore, the expression of miR-363 may contribute to identify patients in need of strategies to select the optimal treatment regimen between chemotherapy and allo-HCST. The AML patients with high miR-363 expression may be preferably recommended for early allo-HSCT.
The possible oncogenic function of miR-363 has been reported previously in T-cell acute lymphoblastic leukemia, multiple myeloma and solid tumors [15, 16]. MiR-363 promotes growth and chemo-resistance in gastric adenocarcinoma by downregulating FBW7 . It was shown that miR-363 is a prognostic marker for hepatocellular carcinoma . However, the function and prognostic role of miR-363 in AML remains unclear. To derive biological insights from AML cases characterized by high miR-363 expression, we identified genes associated with miR-363 expression in AML patients. Interestingly, miR-363 expression was positively correlated with the amounts of leukemogenic transcription factors, including Myb, RUNX3, GATA3, IKZF3, ETS1 and MLLT3. The Myb oncogene, a driver of leukemogenesis, is widely expressed in AML and important for continued proliferation and differentiation blocking activity in AML cells . ETS1 is critical in cell proliferation and differentiation in AML . MLLT3 represents a commonly encountered fusion partner of MLL in translocations t(9;11)(p22;q23), which are related to AML . Notably, we found that the direct target genes (EZH2, KLF6 and PTEN) of miR-363 were downregulated in association with high miR-363 expression. It was shown that loss-of-function mutations of the tumor suppressor gene EZH2 are found in AML . Meanwhile, PTEN plays an essential role in the prevention of leukemogenesis [23, 24]. Indeed, PTEN deletion in hematopoietic cells can induce a myeloproliferative disease within days and transplantable leukemias within weeks. These miR-363 associated genes may participate in the adverse response to chemotherapy in cases highly expressing miR-363. Therefore, the miR-363 related gene-expression profiling signature may support the clinical observation that AML is characterized by the expression of miRNA. However, the mechanisms concerning the regulation of miR-363 expression and subsequent influence of AML treatment outcome remain to be elucidated.
In conclusion, miR-363 levels independently correlate with clinical outcome in a highly heterogeneous cohort of AML cases. MiR-363 expression could greatly contribute to the identification of patients with poor outcome in AML. Expression analysis of miR-363 may be useful to improve the risk stratification of AML patients. Furthermore, allo-HSCT may overcome the unfavorable consequences of high miR-363 expression in AML. Therefore, the expression analysis of miR-363 may help identify cases in need of strategies to select the optimal treatment regimen between chemotherapy and allo-HCST.
acute myeloid leukemia
allogeneic hematopoietic stem cell transplantation
The Cancer Genome Atlas
Paguirigan AL, Smith J, Meshinchi S, Carroll M, Maley C, Radich JP. Single-cell genotyping demonstrates complex clonal diversity in acute myeloid leukemia. Sci Transl Med. 2015;7:281re282.
Dohner H, Estey E, Grimwade D, Amadori S, Appelbaum FR, Buchner T, et al. Diagnosis and management of AML in adults: 2017 ELN recommendations from an international expert panel. Blood. 2017;129:424–47.
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:2209–21.
Yeh CH, Moles R, Nicot C. Clinical significance of microRNAs in chronic and acute human leukemia. Mol Cancer. 2016;15:37.
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.
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: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: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:249–58.
Cancer Genome Atlas Research N, Ley TJ, Miller C, Ding L, Raphael BJ, Mungall AJ, et al. Genomic and epigenomic landscapes of adult de novo acute myeloid leukemia. N Engl J Med. 2013;368:2059–74.
Pattabiraman DR, Gonda TJ. Role and potential for therapeutic targeting of MYB in leukemia. Leukemia. 2013;27:269–77.
Stavast CJ, Leenen PJM, Erkeland SJ. The interplay between critical transcription factors and microRNAs in the control of normal and malignant myelopoiesis. Cancer Lett. 2018;427:28–37.
Stavropoulou V, Kaspar S, Brault L, Sanders MA, Juge S, Morettini S, et al. MLL-AF9 expression in hematopoietic stem cells drives a highly invasive AML expressing EMT-related genes linked to poor outcome. Cancer Cell. 2016;30:43–58.
Klco JM, Spencer DH, Miller CA, Griffith M, Lamprecht TL, O’Laughlin M, et al. Functional heterogeneity of genetically defined subclones in acute myeloid leukemia. Cancer Cell. 2014;25:379–92.
Cornelissen JJ, Versluis J, Passweg JR, van Putten WL, Manz MG, Maertens J, et al. Comparative therapeutic value of post-remission approaches in patients with acute myeloid leukemia aged 40-60 years. Leukemia. 2015;29:1041–50.
Landais S, Landry S, Legault P, Rassart E. Oncogenic potential of the miR-106-363 cluster and its implication in human T-cell leukemia. Cancer Res. 2007;67:5699–707.
Gowda PS, Wildman BJ, Trotter TN, Xu X, Hao X, Hassan MQ, et al. Runx2 suppression by miR-342 and miR-363 inhibits multiple myeloma progression. Mol Cancer Res. 2018;16:1138–48.
Zhang PF, Sheng LL, Wang G, Tian M, Zhu LY, Zhang R, et al. miR-363 promotes proliferation and chemo-resistance of human gastric cancer via targeting of FBW7 ubiquitin ligase expression. Oncotarget. 2016;7:35284–92.
Zhang J, Fan J, Zhou C, Qi Y. miR-363-5p as potential prognostic marker for hepatocellular carcinoma indicated by weighted co-expression network analysis of miRNAs and mRNA. BMC Gastroenterol. 2017;17:81.
Pattabiraman DR, McGirr C, Shakhbazov K, Barbier V, Krishnan K, Mukhopadhyay P, et al. Interaction of c-Myb with p300 is required for the induction of acute myeloid leukemia (AML) by human AML oncogenes. Blood. 2014;123:2682–90.
Lulli V, Romania P, Riccioni R, Boe A, Lo-Coco F, Testa U, et al. Transcriptional silencing of the ETS1 oncogene contributes to human granulocytic differentiation. Haematologica. 2010;95:1633–41.
Ney Garcia DR, Liehr T, Emerenciano M, Meyer C, Marschalek R, Pombo-de-Oliveira Mdo S, et al. Molecular studies reveal a MLL-MLLT3 gene fusion displaced in a case of childhood acute lymphoblastic leukemia with complex karyotype. Cancer Genet. 2015;208:143–7.
Ernst T, Chase AJ, Score J, Hidalgo-Curtis CE, Bryant C, Jones AV, et al. Inactivating mutations of the histone methyltransferase gene EZH2 in myeloid disorders. Nat Genet. 2010;42:722–6.
Yilmaz OH, Valdez R, Theisen BK, Guo W, Ferguson DO, Wu H, et al. Pten dependence distinguishes haematopoietic stem cells from leukaemia-initiating cells. Nature. 2006;441:475–82.
Zhang J, Grindley JC, Yin T, Jayasinghe S, He XC, Ross JT, et al. PTEN maintains haematopoietic stem cells and acts in lineage choice and leukaemia prevention. Nature. 2006;441:518–22.
HZ, NZ and RW designed and performed the computational analyses. TS, YF, YY, QW, JC, SZ, HZ and ZL contributed to statistical analyses. XL, MN and KX designed and wrote the manuscript. All authors read and approved the final manuscript.
The authors thank The Cancer Genome Atlas Network.
The authors declare that they have no competing interests.
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).
Consent for publication
Ethics approval and consent to participate
The research was supported by National Natural Science Foundation of China (81670142, 81772658, 81870163, 81700179); Jiangsu Provincial Key Research and Development Program (BE2017638, BE2017636); Natural Science Foundation of Jiangsu Province (BK20180104, BK20160226); The Foundation of Jiangsu Province Six Talents Peak (2017-WSN-120); Jiangsu Qing Lan Project for Mingshan Niu; Xuzhou Key Research and Application of Basic Research Projects (KC17158); Postgraduate Research & Practice Innovation Program of Jiangsu Province (SJCX18_0704, SJKY19_2110).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.