MiR-340 Is a Biomarker for Selecting Treatment Between Chemotherapy and Allogeneic Transplantation in Acute Myeloid Leukemia

Acute myeloid leukemia (AML) requires refined risk stratification tools to drive decisions concerning effective therapeutic strategies. Here, genome-wide screening was carried out for identifying miRNA molecules capable of predicting treatment outcome in AML patients based on the TCGA dataset. We identified miR-340 as a prognostic factor for selecting treatment between chemotherapy and allogeneic transplantation (allo-HSCT). In multivariable analyses, low miR-340 expression independently predicted reduced OS (HR = 2.07, P = 0.004) and EFS (HR = 1.909, P = 0.01) independent of other well-known prognostic factors. Meanwhile, allo-HSCT overcome deleterious outcomes related to low miR-340. Cases administered allo-HSCT showed markedly improved OS (HR = 0.316, P < 0.0001) and EFS (HR = 0.391, P = 0.002) in comparison with those receiving chemotherapy in the low miR-340 group. Gene expression assessment revealed that elevated miR-340 amounts were negatively correlated with HOXA/HOXB cluster levels, as well as the amounts of the HOX cofactor MEIS1. Strikingly, in silico analysis pointing to HOXA10, HOXB2, and MEIS1 as miR-340 targets. The miR-340 expression may help identify cases requiring strategies for selecting the optimal therapeutic option between chemotherapy and allo-HCST. AML cases showing low miR-340 levels should be strongly considered for early allo-HSCT treatment.


INTRODUCTION
Acute myeloid leukemia (AML), the commonest acute leukemia type affecting adult individuals, results from a fast clonal expansion of cancerous myeloblasts (1,2). AML has a complex and heterogeneous nature (3). Thus, it shows variable outcomes in affected individuals upon administration of current treatments. Cytotoxic chemotherapeutics are considered the mainstay therapeutic option for AML (4,5). Meanwhile, allogeneic hematopoietic stem cell transplantation (allo-HSCT) provides potent antileukemic effects and potential cure in high-risk patients (6,7).
In AML cases, prognosis is associated with multiple intrinsic parameters such as cytogenetic and genetic modifications. For example, NPM1, CEBPA, FLT3, IDH1, IDH2, and TET2 mutations show associations with patient outcome, and are used for AML prognosis (8)(9)(10). Although leukemogenesis is mechanistically well-understood, most AML cases remain uncured. Indeed, no currently available classification system is completely accurate. This suggests an urgent need for identifying new prognostic markers to further ameliorate AML stratification and select rational treatment options.
MicroRNAs (miRNAs) specifically bind mRNA targets and reduce the expression of encoded proteins by translation inhibition (11). Dysregulated miRNAs in AML affects cell growth and hematopoietic differentiation (12,13). Furthermore, aberrant expression of miRNAs is related to patient outcome (14). For instance, elevated miR-181a content predicts favorable outcome in AML with normal cytogenetic features (15). Patients expressing elevated miR-212 amounts tend to show longer survival in all cytogenetic subtypes (16). Meanwhile, overexpression of miR-3151 is related to decreased overall and disease-free survival in cytogenetically normal AML (17). Moreover, combination of miR-3151 with the associated host gene BAALC is more efficient in predicting worse outcome. However, most published reports did not distinguish the effects of chemotherapeutics and allo-HSCT on treatment outcome. It is generally accepted that a marker's value in AML prognosis depends on the therapeutic method administered. Thus, miRNAs could show distinct prognostic values between chemotherapy and allo-HSCT treatment.
Here, we performed genome-wide screening to identify miRNAs with significant prognostic values in AML cases administered chemotherapeutics. We identified miR-340 as a prognostic factor independently of reported clinicomolecular predictive factors in AML. In addition, whether all-HCST overcomes the adverse prognostic impact of miR-340 expression was assessed. To assess molecular mechanisms, the expression levels of genes and miRNAs throughout the genome were determined.

Patients
All the 162 cases assessed were newly diagnosed AML (non-M3 FAB subtypes) based on the WHO classification. RNA-Seq data for the above AML cases were provided by The Cancer Genome Atlas (Project ID: TCGA-LAML) (18). These AML cases were included in a single-institution tissue banking protocol that had approval from the human studies committee at Washington University (WU HSC #01-1014). All samples were collected from patients with previously untreated de novo AML. MicroRNA and mRNA levels were determined before treatment. The risk classification of AML is based on the NCCN recommendation which includes cytogenetics as well as gene mutational status.
Patient treatment followed the NCCN guidelines. Among the 162 cases, 148 (91.4%) received standard induction chemotherapy (3 and 7 days of anthracycline and cytarabine, respectively); the remaining 14 patients received decitabine or azacitidine. Patients with unfavorable risk underwent allo-HSCT if medically suitable for transplantation and there was a matched donor. Some cases with intermediate risk were also administered allo-HSCT during the disease course. A total of 90 cases were solely administered chemotherapy, while 72 also underwent allo-HSCT. Clinical data, such as treatment approach and outcomes, can be found on the TCGA website (Supplementary Table 1).

Gene-Expression Profiling
One hundred and fifty-five of the 162 cases were assessed for both miRNA and mRNA expression levels. Their specimens were employed to identify genes associated with miR-340 expression. In mRNA-seq, genes showing expression levels equal or inferior to a noise threshold of RPKM (Reads per kilobase per million mapped reads) in ≥75% of specimens were excluded. In miRNA-seq, RPM (Reads per million reads) was employed for normalization. Data underwent log2 transformation preanalysis. The mRNA levels were compared between the high and low miRNA groups, with a univariate significance level at 0.005 to identify differentially expressed genes (DEGs). Then, hierarchical clustering using the multiple experiment viewer (MeV) 4.9.0 software was carried out to order gene rows; miRBase Targets V7 and Targetscan Release 7.1 were used for predicting microRNA targets. Gene Ontology enrichment was performed of genes associated with miR-340 using the Database for Annotation, Visualization, and Integrated Discovery (DAVID).

Statistical Analysis
Baseline features were compared between the high and low miRNA groups. Median levels of miR-340 were employed to identify patients with low and high miRNA levels. The miR-340 median values in the chemotherapy and allo-HSCT groups were 82.1 (13.2-227.5) and 75.0 (18.9-286.1), respectively. The Mann-Whitney U-test was employed to assess associations of two continuous variables. Fisher's exact test or the chi-square test was employed for categorical variables. Overall survival (OS) was defined as time between case enrolment and death, censoring live cases at final follow-up. Event-free survival (EFS) was defined as time between case enrolment and induction failure, recurrence, or death. Estimated OS and EFS distributions were determined by Kaplan-Meier curves, with the log-rank test carried out for comparing survival data.
The univariate Cox proportional hazards model was employed for evaluating associations of 340 expression with OS and EFS. The multivariate Cox proportional hazards model was used for identifying factors that affect survival. The factors assessed for model inclusion comprised miR-340 expression amounts; FLT3-ITD, NPM1, DNMT3A, RUNX1, TP53, TET2, MLL-PTD, IDH1/IDH2, and NRAS/KRAS mutation statuses; age; karyotype; and WBC involvement. Parameters with p = 0.10 in univariate analysis were further assessed by multivariate analysis (backward selection), in which p < 0.05 indicated significance. R v3.1.5, GraphPad Prism and SPSS were employed for data analysis, and P < 0.05 indicated statistically significant differences.

Associations of miR-340 Levels With Clinico-Molecular Features
The chemotherapy and allo-HSCT groups were each subdivided into two according to median amounts of miRNAs, respectively. Tables 1, 2 summarize the associations of clinico-genetic features with miR-340 levels. In the chemotherapy group, cases with low miR-340 levels showed a higher rate of first relapse or death within 1 year after enrolment in comparison with those expressing high miR-340 amounts (82.2 vs. 55.6%, P = 0.012). There were no significant associations of other miR-340 expression distributions with clinical features. The above findings suggested that miR-340 expression may predict prognosis independently from known molecular features.

Prognostic Value of miR-340 in Cases Administered Chemotherapy or Allo-HSCT
To obtain prognostic markers for refining AML stratification, a genome-wide screen of miRNAs in AML patients was carried out; miR-340 was detected as a novel molecule predicting prognosis in AML cases treated by chemotherapy. In the chemotherapy group, low miR-340 expressers showed starkly reduced OS (P = 0.0005) and EFS (P = 0.0005) in comparison with high expressers ( Figure 1A). Notably, 5 year OS rates were 35.6 and 5.4% in the high and low miR-340 expression groups, respectively. However, miR-340 levels were not associated with patient outcome in AML Chi-square tests were used for categorical variables.
Frontiers in Oncology | www.frontiersin.org after allo-HCST ( Figure 1B). The above findings indicated that low miR-340 was an unfavorable prognostic factor in AML cases administered chemotherapy.

MiR-340 Is Independently Associated With Clinical Outcome in Patients With AML
To assess whether miR-340 levels independently predict survival in AML, univariable and multivariate cox analyses were carried out. Univariate analysis ( Table 3) revealed low miR-340 had a significant association with reduced OS (HR = 2.393, P = 0.001) and EFS (HR = 2.383, P = 0.001) in patients treated by chemotherapy. In multivariate analysis, miR-340 and many known prognostic factors were evaluated ( Table 3). Low miR-340 expression was independently predictive of reduced OS (HR = 2.07, P = 0.004) and EFS (HR = 1.909, P = 0.01) upon adjustment for age (P < 0.0001) and karyotype risk (P < 0.0001). In patients administered allo-HSCT, univariate analysis revealed TP53-mutant cases displayed adverse OS (P = 0.009); however, the miR-340 expression status was not associated with OS and EFS ( Table 4). Multivariable analysis showed TP53 and FLT3-ITD mutations were independently predictive of reduced OS (P = 0.003 and P = 0.031, respectively). Therefore, the miR-340 expression status did not remain a survival predictor after allo-HSCT.

Allo-HSCT Overcomes the Adverse Prognostic Effects of Low miR-340 Expression in AML
Whether allo-HSCT overcomes the unfavorable effects of low miR-340 expression was assessed. The 162 cases were assigned to two groups based on miR-340 levels. Then, each group was further divided into the chemotherapy and allo-HSCT groups. In cases with low miR-340 expression, significantly longer OS (HR = 0.316, 95% CI 0.167-0.459, P < 0.0001) and EFS (HR = 0.391, 95% CI 0.231-0.622, P = 0.002) were observed after allo-HSCT (Figure 2A). However, there were no differences in OS (P = 0.060) and EFS (P = 0.162) between the allo-HSCT and chemotherapy groups in high miR-340 expressers ( Figure 2B). The above findings indicated miR-340 might constitute a novel molecule for detecting cases requiring strategies to select a rational therapeutic management.

Biological Insights Into miR-340 Function in AML
To gain insights into the biological functions of miR-340, gene expression signatures related to the miRNA expression were assessed in AML patients. As shown in Figure 3, the expression levels of 135 genes showed significant associations with miR-340 expression, with 61 and 74 positively and negatively correlated, respectively. In this study, miR-340 expression was negatively associated with HOXA/HOXB cluster and MEIS1 amounts. The latter genes are important in AML leukemogenesis and self-renewal features in leukemic stem cells (19). We further analyzed the biological significance of these genes by separating AML cases into two groups according to median miR-340 expression. In the low miR-340 expression group, the HOXA1 (P = 0.0112), HOXA2 (P = 0.019), HOXA3 (P < 0.0001), HOXA4 (P = 0.0005), HOXA5 (P < 0.0001), HOXA6 (P < 0.0001), HOXA7 (P < 0.0001), HOXA9 (P = 0.0004), HOXA10 (P < 0.0001), HOXB2 (P = 0.0026), HOXB3 (P < 0.0001), HOXB4 (P < 0.0001), HOXB6 (P < 0.0001), MEIS1 (P = 0.0001), and PRDM16 (P = 0.0002) genes were all expressed at higher levels compared with high miR-340 expression group. Notably, the histone H3K4 methyltransferase PRDM16, a gene that predicts adverse outcome in AML (20), was downregulated in the high miR-340 group. Strikingly, HOXA10, HOXB2, MEIS1 and PRDM16 were identified as miR-340 targets by in silico prediction. These findings corroborated clinical data for AML generated by miRNA analysis. Gene Ontology revealed that genes contributing to biological regulation, cellular and developmental processes, transcription regulation, immune system process, cell apoptosis and proliferation, myeloid cell differentiation, hematopoietic organ development, cellular response to chemical stimulus and embryonic morphogenesis were remarkably overrepresented among DEGs related to miR-340 expression ( Table 5).

DISCUSSION
Identifying new prognostic markers is of great importance in adult AML, because currently available molecular stratification approaches do not completely capture the heterogeneity of such cases (21)(22)(23). The significance of miRNAs as predictive biomarkers in AML patients remains mostly unclear. The present work demonstrated that miR-340 independently predicted prognosis of AML patients administered chemotherapy, providing an important tool for risk-stratifying AML cases. Of greater importance, allo-HSCT overcame miR-340-associated unfavorable outcomes.
Routine therapeutic approaches after AML remission include conventional chemotherapy and allo-HCST (32). However, prognostic biomarkers efficiently guiding treatment selection are lacking. As demonstrated above, low miR-340 expressers treated by allo-HSCT had remarkably prolonged OS and EFS compared with the chemotherapy group. Notably, no benefit of allo-HSCT was found in cases highly expressing miR-340. The above results suggest that allo-HSCT as first line treatment might be of limited value in cases highly expressing miR-340. Thus, miR-340 levels   might help identify individuals requiring strategies for selecting the most suitable therapeutic option between chemotherapy and allo-HCST. Low miR-340 expressing cases may be strongly considered for early allo-HSCT. The tumor suppressor activity of miR-340 has been described in solid malignancies such as breast, colorectal and non-small cell lung cancers, as well as osteosarcoma and glioblastoma (33)(34)(35). However, miR-340's biological role in AML cells is largely undefined. As shown above, miR-340 levels were negatively correlated with HOXA/HOXB cluster and MEIS1 amounts. HOX genes are essential regulators of hematopoietic development and stem cell self-renewal (36). Abnormal overexpression of HOXA and MEIS1 represents a known hallmark of AML (19). Furthermore, artificial overexpression of HOXA7, HOXA9, or HOXA10 in combination with MEIS1 leads to rapid onset AML in animal models (37). High expression of HOXA9 in leukemia blasts predicts adverse outcome in AML. High expression levels of HOXA10 (38,39) and MEIS1 (40) genes also predict poor outcome in AML. Notably, HOXA10, HOXB2, and MEIS1 were predicted in silico to be miR-340 targets in this study. Although they are potential therapeutic targets, developing chemical inhibitors of transcription factor remains challenging. Targeting these genes with synthetic miR-340 compounds may be a potential therapeutic strategy. Taken together, the miR-340 associated gene signature may confirm the notion that AML features miRNA expression changes. However, the mechanisms involved in miR-340 modulation and associated effects on patient outcome in AML after treatment deserve further attention.

CONCLUSIONS
MiR-340 amounts are independently associated with treatment outcome in a highly heterogeneous population of AML patients. MiR-340 amount determination may provide a powerful tool for the identification of a patient subset with adverse outcome in AML, and might help ameliorate stratification risk and inform therapeutic decision making for AML patients. Interestingly, allo-HSCT may overcome miR-340-associated deleterious effects in AML. Therefore, miR-340 expression analysis could help identify patients requiring complex approaches in selecting the most suitable therapeutic option between chemotherapy and allo-HCST.

DATA AVAILABILITY STATEMENT
All datasets supporting the conclusions contained in the present report are included in the manuscript. All clinical data, including treatment approaches and outcomes, are publicly accessible from the TCGA website (https://cancergenome.nih.gov).

AUTHOR CONTRIBUTIONS
MN, NZ, and RW carried out study design and computational analyses. TS, YF, YS, XL, KZ, SZ, and LX participated in statistical analysis. YY and KX contributed to study design and drafted the manuscript. All authors read and approved the final manuscript.