Skip to main content

Expression patterns of microRNAs associated with CML phases and their disease related targets

Abstract

Background

MicroRNAs are important regulators of transcription in hematopoiesis. Their expression deregulations were described in association with pathogenesis of some hematological malignancies. This study provides integrated microRNA expression profiling at different phases of chronic myeloid leukemia (CML) with the aim to identify microRNAs associated with CML pathogenesis. The functions of in silico filtered targets are in this report annotated and discussed in relation to CML pathogenesis.

Results

Using microarrays we identified differential expression profiles of 49 miRNAs in CML patients at diagnosis, in hematological relapse, therapy failure, blast crisis and major molecular response. The expression deregulation of miR-150, miR-20a, miR-17, miR-19a, miR-103, miR-144, miR-155, miR-181a, miR-221 and miR-222 in CML was confirmed by real-time quantitative PCR. In silico analyses identified targeted genes of these miRNAs encoding proteins that are involved in cell cycle and growth regulation as well as several key signaling pathways such as of mitogen activated kinase-like protein (MAPK), epidermal growth factor receptor (EGFR, ERBB), transforming growth factor beta (TGFB1) and tumor protein p53 that are all related to CML. Decreased levels of miR-150 were detected in patients at diagnosis, in blast crisis and 67% of hematological relapses and showed significant negative correlation with miR-150 proved target MYB and with BCR-ABL transcript level.

Conclusions

This study uncovers microRNAs that are potentially involved in CML and the annotated functions of in silico filtered targets of selected miRNAs outline mechanisms whereby microRNAs may be involved in CML pathogenesis.

Introduction

Mammalian microRNAs (miRNA, miR) are short non-coding RNAs that regulate preferentially gene expression by inhibiting translation of specific target mRNAs. MiRNA-mRNA matching is based on imperfect sequence base-pairing with the required complementarity centered over positions 2 - 8 of mRNA's seed sequence [1]. Depending on specific target genes, miRNAs regulate many cellular functions such as developmental timing, signal transduction, apoptosis, cell proliferation and tumorigenesis [25]. Thus, gene expression and role of miRNAs are currently being largely studied in human malignancies and chemical compounds that regulate miRNA levels are potentially very important for developing new treatment strategies in chronic myeloid leukemia (CML). The first miRNA molecules that have been associated with human leukemia pathogenesis were found in chronic lymphocytic leukemia (CLL) [6]. MiR-15 and miR-16 are located in a genomic region that is frequently deleted in CLL, thus the expression of these two miRNAs is downregulated. Other works brought the evidence that many miRNAs are indeed found at chromosomal breakpoints and genomic regions associated with cancer [7, 8].

In CML the following miRNAs were associated with the disease pathogenesis. For instance, the miR-203 was found to be epigenetically silenced in human leukemic Philadelphia chromosome-positive (Ph+) cell lines; this is in line with the observation that BCR-ABL and ABL kinases are miR-203 putative targets [9]. Derivative 9q+ chromosome deletions carrying miR-199b that occurred in some CML patients were associated with miR-199b decrease [10]. Venturini et al. [11] showed miR-17-92 cluster (onkomir-1) to be aberrantly expressed in CD34+ cells of CML patients. Agirre at al. [12] analyzed the expression of 157 miRNAs in mononuclear and CD34+ cells separated from bone marrow of 6 CML patients at diagnosis and found 11 miRNAs (e.g. miR-150, miR-151, miR-221, miR-127, miR-16) aberrantly expressed in CD34+ cells and 53 miRNAs differentially expressed in mononuclear cells (e.g. miR-150, miR-126, miR-221, miR-222, miR-21). Two recent works contributed to the knowledge about expression change in specific microRNAs associated with resistance to imatinib or responsiveness to imatinib after the treatment initiation in CML patients [13, 14]. A group of 19 miRNAs (e.g. miR-191, miR-29a, miR-422b, miR-100, miR-326, miR-26a) were identified as possible predictors for clinical resistance to imatinib in patients with newly diagnosed CML [13]. A relatively rapid increase in the expression of miR-150 and miR-146a and decrease of miR-142-3p and miR-199b-5p in peripheral blood mononuclear cells (PBMCs) of patients newly diagnosed with CML was found two weeks after imatinib initiation [14].

In this study, we used an array platform to characterize differentially expressed miRNAs in peripheral blood total leukocytes of patients at different stages of CML including diagnosis, major molecular response, therapy failure, hematological relapse, accelerated phase and blast crisis with the aim to identify microRNAs associated with pathogenesis of CML. To the best of our knowledge, such integrated microRNA profiling during the course of CML has not yet been performed. Hierarchical clustering analysis based on expression profiles of 49 miRNAs clearly separated patients at diagnosis, hematological relapse and blast crisis from those in major molecular response and therapy failure. We used in silico analyses to better understand the targets of 17 selected miRNAs whose deregulation was confirmed by real-time quantitative PCR (RT-qPCR). Based on our previous results demonstrating that miR-150 downregulation is associated with CML [15], we further validated miR-150 expression in a larger number of patients (n = 70). As MYB represents functionally validated target of miR-150 [16], its gene expression analysis was performed on the same patient cohort. Our data provide significant inverse correlations between miR-150 and MYB expression and BCR-ABL transcript level and indicate that this relationship is potentially important for pathogenesis in CML.

Materials and methods

Patient samples

Twenty four patient samples of total leukocytes from peripheral blood (Table 1) were used to prepare pools representing different CML phases for microarray analysis: diagnosis (n = 5, Dg), major molecular response (n = 5, MMR), therapy failure (n = 5, TF), hematological relapse (n = 5, Hr), and blast crisis (n = 4, BC). Briefly, Dg, Hr and BC contain 100% of Ph+ cells. Therapy failure is defined here as complete hematological response with failure to achieve complete cytogenetic remission (CCgR). Hematological relapse is defined as increased number of WBC (range 14-28*109/L). MMR samples are characterized as BCR-ABL <0.1% (IS). BC samples contain blast cells in peripheral blood from 50% to 79%. Eleven healthy donors of age median 60 (range 45 - 78) and man/woman ratio 3/2 following CML incidence were used to create a control pool.

Table 1 Characteristics of patient samples in the pools

Seventy patient samples of total leukocytes from peripheral blood were used for miR-150 expression validation and MYB expression analyses (Table 2). Of these, 13 represented Dg, 16 = MMR, 14 = TF, 15 = Hr and 12 = AP (accelerated phase) together with BC (AP/BC). Therapy failure is defined here as non CCgR achievement; all patients achieved complete hematological remission and two patients major and minimal cytogenetic response, respectively. Hematological relapse is characterized by increased number of WBC and PLT (median 16*109/L and 448*109/L, respectively).

Table 2 Characteristics of patient samples in groups for miR-150 and MYB expression analysis

The percentage of BCR-ABL transcript level was observed from the routine monitoring using real-time qPCR that is standardized within the frames of international standardization [17]. Mutation analyses were performed by direct sequencing method [18].

All subjects donated their samples with informed consent approved by the Ethic Committee of the Institute of Hematology and Blood Transfusion, Prague.

Sample preparation

Cell pooling was applied for microarray analysis at the aim to reduce individual variability and to find common features of the disease stage. Pooling strategy was performed according to previously described recommendations [19, 20]. Pools consisted of five patient samples; the blast crisis pool contained only four samples due to lack of appropriate material. Samples were selected for pooling according to their similar characteristics listed in Table 1. Each patient contributed to the pool by the same amount of leukocytes (107).

RNA extraction

Two different approaches were initially tested to extract total RNA containing small RNA molecules: acidic phenol-chloroform procedure and miRVana kit (Ambion®, Applied Biosystems, Foster City, CA, USA). The quality and quantity of RNA were evaluated using Agilent 2100 Bioanalyser (Agilent Technologies, Santa Clara, CA, USA) and spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA), respectively. The miRVana kit and phenol-chloroform procedure gave comparable results with respect to RNA quality (RNA integrity numbers 7.8 - 9.1 and 8.5 - 9.0, respectively) and quantity (mean total RNA amounts 5 ug and 3 ug, respectively).

MiRVana kit however gave better 230/280 ratio and therefore was then selected for preparation of samples for microarrays. Acidic phenol-chloroform extraction was used for real-time qPCRs, which is a standardized method in our laboratory for BCR-ABL monitoring in international scale (IS).

Microarray analysis

PIQOR™ miRXplore arrays (Miltenyi Biotech GmbH, Cologne, Germany) were used for miRNA expression profiling and the whole procedure including miRXplore data analysis (control vs. sample) was performed within the genomic facility of the manufacturer. Total RNAs with controlled quality (RIN 8 - 9.2; A260/A280 1.86 - 2.01; A230/A260 1.95 - 2.1) and quantity (1.2 - 4.0 μg) were sent to Miltenyi Biotech laboratory on dry ice. RNA quality and quantity was checked after delivery with the comparable results. Microarray platform contained 872 probes for human miRNAs according to miRBase version 10.1 and an extensive system of controls. Raw data were derived from ImaGene® software (Biodiscovery, El Segundo, USA). Only spots with signal equal to or higher than 50% percentile of the background signal intensities were further analyzed.

The complete microarray data were deposited in Gene Expression Omnibus (GEO) database under the accession number GSE26260 (http://www.ncbi.nlm.nih.gov/geo/).

We applied MultiExperiment Viewer (MeV v4.0 release; http://www.tm4.org/mev) for k-means/medians and hierarchical clustering was performed using average linkage and average dot product metric.

Real-time qPCRs

Real-time qPCR was performed on RotorGene 6000 (Qiagen, San Francisco, CA, USA). The miRNA expression assay kits (Applied Biosystems) specific for selected miRNAs were used to perform reverse-transcriptions and RT-qPCRs. MiR-30c showed stable expression across all the patient and control samples analyzed (stable Ct/ngRNA) and was used as a housekeeping gene for normalization. Relative fold changes of gene expression were assessed using 2-ΔΔCT method. Mean of ΔCT values of 11 healthy donors was used as a calibrator. Results are presented as expression fold change of a patient to a healthy control.

TaqMan Gene Expression Assay (product number Hs00193527; Applied Biosystems) was used for MYB transcript quantification according to the manufacturer's recommendations. The GUS gene was used as the housekeeping gene with the primer set, probe, and protocol adopted from Beillard et al. [21].

Prediction of putative miRNA target genes

The TargetScan Human release 5.1. (http://www.targetscan.org) was used for prediction of miRNA targets. Visualization and Integrated Discovery (DAVID) (http://david.abcc.ncifcrf.gov/) [22, 23] was applied to annotate the biological functions of the predicted targets.

Statistical analyses

Analyses of MYB and miR-150 differential expression between different groups of samples were conducted using Kruskall Wallis's test and Dunn's multiple comparison test. Correlation analyses were calculated using the Spearman's rho correlation test. Statistical analyses and graphs were performed using GraphPad Prison version 4.03 (GraphPad Software, La Jolla, CA, USA).

Results

miRNA expression profiles in CML

Microarray analysis in CML resulted in the detection of 56 differentially expressed miRNAs (samples vs. control). Figure 1 shows three main gene clusters (clusters I, II, III) of altogether 49 miRNAs. A markedly (more than 1.5 fold; others are not indicated here) increased level over the control was found in BC pool for miR-19a, miR-19b, miR-221, miR-126, miR-106a, miR-17, miR-20a and miR-222 that belong to the cluster I. A distant gene cluster II grouped down-regulated miRNAs in BC; more than 1.5 fold change was detected for miR-24, miR-29b, miR-26b, miR-107, miR-103, miR-150, miR-451. The cluster III consisted of miRNAs with increased level in MMR; more than 1.5-fold change was found for miR-663, miR-638 and miR-720.

Figure 1
figure 1

Hierarchical clustering analysis of expression data of 49 microRNAs from microarray analysis. Seven miRNAs are not displayed due to signal at the background level in at least one of the pools analyzed.

A separate cluster was created from MMR and TF pools and was distant from the cluster grouping Dg pool together with Hr and BC pools. MMR and TF pools are represented by the samples from patients during the imatinib treatment with optimal response and failure to achieve CCgR, respectively. The samples from MMR and TF pools were characterized by physiological blood count but different BCR-ABL transcript level and number of Ph+ metaphases (Table 1). The distant cluster of Hr, BC and Dg pools represented CML in progression and at diagnosis, respectively.

Validation of array data and sample pooling

Eight up-regulated (miR-19a, miR-19b, miR-221, miR-222, miR-126, miR-106a, miR-17 and miR-20a) and 3 down-regulated miRNAs (miR-103, miR-150 and miR-451) with more than 2.0-fold change in their expression (see Figure 1) over the control were selected for array data validation and for evaluation of pooling precision by the RT-qPCR. MiRNAs expressions were measured in the pools and in individual samples of each pool. The heatmaps showed comparable results (Figure 2); on comparing array and RT-qPCR data of pools and averaged RT-qPCR data of individuals, hierarchical clustering formed similar gene and sample clusters. We noted also few discrepancies e.g. for miR-126 and miR-451. We evaluated RT-qPCR data of each sample using hierarchical clustering (Figure 3). Seven miRNAs (miR-181a, miR-181b, miR-92a, miR-146a, let7c, miR-144, miR-155) that were not displayed in Figure 1 and 2 due to low signal on the array analysis in at least one of the pools, were included into the RT-qPCR analysis because of strong change in their expressions in BC pool. Three prominent patient clusters and three prominent miRNA clusters were identified (Figure 3). Firstly, a gene cluster distant from the other two consisted of miR-103, miR-150, miR-451 and miR-144. These molecules showed a rather decreased level at Dg, in Hr and BC. The other two closely related clusters consisted of miR-19b, miR-19a, miR-17, miR-20a, miR-92a, miR-106a, miR-222, miR-126, miR-146a, miR-181a, miR-181b, let7c, miR-155 and miR-221 that were up-regulated in BC samples.

Figure 2
figure 2

Hierarchical clustering analysis of miRNA expression data. (A.) microarray - pools, (B.) real-time qPCR - pools and (C.) real-time qPCR - individual samples.

Figure 3
figure 3

Hierarchical clustering analysis of detected levels of 17 selected miRNAs by real-time qPCRs.

The BC samples formed one cluster that was mixed with two samples from hematological relapse (Hr 4, 5), one from diagnosis (Dg 3) and one therapy failure (TF 4). This cluster was related to the cluster that was formed by Hr samples (Hr 1-3), two Dg samples (Dg 1, 4) and one therapy failure (TF 3). The third cluster distant from the other two consisted of all MMR samples, three TFs (TF 1, 2, 5) and was mixed with Dg 2 sample.

Target in silico analyses and functional annotation

The putative target genes were selected for 17 miRNAs with validated expression in different phases of CML (Figure 3). The predicted targets containing highly conserved sites were further studied according to the PCT values [24] from the TargetScan release 5.1 except for miR-106a because of lack data in the database. Because the PCT values are available only in the TargetScan we did not use another databases for miR-106a target prediction to preserve data consistency.

The targets were selected according to PCT equal or higher than 0.1 and 0.5 of miRNAs with low and high number of targets in the database, respectively (see Additional file 1: Table S1). The PCT, ranging between 0 and 1, corresponds to a Bayesian estimate of the probability that a site conserved to a particular branch length is conserved due to miRNA targeting [24]. The 30% of all of the targets with only conserved sites (Additional file 1: Table S1) are putative target genes of more than one of the selected miRNAs.

Functional annotation analysis of predicted targets (Table 3) revealed several biological processes (P < 0.0001). The encoded proteins are involved mainly in the regulation of transcription, intracellular signaling cascades, amino acid phosphorylation, regulation of RNA metabolic processes, regulation of apoptosis, regulation of cell proliferation and protein transport. Several proteins are implicated in hematopoietic or lymphoid organ development (n = 52; e.g. BMI1, WNT3A, MLL5, IL25, CDK6, MYH9, BCL2L11, CRKL, KIT, BCL2, RUNX1, TCF3, PIK3R1, NOTCH2, PKNOX1, SP1, SP3, TGFBR3), regulation of erythrocyte differentiation (ACVR2A, ACVR1B, ETS1, MAFB, SPI1, CDK6, FOXO3, INPP5D, ARNT) and regulation of myeloid cell differentiation (ZFP36, HMGB3, MAFB, KLF10, NDFIP1, SPI1, CDK6, FOXO3, PRDM16, PURB, ARNT, LIF, ACVR1B, ACVR2A, ID2, ETS1, GNAS, INPP5D, RUNX1).

Table 3 Functional annotation of predicted targets

Using KEGG database [25] we analyzed signaling pathways with significant hits (P < 0.0001) for predicted targets involved in endocytosis (hsa04144), pathways in cancer (hsa05200), mTOR signaling pathway (hsa04150), hedgehog signaling pathway (hsa04340), chronic myeloid leukemia (hsa05220), focal adhesion (hsa04510) and Wnt signaling (hsa04310) (Table 3).

Table 4 summarizes predicted targets associated with chronic myeloid leukemia. Most of them are involved in MAPK signaling (BCR, E2F2, E2F3, CBL, RAF1, CRK, CRKL, KRAS, SOS1, MAPK1). TGBR2, SMAD4 and ACVR1B play a role in transforming growth factor beta signaling pathway. Cyclin D1 (CCND1) and cyclin-dependent kinase 6 (CDK6) are important for the cell cycle and in the p53 pathway. Cell cycle is influenced by cyclin-dependent kinase inhibitor 1B (p27). ErbB signaling pathway encompasses PIK3R1 and PIK3R3. RUNX1 known as AML1 contributed to the abnormality in growth inhibition.

Table 4 Target annotation in pathways of chronic myeloid leukemia (hsa05220)

MiR-150 down-regulation and targeted MYB overexpression

Dramatic reduction of miR-150 in BC, at Dg and in Hr and its normal levels in patients under imatinib treatment (MMR and TF) (Figure 3) prompted us to determine this expression pattern on a larger cohort of patients (n = 70; Table 2). Significant down-regulation of the miRNA (p < 0.05) in comparison to healthy controls (n = 11) was confirmed for diagnosis and progressed phases of CML (Figure 4A.). MiR-150 level decreased more than 2-fold in 67% of hematological relapses (n = 10/15). There was no significant change in MMR and TF compared to controls. Among all patient samples analyzed, we found significant inverse correlation of miR-150 expression with BCR-ABL transcript level (p = 0.01; r = -0.501). To test whether miR-150 is regulated by BCR-ABL we have used a Ph+ cell line MOLM-7 and incubated it with two concentrations of imatinib (1uM and 10 uM, Additional file 2: Figure S1A) for total 48 hours. We observed that following reduction of BCR-ABL tyrosine kinase activity (exemplified by decreased intensity of p-CRKL (Additional file 2: Figure S1B)) by imatinib the miR-150 levels were significantly upregulated. This paragraph provides link between levels of one particular microRNA, miR-150, identified by our microarray analysis and CML pathogenesis.

Figure 4
figure 4

Expression analysis of miR-150 (A.) and its target MYB (B.) in different phases of CML in comparison to control (CONT). CONT = 11, Dg = 13 diagnosis, AP/BC = 12 accelerated phase/blast crisis, Hr = 15 hematological relapse, TF = 14 therapy failure (failure to achieve complete cytogenetic remission), MMR = 16 major molecular response. *** P < 0.001; ** P < 0.01; * P < 0.05

MYB is a confirmed target of miR-150 [16]. The expression pattern of MYB during the course of CML has not as yet been reported. We decided to analyze MYB transcript levels in the same cohort of patients. MYB was significantly increased at Dg, in AP/BC and Hr in comparison to controls (p < 0.001), to MMR (p < 0.05) and to TF (p < 0.05). Spearman's rho analysis displayed significant inverse correlation of MYB with miR-150 expression (p = 0.01; r = -0.521) and significant positive correlation between MYB expression and BCR-ABL transcript level (p = 0.01; r = 0.771).

Discussion

Specific microRNAs regulate hematopoietic cell differentiation and development [26]. The main interest is in whether there exists a link between levels of miRNAs and leukemia pathogenesis. The first work dealing with miRNA expression in CML demonstrated enhanced expression of the miR-17-92 cluster in CML CD34+ cells [11]. Other works that reported miRNA aberrant expression in CML appeared very recently. For example, it demonstrated that several miRNAs dysregulated in CML (miR-150, miR-146a, miR-142-3p, miR-199b-5p) were rapidly restored under imatinib treatment [14]. Several miRNAs (e.g. miR-191, miR-29a, miR-422b, miR-100, miR-326, miR-26a) are promising predictors of imatinib resistance in newly diagnosed CML patients [13].

This study investigates microRNA differential expression profiles that were initially analyzed at different stages of CML using microarrays. Pooling of patient samples was applied for microarray analysis to reduce individual variability and to find common features of the disease.

MiRNA array data showed similar expression pattern of 49 miRNAs in imatinib responders with MMR and patients with failure to achieve complete cytogenetic response. As expected, hierarchical clustering assembled the pools of samples at diagnosis, in hematological relapse and blast crisis, while MMR and TF pools formed a separate cluster (Figure 1). Total leukocytes from blast crisis peripheral blood that consisted of more than 50% blasts of each sample in the pool showed the highest number of strongly deregulated miRNAs.

We applied the functional annotation tool DAVID to look for the biological functions of predicted targets with only conserved sites and high PCT values of the 17 miRNAs with real-time qPCR-confirmed up-regulation (miR-19a, miR-19b, miR-17, miR-20a, miR-92a, miR-221, miR-222, miR-126, miR-146a, miR-181a, miR-181b, let7c and miR-155) and down-regulation (miR-103, miR-150, miR-451 and miR-144) in blast crisis (Figure 3). Several targets were involved in the processes that were found to be important in CML; endocytosis [27], mTOR signaling pathway [28, 29], hedgehog signaling [30, 31], focal adhesion [32, 33] and Wnt signaling [34, 35]. We summarized 19 genes with the probability to be targeted by miR-20a, miR-17, miR-19a, miR-103, miR-144, miR-150, miR-155, miR-181a, miR-221 and miR-222. The encoded proteins were annotated in pathways related to the CML (hsa05220). Out of these, 10 targets are involved in MAPK signaling (BCR, E2F2, E2F3, CBL, RAF1, CRK, CRKL, KRAS, SOS1, MAPK1). Interestingly, inhibition of MAPK signaling in Ph+ cell line K562 induced apoptosis [36]. Application of MAPK specific inhibitor U0126 showed synergistic effect with imatinib resulting in CD34+ progenitor reduction in CML [36].

Confirmed increase of miR-19a, miR-19b, miR-17, miR- 20a, miR-92a, miR-106a, miR-221, miR-222, miR-126, miR-146a, miR-181a, miR-181b, let7c and miR-155 was identified in samples of BC pool (Figure 3). This pattern was not found in Dg, Hr, TF or MMR pools. Overexpression of these miRs may be related to the immature character of blasts. Whether the increased level of these miRNAs may contribute to the CML pathogenesis or may simply reflect the stage of the disease is the matter of further investigation. Abnormal expression of onkomir miR-17-92 (miR-17; miR-19a; miR-19b; miR-20a; miR-92a) was described in CML CD34+ cells [11]. Agirre et al. [12] found up-regulated miR-221 and miR-222 in mononuclear cells of CML patients in comparison to healthy controls. MiR-155, miR-106a, miR-146a, miR-181 and miR-126 were reported as deregulated miRNAs in CML [13, 14]. To our knowledge, let7c expression has so far not been described in CML. In this study, our in silico analyses revealed that miR-221 and miR-103 (PCT 0.63 and 0.73, respectively) target PIK3R1. PIK3R3 is predicted to be regulated by miR-19a and miR-181a (PCT 0.92 and 0.81, respectively). PI3K is annotated in ERBB, MAPK and mTOR signaling pathways. KRAS, which is involved in MAPK signaling, is a predicted target of miR-19a (PCT 0.92). MAPK expression may be regulated by onkomirs miR-17 and miR-19a (PCT 0.96 and 0.86, respectively). Interestingly, it was reported that RAS/MAPK signaling may contribute to the survival of BCR-ABL positive cells under imatinib selection pressure [37]. AKT1, a member of the antiapoptotic PI3K pathway, is involved in both, BCR-ABL mediated transformation as well as in response to the BCR-ABL kinase inhibitors. It was shown that the PI3K/AKT/mTOR signaling is activated in imatinib naive cells while under imatinib pressure it may enhance resistance to imatinib [38]. As shown in our real-time qPCR data (Figure 3), the rather decreased levels of miR-181a, miR-221 and miR-19a in some imatinib treated patients, and miR-103 down-regulation in a number of blast crisis, diagnosis and progressed CML may contribute to the increased level of PI3K and thus may be involved in the previously described PI3K/AKT/mTOR signaling activation and in the resistance development in some CML cases. Though no experimental therapy using miRNA modulation has as yet provided significant and curative approach, the knowledge of deregulation of miRNAs specific for CML may facilitate the development of such therapeutic strategies. Several candidate microRNAs (e.g. miR-181a, miR-221, miR-19a, miR-103) regulating expression in CML target important signaling pathways may represent promising candidate targets for CML therapy.

The real-time qPCR validated the down-regulation of miR-150, miR-451, miR-103 and miR-144 overall in individual samples of BC, Hr, Dg pools and in some samples of TF pool (Figure 3). These molecules may be related to the CML pathogenesis and may reflect transformation from chronic to accelerated phases. Agirre et al. found miR-150 downregulation in mononuclear cells and CD34+ cells separated from bone marrow in newly diagnosed CML patients (n = 6) in comparison to healthy donors (n = 6) [12]. MiR-150 was recently described to be downregulated in untreated CML patients [14]. Flamant et al. [14] suggest that miR-150 play a role in leukemic cells and potentially in the more primitive hematopoietic compartment in chronic phase CML patients. This is in line with the knowledge that miR-150 is important in the regulation of hematopoiesis. During normal erythroid differentiation its level is gradually decreased [39], however; it shows the highest expression in mature lymphocytes [40]. Others proved that miR-150 expression increases during B-lymphoid differentiation in contrast to myeloid differentiation. It seems likely that miR-150 regulates the development of other two different blood lineages; B lymphocytes and megakaryocytes [41, 42]. Thus, miR-150 deregulation is found in hematological malignancies; miR-150 is decreased in polycytemia vera reticulocytes [43] and a marked decrease was recently also detected in MDS-del(5q) [44] while, in contrast, a twofold increase was found in CLL lymphocytes [45].

Based on our results [15] and recent results of others we expanded real-time qPCR assays of miR-150 on the larger cohort of CML patients. Decreased level of miR-150 was confirmed in patients at diagnosis, in the majority of patients with hematological relapse and in accelerated phase and blast crisis. Normal miR-150 level was observed in imatinib treated patients with major molecular response and failure to achieve CCgR. Our observations are consistent with the data of Flamant et al. [14] showing rapid increase of miR-150 expression after imatinib treatment initiation in patients with newly diagnosed CML. They further found that low miR-150 expression inversely correlated with white blood count and thus speculated that the level reflected the high leukocyte counts in newly diagnosed CML patients. We showed here a significant inverse correlation of miR-150 expression with BCR-ABL transcript level (p = 0.01; r = -0.501). Non-treated newly diagnosed patients, patients with disease progression and resistant to imatinib showed a high level of BCR-ABL together with high leukocyte count and decreased amount of miR-150. Normal miR-150 level was detected in imatinib responders (MMR) and patients with failure to achieve CCgR (TF) with normal blood count and low BCR-ABL transcript level. As imatinib targets Ph+ cells, the normal level of miR-150 in imatinib treated patients in chronic phase with physiological blood count could be the result of the suppression of leukemic cells and the concomitant recovery of normal hematopoiesis under imatinib treatment. Our in vitro tests showed elevated expression of miR-150 and marked decrease of p-CRKL following imatinib in vitro treatment of Ph+ cell line MOLM-7. These findings suggest a potential functional relationship between miR-150 and BCR-ABL.

Gene expression of MYB in our study showed a significant inverse correlation with miR-150 transcript level (p = 0.01; r = -0,409). MYB is the proven target of miR-150 and encodes a transcriptional factor required for proliferation and survival of normal and leukemic blast cells. A recently published study on a mouse model of blast crisis reported that c-MYB is required for BCR-ABL dependent leukemogenesis [46]. Lidonnici et al. [46] speculated that miR-150 reduction might contribute to the c-MYB upregulation that is likely induced by BCR-ABL, and may be involved in BCR-ABL driven leukemogenesis in CML. Interestingly, we found a significant correlation between MYB expression and BCR-ABL transcript level (p = 0.01; r = 0,782) in CML patients, which is in line with the above described suggestion.

In summary, our data demonstrated that miR-150, miR-20a, miR-17, miR-19a, miR-103, miR-144, miR-155, miR-181a, miR-221 and miR-222 are deregulated in CML. Furthermore, in silico filtering identified targeted genes that are involved in cell cycle, growth inhibition, MAPK, ErBb, transforming growth factor beta and p53 signaling pathways that are reported in CML pathogenesis. MiR-150 expression showed significant negative correlation with its target MYB and with BCR-ABL transcript level. The results of this study outline the mechanisms whereby miRNAs may be implicated in CML pathogenesis. However, if they function in BCR-ABL dependent or independent manner has to be elucidated.

References

  1. Boutla A, Delidakis C, Tabler M: Developmental defects by antisense-mediated inactivation of micro-RNAs 2 and 13 in Drosophila and the identification of putative target genes. Nuc Ac Res. 2003, 31: 4973-4980. 10.1093/nar/gkg707.

    Article  CAS  Google Scholar 

  2. Barbee SA, Estes PS, Cziko AM, Hillebrantd J, Luedeman RA, Coller JM, Johnson N, Howlett IC, Ceng C, Ueda R, Brand AH, Newbury SF, Wilhelm JE, Levine RB, Nakamura A, Parker R, Ramaswami M: Staufen- and FMRP containing neuronal RNPs are structurally and functionally related to somatic P bodies. Neuron. 2006, 52: 997-1009.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Bhattacharyya SN, Habermacher R, Martine U, Closs EI, Filipowicz : Relief of microRNA-mediated translational repression in human cells subjected to stress. Cell. 2006, 125: 1111-1124.

    Article  CAS  PubMed  Google Scholar 

  4. Carthew RW: RNA interference: the fragile X syndrome connection. Curr Biol. 2002, 12: R852-854.

    Article  CAS  PubMed  Google Scholar 

  5. Epanchinstev A, Jung P, Menssen A, Hermeking H: Inducible microRNA expression by an all-in-one episomal vector system. Nuc Ac Res. 34: e119-10.1093/nar/gkl624.

  6. Calin GA, Calin DD, Shimizu M, Bichi R, Zupo S, Noch E, Aldler H, Rattan S, Keating M, Rai K, Rassenti L, Kipps T, Negrini M, Bullrich F, Croce CM: Frequent deletions and down-regulation of micro-RNA genes miR15 and miR16 at 13q14 in chronic lymphocytic leukemia. PNAS. 2002, 26: 15524-15529. 10.1073/pnas.242606799.

    Article  Google Scholar 

  7. Calin GA, Liu CG, Svignani C, Ferracin M, Felli N, Dumitru CD, Shimizu M, Cimmino A, Zupo S, Dono M, Dell' Aquilla ML, Alder H, Rassenti L, Kipps TJ, Bullrich F, Negrini M, Groce CM: MicroRNA profiling reverses distinct signatures in B cell chronic lymphocytic leukemias. PNAS. 2004, 101: 11755-11760.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Zhang L, Huang J, Yang N, Greshock J, Megraw MS, Giannakakis A, Liang S, Naylor TL, Barchetti A, Ward MR, Yao G, Medina A, O'Brien-Jenkins A, Katsaros D, Hatzigeorgiou A, Gimotty PA, Weber BL, Coukos G: Micro-RNAs exhibit high frequency genomic alterations in human cancer. PNAS. 2006, 103: 9136-9141.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Bueno MJ, de Castro IP, de Cedrón MG, Santos J, Calin GA, Cigudosa JC, Croce CM, Piqueras JF, Malumbres M: Genetic and epigenetic silencing of microRNA-203 enhances ABL1 and BCR-ABL1 oncogene expression. Cancer Cell. 2008, 13: 496-506.

    Article  CAS  PubMed  Google Scholar 

  10. Albano F, Anelli L, Zagaria A, Liso V, Rocchi M, Specchia G: MIRN199B downregulation in chronic myeloid leukaemia is associated with deletions on der(9). BJ of Haematology. 2009, 144: 271-273. 10.1111/j.1365-2141.2008.07459.x.

    Article  Google Scholar 

  11. Venturini L, Battmer K, Castoldi M, Schultheis B, Hochhaus A, Muckenthaler MU, Ganser A, Eder M, Scherr M: Expression of the miR-17-92 polycistron in chronic myeloid leukemia (CML) CD34+ cells. Blood. 2007, 109: 4399-405.

    Article  CAS  PubMed  Google Scholar 

  12. Agirre X, Velasco A, San José-Enériz E, Garate L, Bandrés E, Cordeu L, Aparicio O, Saez B, Navarro G, Vilas-Zornoza A, Pérez-Roger I, García-Foncillas J, Torres A, Heiniger A, Calasanz MJ, Fortes P, Román-Gómez J, Prósper F: Down-regulation of hsa-miR-10a in chronic myeloid leukemia CD34+ cells increases USF2-mediated cell growth. Mol Cancer Res. 2008, 6: 1830-1840.

    Article  CAS  PubMed  Google Scholar 

  13. Enériz ESJ, Román-Gómez J, Jiménez-Velasco A, Garate L, Maritn V, Coredu L, Vilas-Zornoza A, Rodríguez-Otero P, Calasanz MJ, Pr'ésper F, Agirre X: MicroRNA expression profiling in imatinib-resistant chronic myeloid leukemia patients without clinically significant ABLI-mutations. Mol Cancer. 2009, 8: 69-72.

    Article  Google Scholar 

  14. Flamant S, Richie W, Guilhot J, Hols J, Bonnet ML, Chomel JC, Guilhot F, Turhan AG, Rasko JEJ: Micro-RNA response to imatinib mesylate in patients with chronic myeloid leukemia. Haematologica. 2010, 95: 1325-1333.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Poláková KM, Lopotová T, Klamová H, Moravcová J: Differential expression of miRNAs during the course of chronic myeloid leukemia. Blood. 2008, 112: s1082-

    Google Scholar 

  16. Barroga CF, Pham H, Kaushansky K: Thrombopoietin regulates c-Myb expression by modulating micro RNA miR-150 and high BIC/miR-155 expression. Exp Hematol. 2008, 36: 1585-1592.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Rulcová J, Zmeková V, Zemanová Z, Klamová H, Moravcová J: The effect of total-ABL, GUS and B2M control genes on BCR-ABL monitoring by real-time RT-PCR. Leuk Res. 2007, 31: 483-491.

    Article  PubMed  Google Scholar 

  18. Polaková KM, Lopotová T, Klamová H, Moravcová J: High-resolution melt curve analysis: Initial screening for mutations in BCR-ABL kinase domain. Leuk Res. 2008, 32: 1236-1243.

    Article  PubMed  Google Scholar 

  19. Kendziorski C, Irizarry RA, Chen KS, Haag JD, Gould MN: On the utility of pooling biological samples in microarray experiments. PNAS. 2005, 102: 4252-4257.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Peng X, Wood CL, Blalock EM, Chen KCh, Landfield PW, Stromberg AJ: Statistical implications of pooling RNA samples for microarray experiments. BMC Bioinf. 2003, 4: 26-35. 10.1186/1471-2105-4-26.

    Article  Google Scholar 

  21. Beillard E, Pallisgaard N, van der Velden VH, Bi W, Dee R, van der Schoot E, Delabesse E, Macintyre E, Gottardi E, Saglio G, Watzinger F, Lion T, van Dongen JJ, Hokland P, Gabert J: Evaluation of candidate control genes for diagnosis and residual disease detection in leukemic patients using 'real-time' quantitative reverse-transcriptase polymerase chain reaction (RQ-PCR) - a Europe against cancer program. Leukemia. 2003, 17: 2474-2486.

    Article  CAS  PubMed  Google Scholar 

  22. Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID Bioinformatics Resources. Nature Protocol. 2009, 4: 44-57.

    Article  CAS  Google Scholar 

  23. Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4: P3-

    Article  PubMed  Google Scholar 

  24. Friedman RC, Farh KKH, Burge ChB, Bartel DP: Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2009, 19: 92-105.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Kanehisa M, Goto S, Furumichi M, Tanabe M, Hirakawa M: KEGG for representation and analysis of molecular networks involving diseases and drugs. Nucleic Acids Res. 2010, 38: D355-D360.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Chen ChZ, Li , Lodish HF, Bartel DP: MicroRNAs modulate hematopoietic lineage differentiation. Science. 2004, 303: 83-86.

    Article  CAS  PubMed  Google Scholar 

  27. Eisendle K, Lang A, Eibl B, Nachbaur D, Glassl H, Fiegl M, Thaler J, Gastl G: Phenotypic and functional deficiencies of leukaemic dendritic cells from patients with chronic myeloid leukaemia. BJ of Haematology. 2003, 120: 63-73. 10.1046/j.1365-2141.2003.03979.x.

    Article  Google Scholar 

  28. Carayol N, Vakana E, Sassano A, Kaur S, Goussetis DJ, Glaser H, Druker BJ, Donato JN, Altman JK, Barr S, Platanias LC: Critical roles for mTORC2- and rapamycin-insensitive mTORC1-complexes in growth and survival of BCR-ABL-expressing leukemic cells. PNAS. 2010, 107: 12469-12474.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  29. Mayerhofer M, Aichleberger KJ, Florian S, Krauth MT, Hauswirth AW, Derdak S, Sperr WR, Esterbauer H, Wagner O, Marosi Ch, Pickl WF, Deininger M, Weisberg E, Druker BJ, Griffin JD, Sillaber Ch, Valent P: Identification of mTOR as a novel bifunctional target in chronic myeloid leukemia: dissection of growth-inhibitory and VEGF-suppressive effects of rapamycin in leukemic cells. FASEB. 2005, 19: 960-962.

    CAS  Google Scholar 

  30. Zhao C, Chen A, Jamieson CH, Fereshteh M, Abrahamsson A, Blum J, Kwon HY, Chute JP, Rizzieri D, Munchhof M, VanArsdale T, Beachy PA, Reya T: Hedgehog signaling is essential for maintenance of cancer stem cells in myeloid leukaemia. Nature. 2009, 458: 776-779.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Dierks C, Beigi R, Guo GR, Zirlik K, Stegert MR, Manley P, Trussell C, Shcmitt-Graeff A, Landwerlin K, Veelken H, Warmuth M: Expansion of Bcr-Abl-positive leukemic stem cells is dependent of Hedgehog pathway activation. Cancer Cell. 2008, 14: 238-249.

    Article  CAS  PubMed  Google Scholar 

  32. Brunton VG, Avizienyte E, Fincham VJ, Serrels B, Metcalf ChA, Saqyer TK, Frame MC: Identification of Src-specific phosphorylation site on focal adhesion kinase: dissection of the role of Src SH2 and catalytic function and their consequences for tumor cell behavior. Cancer Res. 2005, 65: 1335-1342.

    Article  CAS  PubMed  Google Scholar 

  33. Salgia R, Brunkhorst B, Pisick E, Li JL, Lo SH, Chen LB, Griffin JD: Increased tyrosine phosphorylation of focal adhesion proteins in myeloid cell lines expressing p210BCR/ABL. Oncogene. 1995, 11: 1149-1155.

    CAS  PubMed  Google Scholar 

  34. Desphande AJ, Buske Ch: Knocking the Wnt out of the sails of leukemia stem cell development. Cell Stem Cell. 2007, 1: 597-598.

    Article  Google Scholar 

  35. Ress A, Moelling K: Bcr is a negative regulator of the Wnt signalling pathway. EMBO. 2005, 6: 1095-1100. 10.1038/sj.embor.7400536.

    Article  CAS  Google Scholar 

  36. Kang ChD, Yoo SD, Hwanga BW, Kima KW, Kimc DW, Kima ChM, Kima SH, Chunga BS: The inhibition of ERK/MAPK not the activation of JNK/SAPK is primarily required to induce apoptosis in chronic myelogenous leukemic K562 cells. Leuk Res. 2000, 24: 527-534.

    Article  CAS  PubMed  Google Scholar 

  37. Chu S, Holtz M, Gupta M, Bhatia R: BCR/ABL kinase inhibition by imatinib mesylate enhances MAP kinase activity in chronic myelogenous leukemia CD34+ cells. Blood. 2004, 103: 3167-3174.

    Article  CAS  PubMed  Google Scholar 

  38. Burchert A, Wang Y, Cai D, von Bubnoff N, Paschka P, Mueller-Bruesselbach S, Ottmann OG, Duyster J, Hochhasou A, Neubauer A: Compensatory PI3-kinase/Akt/mTor activation regulates resistance development. Leukemia. 2005, 19: 1774-1782.

    Article  CAS  PubMed  Google Scholar 

  39. Bruchova H, Yoon D, Agarwal AM, Mendell J, Prchal JT: Regulated expression of microRNAs in normal and polycythemia vera erythropoiesis. Exp Hematol. 2007, 35: 1657-1667.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Zhou B, Wang S, Mayer C, Bartel DP, Lodish HF: MiR-150, a micro RNA expressed in mature B and T cells, blocs early B cell development when expressed prematurely. Proc Natl Acad Sci. 2007, 104: 1080-7085.

    Google Scholar 

  41. Xiao C, Calado DP, Galler G, Thai TH, Patterson HCh, Wang J, Rajewsky N, Bender TP, Rajewsky K: MiR-150 controls B-cell differentiation by targeting the transcription factor c-Myb. Cell. 2007, 131: 146-159.

    Article  CAS  PubMed  Google Scholar 

  42. Lu J, Guo S, Ebert BL, Zhang H, Peng Y, Bosco J, Pretz J, Schlanger R, Wang YJ, Mak RH, Dombkowski DM, Preffer FI, Scadden DT, Golub TR: MicroRNA-mediated control of cell fate in megakaryocyte-erythrocyte progenitors. Developmental Cell. 2008, 14: 843-853.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  43. Bruchova H, Yoon D, Agarwal AM, Swierczek S, Prchal JT: Erythropoiesis in polycythemia vera is hyper-proliferative and has accelerated maturation. Blood Cells, Molecules, and Diseases. 2009, 43: 81-87.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  44. Hussein K, Theophile K, Busche G, Schlegelberger B, Gohring G, Kreipe H, Bock O: Significant inverse correlation of microRNA-150/MYB and microRNA-222/p27 in myelodysplastic syndrome. Leuk Res. 2010, 34: 328-334.

    Article  CAS  PubMed  Google Scholar 

  45. Fulci V, Chiaretti S, Goldoni M, Azzalin G, Carucci N, Tavolaro S, Castellano L, Magrelli A, Citarella F, Messina M, Maggio R, Peragine N, Santangelo S, Mauro FR, Landgraf P, Tuschl T, Weir DB, Chien M, Russo JJ, Ju J, Sheridan R, Sander Ch, Zavolan M, Guarini A, Foa` R, Macino G: Quantitative technologies establish a novel microRNA profile of chronic lymphocytic leukemia. Blood. 2007, 109: 4944-4951.

    Article  CAS  PubMed  Google Scholar 

  46. Lidonnici MR, Corradini F, Waldron T, Bender TP, Calabretta B: Requirement of c-Myb for p210BCR/ABL-dependent transformation of hematopoietic progenitors and leukemogenesis. Blood. 2008, 111: 4771-4779.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements and Funding

We are grateful to the panel of experts from the CELL study group. We would like to acknowledge Hana Bruchová-Votavová for her valuable advice and critical revision of the manuscript, Bohumír Procházka for statistical analysis, and Kateřina Vlčanová for technical assistance. We are thankful to Prof. Kyra Michalová and Jana Březinová for cytogenetic data. We thank to Silvia Rüberg (Miltenyi Biotech) for technical support. We also thank to Zuzana Ondráčková for westernblott analysis, Markéta Žáčková for the cultivation protocol, and to Jitka Němcová and to Petr Stöckbauer for providing of MOLM-7 cell culture. TS is supported by grants: IGA (NS10310-3/2009), MSMT (NPVII 2B06077), MPO (FR-TI2/509), PB is supported by GAUK (251135 82210).

This work has been supported by GAČR GP301/08/P154 and MZOUHKT2005 (00023736).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Kateřina Machová Poláková.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

KMP - conception and design, in silico analyses, data evaluation and interpretation, manuscript drafting; TL - contribution to the manuscript drafting, real-time qPCR analyses and evaluation; HK - provision of patient samples, clinical data evaluation, critical revisions; PB - in vitro tests, contribution to data evaluation and interpretation, critical revision; MT - clinical data revision; critical revisions of the manuscript; TS - contribution to data evaluation and interpretation, critical revision; JM - critical revisions. All authors read and approved the final manuscript.

Electronic supplementary material

12943_2010_877_MOESM1_ESM.DOC

Additional file 1:Table S1: Number of predicted targets of conserved miRNA families. Table summarizes number of the microRNA targets that were selected according to PCT equal or higher than 0.1 and 0.5. (DOC 34 KB)

12943_2010_877_MOESM2_ESM.PPT

Additional file 2:Figure S1: MiR-150 exression is elevated and BCR-ABL activity is dropped after imatinib treatment in Ph+ MOML-7 cells. Data from in vitro test of miR-150 expression change after Ph+ cell line incubation with imatinib. (A.) MiR-150 expression change after imatinib treatment in Ph+ MOML-7 cells. 3 × 106 MOLM-7 cells were incubated for 24 or 48hrs with or without (white bars; CTRL = control) imatinib. Two different concentrations of imatinib were tested (1 μM - gray bars and 10 μM - black). Cellular RNA was isolated by Trizol (Invitrogen), transcribed using High Capacity cDNA Reverse Transcription Kit (Roche Diagnostics). Real-time qPCR was performed using TaqMan protocol (Roche Diagnostics) and was run on the ABI 7900HT instrument. RNU44 was used as housekeeping gene. Data were evaluated by 2-ΔCt method. The viability of culture with imatinib decreased: 1 μM imatinib- from 94% after 24hrs to 21% after 48hrs; 10 μM imatinib- from 95% after 24hrs to 18% after 48hrs. (B.) The intensity change of p-CRKL after imatinib treatment in Ph+ MOLM-7 cells. The amount of p-CRKL (a client molecule of BCR-ABL tyrosine kinase) was measured by standard western blot analysis (p-CRKL (Tyr207) antibody; Cell Signaling Technology) after 24h culture cultivation with imatinib using both concentrations. Beta-actin (monoclonal antibody Anti-beta-Actin, Sigma) was used as the loading control and was measured by western blot analysis using alkaline phosphatase. The measurement of p-CRKL was not possible to perform in the culture after 48 incubation due to a marked viability decrease and thus to low amount of material. (PPT 222 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Machová Poláková, K., Lopotová, T., Klamová, H. et al. Expression patterns of microRNAs associated with CML phases and their disease related targets. Mol Cancer 10, 41 (2011). https://doi.org/10.1186/1476-4598-10-41

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1476-4598-10-41

Keywords