MicroRNA–mRNA Pairs Associated with Outcome in AML: From In Vitro Cell-Based Studies to AML Patients

Cytarabine is the primary chemotherapeutic agent used for treatment of acute myeloid leukemia (AML). Disease relapse after initial remission remains one of the most pressing therapeutic challenges in the treatment of AML. Relapsed disease is often resistant to cytarabine and subsequent salvage therapy is ineffective. Recent studies have shown that some microRNAs (miRNAs) are associated with prognosis, but have not yet explored the role of miRNAs in cellular response to cytarabine. We identified 20 miRNAs that associate with the in vitro cytarabine chemo-sensitivity or apoptotic response of eight AML cell lines. Out of the 20 miRNAs, data on 18 miRNAs was available in AML patients from The Cancer Genome Atlas database. Our stepwise-integrated analyses (step 1 – miRNA–target mRNA that were significantly correlated in AML patients; step 2 – mRNAs from step 1 with significant association with overall survival (OS)) identified 23 unique miRNA–mRNA pairs predictive of OS in AML patients. As expected HOX genes (HOXA9, HOXB7, and HOXA10) were identified to be regulated by miRs as well as predictive of worse OS. Additionally, miR107-Myb, miR-378-granzyme B involved in granzyme signaling and miR10a-MAP4K4 were identified to be predictive of outcome through integrated analysis. Although additional functional validations to establish clinical/pharmacologic importance of miRNA–mRNA pairs are needed, our results from RNA electrophoretic mobility shift assay confirmed binding of miR-10a, miR-378, and miR-107 with their target genes GALNT1, GZMB, and MYB, respectively. Integration of pathogenic and pharmacologically significant miRNAs and miRNA–mRNA relationships identified in our study opens up opportunities for development of targeted/miRNA-directed therapies.

Introduction

Acute myeloid leukemia (AML) is a hematological malignancy characterized by the presence of immature abnormal myeloid cells in bone marrow. It is a heterogeneous disease with various subtypes classified based on the morphology, immunophenotype, and cytogenetics that are associated with outcome (Döhner et al., 2010). In spite of advances in recent years, 5-year overall survival (OS) is roughly 60% for children and ∼25% for adults (Cancer Facts and Figures, American Cancer Society). Cytarabine (1-β-arabinofuranosylcytosine, ara-C), a nucleoside analog, is the most widely used chemotherapeutic agent used in combination with an anthracycline for treatment of AML. Cytarabine requires intracellular activation to form an active triphosphate metabolite that triggers apoptosis by inhibiting DNA synthesis. Although chemotherapeutic regimens including cytarabine induce complete response in 65–80% of AML patients, the majority of these patients suffer from disease relapse within 2 years of diagnosis (Cros et al., 2004). This can be partly attributed to the development of resistance of leukemic cells to cytarabine-based chemotherapy regimens (Montillo et al., 1998; Estey, 2000). Several factors such as molecular and cytogenetic subtype, differential gene-expression profiles, and epigenetics can account for development of resistance in AML (Sasaki et al., 2008; Peters et al., 2010; Mitra et al., 2011; Cao et al., 2013; Gamazon et al., 2013; Mortland et al., 2013). We have previously identified gene expression signatures in AML patients predicting beneficial and detrimental patterns associated with cytarabine-based response (Lamba et al., 2011). Others have identified gene-expression differences between cytarabine-sensitive and -resistant cell lines in order to understand the molecular mechanisms underlying cytarabine resistance (Abe et al., 2006; Song et al., 2009).

MicroRNAs (miRNAs, miRs) are non-coding RNAs of 22–25 nucleotides that regulate gene expression. miRNAs bind to the complimentary sequence of messenger RNAs (mRNAs). In many cases, this binding suppresses mRNA translation or promotes mRNA degradation, thereby reducing expression at the protein level (Bartel, 2004). It has been shown that miRNAs play an important role in various cancers by regulating genes involved in cell proliferation, differentiation, and apoptosis (Zhang et al., 2010; Paik et al., 2011; Chen et al., 2013; Palma et al., 2014). Likewise, several recent studies have identified various miRNAs that differentiate the disease subgroups, associate with disease development, and associate with clinical prognosis of AML (Dixon-McIver et al., 2008; Jongen-Lavrencic et al., 2008; Marcucci et al., 2011; Starczynowski et al., 2011; Li et al., 2013; Jinlong et al., 2015; Shibayama et al., 2015). However, these valuable contributions have not yet carefully examined the roles of miRNAs in the cellular response of AML to cytarabine. Therefore, in this study, we characterized the association of miRNA expression with apoptotic response to cytarabine to identify candidate miRNAs for further evaluation of their relationship with mRNAs and OS in The Cancer Genome Atlas (TCGA) cohort of 200 AML patients.

Materials and Methods

Figure 1 shows overall study design.

FIGURE 1

FIGURE 1. Overall summary of the study. Study design for identification of microRNAs (miRNA) influencing cytarabine chemo-sensitivity and survival in acute myeloid leukemia (AML) patients.

Cell Culture and Regents

The AML cell lines HL-60, MV-4-11, Kasumi-1, THP-1, AML-193, and KG-1 were obtained from ATCC (Manassas, VA, USA), while ME-1 and MOLM-16 cell lines were obtained from DSMZ (Braunschweig, Germany). All the cell lines were cultured in the media as recommended by the supplier and were maintained in a humidified incubator at 37°C with 5% CO₂. The cells were passaged every 2–3 days in order to maintain them in logarithmic growth phase. Cytarabine and the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) reagent were purchased from Sigma–Aldrich (St. Louis, MO, USA). Stock concentrations for cytarabine (5 mg/ml) were prepared in sterile water and stored in –20°C in aliquots.

Cytotoxicity Assay

Cytarabine cytotoxicity was determined using the MTT assay. Briefly, AML cell lines were plated in 96-well plate at seeding density of 2.5 × 10⁵ cells/ml and incubated at 37°C overnight. After 24 h of recovery time, the cells were exposed to varying concentrations (200, 100, 50, 5, 0.5, 0.01, and 0 μM) of cytarabine. Cell viability was determined 48 h post cytarabine treatment by incubation with MTT reagent followed by measuring the absorbance at 570 nm using Synergy plate reader (BioTek, USA). The percent cell survival at each concentration was calculated using the Gen5^TM Software version 1.11 (Winooski, VT, USA). The area under the survival curve (AUC) was calculated by the trapezoidal method using the GraphPad Prism software version 6 (La Jolla, CA, USA).

Apoptosis Assay

The apoptotic activity of AML cell lines following treatment (48 h) with varying concentration of cytarabine (as indicated above) was determined using the Caspase-Glo^® 3/7 assay as per manufacturer’s instructions (Promega, USA). Forty-eight hours of post cytarabine treatment, luminescence was read using Synergy plate reader (BioTek, USA). The luminescence produced is directly proportional to the caspase activity. The caspase activity at each concentration was normalized to the control and the area under the relative caspase activity curve (AUC) was calculated by the trapezoidal method using the GraphPad Prism software version 6 (La Jolla, CA, USA).

MicroRNA Expression Analysis

For determination of miRNA expression, total RNA was isolated using mirVana^TM miRNA Isolation kit (Life Technologies, USA) as per the manufacturer’s protocol. The RNA quality and concentration was measured using NanoDrop 2000 UV-Vis spectrophotometer (Thermo Scientific, USA). A total of 100 ng of purified total RNA was assayed for determination of 800 human miRNA expression using the nCounter Human v2 miRNA Expression Assay kit (Nanostring Technologies, USA). miRNA expression data normalization was performed using the nSolver^TM Analysis Software (Nanostring Technologies) according to the manufacturer’s instructions. In order to avoid using the miRNAs with a very low expression, we further filtered out the miRNAs with expression counts <30 (two times the mean ± 2SD of negative control value, accounting for the background noise). Total 412 miRNAs with expression counts >30 were evaluated for differential expression between sensitive and resistant cell lines and for their correlation with cytarabine chemosensitivity.

TCGA Data

The clinical outcome data, mRNA expression and miRNA expression data for AML patients were extracted from The Cancer Genome Atlas (TCGA) Data Portal¹ (Cancer Genome Atlas Research Network, 2013). Out of the 200 AML patients in TCGA database, 197 patients had gene expression profiling data available and 187 patients had miRNA expression data available. One hundred and eighty-six patients had both gene expression and miRNA expression data available. Out of the 186 AML patients, 13 patients lacked valid survival information and two patients lacked cytogenetically defined risk information. Thus, data for a total of 186 patients were used to evaluate miRNA–mRNA associations, 173 patients used to evaluate mRNA-OS associations, and miRNA-OS associations (171 for stratified analyses).

Statistical Analysis

For each miRNA, Spearman’s rank-based correlation was used to measure the association with cytarabine treatment response or apoptosis in the eight cell lines. For each miRNA–mRNA pair with predicted binding sites defined by miRBase² (release 21), Spearman’s rank-based correlation was used to evaluate the association of miRNA expression with mRNA expression on the TCGA AML cohort. The p-value for the Spearman statistic was determined by 10,000 permutations. For each miRNA or mRNA, Cox regression [or Jung’s statistic] was used to evaluate the association of expression with OS. FDR was estimated by Pounds and Cheng’s (2006) robust FDR method.

Electrophoretic Mobility Shift Assays

The functional validation for binding efficiencies between selected miRNAs and mRNAs was performed using the electrophoretic mobility shift assays (EMSAs) as described previously (Yu et al., 2015). The binding free energy between the respective mRNA and miRNA pair (demonstrating inverse relationship) was predicted using the RNAhybrid software. The miRNA oligonucleotides were labeled with cy5^TM dye on their 5′ ends. The 2′ O-methyl-modified mRNA oligonucleotides were labeled with IRDye^® 800 (LI-COR Biosciences, USA) dye on their 5′ ends. The labeled oligonucleotides were synthesized by Integrated DNA Technologies (Coralville, IA, USA). RNA EMSA experiment was performed using the LightShift Chemiluminescent RNA EMSA Kit (Thermo Scientific, USA) according to the manufacturer’s protocol. The mRNA oligonucleotide was heated for 10 min at 80°C and then placed on ice in order to relax the secondary structures. In each 20 μl binding reaction, 200 nM miRNA oligonucleotide and/or mRNA oligonucleotide were mixed with RNA EMSA binding buffer and incubated at 25°C for 25 min. The reaction mixtures were separated on a 15% polyacrylamide gel by electrophoresis at 4°C. The binding reactions were electrophoretically transferred onto nylon membrane and the resulting mobility shifts were imaged using Odyssey CLx Infrared System (LI-COR Biosciences, USA).

Results

Cytarabine Chemo-Sensitivity of AML Cell Lines

The AML cell lines showed considerable variability in cytarabine sensitivity as measured by cytotoxicity AUC and apoptosis AUC in the MTT assay (Table 1; Supplementary Figures S1A–H). Based on the cytotoxicity AUC, HL-60, MV-4-11, KG-1, and ME-1 were classified as sensitive (cytarabine AUC < 12000), while MOLM-16, AML-193, Kasumi-1, and THP-1 were classified as resistant cell lines (cytarabine AUC > 12000).

TABLE 1

TABLE 1. Characterization of acute myeloid leukemia (AML) cell lines for cytarabine chemosensitivity.

MicroRNAs Associated with Cytarabine Chemo-Sensitivity in AML Cell Lines

Of the 800 human miRNAs quantitated using the nCounter Nanostring platform we excluded 388 miRNAs due to very low expression and 412 miRNA were analyzed further. Twenty miRNAs were associated with cellular viability (nine miRs) or caspase activation (11 miRs) in AML cell lines post treatment with cytarabine (Table 2). Expression of miR-25-3p, miR-148b-3p, miR-107, miR-374-5p, miR-425-5p were positively associated with AUC for cell survival post-cytarabine treatment and miR-16-5p, miR-24-3p, miR-196a-5p, and miR-155-5p were negatively associated with AUC for cell survival post-cytarabine treatment (p < 0.05). (Selected miRNAs are shown in Supplementary Figure S2.) Expression levels of miR-10a-5p, miR-29a/b-3p, miR-30e-5p, miR-33a-5p, miR-378a/g were positively and expression levels miR-197-3p, miR-27b-3p, miR-324-5p, and miR-421 were negatively associated with AUC for caspase-3/7 activation (apoptosis) post cytarabine treatment (Table 2, p < 0.05). Using Ingenuity pathway analysis tool, the miRNAs that were correlated with cytarabine chemosensitivity were also found to potentially impact important biological process relevant to leukemia/cancer (Supplementary Figure S3).

TABLE 2

TABLE 2. MicroRNAs significantly associated with cytarabine-induced cytotoxicity AUC and cytarabine-induced apoptosis (caspase-3/7 activity).

Pairs of Significantly Correlated mRNAs and miRNAs that Associate with Overall Survival of AML Patients

Of 20 miRNAs identified above, data on 18 were available in AML patients from TCGA database and were tested for associations with risk group and outcome. As shown in Supplementary Table S1, seven of these miRNAs (miR-10a, miR-16, miR-196a, miR-197, miR-421, miR-155, and miR24) demonstrated significant difference in expression levels among the three risk groups. In risk stratified analysis miR107, miR-155, miR-196a, miR-25, and miR29b were associated with worse OS, whereas miR-25 was predictive of better OS in AML patients at p < 0.05.

Using miRBase² (release 21), we determined that 5006 probe sets representing 2830 gene with binding sites for these 18 miRNAs. These 5006 mRNAs and 18 miRNAs belong to 7132 distinct miRNA–mRNA pairs. Using the analysis strategy outlined in Figure 1, we found that 23 of the 7132 miRNA–mRNA pair’s satisfied criteria listed below (Table 3):

(a) significant association between miRNA and target mRNA (p < 0.05; 1532 pairs) and

(b) mRNA expression associated significantly with OS in an unstratified Cox regression model (p < 0.001; FDR < 0.05).

TABLE 3

TABLE 3. miRNA–mRNA pairs predictive of overall survival in AML Patients (data from TCGA).

These 23 pairs included 16 unique genes and 10 unique miRNAs (Table 3 – some mRNAs and some miRNAs belonged to multiple pairs). A positive correlation of mRNA and miRNA was observed for 10 of these pairs and a negative correlation was observed for the other 13 pairs.

Among mRNA-OS associations COL3A1, GALNT1, GALNT7, LTK, MAP4K4, MYB, PAPOLG, RPL35A, TMEM87A, and WDR48 were associated with better OS and HOX family genes (HOXA9, HOXA10, HOXB7), GZMB, SE1L3, and an oncogene PIM1 were associated with inferior outcome (Table 3). Since nine of these genes demonstrated significant association with risk groups we also performed risk-stratified analysis and all but three genes HOXA9, HOXA10, and LTK were significantly associated with OS in risk stratified analysis, indicating that for these genes the observed association with OS might be driven by risk group characteristics. Figure 2 shows the representative correlation plots as well as overall curves for miR10a-GALNT1, miR10a-MAP4K4, miR16-Pim1, miR378-GZMB, and miR107-MYB.

FIGURE 2

FIGURE 2. Representative plots showing correlation between miRNA–target messenger RNA (mRNA) and survival curves of miRNA and mRNA expression with overall survival (OS) in AML patients from The Cancer Genome Atlas (TCGA) database. (A) Correlation plot of miR10a–GALNT1 mRNA levels and Kaplan–Meier survival curves of miR10a and GALNT1 expression with OS. (B) Plots for miR-10a–MAP4K4 pair. (C) Plots for miR-16 and Pim1 pair. (D) Plots for miR-378 and GZMB pair. (E) Plots for miR-107 and Myb pair.

We further utilized Ingenuity pathway analysis tool to map the 16 unique genes identified in integrative analysis and as shown in Figure 3, these genes were associated with cell proliferation, apoptosis, RNA expression, and quantity of hematopoietic progenitor cells, hematological cancer and myeloid leukemia.

FIGURE 3

FIGURE 3. Network analysis. Ingenuity pathway analysis tool was utilized to evaluate the genes identified in miRNA–mRNA-OS analysis. Several genes mapped to biological processes of relevance to hematological malignancies, myeloid leukemia as well as to cell proliferation and apoptosis, RNA expression. Genes associated with better OS are in green and ones with worse OS are in red.

Functional Validation Using RNA Electrophoretic Mobility Shift Assays (RNA EMSA)

To further validate the miRNA–mRNA pairs discovered as above, we performed RNA EMSAs. We first selected miRNA–mRNA pairs from Table 3 that demonstrated significant inverse relationships. Binding free energy for these miRNA–mRNA pairs was calculated using RNAhybrid software. The 3′UTR sequence of mRNA was obtained from the UCSC Genome browser and miRNA sequence was obtained from miRBase software. Of miRNA–mRNAs demonstrating minimum free energy of binding <–24 kcal/mol miR107-MYB; miR378a-GZMB, and miR10a-GALNT1 were validated by RNA EMSAs. EMSA results confirmed binding for miR-107-MYB, miR-10-GALNT1, and miR-378-GZMB and are shown in Figure 4, respectively. As shown in lane 3 in Figure 4, the EMSA results for miR-107, miR-10a, and miR-378 demonstrate binding of this miRNA with their respective mRNA target sequences confirming the thermodynamic stability of these complexes predicted in the in silico analysis. In addition, we observed competition of binding in mRNA–miRNA complexes after adding excess unlabeled specific miRNA probe (Figure 4, lane 4), excess unlabeled mRNA probe (Figure 4, lane 5) but not by adding excess unlabeled non-specific probe (Figure 4, lane 6).

FIGURE 4

FIGURE 4. Validation of binding interaction between miRNA–mRNA by RNA electrophoretic mobility shift assays (RNA EMSAs): (A) miR107–MYB; (B) miR10a–GALNT1 and (C) miR378–GZMB. RNA EMSAs were performed using cy5-labeled respective miRNA oligonucleotide and 2′-O-methyl modified and IRD-800 labeled mRNA oligonucleotides. Lanes 1 and 2 show the mobility of the respective labeled mRNA or miRNA oligonucleotide. Lane 3 shows the mobility of the labeled miRNA oligonucleotide with its target mRNA oligonucleotide. Lanes 4 and 6 show the mobility of labeled mRNA oligonucleotide in presence of excess unlabeled specific competitor miRNA oligonucleotide and excess unlabeled non-specific competitor (NC). Lane 5 shows the mobility of labeled mRNA oligonucleotide in presence of excess unlabeled specific competitor mRNA oligonucleotide.

Discussion

Acute myeloid leukemia is a heterogeneous disease with dismal outcome. Additional complexity is added by very heterogeneous nature of AML with cytogenetic abnormalities used for risk classification in AML patients. Although cytarabine has been the backbone of AML chemotherapy for more than 50 years, there are still gaps in our understanding of the molecular mechanisms contributing to development of drug resistance in AML. Although advances in supportive care have improved, the treatment strategies have not changed much with cytarabine being still the main player. Thus, understanding the molecular mechanisms underlying cytarabine resistance will be of great interest in developing predictive models of outcome as well as for developing novel therapeutic strategies. Recent research has shown the significant role of miRNAs in normal hematopoiesis (Chen et al., 2004; Garzon and Croce, 2008; O’Connell et al., 2008) as well as miRNA deregulation in AML. miRNAs (miR-125, miR-146, miR-142, miR-155, miR-29, miR-181, let-7a, etc.) of potential prognostic significance have been identified (Khalaj et al., 2014). Recent data supports putative diagnostic role of miR-155 in all AML subgroups and miR-196b within M4-5 subgroups (Yan et al., 2015). miR-9 has been also implicated in promoting proliferation of leukemic cells in normal karyotype AML by targeting Hes1 (Tian et al., 2015). However, the role of miRNAs in development of resistance to cytarabine and thus inferior clinical outcome has not been investigated in detail. In this study, we report results of our genome-wide evaluation to identify miRNAs associated with cytarabine sensitivity in eight AML cell lines as well as their impact on clinical outcome in AML patients (Figure 1 outlines overall study design and results).

We screened expression levels of 800 human miRNAs in eight AML cell lines and after filtering out miRNAs with very low expression, evaluated 412 miRNAs for association with cytarabine chemo-sensitivity, measured as cell viability and apoptosis induction following cytarabine treatment. Twenty unique miRNAs were predictive of cytarabine chemo-sensitivity in AML cell lines, and 18 of these were further evaluated in AML patients from TCGA database. Seven of these were differentially expressed among AML risk groups (favorable, intermediate and poor; p < 0.05) and after risk stratification five miRNAs (miR-107, miR-155, miR-25, miR-29b, and miR-196a) were associated with OS (Supplementary Table S1). All but miR-25 were associated with worse OS in AML patients.

Among the miRNAs associated with OS, miR-155 is located in a non-coding RNA transcript cluster called B-cell integration cluster (BIC), which has been shown to cooperate with c-Myc (Clurman and Hayward, 1989; Tam, 2001). miR-155 is considered as an oncomiR with implications in pathogenesis of AML (O’Connell et al., 2008; Palma et al., 2014); it has been associated with SHIP1 (negatively regulator of PI3K/AKT pathway) and CEBP-β (Gorgoni et al., 2002; O’Connell et al., 2009). Our results show that miR-155 is associated with cytarabine sensitivity which is not in consensus with its association with inferior outcome, thereby indicating that miR-155 (which is also differentially expressed among risk groups) might have significant impact on disease pathogenesis but might not be impacting drug response. miR-29 family members are regulators of myeloid differentiation and have been shown to be deregulated in AML (Garzon et al., 2008, 2009a; Han et al., 2010). MiR29a/29b have been associated with expression levels of oncogenes MCL1, CDK6, IGFR, and JAK2 (Mott et al., 2007; Garzon et al., 2009a) as well as have been shown to target DNA modifying genes DNMTs and TET2 (Garzon et al., 2009b; Cheng et al., 2013).

In step-wise integrated analysis we identified 23 miRNA–mRNA pairs predictive of survival in AML patients and these pairs were mapped to 16 unique mRNAs.

We further validated miR107-MYB, miR-10a-GALNT1, and miR378-GZMB miRNA–mRNA pairs using electrophoretic mobility shift assays (Figure 4), which confirmed the binding of these oligos as supported by in silico analysis indicating miRNAs in regulating gene expression of these target genes by binding to specific seed sequences.

As expected HOX genes (HOXA9, HOXA10, and HOXB7) were predictive of worse outcome. miR-196a-1 which was associated cytarabine in vitro sensitivity in AML cell lines was positively correlated in expression with HOXA9 and HOXB7 as well as with AML risk groups. miR-196a-1 gene co-localizes with HOXB gene cluster between HOXB8 and HOXB13, positive correlation observed in AML patients between HOXB7 and mir-196a might be due to transcriptional co-regulation. miR-196a has been previously shown to be positively correlated with several HOX family members including HOXB7 and HOXA9 (Debernardi et al., 2007). In addition to miR-196a, we observed significant negative correlation between miR-16-HOXA10, miR-421-HOXA9, and miR-30e-HOXB7.

Among other miR target genes that were associated with worse OS were PIM-1 (pim1 oncogene), GZMB (granzyme B), and SEL1L3 (Sel-1 Suppressor Of Lin-12-Like 3). Pim1 is a serine/threonine protein kinase that has role in cell survival and cell proliferation. HOXA9 is transcriptional activator of Pim-1, which is further involved in regulation of MYC transcriptional activity, regulation of cell cycle progression, and phosphorylation and inhibition of proapoptotic proteins (BAD, MAP3K5, and FOXO3) thereby by contributing to its oncogenic activity. Pim-1 is also involved in inactivating MAP3K5 by phosphorylation thereby inhibiting MAP3K5-mediated phosphorylation of JNK and JNK/p38MAPK subsequently reducing caspase-3 activation and cell apoptosis. Pim-1 seems like a potential target for drug discovery, in fact in pediatric preclinical models, Pim1 inhibitor SGI-1776 has been shown to induce complete response to subcutaneous MV4:11 leukemia (Batra et al., 2012) as well as inhibit proliferation in other malignancies such as CLL, B cell lymphoma, multiple myeloma, etc. (Cervantes-Gomez et al., 2013, 2015; Yang et al., 2013).

GZMB was negatively regulated by miR-378 and higher expression was associated with worse OS in AML. GZMB is a key player in Granzyme signaling pathway, which is a lymphocyte granular serine protease that cleaves its substrates at Asp residues. GZMB is expressed in cytotoxic T lymphocytes (CTL) and NK cells and is primary mediator of apoptosis by CTL in cell-mediated immune response. GZMB seems to play critical role antibody –dependent cellular cytotoxicity (Elavazhagan et al., 2015).

Among miR target genes that were associated with good response were family members of Polypeptide N-Acetylgalactosaminyltransferases (GALNT1 and GALNT7), MAP4K4, TMEM87A, and COL3A1. GALNT1 was correlated negatively in expression with miR-10a, miR-16-2 and positively with miR-30e, which also demonstrated positive correlation with GALNT7 expression. MAP4K4 belongs to serine/threonine protein kinase family and has been specifically implicated in activation of MAPK8/JNK pathway.

Both MYB and TMEM87 were inversely associated with miR-107, which was associated with cytarabine resistance in cell lines as well as worse OS in AML patients (Tables 2 and 3). Oncogenic role of miR-107 in regulating tumor invasion and metastasis in gastric cancer by targeting DICER1 (Li et al., 2011) and in colorectal cancer by targeting metastasis suppressors death-associated protein kinase (DAPK) and Krüppel-like factor 4 (KLF4; Chen et al., 2012) has been proposed. In AML patients, we observed negative correlation of DICER with miR-107 (p < 0.01) although DICER expression was not predictive of outcome. Although TMEM87A a transmembrane protein has not been well studied, MYB, a V-Myb Avian Myeloblastosis Viral Oncogene Homolog has been implicated in leukemogenesis. Myb is reported to be overexpressed in AML and results from recent studies shows its potential role in interplay between C/EBPα activity for transcriptional regulation of FLT3 expression (Volpe et al., 2013). Recent report in luminal breast cancer demonstrated for the first time the potential tumor suppressor role of c-Myb gene, (Thorner et al., 2010), which is in concordance with TCGA results with Myb expression associated with better OS thereby warranting further investigation of Myb gene on its impact on treatment outcome in AML.

In summary, although several studies have established prognostic significance of miRNAs (miR-155, miR-29, miR-16, etc.) in AML, role of miRNAs in cytarabine chemosensitivity and development of resistance as a contributor of inferior outcome has not been well studied. In this report, we performed genome-wide miRNA profiling of eight AML cell lines and identified 20 miRNAs predictive of differential AML in vitro chemo-sensitivity (by measuring both cytarabine induced cell death and apoptosis). These were further investigated in AML patients using data from TCGA database, and in an integrated three way analyses of miRNA–target mRNA pairs with significant association and OS we identified 23 miRNA–mRNA-OS pairs of therapeutic importance in AML patients. Although additional functional validation studies to establish clinical/pharmacological importance of miRNA–mRNA pairs are needed, our preliminary data on RNA EMSAs confirmed binding of miR-107-MYB, miR10a-GALNT1, and miR-378-GZMB. Integration of pathogenic and pharmacologically significant miRNAs and miRNA–mRNA relationships opens up opportunities for development of targeted/miRNA-directed therapies.

Author Contributions

JL, VL, and NB were involved in designing the study; NB, LC, and MS performed the experiments; NB, JL, XC, and SP performed statistical analysis. All authors contributed to manuscript writing. All authors read and approved the final manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgment

This work is supported by funds from NIH, National Cancer Institute funding R01-CA132946 (Lamba and Pounds).

Supplementary Material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fphar.2015.00324

TABLE S1 | Association of MicroRNA with OS in AML Patients.

FIGURE S1 | Characterization of acute myeloid leukemia (AML) Cell lines for Cytarabine chemosensitivity. AML cell lines (n = 8) were cultured in the respective media and standard culturing conditions. (A–H) Cytarabine-induced cytotoxicity was determined by MTT assay after 48 h of drug treatment and area under the survival curve (AUC) was calculated using the cell survival data. Cytarabine-induced apoptosis was determined using caspase-3/7 assay after 48 h of drug treatment and relative caspase AUC was calculated.

FIGURE S2 | MicroRNA (miRNA) and cytarabine cell cytotoxicity. Correlation plots showing association of selective miRNAs with cell survival AUC post cytarabine treatment of eight AML cell lines.

FIGURE S3 | Network of miRNAs and important biological processes. miRNAs predictive of in vitro cytarabine chemosensitivity in cell lines as well as with overall survival (OS) in AML patients were analyzed using Ingenuity pathway analysis tools and pathways identified are depicted in the network.

Footnotes

1. ^ cancergenome.nih.gov/
2. ^ www.mirbase.org

References

Abe, S., Funato, T., Takahashi, S., Yokoyama, H., Yamamoto, J., Tomiya, Y., et al. (2006). Increased expression of insulin-like growth factor i is associated with Ara-C resistance in leukemia. Tohoku J. Exp. Med. 209, 217–228. doi: 10.1620/tjem.209.217

PubMed Abstract | CrossRef Full Text | Google Scholar

Bartel, D. P. (2004). MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 116, 281–297. doi: 10.1016/S0092-8674(04)00045-5

CrossRef Full Text | Google Scholar

Batra, V., Maris, J. M., Kang, M. H., Reynolds, C. P., Houghton, P. J., Alexander, D., et al. (2012). Initial testing (stage 1) of SGI-1776, a PIM1 kinase inhibitor, by the pediatric preclinical testing program. Pediatr. Blood Cancer 59, 749–752. doi: 10.1002/pbc.23364

PubMed Abstract | CrossRef Full Text | Google Scholar

Cancer Genome Atlas Research Network (2013). Genomic and epigenomic landscapes of adult de novo acute myeloid leukemia. N. Engl. J. Med. 368, 2059–2074. doi: 10.1056/NEJMoa1301689

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, X., Mitra, A. K., Pounds, S., Crews, K. R., Gandhi, V., Plunkett, W., et al. (2013). RRM1 and RRM2 pharmacogenetics: association with phenotypes in HapMap cell lines and acute myeloid leukemia patients. Pharmacogenomics 14, 1449–1466. doi: 10.2217/pgs.13.131

PubMed Abstract | CrossRef Full Text | Google Scholar

Cervantes-Gomez, F., Chen, L. S., Orlowski, R. Z., and Gandhi, V. (2013). Biological effects of the Pim kinase inhibitor, SGI-1776, in multiple myeloma. Clin. Lymphoma Myeloma Leuk. 13(Suppl. 2), S317–S329. doi: 10.1016/j.clml.2013.05.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Cervantes-Gomez, F., Lavergne, B., Keating, M. J., Wierda, W. G., and Gandhi, V. (2015). Combination of Pim kinase inhibitors and Bcl-2 antagonists in chronic lymphocytic leukemia cells. Leuk. Lymphoma 1–9. doi: 10.3109/10428194.2015.1063141 [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, C. Z., Li, L., Lodish, H. F., and Bartel, D. P. (2004). MicroRNAs modulate hematopoietic lineage differentiation. Science 303, 83–86. doi: 10.1126/science.1091903

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, H. Y., Lin, Y. M., Chung, H. C., Lang, Y. D., Lin, C. J., Huang, J., et al. (2012). miR-103/107 promote metastasis of colorectal cancer by targeting the metastasis suppressors DAPK and KLF4. Cancer Res. 72, 3631–3641. doi: 10.1158/0008-5472.CAN-12-0667

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Jacamo, R., Konopleva, M., Garzon, R., Croce, C., and Andreeff, M. (2013). CXCR4 downregulation of let-7a drives chemoresistance in acute myeloid leukemia. J. Clin. Invest. 123, 2395–2407. doi: 10.1172/JCI66553

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, J., Guo, S., Chen, S., Mastriano, S. J., Liu, C., D’alessio, A. C., et al. (2013). An extensive network of TET2-targeting MicroRNAs regulates malignant hematopoiesis. Cell Rep. 5, 471–481. doi: 10.1016/j.celrep.2013.08.050

PubMed Abstract | CrossRef Full Text | Google Scholar

Clurman, B. E., and Hayward, W. S. (1989). Multiple proto-oncogene activations in avian leukosis virus-induced lymphomas: evidence for stage-specific events. Mol. Cell. Biol. 9, 2657–2664. doi: 10.1128/MCB.9.6.2657

PubMed Abstract | CrossRef Full Text | Google Scholar

Cros, E., Jordheim, L., Dumontet, C., and Galmarini, C. M. (2004). Problems related to resistance to cytarabine in acute myeloid leukemia. Leuk. Lymphoma 45, 1123–1132. doi: 10.1080/1042819032000159861

CrossRef Full Text | Google Scholar

Debernardi, S., Skoulakis, S., Molloy, G., Chaplin, T., Dixon-McIver, A., and Young, B. D. (2007). MicroRNA miR-181a correlates with morphological sub-class of acute myeloid leukaemia and the expression of its target genes in global genome-wide analysis. Leukemia 21, 912–916. doi: 10.1038/sj.leu.2404605

PubMed Abstract | CrossRef Full Text | Google Scholar

Dixon-McIver, A., East, P., Mein, C. A., Cazier, J. B., Molloy, G., Chaplin, T., et al. (2008). Distinctive patterns of microRNA expression associated with karyotype in acute myeloid leukaemia. PLoS ONE 3:e2141. doi: 10.1371/journal.pone.0002141

PubMed Abstract | CrossRef Full Text | Google Scholar

Döhner, H., Estey, E. H., Amadori, S., Appelbaum, F. R., Büchner, T., Burnett, A. K., et al. (2010). Diagnosis and management of acute myeloid leukemia in adults: recommendations from an international expert panel, on behalf of the European LeukemiaNet. Blood 115, 453–474. doi: 10.1182/blood-2009-07-235358

PubMed Abstract | CrossRef Full Text | Google Scholar

Elavazhagan, S., Fatehchand, K., Santhanam, V., Fang, H., Ren, L., Gautam, S., et al. (2015). Granzyme B expression is enhanced in human monocytes by TLR8 agonists and contributes to antibody-dependent cellular cytotoxicity. J. Immunol. 194, 2786–2795. doi: 10.4049/jimmunol.1402316

PubMed Abstract | CrossRef Full Text | Google Scholar

Estey, E. H. (2000). Treatment of relapsed and refractory acute myelogenous leukemia. Leukemia 14, 476–479. doi: 10.1038/sj.leu.2401568

CrossRef Full Text | Google Scholar

Gamazon, E. R., Lamba, J. K., Pounds, S., Stark, A. L., Wheeler, H. E., Cao, X., et al. (2013). Comprehensive genetic analysis of cytarabine sensitivity in a cell-based model identifies polymorphisms associated with outcome in AML patients. Blood 121, 4366–4376. doi: 10.1182/blood-2012-10-464149

PubMed Abstract | CrossRef Full Text | Google Scholar

Garzon, R., and Croce, C. M. (2008). MicroRNAs in normal and malignant hematopoiesis. Curr. Opin. Hematol. 15, 352–358. doi: 10.1097/MOH.0b013e328303e15d

PubMed Abstract | CrossRef Full Text | Google Scholar

Garzon, R., Heaphy, C. E., Havelange, V., Fabbri, M., Volinia, S., Tsao, T., et al. (2009a). MicroRNA 29b functions in acute myeloid leukemia. Blood 114, 5331–5341. doi: 10.1182/blood-2009-03-211938

PubMed Abstract | CrossRef Full Text | Google Scholar

Garzon, R., Liu, S., Fabbri, M., Liu, Z., Heaphy, C. E., Callegari, E., et al. (2009b). MicroRNA-29b induces global DNA hypomethylation and tumor suppressor gene reexpression in acute myeloid leukemia by targeting directly DNMT3A and 3B and indirectly DNMT1. Blood 113, 6411–6418. doi: 10.1182/blood-2008-07-170589

PubMed Abstract | CrossRef Full Text | Google Scholar

Garzon, R., Volinia, S., Liu, C. G., Fernandez-Cymering, C., Palumbo, T., Pichiorri, F., et al. (2008). MicroRNA signatures associated with cytogenetics and prognosis in acute myeloid leukemia. Blood 111, 3183–3189. doi: 10.1182/blood-2007-07-098749

PubMed Abstract | CrossRef Full Text | Google Scholar

Gorgoni, B., Maritano, D., Marthyn, P., Righi, M., and Poli, V. (2002). C/EBP beta gene inactivation causes both impaired and enhanced gene expression and inverse regulation of IL-12 p40 and p35 mRNAs in macrophages. J. Immunol. 168, 4055–4062. doi: 10.4049/jimmunol.168.8.4055

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, Y. C., Park, C. Y., Bhagat, G., Zhang, J., Wang, Y., Fan, J. B., et al. (2010). microRNA-29a induces aberrant self-renewal capacity in hematopoietic progenitors, biased myeloid development, and acute myeloid leukemia. J. Exp. Med. 207, 475–489. doi: 10.1084/jem.20090831

PubMed Abstract | CrossRef Full Text | Google Scholar

Jinlong, S., Lin, F., Yonghui, L., Li, Y., and Weidong, W. (2015). Identification of let-7a-2-3p or/and miR-188-5p as prognostic biomarkers in cytogenetically normal acute myeloid leukemia. PLoS ONE 10:e0118099. doi: 10.1371/journal.pone.0118099

PubMed Abstract | CrossRef Full Text | Google Scholar

Jongen-Lavrencic, M., Sun, S. M., Dijkstra, M. K., Valk, P. J., and Lowenberg, B. (2008). MicroRNA expression profiling in relation to the genetic heterogeneity of acute myeloid leukemia. Blood 111, 5078–5085. doi: 10.1182/blood-2008-01-133355

PubMed Abstract | CrossRef Full Text | Google Scholar

Khalaj, M., Tavakkoli, M., Stranahan, A. W., and Park, C. Y. (2014). Pathogenic microRNA’s in myeloid malignancies. Front. Genet. 5:361. doi: 10.3389/fgene.2014.00361

CrossRef Full Text | Google Scholar

Lamba, J. K., Crews, K. R., Pounds, S. B., Cao, X., Gandhi, V., Plunkett, W., et al. (2011). Identification of predictive markers of cytarabine response in AML by integrative analysis of gene-expression profiles with multiple phenotypes. Pharmacogenomics 12, 327–339. doi: 10.2217/pgs.10.191

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., Zhang, Y., Shi, Y., Dong, G., Liang, J., Han, Y., et al. (2011). MicroRNA-107, an oncogene microRNA that regulates tumour invasion and metastasis by targeting DICER1 in gastric cancer. J. Cell Mol. Med. 15, 1887–1895. doi: 10.1111/j.1582-4934.2010.01194.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Gao, L., Luo, X., Wang, L., Gao, X., Wang, W., et al. (2013). Epigenetic silencing of microRNA-193a contributes to leukemogenesis in t(8;21) acute myeloid leukemia by activating the PTEN/PI3K signal pathway. Blood 121, 499–509. doi: 10.1182/blood-2012-07-444729

PubMed Abstract | CrossRef Full Text | Google Scholar

Marcucci, G., Haferlach, T., and Dohner, H. (2011). Molecular genetics of adult acute myeloid leukemia: prognostic and therapeutic implications. J. Clin. Oncol. 29, 475–486. doi: 10.1200/JCO.2010.30.2554

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitra, A. K., Crews, K. R., Pounds, S., Cao, X., Feldberg, T., Ghodke, Y., et al. (2011). Genetic variants in cytosolic 5′-nucleotidase II are associated with its expression and cytarabine sensitivity in HapMap cell lines and in patients with acute myeloid leukemia. J. Pharmacol. Exp. Ther. 339, 9–23. doi: 10.1124/jpet.111.182873

PubMed Abstract | CrossRef Full Text | Google Scholar

Montillo, M., Mirto, S., Petti, M. C., Latagliata, R., Magrin, S., Pinto, A., et al. (1998). Fludarabine, cytarabine, and G-CSF (FLAG) for the treatment of poor risk acute myeloid leukemia. Am. J. Hematol. 58, 105–109. doi: 10.1002/(SICI)1096-8652(199806)58:2<105::AID-AJH3>3.0.CO;2-W

CrossRef Full Text | Google Scholar

Mortland, L., Alonzo, T. A., Walter, R. B., Gerbing, R. B., Mitra, A. K., Pollard, J. A., et al. (2013). Clinical significance of CD33 nonsynonymous single-nucleotide polymorphisms in pediatric patients with acute myeloid leukemia treated with gemtuzumab-ozogamicin-containing chemotherapy. Clin. Cancer Res. 19, 1620–1627. doi: 10.1158/1078-0432.CCR-12-3115

PubMed Abstract | CrossRef Full Text | Google Scholar

Mott, J. L., Kobayashi, S., Bronk, S. F., and Gores, G. J. (2007). mir-29 regulates Mcl-1 protein expression and apoptosis. Oncogene 26, 6133–6140. doi: 10.1038/sj.onc.1210436

PubMed Abstract | CrossRef Full Text | Google Scholar

O’Connell, R. M., Chaudhuri, A. A., Rao, D. S., and Baltimore, D. (2009). Inositol phosphatase SHIP1 is a primary target of miR-155. Proc. Natl. Acad. Sci. U.S.A. 106, 7113–7118. doi: 10.1073/pnas.0902636106

PubMed Abstract | CrossRef Full Text | Google Scholar

O’Connell, R. M., Rao, D. S., Chaudhuri, A. A., Boldin, M. P., Taganov, K. D., Nicoll, J., et al. (2008). Sustained expression of microRNA-155 in hematopoietic stem cells causes a myeloproliferative disorder. J. Exp. Med. 205, 585–594. doi: 10.1084/jem.20072108

PubMed Abstract | CrossRef Full Text | Google Scholar

Paik, J. H., Jang, J. Y., Jeon, Y. K., Kim, W. Y., Kim, T. M., Heo, D. S., et al. (2011). MicroRNA-146a downregulates NFκB activity via targeting TRAF6 and functions as a tumor suppressor having strong prognostic implications in NK/T cell lymphoma. Clin. Cancer Res. 17, 4761–4771. doi: 10.1158/1078-0432.CCR-11-0494

PubMed Abstract | CrossRef Full Text | Google Scholar

Palma, C. A., Al Sheikha, D., Lim, T. K., Bryant, A., Vu, T. T., Jayaswal, V., et al. (2014). MicroRNA-155 as an inducer of apoptosis and cell differentiation in Acute Myeloid Leukaemia. Mol. Cancer 13:79. doi: 10.1186/1476-4598-13-79

PubMed Abstract | CrossRef Full Text | Google Scholar

Peters, G. J., Hodzic, J., Ortega, B., Giovannetti, E., Adema, A. D., Broekhuizen, R., et al. (2010). Methylation specific PCR to characterize methylation of the promoter of deoxycytidine kinase. Nucleosides Nucleotides Nucleic Acids 29, 408–413. doi: 10.1080/15257771003730078

PubMed Abstract | CrossRef Full Text | Google Scholar

Pounds, S., and Cheng, C. (2006). Robust estimation of the false discovery rate. Bioinformatics 22, 1979–1987. doi: 10.1093/bioinformatics/btl328

PubMed Abstract | CrossRef Full Text | Google Scholar

Sasaki, S., Kobunai, T., Kitayama, J., and Nagawa, H. (2008). DNA methylation and sensitivity to antimetabolites in cancer cell lines. Oncol. Rep. 19, 407–412. doi: 10.3892/or.19.2.407

CrossRef Full Text | Google Scholar

Shibayama, Y., Kondo, T., Ohya, H., Fujisawa, S., Teshima, T., and Iseki, K. (2015). Upregulation of microRNA-126-5p is associated with drug resistance to cytarabine and poor prognosis in AML patients. Oncol. Rep. 33, 2176–2182. doi: 10.3892/or.2015.3839

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, J. H., Kim, S. H., Kweon, S. H., Lee, T. H., Kim, H. J., and Kim, T. S. (2009). Defective expression of deoxycytidine kinase in cytarabine-resistant acute myeloid leukemia cells. Int. J. Oncol. 34, 1165–1171.

PubMed Abstract | Google Scholar

Starczynowski, D. T., Morin, R., Mcpherson, A., Lam, J., Chari, R., Wegrzyn, J., et al. (2011). Genome-wide identification of human microRNAs located in leukemia-associated genomic alterations. Blood 117, 595–607. doi: 10.1182/blood-2010-03-277012

PubMed Abstract | CrossRef Full Text | Google Scholar

Tam, W. (2001). Identification and characterization of human BIC, a gene on chromosome 21 that encodes a noncoding RNA. Gene 274, 157–167. doi: 10.1016/S0378-1119(01)00612-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Thorner, A. R., Parker, J. S., Hoadley, K. A., and Perou, C. M. (2010). Potential tumor suppressor role for the c-Myb oncogene in luminal breast cancer. PLoS ONE 5:e13073. doi: 10.1371/journal.pone.0013073

PubMed Abstract | CrossRef Full Text | Google Scholar

Tian, C., You, M. J., Yu, Y., Zhu, L., Zheng, G., and Zhang, Y. (2015). MicroRNA-9 promotes proliferation of leukemia cells in adult CD34-positive acute myeloid leukemia with normal karyotype by downregulation of Hes1. Tumour Biol. doi: 10.1007/s13277-015-4581-x [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Volpe, G., Walton, D. S., Del Pozzo, W., Garcia, P., Dasse, E., O’neill, L. P., et al. (2013). C/EBPalpha and MYB regulate FLT3 expression in AML. Leukemia 27, 1487–1496. doi: 10.1038/leu.2013.23

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, W., Xu, L., Sun, Z., Lin, Y., Zhang, W., Chen, J., et al. (2015). MicroRNA biomarker identification for pediatric acute myeloid leukemia based on a novel bioinformatics model. Oncotarget 6, 26424–26436. doi: 10.18632/oncotarget.4459

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Q., Chen, L. S., Neelapu, S. S., and Gandhi, V. (2013). Combination of Pim kinase inhibitor SGI-1776 and bendamustine in B-cell lymphoma. Clin. Lymphoma Myeloma Leuk. 13(Suppl. 2), S355–S362. doi: 10.1016/j.clml.2013.05.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, D., Green, B., Marrone, A., Guo, Y., Kadlubar, S., Lin, D., et al. (2015). Suppression of CYP2C9 by microRNA hsa-miR-128-3p in human liver cells and association with hepatocellular carcinoma. Sci. Rep. 5:8534. doi: 10.1038/srep08534

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X. J., Ye, H., Zeng, C. W., He, B., Zhang, H., and Chen, Y. Q. (2010). Dysregulation of miR-15a and miR-214 in human pancreatic cancer. J. Hematol. Oncol. 3:46. doi: 10.1186/1756-8722-3-46

PubMed Abstract | CrossRef Full Text | Google Scholar