Overexpression of primary microRNA 221/222 in acute myeloid leukemia
© Rommer et al.; licensee BioMed Central Ltd. 2013
Received: 7 February 2013
Accepted: 10 July 2013
Published: 29 July 2013
Acute myeloid leukemia (AML) is a hematopoietic malignancy with a dismal outcome in the majority of cases. A detailed understanding of the genetic alterations and gene expression changes that contribute to its pathogenesis is important to improve prognostication, disease monitoring, and therapy. In this context, leukemia-associated misexpression of microRNAs (miRNAs) has been studied, but no coherent picture has emerged yet, thus warranting further investigations.
The expression of 636 human miRNAs was compared between samples from 52 patients with AML and 13 healthy individuals by highly specific locked nucleic acid (LNA) based microarray technology. The levels of individual mature miRNAs and of primary miRNAs (pri-miRs) were determined by quantitative reverse transcriptase (qRT) PCR. Transfections and infections of human cell lines were performed using standard procedures.
64 miRNAs were significantly differentially expressed between AML and controls. Further studies on the clustered miRNAs 221 and 222, already known to act as oncogenes in other tumor types, revealed a deficiency of human myeloid cell lines to process vector derived precursor transcripts. Moreover, endogenous pri-miR-221/222 was overexpressed to a substantially higher extent than its mature products in most primary AML samples, indicating that its transcription was enhanced, but processing was rate limiting, in these cells. Comparison of samples from the times of diagnosis, remission, and relapse of AML demonstrated that pri-miR-221/222 levels faithfully reflected the stage of disease.
Expression of some miRNAs is strongly regulated at the posttranscriptional level in AML. Pri-miR-221/222 represents a novel molecular marker and putative oncogene in this disease.
KeywordsAML miR-221 pri-miRNA lncRNA
Acute myeloid leukemia (AML) is a frequently fatal malignant disease of hematopoietic stem and progenitor cells (HSPCs). Prognostic factors include patient age, antecedent hematological disease, preceding cytotoxic treatments for a primary disorder, and the presence of specific cytogenetic, molecular, and epigenetic aberrations [1–7]. Identification and investigation of such somatic genetic alterations has enhanced our understanding of disease biology, augments diagnosis and prognostication, and may aid monitoring of the course of disease [1, 3, 7, 8]. Although cytogenetics and molecular genetics have already facilitated great progress in these respects, novel technologies like next generation sequencing and biological breakthroughs like the identification of microRNAs (miRNAs) and long noncoding RNAs (lncRNAs) as entirely novel gene classes still provide fundamental additional insights. In this context, the expression, functions, and potential prognostic value of miRNAs in AML have been studied by a number of research groups over the past few years [9–15]. In contrast, little is known about the expression and roles of lncRNAs in this disease.
miRNAs are small (~22 nucleotide, nt) RNA molecules which typically regulate multiple target genes at the levels of mRNA stability and translation efficiency [16, 17]. They are excised through sequential processing steps from larger, usually polymerase II transcribed, RNA molecules termed primary miRNAs (pri-miRNAs). One pri-miRNA molecule may give rise to one or several miRNA species, with the respective genomic region referred to as a miRNA cluster in the latter case . Classical pri-miRNA processing consists of two endonucleolytic cleavage steps: the first one is carried out in the nucleus by the RNase Drosha and gives rise to ~70 nt long, hairpin shaped precursor (pre-) miRNA molecules. These are exported to the cytoplasm and further processed by another RNase, Dicer, to yield mature miRNAs [18–21]. After incorporation into the RNA induced silencing complex (RISC), miRNAs recognize their target mRNAs through a 7–8 nt seed sequence and repress their translation and stability [16, 17, 20, 21]. In addition to this well described way of action of mature miRNAs, biochemical and biological functions have recently been ascribed to some pri-miRNAs [22, 23], thus greatly expanding the potential biological relevance of miRNA genes.
Large scale expression analyses of mature miRNAs revealed specific miRNA patterns that were associated with recurrent genetic abnormalities in AML [10–12, 24]. Certain miRNAs were deregulated in AML compared to healthy bone marrow (BM), peripheral blood (PB), or CD34 positive (CD34+) HSPCs [10–12, 24, 25], and miRNA signatures predictive of survival have been identified [9, 10]. However, results from different research groups exhibited only limited concordance, most likely due to differences between patient collectives, types of healthy control, and methodologies .
In this study, we re-addressed the question which miRNAs are differentially expressed between AML and healthy controls, and therefore might contribute to leukemogenesis. Using highly specific locked nucleic acid (LNA) oligonucleotide arrays capable of detecting 559 annotated and 77 proprietary human miRNAs, 64 miRNAs were found to be significantly misexpressed in AML. Further studies on the clustered miRNAs 221 and 222 revealed that pri-miR221/222 was overexpressed to a much higher extent than its mature products, which is best explained by an increased rate of transcription on the background of a limited processing capacity. Because pri-miR-221/222 is strongly overexpressed in a large proportion of patients with AML, it may represent a novel molecular marker and a putative oncogene in this disease.
Primary samples from patients and controls
This study was approved by the Ethics Committee of the Medical University of Vienna (EK no 609/2011). In some cases, archival samples were used. Written informed consent was obtained prior to sample collection from all other subjects. Data were analyzed anonymously.
For miRNA microarray hybridization, 52 AML PB samples and 13 control samples were used (Additional file 1: Tables S1A and B). The AML samples had been collected at the time of diagnosis, strongly enriched for leukemic blasts using Ficoll, and vitally frozen. Of the healthy control samples, three were CD34+ HPSCs enriched from bone marrow (BM) and five were BM mononuclear cells (MNCs) (all purchased from Lonza, Basel, Switzerland). Another five were PB samples from healthy volunteers, and were subjected to erythrocyte lysis prior to RNA isolation (Additional file 1: Table S1A). Twenty-two and 9 of these AML and control samples, respectively, were also used for qRT-PCR for mature miR-221 and 222. Characteristics of 27 additional diagnostic AML samples used to determine the ratio between pri-miR-221/222 and mature miR-221, as well as the expression of pri-miR-221/222 and other pri-miRNAs are summarized in Additional file 1: Table S1C. Eight additional control samples used for qRT-PCR experiments are included in Additional file 1: Table S1A. To investigate pri-miR-221/222 expression during the course of disease, four paired samples from the time of diagnosis and remission and three paired samples from the time of diagnosis and relapse were analyzed (Additional file 1: Table S1D).
Production and hybridization of miRNA microarrays, and primary data analysis
miRNA microarrays were produced as previously described [26, 27]. Briefly, LNA modified oligonucleotide probes specific for the 559 human miRNAs contained in miRBase version 9.2 (http://www.mirbase.org/) as well as probes for 77 proprietary miRPlus sequences (Exiqon, Vedbaek, Denmark) were spotted onto Hisens epoxy-coated glass slides (Schott Nexterion, Louisville, KY, USA) in eight replicates.
RNA was extracted with Trizol (Life Technologies, Carlsbad, CA, USA), subjected to quality control using a Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA), and labelled with Hy3 using the miRCURY™ LNA microRNA array labelling kit (Exiqon). A Hy5 labelled mix of all samples was included in each hybridization as a common reference. Arrays were hybridized overnight at 60°C in microarray hybridization chambers (Corning, Corning, NY, USA). Hybridization and wash buffers from the miRCURY LNA microRNA Array Kit (Exiqon) were used according to the manufacturer’s instructions. Arrays were scanned with a GenePix 4100A Microarray Scanner and evaluated with GenpixPro 5.1 software (Molecular Devices, Sunnyvale, CA, USA). Primary data analysis was performed using ArrayNorm . Features were filtered for low quality spots, the local background was subtracted, spots were normalized to the global mean, and the ratios between the sample of interest and the common reference were log2 transformed.
Cell lines, miRNA expression vectors, stable transfections, and infections
KG1 and KG1a cells  were obtained from the German Collection of Microorganisms and Cell Cultures, Braunschweig, Germany. HNT34  and MPD  cells were kindly provided by Dr. Hiroyuki Hamaguchi, Musashino Red Cross Hospital, Tokyo, Japan, and Dr. Cassandra Paul, Wright State University, Dayton. Ohio, USA, respectively. HL60 , HEL , U937 , HeLa , MCF7 , and 293 T  cells were obtained from cell culture collections of the Medical University of Vienna. All cell lines were regularly tested for mycoplasma contamination.
The human myeloid cell lines KG1, HL60, HEL, U937, MPD, and HNT34 were cultured in Roswell Park Memorial Institute (RPMI) 1640 medium supplemented with 10%Fetal Bovine Serum (FBS) and 1%Penicillin-Streptomycin-Glutamine (PSG; all from Life Technologies, Carlsbad, CA, USA) in a humidified incubator at 37°C and 5%CO2. The same media were used for KG1a cells, except that they contained 20%FBS. The adherent cell lines HeLa, MCF7, and 293 T were maintained in Dulbecco’s Modified Eagle Medium (DMEM; Life Technologies) with 10%FBS and 1%PSG.
Plasmid miR-Vec-221/222, which contains the miR-221/222 cluster along with a blasticidin resistence gene in the pMSCV backbone, as well as the corresponding empty vector (miR-Vec) were kindly provided by the Agami lab . They were stably transfected into HL60 cells by electroporation and transfectants were selected using 2 μg/ml Blasticidin (Invivogen, San Diego, CA, USA).
pEZX-MR03-miR-221 (Homo sapiens microRNA miR-221 stem-loop expression clone, #HMIR0369), an HIV based lentiviral vector containing the miR-221 precursor, and pEZX-MR01-control (#CMIR0001-MR01), which contains a scrambled sequence instead, were obtained from GeneCopoeia (Rockville, Maryland, USA). They were infected into the myeloid cell lines HL60, KG-1, and KG-1a using standard procedures . Infected cells were sorted for GFP positivity on a FACS Aria (BD Biosciences, NJ, USA).
RNA extraction, cDNA synthesis, and quantitative reverse transcriptase PCR (qRT-PCR)
Total RNA was isolated from primary samples and cell lines using Trizol (Life Technologies). For detection of pri-miRNAs, RNA was treated with DNase I and converted to cDNA with M-MLV reverse transcriptase primed by random hexamer oligonucleotides (all reagents from Life Technologies). qRT-PCR was performed in an ABI Step One Plus sequence detection system (Applied Biosystems, Life Technologies) using the Mesa Green qPCR Master Mix Plus (Eurogentec, Liège, Belgium) and the primers listed in Additional file 2: Table S2A (synthesized by MWG Eurofins, Ebersberg, Germany). All primer pairs were subjected to standard curve analysis and yielded slopes between -3.0 and -3.5, indicating optimal or near optimal amplification efficiencies. Levels of mature miR-221 and miR-222 and of RNU6B were measured using Taqman assays (ID000524, hsa-miR-221 TaqMan Assay; ID002276, hsa-miR-222 TaqMan Assay; ID001093, mature miR control-RNU6B TaqMan Assay; Applied Biosystems). qRT-PCR reactions were performed in triplicate. The relative expression of pri-miRNAs compared to the housekeeping gene beta-2-microglobulin, and of mature miRNAs relative to RNU6B, were calculated according to the ΔΔCt method .
RNA deep sequencing analysis and prediction of putative pri-miR 221/222 transcripts
RNA deep sequencing optimized for detection of lncRNAs was performed as described . In short, total RNA was extracted from human Hs27 foreskin fibroblasts with TRIreagent (Sigma-Aldrich, Seelze, Germany), treated with DNaseI (DNA-free kit, Ambion, Life Technologies), depleted of ribosomal RNA using RiboMinus Transcriptome Isolation Kit Human/Mouse (Life Technologies) and Ribo-Zero rRNA Removal Kit Human/Mouse/Rat (Epicentre Biotechnologies, Madison, WI, USA), and fragmented by hydrolysis. Double stranded cDNA was generated with SuperScript II Reverse Transcriptase (Life Technologies). The RNA-Seq library was prepared using the Chip-Seq DNA Sample Prep Kit (Illumina, San Diego, CA, USA), and sequenced using an Illumina Genome Analyzer II and Illumina HiSeq 2000 systems. 36 bp and 51 bp single end reads were aligned to human genome build hg18 using Bowtie (http://bowtie-bio.sourceforge.net/index.shtml) and visualized on the University of California Santa Cruz (UCSC) genome browser (http://genome.ucsc.edu/) (Vlatkovic IM, manuscript in preparation), together with publically available global run-on sequencing (GRO-Seq) data , ENCODE histone modification Chip-Seq data generated by the BROAD Institute, and RefSeq Genes.
Cloning of reporter vectors, transient transfections, and luciferase reporter assays
Two fragments surrounding the transcription start site predicted for the putative 28.2 kb pri-miR-221/222 transcript were amplified from PB leukocyte genomic DNA using the primers shown in Additional file 2: Table S2B and Phusion High Fidelity Polymerase (New England Biolabs, Ipswich, MA, USA). The resulting PCR products were ligated into the pGL3-Promotor (pGL3-P) vector (Promega, Madison, WI, USA) using the SacI and XhoI restriction sites engineered onto the PCR primers to generate reporter vectors pGL3-P(-1874/+45) and pGL3-P(+17/+1952).
Subconfluent cultures of 293 T, HeLa and MCF7 cells growing in 24-well plates were transfected with 1 μg of the indicated pGL3-P vector derivative and 30 ng of the renilla luciferase vector pGL 4.70 (Promega, Madison, WI, USA), using JetPei Transfection Reagent (Polyplus Transfections, Illkirch, France) according to the manufacturer’s instructions. 48 h later, luciferase activities were measured using a Tristar LB941 (Berthold Technologies, Bad Wildbad, Germany) and the Dual-Luciferase Reporter Assay (Promega, Madison, WI, USA). Firefly luciferase activity was normalized to Renilla luciferase activity to control for transfection efficiency.
For statistical analyses of miRNA microarray data, only features with a valid human miRNA annotation, and which were detectable in at least half of the relevant samples (i.e., 27, 2, 3, and 3 of the AML, control CD34+, control BM, and control PB samples, respectively) were considered. One-way ANOVA was performed to test if a miRNA was differentially expressed in one of the classes, and a moderated t-test (R/Biconductor package limma) was used to test if a miRNA was differentially expressed between two respective classes. Resulting p-values were adjusted for multiple hypothesis testing based on the false discovery rate (FDR) .
Spearman’s rank correlation was used to determine associations between miRNA expression as a continuous variable and the clinical parameters age, white blood cell count, percentage of blasts, lactate dehydrogenase (LDH), and cytogenetic risk. The point biserial correlation coefficient was used to show associations between miRNA expression and the dichotomous parameters gender and achievement of complete remission (CR), and significance was probed using the Wilcoxon rank-sum test. Associations between miRNA expression levels and the FAB type were given by the eta coefficient and significance was calculated by one-way ANOVA. All p-values were adjusted for multiple hypothesis testing based on the FDR . All calculations were performed using R.
For experiments using cultured cells at least three independent biological replicates were performed. Results are expressed as means ± standard error of the mean (SEM), or, where indicated, one representative experiment with standard deviations (SDs) from technical replicates is shown. Student’s two-tailed t-test at a significance level of 0.05 was employed to probe differences between groups for statistical significance.
Microarray analyses identify miRNAs deregulated in AML compared to healthy controls
To compare miRNA expression profiles between AML and normal controls, RNA from 52 AML samples (Additional file 1: Table S1B; enriched for leukemic blasts through Ficoll purification) and 13 healthy donors (5 peripheral blood (PB), 5 bone marrow (BM), and 3 BM CD34+; Additional file 1: Table S1A) was hybridized to microarrays containing LNA probes specifically detecting 559 human miRNAs deposited in miRBase 9.2 as well as 77 proprietary miRPlus sequences. Only miRNAs that were expressed in more than half of each the AML, healthy PB, healthy BM, and healthy CD34+ samples were considered for the respective comparisons. At a false discovery rate (FDR) of <5%, 3 miRNAs were up- and 10 downregulated in AML versus control CD34+ cells. 23 miRNAs were up- and 22 downregulated in AML compared to healthy BM, and 17 miRNAs were expressed at higher and 19 at lower levels in AML than in normal PB. 22 of these miRNAs were deregulated compared to two out of the three control tissues, and only four showed significantly altered expression with respect to all three types of control (Additional file 3: Table S3).
miRNA expression patterns were also tested for associations with preclinical and clinical parameters, namely, sex, age, French-Amercian-British (FAB) type, white blood cell count (WBC), blast percentage, lactate dehydrogenase (LDH) levels, cytogenetic risk, and achievement of complete remission (CR). At an FDR of <5%, significant associations were found between miRNA levels and blast percentage, cytogenetic risk, and sex (Additional file 4: Table S4).
Myeloid leukemia cell lines and primary AML samples display a limited pri-miR-221/222 processing capacity
To study the biological effects of miR-221 and miR-222 in human myeloid cells, we introduced suitable expression vectors into the cell lines HL60, KG1, and KG1a. However, neither blasticidin resistant HL60 clones obtained after transfection with miR-Vec-221/222 , nor GFP positive HL60, KG1, or KG1a cell populations sorted after infection with pEZX-MR03-miR-221 exhibited increased expression of mature miR-221 or miR-222 (Additional file 5: Figure S1A, and data not shown). In contrast, miR-221 was clearly elevated after transient transfection of HeLa cells with pEZX-MR03-miR-221 (Additional file 5: Figure S1B), demonstrating that this vector was functional. Furthermore, qPCR showed that, other than the mature miR-221, its precursor was indeed overexpressed in the pEZX-MR03-miR-221 infected myeloid cell lines (Additional file 5: Figure S1C), indicating that these cell lines are able to transcribe pre-miR-221 from the pEZX-MR03 vector, but process it inefficiently, possibly due to a saturation effect.
pri-miR-221/222 is universally overexpressed in AML
pri-miRNAs are not generally overexpressed in AML
The pri-miR-221/222 transcript overexpressed in AML is at least 23 kb long
To investigate whether any of these transcripts was expressed in AML, we designed primer pairs for regions near the 5’-ends of each of them, and another pair for a region 3’ of miR-221 (Figure 5A). qRT-PCR on 3 AML and 2 healthy control PB samples showed that the amplicons near the predicted 5’-ends of the 5.6 and 28.2 kb transcripts were strongly elevated in AML, while the amplicon near the 5’-end of the putative 108.5 kb transcript and the 3’ amplicon displayed inconsistent behaviour (Figure 5B-F). Even though these experiments do not formally prove that the amplicons exhibiting similar expression patterns belong to the same transcript, this interpretation is supported by the RNA sequencing data, which strongly suggest the expression of one continuous RNA from this region, and by the fact that according to the UCSC genome browser no other known genes or transcripts are located there. Based on the positions of the primers used (Figure 5A), our results indicate that the pri-miR-221/222 transcript overexpressed in AML starts no more than 2 kb downstream of the 5’-end of the putative 28.2 kb transcript and ends less than 1.5 kb downstream of the miR-221 sequence, making for a total length of at least 23 kb.
To confirm that the region surrounding the presumed pri-miR-221/222 transcriptional start site (TSS) does indeed have transcription activating potential, two fragments encompassing the regions from -1874 to +45 and from +17 to +1952 relative to the estimated start site of the predicted 28.2 kb transcript were cloned into the luciferase reporter vector pGL3-Promoter to yield plasmids pGL3-P(-1874/+45) and pGL3-P(+17/+1952). Because human hematopoietic cell lines are notoriously difficult to transfect, luciferase assays were performed in readily transfectable adherent cell lines. The luciferase activity of pGL3-P(+17/+1952) was approximately 13-fold higher than that of empty pGL3-P in MCF7 cells, while the upstream fragment increased the transcriptional activity of the reporter only about twofold (Figure 5G). Comparable results were obtained in 293 T and HeLa cells (Additional file 6: Figure S2). These data confirm that a region in the immediate vicinity of the predicted pri-miR-221/222 TSS can indeed promote transcription in human cells.
pri-miR-221/222 levels reflect the disease stage in AML
To identify miRNAs potentially contributing to leukemogenesis, the miRNA expression profiles of 52 AML samples were compared to those of hematopoietic cells from 13 healthy donors. Using microarrays containing highly specific LNA probes for over 600 human miRNAs - a substantially larger number than in most previous studies - 64 miRNAs were found to be misexpressed in AML in a statistically significant manner. Interestingly, subsequent studies on miR-221, one of the miRNAs most consistently overexpressed in AML, revealed a limited capacity of malignant myeloid cells to process its precursor forms. Firstly, only precursor, but not mature miR-221 levels were increased after introduction of appropriate expression vectors into several different human hematopoietic cell lines, and secondly, primary AML cells overexpressed pri-miR-221/222 to a much higher extent than its mature products. Because processing of vector-borne miRNA precursors, like that of endogenous pri-miRNAs, requires consecutive cleavage by Drosha and Dicer, the most plausible explanation of these results is a limited pri-miR-221/222 processing capacity at the Drosha level. pri-miR-221/222 overexpression was present in a high proportion of the patients investigated in this study (84%; 31/37 patients in the experiment shown in Figure 3, using the mean plus 3 standard deviations of the respective healthy controls as a cutoff for overexpression), so that the conclusion appears justified that it is a common phenomenon in AML. Underscoring the association between pri-miR-221/222 overexpression and AML, analysis of samples from the times of diagnosis, remission, and relapse showed that pri-miR-221/222 levels faithfully reflected the stage of disease.
As explanations for the elevated pri-miR-221/222 levels present in AML, pathologically increased transcription on the background of a physiologically limited processing capacity, or an acquired processing deficiency on the background of unaltered transcriptional activity may be considered. In disfavour of the latter possibility, a processing defect at the Drosha level does not necessarily lead to an accumulation of pri-miR molecules [52–56]. Also, it could not easily be reconciled with our own and published [12, 24, 25, 44] observations that mature miR-221 and 222 are not down-, but upregulated in AML. Increased transcription of pri-miR-221/222, in concert with a limited capacity of hematopoietic cells to process this molecule, therefore likely accounts for its induction in AML. This has not been proven, however, and will be the subject of future investigations.
A role for restricted miRNA processing capacities in tumor cells as well as in developmentally immature cells has been previously suggested by Thomson et al. . These authors found that mature let-7 g was almost undetectable in embryonic stem (ES) cells and P19 teratocarinoma cells, although the corresponding pri-miRNA was as abundant as in E14.5 embryos and differentiated P19 cells, which expressed high levels of let-7 g. Similarly, decreased expression of 22 miRNAs in a panel of human tumors of various origins was not accompanied by a reduction in the corresponding pri-miRNAs, indicating that it was due to an acquired processing defect rather than to transcriptional downregulation . In contrast, our data show that in AML transcriptional upregulation on the background of a probably a priori limited processing capacity causes an elevation of pri-miR-221/222 levels. Despite of these differences, both studies concur to suggest a role for restricted miRNA processing in tumorigenesis.
Since pri-miR-221/222 is upregulated in a large majority of AML, either this overexpression per se or an event leading to it should be essential for leukemogenesis. In principle, elevation of pri-miR-221/222 levels could be a bystander effect of the induction of a nearby oncogene. However, according to human genome assembly 19 (http://genome.ucsc.edu/) no known genes are located within a genomic region that extends from 0.7 Mb centromeric to 0.5 Mb telomeric of the mature miR-221 sequence. Also, co-activation of pri-miR-221/222 with an adjacent oncogene would not be able to explain why other pri-miRNAs are also overexpressed in AML. Alternatively, a transcription factor whose deregulation contributes to misexpression of cancer genes could also activate pri-miR-221/222 and other pri-miRNAs in trans. The possible identity of such a transcription factor is presently obscure. Thirdly, an acquired pri-miRNA processing deficiency might contribute to leukemogenesis by diminishing the production of mature miRNAs. This appears unlikely for the reasons outlined above. Fourth, inhibition of pri-miRNA processing could be a side effect of mutations whose primary oncogenic effect is unrelated to miRNA biogenesis. For example, mutations inactivating the transcriptional competence of p53 also interfered with miRNA processing . However, this defect was not accompanied by increased pri-miRNA levels , and furthermore p53 mutations are infrequent in AML , so that such a mechanism remains speculative at present.
Thus, even though a number of possibilities can be conceived to explain pri-miR-221/222 overexpression as a bystander effect of other oncogenic events, at present none of these is sufficiently plausible to discredit the possibility that pri-miR-221/222 may itself act as an oncogene. lncRNAs have recently received great attention as possible cancer genes with diverse ways of action, e.g. the regulation of chromatin and transcription , and pri-miRNAs may represent a novel type of lncRNA. In a recent study, the primary transcripts of the let-7 family of miRNAs (pri-let-7) were reported to be able to regulate gene expression , thus ascribing a biochemical function to this type of molecule. Furthermore, some lncRNAs have been reported to serve as precursors for small RNA species in addition to their well described functions as macro ncRNAs [23, 60]. Of particular interest, H19, an imprinted RNA that is overexpressed in many tumors and displays oncogenic functions [61, 62], has been found to host a miRNA . It can therefore be considered as an example of a pri-miRNA that at the same time acts as an oncogenic lncRNA.
In this study we identified 64 miRNAs that were differentially expressed between AML and healthy controls, and may therefore contribute to leukemogenesis. Moreover, we found that pri-miR-221/222 was strongly overexpressed, but inefficiently processed, in a large proportion of AML, indicating that it may represent a novel molecular marker and a putative oncogene in this disease.
Acute myeloid leukemia
Human genome assembly
- FAB classification:
French American British classification
Fetal bovine serum
False discovery rate
Global run-on sequencing
Hematopoietic stem and progenitor cell
Locked nucleic acid
Long noncoding RNA
Normal bone marrow
Normal peripheral blood
Quantitative reverse transcriptase polymerase chain reaction
RNA induced silencing complex
Standard error of the mean
Transcriptional start site
University of California Santa Cruz
White blood cell count.
This work was funded by the Austrian National Bank (OeNB), grant no 12505, the Austrian Science Foundation (FWF), grants no P19795 and P20920 to R.W., and a Clinical Research Grant of the Austrian Society of Hematology and Oncology (OeGHO). None of these funding bodies had any role in the design of the study, the collection, analysis, and interpretation of data, the writing of the manuscript, or the decision to submit it for publication.
miR-Vec-221 and empty miR-Vec were kindly provided by the Agami lab . Michael Wagner and Holger Daims of the University of Vienna kindly allowed access to their microarray scanner. Karin Lind and Trang Le are gratefully acknowledged for administrating patient samples and data. Finally, we thank Denise Barlow (CeMM, Vienna, Austria) for the unpublished RNA-seq data.
- Mrozek K, Marcucci G, Paschka P, Whitman S, Bloomfield C: Clinical relevance of mutations and gene-expression changes in adult acute myeloid leukemia with normal cytogenetics: are we ready for a prognostically prioritized molecular classification?. Blood. 2007, 109 (2): 431-448. 10.1182/blood-2006-06-001149.View ArticlePubMedPubMed CentralGoogle Scholar
- Lowenberg B: Acute myeloid leukemia: the challenge of capturing disease variety. Hematology Am Soc Hematol Educ Program. 2008, 1-11.Google Scholar
- Grimwade D, Walker H, Oliver F, Wheatley K, Harrison C, Harrison G, Rees J, Hann I, Stevens R, Burnett A, et al: The importance of diagnostic cytogenetics on outcome in AML: analysis of 1,612 patients entered into the MRC AML 10 trial. The Medical Research Council Adult and Children’s leukaemia working parties. Blood. 1998, 92 (7): 2322-2333.PubMedGoogle Scholar
- Schoch C, Kern W, Schnittger S, Buchner T, Hiddemann W, Haferlach T: The influence of age on prognosis of de novo acute myeloid leukemia differs according to cytogenetic subgroups. Haematologica. 2004, 89 (9): 1082-1090.PubMedGoogle Scholar
- Deneberg S, Guardiola P, Lennartsson A, Qu Y, Gaidzik V, Blanchet O, Karimi M, Bengtzen S, Nahi H, Uggla B, et al: Prognostic DNA methylation patterns in cytogenetically normal acute myeloid leukemia are predefined by stem cell chromatin marks. Blood. 2011, 118 (20): 5573-5582. 10.1182/blood-2011-01-332353.View ArticlePubMedGoogle Scholar
- Sill H, Olipitz W, Zebisch A, Schulz E, Wolfler A: Therapy-related myeloid neoplasms: pathobiology and clinical characteristics. Br J Pharmacol. 2011, 162 (4): 792-805. 10.1111/j.1476-5381.2010.01100.x.View ArticlePubMedPubMed CentralGoogle Scholar
- Dohner H, Estey E, Amadori S, Appelbaum F, Buchner T, Burnett A, Dombret H, Fenaux P, Grimwade D, Larson R, et al: Diagnosis and management of acute myeloid leukemia in adults: recommendations from an international expert panel, on behalf of the European LeukemiaNet. Blood. 2010, 115 (3): 453-474. 10.1182/blood-2009-07-235358.View ArticlePubMedGoogle Scholar
- Kern W, Haferlach C, Haferlach T, Schnittger S: Monitoring of minimal residual disease in acute myeloid leukemia. Cancer. 2008, 112 (1): 4-16. 10.1002/cncr.23128.View ArticlePubMedGoogle Scholar
- Marcucci G, Radmacher M, Maharry K, Mrozek K, Ruppert A, Paschka P, Vukosavljevic T, Whitman S, Baldus C, Langer C, et al: MicroRNA expression in cytogenetically normal acute myeloid leukemia. N Engl J Med. 2008, 358 (18): 1919-1928. 10.1056/NEJMoa074256.View ArticlePubMedGoogle Scholar
- Garzon R, Volinia S, Liu C, Fernandez-Cymering C, Palumbo T, Pichiorri F, Fabbri M, Coombes K, Alder H, Nakamura T, et al: MicroRNA signatures associated with cytogenetics and prognosis in acute myeloid leukemia. Blood. 2008, 111 (6): 3183-3189. 10.1182/blood-2007-07-098749.View ArticlePubMedPubMed CentralGoogle Scholar
- Jongen-Lavrencic M, Sun S, Dijkstra M, Valk P, Lowenberg B: MicroRNA expression profiling in relation to the genetic heterogeneity of acute myeloid leukemia. Blood. 2008, 111 (10): 5078-5085. 10.1182/blood-2008-01-133355.View ArticlePubMedGoogle Scholar
- Dixon-McIver A, East P, Mein C, Cazier J, Molloy G, Chaplin T, Andrew Lister T, Young B, Debernardi S: Distinctive patterns of microRNA expression associated with karyotype in acute myeloid leukaemia. PLoS One. 2008, 3 (5): e2141-10.1371/journal.pone.0002141.View ArticlePubMedPubMed CentralGoogle Scholar
- Wieser R, Scheideler M, Hackl H, Engelmann M, Schneckenleithner C, Hiden K, Papak C, Trajanoski Z, Sill H, Fonatsch C: microRNAs in acute myeloid leukemia: expression patterns, correlations with genetic and clinical parameters, and prognostic significance. Gene Chromosome Canc. 2010, 49 (3): 193-203.Google Scholar
- Bryant A, Lutherborrow M, Ma D: The clinicopathological relevance of microRNA in normal and malignant haematopoiesis. Pathology. 2009, 41 (3): 204-213. 10.1080/00313020902756287.View ArticlePubMedGoogle Scholar
- Pelosi E, Labbaye C, Testa U: MicroRNAs in normal and malignant myelopoiesis. Leuk Res. 2009, 33: 1584-1593. 10.1016/j.leukres.2009.04.039.View ArticlePubMedGoogle Scholar
- Pillai R, Bhattacharyya S, Filipowicz W: Repression of protein synthesis by miRNAs: how many mechanisms?. Trends Cell Biol. 2007, 17 (3): 118-126. 10.1016/j.tcb.2006.12.007.View ArticlePubMedGoogle Scholar
- Breving K, Esquela-Kerscher A: The complexities of microRNA regulation: mirandering around the rules. Int J Biochem Cell Biol. 2010, 42 (8): 1316-1329. 10.1016/j.biocel.2009.09.016.View ArticlePubMedGoogle Scholar
- Kim V, Nam J: Genomics of microRNA. Trends Genet. 2006, 22 (3): 165-173. 10.1016/j.tig.2006.01.003.View ArticlePubMedGoogle Scholar
- Newman M, Hammond S: Emerging paradigms of regulated microRNA processing. Genes Dev. 2010, 24 (11): 1086-1092. 10.1101/gad.1919710.View ArticlePubMedPubMed CentralGoogle Scholar
- Winter J, Jung S, Keller S, Gregory R, Diederichs S: Many roads to maturity: microRNA biogenesis pathways and their regulation. Nat Cell Biol. 2009, 11 (3): 228-234. 10.1038/ncb0309-228.View ArticlePubMedGoogle Scholar
- Krol J, Loedige I, Filipowicz W: The widespread regulation of microRNA biogenesis, function and decay. Nat Rev Genet. 2010, 11 (9): 597-610.PubMedGoogle Scholar
- Trujillo R, Yue S, Tang Y, O’Gorman W, Chen C: The potential functions of primary microRNAs in target recognition and repression. EMBO J. 2010, 29 (19): 3272-3285. 10.1038/emboj.2010.208.View ArticlePubMedPubMed CentralGoogle Scholar
- Koerner M, Pauler F, Huang R, Barlow D: The function of non-coding RNAs in genomic imprinting. Development. 2009, 136 (11): 1771-1783. 10.1242/dev.030403.View ArticlePubMedPubMed CentralGoogle Scholar
- Li Z, Lu J, Sun M, Mi S, Zhang H, Luo R, Chen P, Wang Y, Yan M, Qian Z, et al: Distinct microRNA expression profiles in acute myeloid leukemia with common translocations. Proc Natl Acad Sci USA. 2008, 105 (40): 15535-15540. 10.1073/pnas.0808266105.View ArticlePubMedPubMed CentralGoogle Scholar
- Isken F, Steffen B, Merk S, Dugas M, Markus B, Tidow N, Zuhlsdorf M, Illmer T, Thiede C, Berdel W, et al: Identification of acute myeloid leukaemia associated microRNA expression patterns. Br J Haematol. 2008, 140 (2): 153-161.View ArticlePubMedGoogle Scholar
- Merkel O, Hamacher F, Laimer D, Sifft E, Trajanoski Z, Scheideler M, Egger G, Hassler M, Thallinger C, Schmatz A, et al: Identification of differential and functionally active miRNAs in both anaplastic lymphoma kinase (ALK) + and ALK- anaplastic large-cell lymphoma. Proc Natl Acad Sci USA. 2010, 107 (37): 16228-16233. 10.1073/pnas.1009719107.View ArticlePubMedPubMed CentralGoogle Scholar
- Hackl M, Brunner S, Fortschegger K, Schreiner C, Micutkova L, Muck C, Laschober G, Lepperdinger G, Sampson N, Berger P, et al: miR-17, miR-19b, miR-20a, and miR-106a are down-regulated in human aging. Aging Cell. 2010, 9 (2): 291-296. 10.1111/j.1474-9726.2010.00549.x.View ArticlePubMedPubMed CentralGoogle Scholar
- Pieler R, Sanchez-Cabo F, Hackl H, Thallinger G, Trajanoski Z: ArrayNorm: comprehensive normalization and analysis of microarray data. Bioinformatics. 2004, 20 (12): 1971-1973. 10.1093/bioinformatics/bth174.View ArticlePubMedGoogle Scholar
- Koeffler H, Golde D: Acute myelogenous leukemia: a human cell line responsive to colony-stimulating activity. Science. 1978, 200 (4346): 1153-1154. 10.1126/science.306682.View ArticlePubMedGoogle Scholar
- Hamaguchi H, Suzukawa K, Nagata K, Yamamoto K, Yagasaki F, Morishita K: Establishment of a novel human myeloid leukemia cell line (HNT-34) with t(3;3)(q21;q26), t(9;22)(q34;q11) and the expression of EVI1 gene, p210 and p190 BCR/ABL chimaeric transcripts from a patient with AML after MDS with 3q21q26 syndrome. Br J Haematol. 1997, 98: 399-407. 10.1046/j.1365-2141.1997.2143029.x.View ArticlePubMedGoogle Scholar
- Paul C, Aly E, Lehman J, Page S, Gomez-Cambronero J, Ackerman S, Baumann M: Human cell line that differentiates to all myeloid lineages and expresses neutrophil secondary granule genes. Exp Hematol. 2000, 28 (12): 1373-1380. 10.1016/S0301-472X(00)00552-X.View ArticlePubMedGoogle Scholar
- Collins S, Gallo R, Gallagher R: Continuous growth and differentiation of human myeloid leukaemic cells in suspension culture. Nature. 1977, 270 (5635): 347-349. 10.1038/270347a0.View ArticlePubMedGoogle Scholar
- Martin P, Papayannopoulou T: HEL cells: a new human erythroleukemia cell line with spontaneous and induced globin expression. Science. 1982, 216 (4551): 1233-1235. 10.1126/science.6177045.View ArticlePubMedGoogle Scholar
- Sundstrom C, Nilsson K: Establishment and characterization of a human histiocytic lymphoma cell line (U-937). Int J Cancer. 1976, 17 (5): 565-577. 10.1002/ijc.2910170504.View ArticlePubMedGoogle Scholar
- Scherer W, Syverton J, Gey G: Studies on the propagation in vitro of poliomyelitis viruses. IV. Viral multiplication in a stable strain of human malignant epithelial cells (strain HeLa) derived from an epidermoid carcinoma of the cervix. J Exp Med. 1953, 97 (5): 695-710. 10.1084/jem.97.5.695.View ArticlePubMedPubMed CentralGoogle Scholar
- Soule H, Vazguez J, Long A, Albert S, Brennan M: A human cell line from a pleural effusion derived from a breast carcinoma. J Natl Cancer Inst. 1973, 51 (5): 1409-1416.PubMedGoogle Scholar
- Graham F, Smiley J, Russell W, Nairn R: Characteristics of a human cell line transformed by DNA from human adenovirus type 5. J Gen Virol. 1977, 36 (1): 59-74. 10.1099/0022-1317-36-1-59.View ArticlePubMedGoogle Scholar
- Le Sage C, Nagel R, Egan D, Schrier M, Mesman E, Mangiola A, Anile C, Maira G, Mercatelli N, Ciafre S, et al: Regulation of the p27(Kip1) tumor suppressor by miR-221 and miR-222 promotes cancer cell proliferation. EMBO J. 2007, 26 (15): 3699-3708. 10.1038/sj.emboj.7601790.View ArticlePubMedPubMed CentralGoogle Scholar
- Platzer B, Jorgl A, Taschner S, Hocher B, Strobl H: RelB regulates human dendritic cell subset development by promoting monocyte intermediates. Blood. 2004, 104 (12): 3655-3663. 10.1182/blood-2004-02-0412.View ArticlePubMedGoogle Scholar
- Livak K, Schmittgen T: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.View ArticlePubMedGoogle Scholar
- Huang R, Jaritz M, Guenzl P, Vlatkovic I, Sommer A, Tamir I, Marks H, Klampfl T, Kralovics R, Stunnenberg H, et al: An RNA-Seq strategy to detect the complete coding and non-coding transcriptome including full-length imprinted macro ncRNAs. PLoS One. 2011, 6 (11): e27288-10.1371/journal.pone.0027288.View ArticlePubMedPubMed CentralGoogle Scholar
- Core L, Waterfall J, Lis J: Nascent RNA sequencing reveals widespread pausing and divergent initiation at human promoters. Science. 2008, 322 (5909): 1845-1848. 10.1126/science.1162228.View ArticlePubMedPubMed CentralGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc B. 1995, 57 (1): 289-300.Google Scholar
- Cammarata G, Augugliaro L, Salemi D, Agueli C, La Rosa M, Dagnino L, Civiletto G, Messana F, Marfia A, Bica M, et al: Differential expression of specific microRNA and their targets in acute myeloid leukemia. Am J Hematol. 2010, 85 (5): 331-339.PubMedGoogle Scholar
- Ciafre S, Galardi S, Mangiola A, Ferracin M, Liu C, Sabatino G, Negrini M, Maira G, Croce C, Farace M: Extensive modulation of a set of microRNAs in primary glioblastoma. Biochem Biophys Res Commun. 2005, 334 (4): 1351-1358. 10.1016/j.bbrc.2005.07.030.View ArticlePubMedGoogle Scholar
- Mercatelli N, Coppola V, Bonci D, Miele F, Costantini A, Guadagnoli M, Bonanno E, Muto G, Frajese G, De Maria R, et al: The inhibition of the highly expressed miR-221 and miR-222 impairs the growth of prostate carcinoma xenografts in mice. PLoS One. 2008, 3 (12): e4029-10.1371/journal.pone.0004029.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee E, Gusev Y, Jiang J, Nuovo G, Lerner M, Frankel W, Morgan D, Postier R, Brackett D, Schmittgen T: Expression profiling identifies microRNA signature in pancreatic cancer. Int J Cancer. 2007, 120 (5): 1046-1054.View ArticlePubMedPubMed CentralGoogle Scholar
- Fornari F, Gramantieri L, Ferracin M, Veronese A, Sabbioni S, Calin G, Grazi G, Giovannini C, Croce C, Bolondi L, et al: MiR-221 controls CDKN1C/p57 and CDKN1B/p27 expression in human hepatocellular carcinoma. Oncogene. 2008, 27 (43): 5651-5661. 10.1038/onc.2008.178.View ArticlePubMedGoogle Scholar
- Miller T, Ghoshal K, Ramaswamy B, Roy S, Datta J, Shapiro C, Jacob S, Majumder S: MicroRNA-221/222 confers tamoxifen resistance in breast cancer by targeting p27Kip1. J Biol Chem. 2008, 283 (44): 29897-29903. 10.1074/jbc.M804612200.View ArticlePubMedPubMed CentralGoogle Scholar
- Di Leva G, Gasparini P, Piovan C, Ngankeu A, Garofalo M, Taccioli C, Iorio M, Li M, Volinia S, Alder H, et al: MicroRNA cluster 221–222 and estrogen receptor alpha interactions in breast cancer. J Natl Cancer Inst. 2010, 102 (10): 706-721. 10.1093/jnci/djq102.View ArticlePubMedPubMed CentralGoogle Scholar
- Bernstein B, Kamal M, Lindblad-Toh K, Bekiranov S, Bailey D, Huebert D, McMahon S, Karlsson E, Kulbokas E, Gingeras T, et al: Genomic maps and comparative analysis of histone modifications in human and mouse. Cell. 2005, 120 (2): 169-181. 10.1016/j.cell.2005.01.001.View ArticlePubMedGoogle Scholar
- Noh H, Hong S, Dong Z, Pan Z, Jing Q, Huang S: Impaired MicroRNA processing facilitates breast cancer cell invasion by upregulating Urokinase-type plasminogen activator expression. Gene Canc. 2011, 2 (2): 140-150. 10.1177/1947601911408888.View ArticleGoogle Scholar
- Thomson J, Newman M, Parker J, Morin-Kensicki E, Wright T, Hammond S: Extensive post-transcriptional regulation of microRNAs and its implications for cancer. Genes Dev. 2006, 20 (16): 2202-2207. 10.1101/gad.1444406.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee Y, Ahn C, Han J, Choi H, Kim J, Yim J, Lee J, Provost P, Radmark O, Kim S, et al: The nuclear RNase III Drosha initiates microRNA processing. Nature. 2003, 425 (6956): 415-419. 10.1038/nature01957.View ArticlePubMedGoogle Scholar
- Landthaler M, Yalcin A, Tuschl T: The human DiGeorge syndrome critical region gene 8 and Its D. melanogaster homolog are required for miRNA biogenesis. Curr Biol. 2004, 14 (23): 2162-2167. 10.1016/j.cub.2004.11.001.View ArticlePubMedGoogle Scholar
- Nakamura T, Canaani E, Croce C: Oncogenic All1 fusion proteins target Drosha-mediated microRNA processing. Proc Natl Acad Sci USA. 2007, 104 (26): 10980-10985. 10.1073/pnas.0704559104.View ArticlePubMedPubMed CentralGoogle Scholar
- Suzuki H, Yamagata K, Sugimoto K, Iwamoto T, Kato S, Miyazono K: Modulation of microRNA processing by p53. Nature. 2009, 460 (7254): 529-533. 10.1038/nature08199.View ArticlePubMedGoogle Scholar
- Haferlach C, Dicker F, Herholz H, Schnittger S, Kern W, Haferlach T: Mutations of the TP53 gene in acute myeloid leukemia are strongly associated with a complex aberrant karyotype. Leukemia. 2008, 22 (8): 1539-1541. 10.1038/leu.2008.143.View ArticlePubMedGoogle Scholar
- Huarte M, Rinn J: Large non-coding RNAs: missing links in cancer?. Hum Mol Genet. 2010, 19 (R2): R152-R161. 10.1093/hmg/ddq353.View ArticlePubMedPubMed CentralGoogle Scholar
- Wilusz J, Sunwoo H, Spector D: Long noncoding RNAs: functional surprises from the RNA world. Genes Dev. 2009, 23 (13): 1494-1504. 10.1101/gad.1800909.View ArticlePubMedPubMed CentralGoogle Scholar
- Looijenga L, Verkerk A, De Groot N, Hochberg A, Oosterhuis J: H19 in normal development and neoplasia. Mol Reprod Dev. 1997, 46 (3): 419-439. 10.1002/(SICI)1098-2795(199703)46:3<419::AID-MRD22>3.0.CO;2-S.View ArticlePubMedGoogle Scholar
- Matouk I, DeGroot N, Mezan S, Ayesh S, Abu-lail R, Hochberg A, Galun E: The H19 non-coding RNA is essential for human tumor growth. PLoS One. 2007, 2 (9): e845-10.1371/journal.pone.0000845.View ArticlePubMedPubMed CentralGoogle Scholar
- Cai X, Cullen B: The imprinted H19 noncoding RNA is a primary microRNA precursor. RNA. 2007, 13 (3): 313-316. 10.1261/rna.351707.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/13/364/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.