- Research article
- Open Access
- Open Peer Review
Selection of three miRNA signatures with prognostic value in non-M3 acute myeloid leukemia
BMC Cancervolume 19, Article number: 109 (2019)
MiRNAs that are potential biomarkers for predicting prognosis for acute myeloid leukemia (AML) have been identified. However, comprehensive analyses investigating the association between miRNA expression profiles and AML survival remain relatively deficient.
In the present study, we performed multivariate Cox’s analysis and principal component analysis (PCA) using data from The Cancer Genome Atlas (TCGA) to identify potential molecular signatures for predicting non-M3 AML prognosis.
We found that patients who were still living were significantly younger at diagnosis than those who had died (P = 0.001). In addition, there was a marked difference in living status among different risk category groups (P = 0.022). A multivariate Cox model suggested that three miRNAs were potential biomarkers of non-M3 AML prognosis, including miR-181a-2, miR-25 and miR-362. Subsequently, PCA analyses were conducted to comprehensively represent the expression levels of these three miRNAs in each patient with a PCA value. According to the log-rank test, AML outcome for patients with lower PCA values was significantly different from those with higher PCA values (P < 0.001). Further bioinformatic analysis revealed the biological functions of the selected miRNAs.
We conducted a comprehensive analysis of TCGA non-M3 AML data, identifying three miRNAs that are significantly correlated with AML survival. PCA values for the identified miRNAs are valuable for predicting AML prognosis.
Acute myeloid leukemia (AML) is the most common acute leukemia in adult patients and is an aggressive haematological malignancy. It is a haematological malignancy that is promoted by various factors, including environmental exposure and several abnormality in cellular and molecular level . Despite decades of effort and achievements, the detailed mechanism of AML remains unclear. Studies have demonstrated that several biological factors may be involved in the pathogenesis and clinical outcomes of AML. As determining patient prognosis is important for individualized treatment, it is essential to comprehensively investigate biological factors that may be associated with AML outcomes and to select valuable biomarkers for prognosis prediction.
MicroRNAs (miRNAs), small noncoding RNA species approximately 22 nucleotides in length, have been demonstrated to play a predominant role in translational and post-transcriptional regulation. miRNAs participate in a wide range of biological processes, and their expression levels are frequently abnormal in human cancers . MiRNAs are involved in the pathogenesis of AML and influence the differentiation, proliferation of haematopoietic cells and clinical treatment response [3,4,5]. Expression levels of miRNA have been identified as potential biomarkers for predicting prognosis in various human cancers, including AML [6,7,8,9]. However, most studies have focused on one or few miRNAs, while synthesizing several miRNA signatures to comprehensively predict AML outcome has yet to be done.
In most cases of AML, patients at the same clinical stage may experience completely different clinical outcomes. Therefore, there is a critical need to explore novel biomarkers that facilitate the clinical prediction of AML prognosis. In 2013, The Cancer Genome Atlas (TCGA) Research Network reported an analysis of the genomes of 200 AML cases, along with RNA and microRNA sequencing data. They found that miRNA signatures were associated with unfavourable, intermediate, and favourable cytogenetic risk categories . Herein, we performed comprehensive analysis using data from TCGA, which contains200 AML patients and expression data for705 miRNAs. We sought to identify potential molecular signatures to predict non-M3 AML prognosis. In addition, bioinformatics analyses were performed to evaluate functions of the identified miRNAs. Our results aid in further characterizing mechanisms of non-M3 AML pathogenesis.
AML dataset in the TCGA database
MiRNA-seq data for 705 miRNAs and clinical data for AML patients were downloaded from the TCGA data portal (https://portal.gdc.cancer.gov/), an interactive data system for researchers to query, download, upload, and analyse harmonized cancer genomic data sets. Since M3 subtype has different biologic features and favourable outcome from other AML subtypes, we only included 179 non-M3 patients in the present analysis. Patients were treated in accordance with National Comprehensive Cancer Network (NCCN) guidelines (www.nccn.org). Many unfavourable risk and intermediate risk patients underwent allogeneic transplantation at some point in their disease course. RNA was extracted from whole blood, and patient samples were collected prior to treatment. Patients without miRNA expression data or clinical information were excluded. As a result, there were 147 non-M3 AML patients included in our study. Clinical demographics (e.g., gender, age, risk category) and survival data (overall and event-free survival) for these patients were also obtained from the TCGA dataset.
Multivariate survival analysis
We used clinical information, including age, gender, white blood cell (WBC) count to build a multivariate Cox proportional hazard ratio model. Clinical parameters were selected as significant using 0.05 as the cut off. miRNA with missing expression data > 75% were deleted for the following analyses. A multivariate Cox model was used to evaluate the association between expression level of each miRNA and prognosis of a patient. Subsequently, false discovery rate (FDR) adjustment was conducted, and significant miRNAs were filtered out with a cut off of 0.05. As a result, we screened out three miRNAs that were significantly related to survival time in AML as potential biomarkers of prognosis. Subsequently, One-way ANOVA test was performed to investigate whether there was significant difference of miRNA expression level among different sub-groups of AML.
Analyses of PCA and the association between miRNA expression and AML survival
To determine whether selected miRNAs were reliable indices for clinical prediction of AML prognosis, we conducted principal component analysis (PCA) to obtain a value for each patient based on the expression level of these miRNAs. In addition, the median of all the PCA values was used as a cut off to separate patients into one of two groups. We performed the log-rank test to investigate whether there was a significant difference between the survival of the two groups, and we used Kaplan-Meier plots to visualize survival curves.
Hierarchical clustering analysis
Patients were divided into two groups according to their survival status (i.e., non-censored patients were recognized as the poor-outcome group, while censored patients denoted the good-outcome group). We performed unsupervised hierarchical clustering analysis using the Cluster/Tree View program to give a general overview of the expression of the selected miRNAs in the good-outcome and poor-outcome groups. A heat map was created with expression levels of the three miRNAs and grouping information from AML patients to visualize the relationship between miRNA expression level and overall survival (OS).
Integrated gene expression
We also downloaded RNA-sequencing data for AML patients from the TCGA database. A total of 135 cases with both coding gene expression levels and miRNA expression data were included in subsequent correlation analyses. We performed a Pearson correlation for each miRNA and all the coding genes and then selected the top 100 coding genes with the highest Pearson correlation coefficient. The top 100 genes were considered potential targets of the miRNAs because their expression levels were most correlated with those of the miRNAs.
Gene ontology (GO) analysis
We used the online software DAVID (Database for Annotation, Visualization and Integrated Discovery, http://david.abcc.ncifcrf.gov/home.jsp) to analyse the molecular functions and biological pathways of the coding genes most correlated with the selected miRNAs (top100 of Pearson correlation coefficient). Names of the correlated coding genes were input into DAVID, and the results showed an annotation of input genes utilizing Gene Ontology (GO) terms (http://www.geneontology.org/).
Association between clinical factors and AML survival
Relative data were downloaded from TCGA database. Patients without miRNA expression data or clinical information were excluded. As a result, 147 non-M3 AML patients were included in our study. Among them, 96 were non-censored while the other 51 patients were censored. We subsequently investigated the association between patient clinical characteristics and survival status. As shown in Table 1, patients who were still living were significantly younger at diagnosis than those who were deceased (49.2 ± 14.4 VS. 58.2 ± 16.2, P = 0.001). In addition, there was a remarkable difference in living status among different groups of risk category (P = 0.022). However, we did not find any differences in the proportion of living patients between males and females (P = 0.932). Furthermore, the distribution of survivors among different FAB classifications of AML were not significantly different (P = 0.767).
Selection of miRNAs significantly associated with AML prognosis
After conducting multivariate Cox analyses of expression levels of all the miRNAs and clinical OS of AML cases, we selected three miRNAs as biomarkers of AML prognosis because their expression levels were significantly associated with OS time and survival status, with an FDR correlation less than 0.05. As shown in Table 2, higher levels of hsa.mir.362 indicated poorer prognosis (i.e. hazard ratio > 1), while higher levels of hsa.mir.25 and hsa.mir.181a.2 were related with better outcome (i.e. hazard ratio < 1). The workflow of carrying out the 3 miRNA signature was provided in Fig. 1 and a list of the top 50 miRNAs in the multivariate analysis after adjustment by FDR is shown in Additional file 1: Table S1.
In addition, we also found that expression of these miRNA was remarkable different among different AML subgroup (shown in Table 3 and Fig. 2), indicating a potential biological role of these selected miRNAs in development of AML.
MiRNA biomarkers and their association with AML survival status
After selection of the three miRNAs, unsupervised hierarchical clustering analysis was conducted based on expression patterns for the AML patients. We found that samples could be separated into two groups (i.e. good prognosis vs. poor prognosis) according to expression levels of the 3 miRNAs (Fig. 3).
Subsequently, we performed PCA analyses to comprehensively represent the expression levels of the three miRNAs of each patient with a PCA value. We subsequently separated all patients into one of two groups with a cut off set as the median of all PCA values. According to log-rank test, we found that the AML outcome for patients with higher PCA values was significantly different from those with lower PCA values (P < 0.001). OS of the two groups is graphical displayed by a Kaplan-Meier curve (Fig. 4). AML patients with lower PCA values, which were derived from comprehensive expression levels of the three selected miRNAs, exhibit significantly poorer prognosis than patients with lower PCA values. Furthermore, difference of PCA values between survival AML patients and dead patients was shown in Fig. 5 (P = 0.024), supporting a reliable role of PCA value in predicting AML prognosis.
Correlated genes and bioinformatic analyses
Pearson correlation analysis was performed using TCGA coding gene expression data as well as expression levels of the three miRNAs to identify the most correlated coding genes and understand the potential biological functions of these miRNAs. Top 5 GO terms (molecular function) for the top 100 correlated coding genes for each miRNA are shown in Fig. 6. According to these bioinformatic analyses, we inferred that these miRNAs may be involved in AML prognosis via the regulation of various biological process, including RNA binding (miR-181a-2 and miR-25), chromatin binding (miR-181a-2 and miR-25), protein binding (miR-181a-2, miR-25 and miR-362) etc.. Additional file 2: Table S2 lists the top 10 coding genes correlated with each miRNA. In order to test the reliability of relationship between miRNAs and correlated coding genes (i.e. top coding genes with expression level correlated with each miRNAs), we conducted bioinformatics analyses using Targetscan (http://www.targetscan.org/vert_72/). Finally, we confirmed that many of the TOP 10 genes have binding site for the correlated miRNAs. For instance, a total of 5 genes (i.e. FCF1, PAIP1, MTPAP, ZNF124 and RNASEH2B) listed in Additional file 2: Table S2 have binding site of miR-25.
Furthermore, enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of correlated coding genes are shown in Table 4. It should be noted that correlated genes of miR-181a-2 and miR-25 enriched in the same pathways (i.e. transcriptional misregulation in cancer, glutamatergic synapse and microRNAs in cancer).
In past decades, numerous studies have investigated biomarkers for the prognosis of cancers. Among various kinds of biomarkers, miRNAs have attracted increasing attention due to their stability in tumour tissue and the ease with which they can be detected in peripheral blood [11,12,13]. The pathogenesis and prognosis of many cancers has implicated the deregulation of miRNAs, which play vital roles in regulating the expression of coding genes involved in a wide range of biological processes [14, 15]. Therefore, several studies have investigated the clinical value of miRNAs for the prediction of cancer outcomes [9, 16, 17]. Koolivand et al. observed that miR-155 might function as an oncomir in AML and can be a prognosis biomarker for AML , while Xu H et al. have reported role for miR-135a in predicting poor prognosis in AML . We inferred from these studies that miRNAs can be valuable biomarkers for the prediction of AML prognosis. However, most current studies have focused on single miRNAs, while systematic analyses of association between whole miRNA expression profiles and AML outcomes remain relatively rare.
In the present study, we obtained data for miRNAs and coding gene expression microarrays in AML patients from the TCGA database. Clinical demographics, including age, gender, risk category, FAB morphology category, as well as OS time and event status, were also downloaded. Analyses showed that patient age and risk category were significantly related to survival status. Previous studies have demonstrated that elderly patients have poorer AML outcomes, partly due to their poorer tolerance to chemotherapy treatment and their reluctance to accept intensive regimens [20, 21]. Our results also showed that the percentage of living patients was remarkably lower in the older group than in the younger group. In addition, we found that the frequency distribution of living patients was significantly different among groups of AML risk categories, which were mainly distinguished by genetic abnormalities . Our results are in accordance with the current consensus that risk stratification of AML is an important prognostic factor [23, 24].
Subsequently, three miRNAs were identified to be significantly associated with AML survival. We subsequently performed PCA analyses to extract data for the expression levels of the three miRNAs, representing them as one PCA value for one each patient. Subsequently, analysis suggested that PCA values can separate AML into two groups with significantly different prognoses. Compared with most studies that identified single miRNAs associated with AML outcome [5, 25], our results are more comprehensive and illustrate the clinical value of a panel of three miRNAs that predict AML survival. In further studies, a scoring system will be constructed from expression level of these three miRNAs to predict prognosis of AML in a quantitative manner.
Recently, Chuang MK identified a 3-microRNA scoring system for prognostication of AML , and their list of miRNA signatures (i.e., miR-9, miR-155, miR-203) was different from ours, which were derived from the TCGA study . We speculated this difference was partly due to different populations (i.e., Chuang’s study was conducted in a Taiwanese population). Furthermore, differential analysis methods could also be a reason for the discrepancy.
In addition to this, Gao HY et al. have reported another work identifying some miRNAs associated with AML prognosis, using TCGA data. We have compared our present result with Gao HY’s report  and found that only miR-362was in common. However, patients of all AML subtypes were included in Gao’s study, while only non-M3 patients were included in ours because AML-M3 patients have remarkable different biologic features and good outcome from other AML subtypes. Besides, in Gao’ study, they drawn KM curve for each miRNA and miRNAs exhibiting significant differences were selected as candidate factors, while we performed multivariate Cox’s analysis to identify potential miRNAs associated with AML survival. Therefore, we inferred that different study populations as well as different analyse method may be the reason for differences between Gao’s study and ours. On the other hand, we have noticed that miR-362 was identified in both our and Gao’s study, with the same variation trend (i.e. higher expression level predicted poorer prognosis). After a literature reviewing, we found Jin J et al. have also reported high miR-362 expression is associated with poorer OS in AML . There observations provided strong evidence for the outcome predicting value of miR-362 in AML.
With respect to the other miRNAs identified to be valuable predictors of AML outcome, their biological role in cancer prognosis was also investigated through a literature review. As for miR-181a-2, a multiple centre study have demonstrated that higher miR-181a expression was significantly associated with better outcome of AML patients . In xenograft mouse models, ectopic miR-181a expression have also shown to inhibits tumor growth . This miRNA also plays determinant role in the sensitization of leukemic resistant cells to DNR and NK cell . Besides, miR-181a has shown to sensitize resistant leukemia HL-60/Ara-C cells to Ara-C by inducing apoptosis . While for miR-25, one study conducted by Wang Y et al. has reported that it was significantly associated with overall survival of 53 AML patients and its expression was associated with a good outcome . These studies provide further evidence for biological role of miR-181a-2 and miR-25 in AML progression.
Subsequently, bioinformatic analyses were performed to further reveal potential molecular functions and biological pathways for these three miRNAs. Firstly, we tested the reliability of relationship between miRNAs and correlated coding genes by conducting bioinformatics analyses using Targetscan, and further biological experiment will be conducted to confirm the association between other miRNAs and coding genes correlated with them. Then GO analysis revealed molecular activities of the genes most correlated with each miRNA. We found the correlated coding genes for the three miRNAs were relatively concentrated in RNA binding, chromatin binding and protein binding. Therefore, we can infer the biological function of the selected three miRNAs more rationally. In addition, pathway analyses also revealed that correlated coding gene expression for the identified miRNAs participated in transcriptional misregulation in cancer. Some of the potentially related pathways identified were interesting, such as mineral absorption. Previous studies have already shown that abnormalities in mineral homeostasis are involved with acute leukemia , but a detailed mechanism for this has not been demonstrated. It should be also noted that correlated genes of miR-181a-2 and miR-25 enriched in the same pathways, indicating a similar mechanism of the two miRNAs in pathogenesis of AML. Recently, Chang H et al. identified robust survival sub-pathways in AML using sample-matched miRNA and gene expression data sets from the TCGA database. Our results are different from their study partly because of the divergent study design, and the pathways we revealed may be a valuable clue for AML exploration. The results of the bioinformatics analyses provided us potential clues for further investigation of the biological function of the three identified miRNAs.
Although our present study have identified potential target genes and related pathways of these 3 miRNAs by bioinformatic approaches, in-depth mechanical approaches are still needed to verify our bioinformatic result. This is a limitation of our study. However, some of the top correlated coding genes have already been reported to be mutated in AML (e.g. CCND2, STMN1 and TFEB), which provided indirect evidence for our prediction that the selected 3 miRNAs may participate in AML by targeting these “top genes” [34,35,36]. In addition, the case numbers of TCGA data downloaded were too small to be split into training and validation cohorts. This is a limitation of our present study. In the future, we will collect more AML samples to create a validation cohort for our results.
In conclusion, the present study identified three miRNAs through comprehensive analysis of TCGA AML data that are significantly related with AML survival. PCA values of the identified miRNAs play important roles assessing the predictive values of these miRNAs for AML prognosis. Further studies focused on a scoring system to predict AML survival should be conducted.
Acute myeloid leukemia
Database for annotation, visualization and integrated discovery
False discovery rate
Kyoto Encyclopedia of Genes and Genomes
National Comprehensive Cancer Network
Principal component analysis
The Cancer Genome Atlas
White blood cell
Prada-Arismendy J, Arroyave JC, Rothlisberger S. Molecular biomarkers in acute myeloid leukemia. Blood Rev. 2017;31(1):63–76.
Sole C, et al. miRNAs in B-cell lymphoma: molecular mechanisms and biomarker potential. Cancer Lett. 2017;405:79–89.
Wallace JA, O’Connell RM. MicroRNAs and acute myeloid leukemia: therapeutic implications and emerging concepts. Blood. 2017;130(11):1290–301.
Khalaj M, et al. miR-99 regulates normal and malignant hematopoietic stem cell self-renewal. J Exp Med. 2017.
Favreau AJ, et al. miR-199b, a novel tumor suppressor miRNA in acute myeloid leukemia with prognostic implications. Exp Hematol Oncol. 2015;5:4.
Chuang MK, et al. A 3-microRNA scoring system for prognostication in de novo acute myeloid leukemia patients. Leukemia. 2015;29(5):1051–9.
Lubov J, et al. Meta-analysis of microRNAs expression in head and neck cancer: uncovering association with outcome and mechanisms. Oncotarget. 2017;8(33):55511–24.
Yuan GQ, et al. A 4-miRNAs signature predicts survival in glioblastoma multiforme patients. Cancer Biomark. 2017.
Wu HH, Lin WC, Tsai KW. Advances in molecular biomarkers for gastric cancer: miRNAs as emerging novel cancer markers. Expert Rev Mol Med. 2014;16:e1.
Cancer Genome Atlas Research, N, et al. Genomic and epigenomic landscapes of adult de novo acute myeloid leukemia. N Engl J Med. 2013;368(22):2059–74.
Souza MF, et al. Circulating mRNAs and miRNAs as candidate markers for the diagnosis and prognosis of prostate cancer. PLoS One. 2017;12(9):e0184094.
Elhamamsy AR, et al. Circulating miR-92a, miR-143 and miR-342 in plasma are novel potential biomarkers for acute myeloid leukemia. Int J Mol Cell Med. 2017;6(2):77–86.
Hamam R, et al. Circulating microRNAs in breast cancer: novel diagnostic and prognostic biomarkers. Cell Death Dis. 2017;8(9):e3045.
Wang Y, et al. microRNA-1297 inhibits the growth and metastasis of colorectal Cancer by suppressing cyclin D2 expression. DNA Cell Biol. 2017;36(11):991–9.
Alam F, et al. The role of p53-microRNA 200-Moesin axis in invasion and drug resistance of breast cancer cells. Tumour Biol. 2017;39(9):1010428317714634.
Izquierdo L, et al. Prognostic microRNAs in upper tract urothelial carcinoma: multicenter and international validation study. Oncotarget. 2017;8(31):51522–9.
Bailey ST, Westerling T, Brown M. Loss of estrogen-regulated microRNA expression increases HER2 signaling and is prognostic of poor outcome in luminal breast cancer. Cancer Res. 2015;75(2):436–45.
Koolivand M, et al. Alleviating the progression of acute myeloid leukemia (AML) by sulforaphane through controlling miR-155 levels. Mol Biol Rep. 2018.
Xu H, Wen Q. Downregulation of miR135a predicts poor prognosis in acute myeloid leukemia and regulates leukemia progression via modulating HOXA10 expression. Mol Med Rep. 2018;18(1):1134–40.
Thein MS, et al. Outcome of older patients with acute myeloid leukemia: an analysis of SEER data over 3 decades. Cancer. 2013;119(15):2720–7.
Heiblig M, et al. Acute myeloid leukemia in the elderly (age 70 yr or older): long-term survivors. Eur J Haematol. 2017;98(2):134–41.
Dohner H, et al. Diagnosis and management of AML in adults: 2017 ELN recommendations from an international expert panel. Blood. 2017;129(4):424–47.
Estey EH. Acute myeloid leukemia: 2014 update on risk-stratification and management. Am J Hematol. 2014;89(11):1063–81.
O’Donnell MR. Risk stratification and emerging treatment strategies in acute myeloid leukemia. J Natl Compr Cancer Netw. 2013;11(5 Suppl):667–9.
Butrym A, et al. Low expression of microRNA-204 (miR-204) is associated with poor clinical outcome of acute myeloid leukemia (AML) patients. J Exp Clin Cancer Res. 2015;34:68.
Gao HY, et al. Screening of prognostic risk microRNAs for acute myeloid leukemia. Hematology. 2018:1–9.
Ma QL, et al. MiR-362-5p as a novel prognostic predictor of cytogenetically normal acute myeloid leukemia. J Transl Med. 2018;16(1):68.
Schwind S, 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.
Hickey CJ, et al. Lenalidomide-mediated enhanced translation of C/EBPalpha-p30 protein up-regulates expression of the antileukemic microRNA-181a in acute myeloid leukemia. Blood. 2013;121(1):159–69.
Nanbakhsh A, et al. miR-181a modulates acute myeloid leukemia susceptibility to natural killer cells. Oncoimmunology. 2015;4(12):e996475.
Bai H, et al. miR-181a sensitizes resistant leukaemia HL-60/Ara-C cells to Ara-C by inducing apoptosis. J Cancer Res Clin Oncol. 2012;138(4):595–602.
Wang Y, et al. MicroRNAs expression signatures are associated with lineage and survival in acute leukemias. Blood Cells Mol Dis. 2010;44(3):191–7.
Atkinson SA, et al. Bone and mineral abnormalities in childhood acute lymphoblastic leukemia: influence of disease, drugs and nutrition. Int J Cancer Suppl. 1998;11:35–9.
Eisfeld AK, et al. Mutations in the CCND1 and CCND2 genes are frequent events in adult patients with t(8;21)(q22;q22) acute myeloid leukemia. Leukemia. 2017;31(6):1278–85.
Machado-Neto JA, et al. Stathmin 1 is involved in the highly proliferative phenotype of high-risk myelodysplastic syndromes and acute leukemia cells. Leuk Res. 2014;38(2):251–7.
Huan C, et al. Renal carcinoma-associated transcription factors TFE3 and TFEB are leukemia inhibitory factor-responsive transcription activators of E-cadherin. J Biol Chem. 2005;280(34):30225–35.
National Natural Science Foundation of China (81602913, 81670155) have supported design and data collection of the study; Natural Science Foundation of Jiangsu Province (BK20150085) have supported data collection; Nanjing Medical Science and Technology Development Project (YKK15140) have supported data analysis; Nanjing Medical Science and Technology Development Foundation (QRX17164) and Nanjing Science and Technology Development Foundation (201503043) have supported interpretation of data and writing of the manuscript; Technology Development Foundation of Nanjing Medical University (2014NJMUZD079) have also supported writing of the manuscript.
Availability of data and materials
Our present study obtained open data from the TCGA database. (https://cancergenome.nih.gov/).
Ethics approval and consent to participate
The study was approved by the Ethics Committee of Children’s Hospital of Nanjing Medical University and conformed to the provisions of the Declaration of Helsinki. Data of our present study was downloaded from an open database TCGA, so there was no informed consent from participants.
Consent for publication
The authors have no conflict of interest to report.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Yao Xue, Yuqiu Ge and Meiyun Kang are co-first authors.
Table S1. Title of data: List of top 50 miRNAs in the multivariate analysis after adjustment by FDR. Description of data: This Table listed TOP 50 miRNAs after multivariate analysis and the top 3 miRNAs were significantly associated with AML prognosis. (DOCX 18 kb)
Table S2. Title of data: Top ten coding genes correlated with each selected miRNAs. Description of data: This table listed top 10 coding genes associated with the 3 selected miRNAs. This correlation was analyzed by Pearson correlation and the expression level of top 10 coding genes was associated with 3 miRNAs with the highest Pearson correlation coefficient. (DOCX 16 kb)