- Open Access
Hypermethylation and down-regulation of DLEU2 in paediatric acute myeloid leukaemia independent of embedded tumour suppressor miR-15a/16-1
Molecular Cancer volume 13, Article number: 123 (2014)
Acute Myeloid Leukaemia (AML) is a highly heterogeneous disease. Studies in adult AML have identified epigenetic changes, specifically DNA methylation, associated with leukaemia subtype, age of onset and patient survival which highlights this heterogeneity. However, only limited DNA methylation studies have elucidated any associations in paediatric AML.
We interrogated DNA methylation on a cohort of paediatric AML FAB subtype M5 patients using the Illumina HumanMethylation450 (HM450) BeadChip, identifying a number of target genes with p <0.01 and Δβ >0.4 between leukaemic and matched remission (n = 20 primary leukaemic, n = 13 matched remission). Amongst those genes identified, we interrogate DLEU2 methylation using locus-specific SEQUENOM MassARRAY® EpiTYPER® and an increased validation cohort (n = 28 primary leukaemic, n = 14 matched remission, n = 17 additional non-leukaemic and cell lines). Following methylation analysis, expression studies were undertaken utilising the same patient samples for singleplex TaqMan gene and miRNA assays and relative expression comparisons.
We identified differential DNA methylation at the DLEU2 locus, encompassing the tumour suppressor microRNA miR-15a/16-1 cluster. A number of HM450 probes spanning the DLEU2/Alt1 Transcriptional Start Site showed increased levels of methylation in leukaemia (average over all probes >60%) compared to disease-free haematopoietic cells and patient remission samples (<24%) (p < 0.001). Interestingly, DLEU2 mRNA down-regulation in leukaemic patients (p < 0.05) was independent of the embedded mature miR-15a/16-1 expression. To assess prognostic significance of DLEU2 DNA methylation, we stratified paediatric AML patients by their methylation status. A subset of patients recorded methylation values for DLEU2 akin to non-leukaemic specimens, specifically patients with sole trisomy 8 and/or chromosome 11 abnormalities. These patients also showed similar miR-15a/16-1 expression to non-leukaemic samples, and potential improved disease prognosis.
The DLEU2 locus and embedded miRNA cluster miR-15a/16-1 is commonly deleted in adult cancers and shown to induce leukaemogenesis, however in paediatric AML we found the region to be transcriptionally repressed. In combination, our data highlights the utility of interrogating DNA methylation and microRNA in combination with underlying genetic status to provide novel insights into AML biology.
Acute myeloid leukaemia (AML) is the third most common form of leukaemia in children, typically characterised by the rapid proliferation of primitive haematopoietic myeloid progenitor cells. Paediatric AML is a highly heterogeneous disease, which presents a major barrier towards the development of accurate disease classification, risk stratification and targeted therapies within the clinic. The French-American-British (FAB) and more recently World Health Organisation (WHO) classifications of leukaemia take into account cell morphology, cytogenetic aberrations and common genetic lesions. However, not all patients fall into these well-defined categories. Additionally, the recurrent chromosomal and genetic lesions frequently found in AML fail to induce leukaemogenesis and do not explain the recognised clinical heterogeneity[4, 5].
One of the hallmarks of nearly all human cancers is the disruption of the epigenetic profile, including gross aberrations in DNA methylation. Increasing evidence in adult AML has indicated that epigenetic events play critical roles in the onset, progression, and outcome of AML and may help tailor disease treatment. However, the need for similar elucidations in childhood disease is paramount. Aberrant methylation of cytosine residues at palindromic CpG sites (often clustered in dense CpG ‘islands’) near gene promoter regions is widely studied in carcinogenesis and haematological malignancies[6, 7]. It is now well established that elevated DNA methylation is an important mechanism of gene transcriptional inactivation[8, 9] and genes such as ESR1, IGSF4 and CDKN2B/p15 are epigenetically silenced in adult leukaemia. Previous studies have subdivided adult AML into 16 epigenetic sub-groups based on DNA methylation signatures, correlating with patient clinical outcome and distinct from both normal haematopoietic cells and normal stages of myeloid differentiation. Despite such emerging findings in an adult context, the utility of individual DNA methylation disruptions in paediatric AML has yet to be fully evaluated.
MicroRNA (miRNA) represent an alternative epigenetic regulator, having been implicated in the regulation of critical gene expression networks in plants and animals. The role of miRNA in haematopoiesis, cancer and disease is also beginning to be appreciated[10, 11]. The global influence of individual miRNA on the genome is difficult to dissect, as miRNA can modulate the expression of hundreds of genes, and each gene can harbour binding sites for several miRNA. Human miRNA are initially transcribed (pri-miRNA), and processed by several complexes to form a 70 bp hairpin-loop (pre-miRNA). After successive enzymatic steps, a miRNA:miRNA* complementary duplex is formed where the ‘functional’ strand is combined with RISC (RNA Induced Silencing Complex) and Argonaute proteins to guide, and inhibit, specific target messenger RNA (mRNA) through base pair recognition[14, 15]. However, the miRNA transcriptome is becoming increasingly complex, emphasised by Next Generation Sequencing (NGS) technologies. NGS has highlighted that alternate miRNA* transcripts, as well as miRNA sequence variants (isomiRs) may play a biological role, similar to their canonical miRNA relatives[17, 18].
Links between miRNA deregulation and cancer diagnosis were first identified in adult Chronic Lymphoblastic Leukaemia (CLL), where the loss or down-regulation of tumour-suppressing miRNA cluster miR-15a/16-1 directly caused leukaemic transformation[19, 20]. At present, no such association has been identified for childhood leukaemia. The expression of paediatric disease-associated miRNA has to date only identified a distinction between leukaemia of different lineages and the differentiation of rearranged AMLs within a limited number of cytogenetic subtypes[21, 22]. Paediatric MLL can be distinguished from others by differentially expressed miR-126, miR-146a, miR-181a/b/d, miR-100, miR-21, miR-196a/b, miR-29 and miR-125b. However concordance among studies is often low and the mechanism of deregulation is often unknown[22, 23].
Genes encoding miRNA can be regulated epigenetically in a similar manner to protein coding genes. Studies have demonstrated epigenetically regulated miRNA in adult AML, including hypermethylation and down-regulation of miR-124a and associated deregulation of target mRNA EVI1, CEBPA and CDK6 independent of diagnostic cytogenetic subtype (reviewed in). Additionally, miR-193a targeting KIT, and miR-14b targeting CREB have been identified in adult investigations as specifically controlled by DNA methylation (reviewed in). However, the identification of DNA methylation and miRNA expression connections in paediatric leukaemia is lacking.
Paediatric AML has distinct cytogenetic and clinical features relative to their adult counterparts[5, 21, 24–26]. Therefore, there is a critical need to improve our understanding of the biology of childhood leukaemia as separate entities, distinct from adult disease. Cognisant of this, we aimed to identify differential DNA methylation within paediatric AML on a genome-scale using defined clinical subtypes and age-matched controls. We identified a number of significantly altered DNA methylation loci, with associated gene and miRNA expression change, between paediatric AML and non-leukaemic counterparts. Specifically we describe here the epigenetic deregulation of DLEU2, which has associated alterations in downstream miR-15a/16-1 miRNA cluster expression.
Results and discussion
The DLEU2 gene is specifically hypermethylated and repressed in paediatric AML subtype M5
The FAB subtype M5 (monocytic/blastic leukaemia) is a distinct subtype with characteristic chromosomal abnormalities including t(8; 16), +8 and various translocations involving 11q23 and the MLL locus such as t(9;11), t(10;11), t(11;19) and others. AML subtype M5 also has a high proportion of cytogenetically normal (CN-AML) patients, and those with complex karyotypes. A combination of these factors add to the overall unfavourable outcome of paediatric M5 diagnosis. Genome-scale methylation profiling of AML M5 bone marrow samples identified 3,352 significantly differentially methylated probes (DMPs) between paediatric AML FAB M5 (n = 20) and matching non-leukaemic (n = 17) samples. Applying more stringent feature selection criteria of an adjusted p-value <0.01 and ∆β of >0.4 reduced the number of DMPs to 137 (Additional file1).
The list of DMPs included several localising to the long non-coding RNA DLEU2 and the embedded miRNA cluster miR-15a/16-1[19, 20], previously implicated in adult leukaemic[19, 20, 30]. To date, disruption of this region has not been observed in paediatric cancers and as such we chose to focus on DLEU2 in subsequent analysis. A total of three DLEU2 DMPs had an adjusted p-value <0.01 and ∆β of >0.4 (cg05394800, cg20529344, cg12883980). These probes were located within three CpG islands (Chr 13: 50, 690,000-50,708,000: UCSC human hg19 assembly) at the DLEU2/Alt1 transcriptional start site (TSS) and ‘north shore’ (Figure 1), a region up to 2 kb upstream from the DLEU2/Alt1 TSS CpG island under investigation. Henceforth this will be referred to as the DLEU2 promoter.
A divergent DNA methylation profile was observed between all cases of paediatric AML with elevated methylation (mean 64% (48-80% CI) across significant probes and non-leukaemic samples (23% methylation; 7-39% CI of mean) (p < 0.001 for all significant probes) (Figure 1; Additional file2). DNA methylation levels at CpG sites interrogated by cg05394800, cg20529344, cg12883980 were confirmed using SEQUENOM EpiTYPER MassArray (Additional file3 and Additional file4) in the discovery sample set of 20 leukaemic patients and 17 non-leukaemic controls, as well as a validation set of paediatric AML M5 patients with heterogeneous cytogenetic diagnoses (n = 19) (CN-AML (n = 5), MLL (n = 7), RUNX1 (n = 2), WT1 (n = 1)) (Additional file5). This analysis confirmed specific DLEU2 promoter hypermethylation in paediatric AML patients (p < 0.001) (Additional file4), with additional genes in this region, including DLEU1, hsa-miR-15a/16-1 microRNA cluster, TRIM13, KCNRG and hsa-miR-3613 showing no differences in methylation between cases and controls (Additional file6). This highlights the potential regional specificity of the observed DLEU2 methylation change in association with AML.
We observed a significant down-regulation of DLEU2 gene expression in paediatric AML (0.07 Fold Change (FC); p = 0.014), and a significant inverse correlation between promoter DNA methylation and gene expression levels (p = 0.0001, Additional file7). Recent studies in adult CLL have also identified a negative correlation between DLEU2 promoter methylation (DLEU2/Alt1) and gene expression. Interestingly there was no change in gene expression for any other genes in this region (TRIM13 = 1.75 FC; DLEU1 = 1.05 FC. Figure 2A).
DLEU2 and embedded miR-15a/16-1 are regulated independently in paediatric AML
The miR-15a/16-1 cluster has been described to have potent tumour suppressor activity, targeting numerous oncogenic and cell cycle regulatory genes[19, 33]. The cluster is embedded within intron 4 of DLEU2, and has been speculated that expression is driven by the DLEU2 promoter[20, 30]. However, we found no correlation between DLEU2 expression and miR-15a/16-1 expression in paediatric AML, nor down-regulation of the miR-15a/16-1 miRNA cluster in relation to increasing DLEU2 promoter DNA methylation (Figure 2B). In contrast to previous reports for adult leukaemia, no significant change in mature miR-16 expression was observed between paediatric AML and control samples, a result we have reported elsewhere. These observations are independent of the homologous miR-16 cluster embedded in SMC4 on chromosome 3q26, which shows no significant expression or DNA methylation changes in association with AML (Additional file8). Similar results have been reported for adult CLL. Taken together this data suggests an alternate mode of regulation for the miR-15a/16-1 cluster, outside of the DLEU2 promoter region.
Interestingly, recent research has indicated the processing mechanisms of miRNA may be affected by cancer, such that mature miRNA expression becomes disassociated from precursor miRNA levels, and also from the levels of the host gene. We found the miR-15a/16-1 primary precursor transcript is in fact expressed up to three-fold higher in AML patients compared to non-leukaemic counterparts (Figure 2B). The increase in primary transcript appears to correspond to an increase in mature miR-16-1* (2.52 FC) and miR-15a* (2.2 FC) expression, and a moderate increase in miR-15a (1.5 FC). However, individual patient samples show differential degrees of over-expression of one, or all, of the mature species despite the common precursor (Additional file9).
Canonical miR-15a/16-1 microRNA species target numerous oncogenic and integral cell cycle regulatory genes such as BCL-2, MCL-1, CCND1, CDK6, BMI-1, RASSF5, IGSF4, c-MYB and WNT3A[33, 36–42]. To identify the potential significance of increasing alternate miR-15a/16-1 transcripts over their canonical counterparts in paediatric AML, we investigated miRNA target genes using prediction tools, and undertook gene ontology analysis. We found miR-16-1* is three-fold enriched for targeting components of the RNA processing and splicing machinery such as hnRNP regulators, SRRM1 and NUDT2. Down-regulation of these genes by increased expression of a targeting miRNA may be contributing to miRNA/host gene disassociation in DLEU2/miRNA-15a/16-1 expression (Additional file10). One target of miR-15a* was found to be ADORA2A, an important G-protein receptor critical in tissue-specific and systemic inflammatory responses[43, 44]. We have identified the down-regulation of ADORA2A expression within paediatric AML, not mediated through DNA methylation, but potentially through the up-regulation of the miR-15a/16-1 cluster (data not shown). Moreover, miR-15a* and miR-16-1* are 60–90 fold enriched for targeting genes involved in the intrinsic apoptotic pathway such as BCL-2 L11, PPIF and DNM1L (Additional file10). The inhibition of these genes to act upon the cell, mediated through targeting miRNA, can potentially encourage continued cell growth and perpetuation of AML, through negating apoptotic signalling. The over expression of alternate mature miRNA in paediatric AML may also have consequences for integral cell cycle pathways, but through different mechanisms to their canonical counterparts.
DLEU2 interrogation identifies a novel subclass of paediatric AML
In addition to assessing the association of DLEU2 promoter DNA methylation and miR-15a/16-1 miRNA expression with paediatric AML, we also assessed correlations with distinct clinical and diagnostic variables. We found no significant association between DLEU2 methylation/expression with sex, age of disease onset or relapse status nor common gene abnormalities such as FLT3 or MLL, as has been found in a number of other association studies[6, 45]. DLEU2 methylation clustering based on diagnostic FAB subtype revealed a wide range of values (M5a: 16-91%; M5b: 48-90%; ‘Other’ subtypes: 51-93%. Additional file11), indicating mean methylation in isolation does not stratify paediatric AML subtype. Additionally, we found there to be no correlation between percentage of leukaemic blasts at patient diagnosis and the DNA methylation status at the DLEU2 promoter region (data not shown).
Previous adult leukaemia studies have identified novel methylation subgroups, independent of cytogenetic profile[4, 46]. We identified a subset of patients displaying low DLEU2 promoter methylation comprising patients with trisomy 8 and/or chromosome 11q abnormalities, independent of MLL rearrangements (separated into a t(11)/+8 subgroup - Figure 3A, Additional file11). This delineation identified t(11)/+8 subgroup as having comparable DLEU2 methylation levels to non-leukaemic specimens (Figure 3B), and significantly lower methylation compared to the remaining paediatric AML samples (M5a/b p < 0.05; M1/M2/M4 p < 0.0001 Figure 3B). Classifying t(11)/+8 patients as a separate subgroup further delineated M5a, M5b and M1/M2/M4 AML patients from non-leukaemic counterparts according to average methylation status (M5a MM = 0.81, M5b MM = 0.78, M1/M2/M4 MM = 0.87. p < 0.0001. Figure 3B).
Cases carrying 11q rearrangements are the most heterogeneous of paediatric AML[27, 29, 47], linked to 50–104 translocation fusion partners to date[29, 46, 48]. Trisomy 8 is a frequently reported aberration in adult and paediatric AML. However little is known about the gain of chromosome 8 in isolation and its relationship to disease onset. As such this has been speculated to be a disease modulating secondary event. Lower DNA methylation in the t(11)/+8 subgroup may potentially confer a better prognosis, as it is well documented that hypermethylation of the DLEU2 region is involved in leukaemic transformation in adults[19, 20, 30]. Analysis of the t(11)/+8 subgroup in isolation also revealed no significant miRNA expression differences relative to non-leukaemic samples (Figure 4A), and it is well documented that miR-15a/16-1 abnormalities are also associated with leukaemic transformation in adults[19, 20]. We additionally identified a trend towards decreased risk of relapse in t(11)/+8 subgroup compared to other subtypes, and a trend towards better survival outcomes (Additional file12). These analyses combined may elucidate a connection between DLEU2 promoter DNA methylation (Figure 3B), miRNA expression (Figure 4A) and prognostic outcomes for this subgroup of paediatric AML (Additional file12).
The classification of t(11)/+8 cases as an independent subgroup revealed a significant increase in pri-miR-15a expression in traditional FAB subtype M5b from non-leukaemic (20.29-fold, p < 0.001; Figure 4C) and from M5a/M1/M2/M4 (Figure 4B). The increase translates to a >2-fold increase in miR-15a, 4.87-fold increase in miR-15a*, and 9.86-fold increase in miR-16-1* expression. Based on DLEU2 DNA methylation status, in conjunction with miRNA expression analysis, we speculate that pri-miR-15a expression alone may be a useful biomarker to distinguish M5b FAB subtype from all other AML subtypes.
Previous research within paediatric AML has shown that well-defined cytogenetic subgroupings exhibit a wide range of genomic and epigenomic heterogeneity. Linking specific epigenetic features to clinical parameters has the potential to identify pathological drivers of disease and develop enhanced molecular approaches for diagnosis, prognosis and refinement of treatment. Profiling epigenetic regulators has provided clinically relevant biomarkers for adult cancers; however generally the same cannot be said of childhood cancers. Our analysis has identified hypermethylation induced down-regulation of the DLEU2 gene in paediatric AML. The related expression changes of the embedded miR-15a/16-1 microRNA cluster in paediatric AML has the potential to contribute to leukaemic transformation, with a switch from canonical to alternate miRNA family members, which in turn may modulate the expression of downstream regulatory genes. Treating paediatric AML model systems with epigenetic modifying drugs will allow a more comprehensive analysis towards the clinical applications and potential patient-focussed therapeutic interventions for children. Our results highlight the need for further specific interrogation of paediatric AML subtypes as distinctive biological entities, separate from adult disease. Further studies utilising larger patient cohorts are required to explore the complex interplay between the epigenetic regulation of genes harbouring microRNA, taking into account alternate miRNA transcript expression.
This study was approved by the Royal Children’s Hospital (RCH), Melbourne, Ethics Committee (HREC reference #27138E). Samples used consisted of snap frozen bone marrow specimens and archived bone marrow biopsies taken at diagnosis, patient remission/follow-up or at relapse from paediatric acute myeloid leukaemia (AML) cases. The diagnosis of AML was established according to the criteria of the French–American–British (FAB) classification by standard morphological and cytological methods. The evaluation for mutations in MLL, FLT3, RUNX1 and WT1 genes were assessed in a small number of patients as a part of the initial assessment. No patient recorded any chromosome 13 deletions or abnormalities listed in their clinical diagnoses. The median percentage of leukaemic blasts at patient diagnosis is 88% (67–93 95% CI).
Patient samples used in our study consisted of archived, air-dried bone marrow smear slides. The utility of these samples for DNA methylation and miRNA expression analysis has been outlined previously[34, 50]. Cryogenically frozen patient bone marrows were also used where available. All patients were <18 years of age. We chose to focus specifically on the FAB subtype M5 (a and b) as it is a common subtype. Our cohort included comparable ratios between males and females, and also similar numbers of children that relapsed or not. (Additional file5). Control samples from bone marrow of unrelated and unaffected children were analysed in parallel, as well as multiple cell lines, including adult leukaemia (REH, CCRF-CEM), paediatric AML (Kasumi-1, THP-1, MV-4-11 and AML-193) obtained from American Type Culture Collection (ATCC) and subjected to characterisation using the ATCC Proficiency Standard program. Fluorescent Activated Cell Sorted (FACS) isolated haematopoietic progenitor cell populations from unrelated and unaffected paediatric donors (CD19+, CD33+, CD34+ and CD45+ populations, from herein known as CD sorted cells), were also used in this study. Our ‘non-leukaemic’ group consisted of all individual CD sorted cell populations as well as non-leukaemic remission and follow-up slides and bone marrow from paediatric patients. We compared the DLEU2/Alt1 DNA methylation for all ‘non-leukaemic’ specimens, and found no significant differences across these non-leukaemic samples (Additional file13).
DNA extraction, quality control and methylation analysis
Genomic DNA was extracted using the phenol/chloroform method and bisulphite converted in accordance with manufacturer’s protocols, as reported previously. DNA extracts were checked for quality and quantity using a NanoDrop® ND-1000 spectrophotometer (Thermo Fisher Scientific Inc., Scoresby, Victoria, Australia). All DNA was stored at -80°C.
Bisulphite conversion of genomic DNA was performed using the MethylEasy Xceed Bisulfite Modification Kit (Human Genetic Signatures, Sydney, AUST). The converted samples were processed by the Australian Genome Research Facility (AGRF, Melbourne, Australia) and analysed using the Illumina HumanMethylation450 (HM450) Bead Chip arrays according to manufacturer’s protocol. Illumina Genome Studio software was used to extract the raw M-values and probe intensities for downstream processing. Samples with a probe p-detection value of p <0.05 were retained. Data were quantile normalized using lumi and analysed using LIMMA. All analysis was performed using the R statistical software package on M-values that were converted to β-values (0 = unmethylated, 1 = fully methylated) for reporting and biological interpretation. To eliminate sex bias, probes hybridizing to sex chromosomes were filtered out, leaving 366,553 probes common to all samples in the final dataset. Further filtering was based on the degree of difference in β between sample groups, indicated as ∆β calculated using the formula:
SEQUENOM MassARRAY® EpiTYPER® was used to measure locus-specific methylation. Sixty-one samples were analysed, including the 42 analysed on the HM450 platform (Additional file5), an additional 3 cell lines, 8 leukaemic and 8 non-leukaemic whole bone marrow aspirate samples. Primers for analysis were designed using SEQUENOM EpiDesigner software (http://www.epidesigner.com) and the sequences for these are listed in Additional file14. Gene ontology and pathway analysis of genes associated with significantly altered DNA methylation probes were analysed through the use of IPA (Ingenuity® Systems, http://www.ingenuity.com) and GOrilla.
RNA extraction, quality control and expression analysis
Small RNA extraction was performed using TRIzol® (Ambion®) for patient bone marrow samples, mononuclear cells and cell lines in accordance with the manufacturer’s instructions, or the High Pure miRNA Isolation Kit (Roche) for archived slide samples as previously described. Before RNA extraction of fresh bone marrow aspirates, samples were processed using Ficoll-Paque™ (GE Healthcare, Piscataway USA) to isolate the mononuclear cell population. This was immediately cryo-frozen or stored in RNAlater® (Ambion® by Life Technologies, Mulgrave, Victoria Australia) for later extraction. The concentration and purity of all RNA samples was assessed using the NanoDrop® ND-1000 spectrophotometer (Thermo Fisher Scientific Inc., Scorsby, Victoria, Australia). All RNA was stored at -80°C.
Seventy samples were used in the interrogation of miRNA expression, 64 of which were also used for HM450 and SEQUENOM methylation analysis (Additional file5). TaqMan® MicroRNA Reverse Transcription kit and singleplex TaqMan ® microRNA Assays (Applied Biosystems, Life Technologies) (assays listed in Additional file15) were utilised according to the manufacturer’s instructions before routine quantitative real-time PCR (qRT-PCR) was performed using Applied Biosystems 7300 Sequence Detection System. A subset of 29 high quality samples which included all cell lines as previously described, CD19+, CD33+, CD34+, CD45+ as well as 8 leukaemic and 9 non-leukaemic patient whole bone marrows, were used for gene and pri-miR expression analysis. This was due to longer RNA species often being degraded in archived specimens[59, 60]. The SuperScript® VILO™ cDNA Synthesis Kit (Life Technologies) and TaqMan® gene expression assays were used for gene/pri-miR expression as per manufacturer’s instructions (Additional file15 and Additional file16) before performing qRT-PCR analysis according to the manufacturer’s instructions. All qRT-PCR samples were analysed in duplicate.
Expression data was normalized using the synthetic small nucleolar RNA, C/D Box 44 (RNU44) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) for gene/pri-miR expression; and RNU44 and hsa-miR-26b for miRNA analysis, previously identified as having the greatest stability in our samples. Data were analysed using DataAssist™ software (v.3.0, Applied Biosystems, Life Technologies: http://www.lifetechnologies.com/au/en/home/technical-resources/software-downloads/dataassist-software), and all p-values were adjusted using Benjamini-Hochberg False Discovery Rate.
All analyses used the combined average of non-leukaemic primary bone marrow tissue and CD sorted non-leukaemic cell expression values as the Reference Group (calibrator). Fold Change (FC) was calculated using the Livak method of 2-∆∆Ct plotting fold change and associated p-values. The fold change reported here is the difference of the means of each group, and a fold change of >2 between disease and non-leukaemic groups were considered noteworthy. microRNA gene target prediction was assessed using microRNA.org[62, 63] and miRWalk. The top 500 target genes for each miRNA (Additional file17) were used for gene ontology/pathway analysis (Additional file10) by IPA (Ingenuity® Systems, http://www.ingenuity.com) and GOrilla.
Where groups were compared with more than one variable measured (non-repeated measures), a Two-Way ANOVA with Bonferroni Post-testing was performed. Where groups were compared with only one variable measured, a Kruskal-Wallis test was performed with Dunns Post-testing. These analyses were undertaken on normalized data. Correlation assessments utilised Spearman Rank tests. Two-tailed p-values were utilised for assessment of statistical significance, and results with p-values p < 0.05 were considered to be statistically significant. Survival and relapse curves were estimated by the Kaplan-Meier method. Statistical analyses were performed with GraphPad PRISM (Version 5) for Windows, (GraphPad Software, USA, http://www.graphpad.com). CI refers to 95% Confidence Intervals.
Acute myeloid leukaemia
Cytogenetically normal acute myeloid leukaemia
Deleted in lymphocytic leukaemia 2
Differentially methylated probe
French-American-British leukaemia classification method
Illumina humanmethylation450 bead chip array
Transcriptional start site.
Faderl S, Talpaz M, Estrov Z, O’Brien S, Kurzrock R, Kantarjian HM: The biology of chronic myeloid leukemia. N Engl J Med. 1999, 341: 164-172. 10.1056/NEJM199907153410306.
Bennett JM, Catovsky D, Daniel M-T, Flandrin G, Galton DAG, Gralnick HR, Sultan C: Proposals for the classification of the acute leukaemias. Br J Haematol. 1976, 33: 451-458. 10.1111/j.1365-2141.1976.tb03563.x.
Vardiman JW, Thiele J, Arber DA, Brunning RD, Borowitz MJ, Porwit A, Harris NL, Le Beau MM, Hellstrom-Lindberg E, Tefferi A, Bloomfield CD: The 2008 revision of the World Health Organization (WHO) classification of myeloid neoplasms and acute leukemia: rationale and important changes. Blood. 2009, 114: 937-951. 10.1182/blood-2009-03-209262.
Figueroa ME, Lugthart S, Li Y, Erpelinck-Verschueren C, Deng X, Christos PJ, Schifano E, Booth J, van Putten W, Skrabanek L, Campagne F, Mazumdar M, Greally JM, Valk PJM, Lowenberg B, Delwel R, Melnick AM: DNA methylation signatures identify biologically distinct subtypes in acute myeloid leukemia. Cancer Cell. 2010, 17: 13-27. 10.1016/j.ccr.2009.11.020.
Pui C-H, Carroll WL, Meshinchi S, Arceci RJ: Biology, risk stratification, and therapy of pediatric acute leukemias: an update. J Clin Oncol. 2011, 29: 551-565. 10.1200/JCO.2010.30.7405.
Alvarez S, Suela J, Valencia A, Fernandez A, Wunderlich M, Agirre X, Prosper F, Martin-Subero JI, Maiques A, Acquadro F, Rodriguez Perales S, Calasanz MJ, Roman-Gomez J, Siebert R, Mulloy JC, Cervera J, Sanz MA, Esteller M, Cigudosa JC: DNA methylation profiles and their relationship with cytogenetic status in adult acute myeloid leukemia. PLoS One. 2010, 5: e12197-10.1371/journal.pone.0012197.
Gardiner-Garden M, Frommer M: CpG islands in vertebrate genomes. J Mol Biol. 1987, 196: 261-282. 10.1016/0022-2836(87)90689-9.
Galm O, Wilop S, Lüders C, Jost E, Gehbauer G, Herman JG, Osieka R: Clinical implications of aberrant DNA methylation patterns in acute myelogenous leukemia. Ann Hematol. 2005, 84: 39-46. 10.1007/s00277-005-0005-0.
Bird AP: DNA methylation patterns and epigenetic memory. Genes Dev. 2002, 16: 6-21. 10.1101/gad.947102.
Bissels U, Bosio A, Wagner W: MicroRNAs are shaping the hematopoietic landscape. Haematologica. 2012, 97: 160-167. 10.3324/haematol.2011.051730.
Calin GA, Croce CM: MicroRNA signatures in human cancers. Nat Rev Cancer. 2006, 6: 857-866. 10.1038/nrc1997.
Krek A, Grun D, Poy MN, Wolf R, Rosenberg L, Epstein EJ, MacMenamin P, da Piedade I, Gunsalus KC, Stoffel M, Rajewsky N: Comninatorial microRNA target predictions. Nat Genet. 2005, 37: 495-500. 10.1038/ng1536.
Calin GA: MicroRNAs and cancer: what we know and what we still have to learn. Genome Med. 2009, 1: 78-10.1186/gm78.
Tsai LM, Yu D: MicroRNAs in common diseases and potential therapeutic applications. Clin Exp Pharmacol Physiol. 2010, 37: 102-107. 10.1111/j.1440-1681.2009.05269.x.
Zhao G, Yu D, Weiss MJ: MicroRNAs in erythropoiesis. Curr Opin Hematol. 2010, 17: 155-162.
Marti E, Pantano L, Banez-Coronel M, Llorens F, Minones-Moyano E, Porta S, Sumoy L, Ferrer I, Estivill X: A myriad of miRNA variants in control and Huntington’s disease brain regions detected by massively parallel sequencing. Nucleic Acids Res. 2010, 38: 7219-7235. 10.1093/nar/gkq575.
Llorens F, Bañez-Coronel M, Pantano L, Antonio del Río J, Ferrer I, Estivill X, Martí E: A highly expressed miR-101 isomiR is a functional silencing small RNA. BMC Genomics. 2013, 14: 104-10.1186/1471-2164-14-104.
Neilsen CT, Goodall GJ, Bracken CP: IsomiRs – the overlooked repertoire in the dynamic microRNAome. Trends Genet. 2012, 28: 544-549. 10.1016/j.tig.2012.07.005.
Calin GA, Dumitru CD, Shimizu M, Bichi R, Zupo S, Noch E, Alder H, Rattan S, Keating MJ, Rai K, Rassenti L, Kipps TJ, Negrini M, Bullrich F, Croce CM: Frequent deletions and down-regulation of micro-RNA genes miR15 and miR16 at 13q14 in chronic lymphocytic leukemia. Proc Natl Acad Sci U S A. 2002, 99: 15524-15529. 10.1073/pnas.242606799.
Klein U, Lia M, Crespo M, Siegel R, Shen Q, Mo T, Ambesi-Impiombato A, Califano A, Migliazza A, Bhagat G, Dalla-Favera R: The DLEU2/miR-15a/16-1 cluster controls B cell proliferation and its deletion leads to chronic lymphocytic leukemia. Cancer Cell. 2010, 17: 28-40. 10.1016/j.ccr.2009.11.019.
Daschkey S, Rottgers S, Giri A, Bradtke J, Teigler-Schlegel A, Meister G, Borkhardt A, Landgraf P: MicroRNAs distinguish cytogenetic subgroups in pediatric AML and contribute to complex regulatory networks in AML-relevant pathways. PLoS One. 2013, 8: e56334-10.1371/journal.pone.0056334.
Agirre X, Martinez-Climent JA, Odero MD, Prosper F: Epigenetic regulation of miRNA genes in acute leukemia. Leukemia. 2012, 26: 395-403. 10.1038/leu.2011.344.
Chatterton Z, Morenos L, Saffery R, Craig JM, Ashley D, Wong NC: DNA methylation and miRNA expression profiling in childhood B cell acute lymphoblastic leukemia. Epigenomics. 2010, 2: 697-708. 10.2217/epi.10.39.
Meshinchi S, Alonzo TA, Stirewalt DL, Zwaan M, Zimmerman M, Reinhardt D, Kaspers GJL, Heerema NA, Gerbing R, Lange BJ, Radich JP: Clinical implications of FLT3 mutations in pediatric AML. Blood. 2006, 108: 3654-3661. 10.1182/blood-2006-03-009233.
Juhl-Christensen C, Ommen HB, Aggerholm A, Lausen B, Kjeldsen E, Hasle H, Hokland P: Genetic and epigenetic similarities and differences between childhood and adult AML. Pediatr Blood Cancer. 2012, 58: 525-531. 10.1002/pbc.23397.
Liang D-C, Liu H-C, Yang C-P, Jaing T-H, Hung I-J, Yeh T-C, Chen S-H, Hou J-Y, Huang Y-J, Shih Y-S, Huang Y-H, Lin T-H, Shih L-Y: Cooperating gene mutations in childhood acute myeloid leukemia with special reference on mutations of ASXL1, TET2, IDH1, IDH2, and DNMT3A. Blood. 2013, 121: 2988-2995. 10.1182/blood-2012-06-436782.
Tallman MS, Kim HT, Paietta E, Bennett JM, Dewald G, Cassileth PA, Wiernik PH, Rowe JM: Acute Monocytic Leukemia (French-American-British classification M5) does not have a worse prognosis than other subtypes of Acute Myeloid Leukemia: a report from the eastern cooperative oncology group. J Clin Oncol. 2004, 22: 1276-1286. 10.1200/JCO.2004.08.060.
Bullinger L, Döhner K, Bair E, Fröhling S, Schlenk RF, Tibshirani R, Döhner H, Pollack JR: Use of gene-expression profiling to identify prognostic subclasses in adult acute myeloid leukemia. N Engl J Med. 2004, 350: 1605-1616. 10.1056/NEJMoa031046.
Schoch C, Schnittger S, Klaus M, Kern W, Hiddemann W, Haferlach T: AML with 11q23/MLL abnormalities as defined by the WHO classification: incidence, partner chromosomes, FAB subtype, age distribution, and prognostic impact in an unselected series of 1897 cytogenetically analyzed AML cases. Blood. 2003, 102: 2395-2402. 10.1182/blood-2003-02-0434.
Parker H, Rose-Zerilli MJJ, Parker A, Chaplin T, Wade R, Gardiner A, Griffiths M, Collins A, Young BD, Oscier DG, Strefford JC: 13q deletion anatomy and disease progression in patients with chronic lymphocytic leukemia. Leukemia. 2011, 25: 489-497. 10.1038/leu.2010.288.
Bibikova M, Barnes BT, Chan T, Ho V, Klotzle B, Le JM, Delano D, Zhang L, Schroth GP, Gunderson KL, Fan J-B, Shen R: High density DNA methylation array with single CpG site resolution. Genomics. 2011, 98: 288-295. 10.1016/j.ygeno.2011.07.007.
Sarova I, Brezinova J, Zemanova Z, Bystricka D, Krejcik Z, Soukup P, Vydra J, Cermak J, Jonasova A, Michalova K: Characterization of chromosome 11 breakpoints and the areas of deletion and amplification in patients with newly diagnosed acute myeloid leukemia. Genes Chromosomes Cancer. 2013, 52: 619-635.
Calin GA, Cimmino A, Fabbri M, Ferracin M, Wojcik SE, Shimizu M, Taccioli C, Zanesi N, Garzon R, Aqeilan RI, Alder H, Volinia S, Rassenti L, Liu X, Liu C-G, Kipps TJ, Negrini M, Croce CM: MiR-15a and miR-16-1 cluster functions in human leukemia. Proc Natl Acad Sci U S A. 2008, 105: 5166-5171. 10.1073/pnas.0800121105.
Morenos L, Saffery R, Mechinaud F, Ashley D, Elwood N, Craig JM, Wong NC: Evaluation of MicroRNA expression in patient bone marrow aspirate slides. PLoS One. 2012, 7: e42951-10.1371/journal.pone.0042951.
Monteys AM, Spengler RM, Wan J, Tecedor L, Lennox KA, Xing Y, Davidson BL: Structure and activity of putative intronic miRNA promoters. RNA. 2010, 16: 495-505. 10.1261/rna.1731910.
Corcoran MM, Hammarsund M, Zhu C, Lerner M, Kapanadze B, Wilson B, Larsson C, Forsberg L, Ibbotson RE, Einhorn S, Oscier DG, Grander D, Sangfelt O: DLEU2 encodes an antisense RNA for the putative bicistronic RFP2/LEU5 gene in humans and mouse. Genes Chromosomes Cancer. 2004, 40: 285-297. 10.1002/gcc.20046.
Liu Q, Fu H, Sun F, Zhang H, Tie Y, Zhu J, Xing R, Sun Z, Zheng X: miR-16 family induces cell cycle arrest by regulating multiple cell cycle genes. Nucleic Acids Res. 2008, 36: 5391-5404. 10.1093/nar/gkn522.
Aqeilan RI, Calin GA, Croce CM: miR-15a and miR-16-1 in cancer: discovery, function and future perspectives. Cell Death Differ. 2010, 17: 215-220. 10.1038/cdd.2009.69.
Cittelly DM, Das PM, Salvo VA, Fonseca JP, Burow ME, Jones FE: Oncogenic HER2D16 suppresses miR-15a/16 and deregulates BCL-2 to promote endocrine resistance of breast tumors. Carcinogenesis. 2010, 31: 2049-2057. 10.1093/carcin/bgq192.
Hanlon K, Rudin CE, Harries LW: Investigating the targets of MIR-15a and MIR-16-1 in patients with Chronic Lymphocytic Leukemia (CLL). PLoS One. 2009, 4: e7169-10.1371/journal.pone.0007169.
Roccaro AM, Sacco A, Thompson B, Leleu X, Azab AK, Azab F, Runnels J, Jia X, Ngo HT, Melhem MR, Lin CP, Ribatti D, Rollins BJ, Witzig TE, Anderson KC, Ghobrial IM: MicroRNAs 15a and 16 regulate tumor proliferation in multiple myeloma. Blood. 2009, 113: 6669-6680. 10.1182/blood-2009-01-198408.
Zhao H, Kalota A, Jin S, Gewirtz AM: The c-myb proto-oncogene and microRNA-15a comprise an active autoregulatory feedback loop in human hematopoietic cells. Blood. 2009, 113: 505-516. 10.1182/blood-2008-01-136218.
Reutershan J, Cagnina RE, Chang D, Linden J, Ley K: Therapeutic anti-inflammatory effects of myeloid cell adenosine receptor A2a stimulation in lipopolysaccharide-induced lung injury. J Immunol. 2007, 179: 1254-1263. 10.4049/jimmunol.179.2.1254.
Wolff F, Loipetzberger A, Gruber W, Esterbauer H, Aberger F, Frischauf AM: Imiquimod directly inhibits Hedgehog signalling by stimulating adenosine receptor/protein kinase A-mediated GLI phosphorylation. Oncogene. 2013, 32: 5574-5581. 10.1038/onc.2013.343.
Andersson A, Ritz C, Lindgren D, Eden P, Lassen C, Heldrup J, Olofsson T, Rade J, Fontes M, Porwit-MacDonald A, Behrendtz M, Hoglund M, Johansson B, Fioretos T: Microarray-based classification of a consecutive series of 121 childhood acute leukemias: prediction of leukemic and genetic subtype as well as of minimal residual disease status. Leukemia. 2007, 21: 1198-1203. 10.1038/sj.leu.2404688.
Neff T, Armstrong SA: Recent progress toward epigenetic therapies: the example of mixed lineage leukemia. Blood. 2013, 121: 4847-4853. 10.1182/blood-2013-02-474833.
Haferlach T, Schoch C, Schnittger S, Kern W, Loffler H, Hiddemann W: Distinct genetic patterns can be identified in acute monoblastic and acute monocytic leukaemia (FAB AML M5a and M5b): a study of 124 patients. Br J Haematol. 2002, 118: 426-431. 10.1046/j.1365-2141.2002.03599.x.
Meyer C, Kowarz E, Hofmann J, Renneville A, Zuna J, Trka J, Abdelali RB, Macintyre E, De Braekeleer E, De Braekeleer M, Delabesse E, de Oliveira MP, Cave H, Clappier E, van Dongen JJM, Balgobind BV, van den Heuvel-Eibrink MM, Beverloo HB, Panzer-Grumayer R, Teigler-Schlegel A, Harbott J, Kjeldsen E, Schnittger S, Koehl U, Gruhn B, Heidenreich O, Chan LC, Yip SF, Krzywinski M, Eckert C: New insights to the MLL recombinome of acute leukemias. Leukemia. 2009, 23: 1490-1499. 10.1038/leu.2009.33.
Schoch C, Kohlmann A, Dugas M, Kern W, Schnittger S, Haferlach T: Impact of trisomy 8 on expression of genes located on chromosome 8 in different AML Subgroups. Genes Chromosomes Cancer. 2006, 45: 1164-1168. 10.1002/gcc.20380.
Wong NC, Ashley D, Chatterton Z, Parkinson-Bates M, Ng HK, Halemba MS, Kowalczyk A, Bedo J, Wang Q, Bell K, Algar E, Craig JM, Saffery R: A distinct DNA methylation signature defines pediatric pre-B cell acute lymphoblastic leukemia. Epigenetics. 2012, 7: 535-541. 10.4161/epi.20193.
Rosenfeld C, Goutner A, Choquet C, Venuat AM, Kayibanda B, Pico JL: Phenotypic characterization of a unique non-T, non-B acute lymphoblastic leukaemia cell line. Nature. 1977, 267: 841-843. 10.1038/267841a0.
Kaplan J, Mastrangelo R, Peterson WDJ: Childhood lymphoblastic lymphoma, a cancer of Thymus-derived lymphocytes. Cancer Res. 1974, 34: 521-525.
Asou H, Tashiro S, Hamamoto K, Otsuji A, Kita K, Kamada N: Establishment of a human acute myeloid leukemia cell line (Kasumi-1) with 8;21 chromosome translocation. Blood. 1991, 77: 2031-2036.
Tsuchiya S, Yamabe M, Yamaguchi Y, Kobayashi Y, Konno T, Tada K: Establishment and characterization of a human acute monocytic leukemia cell line (THP-1). Int J Cancer. 1980, 26: 171-176. 10.1002/ijc.2910260208.
Lange B, Valtieri M, Santoli D, Caracciolo D, Mavilio F, Gemperlein I, Griffin C, Emanuel B, Finan J, Nowell P: Growth factor requirements of childhood acute leukemia: establishment of GM-CSF-dependent cell lines. Blood. 1987, 70: 192-199.
Du P, Kibbe WA, Lin SM: lumi: a pipeline for processing Illumina microarray. Bioinformatics. 2008, 24: 1547-1548. 10.1093/bioinformatics/btn224.
Smyth GK: limma: Linear Models for Microarray Data. 2005, New York, USA: Springer
Eden E, Navon R, Steinfeld I, Lipson D, Yakhini Z: GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinformatics. 2009, 10: 48-10.1186/1471-2105-10-48.
Li J, Smyth P, Flavin R, Cahill S, Denning K, Aherne S, Guenther SM, O’Leary JJ, Sheils O: Comparison of miRNA expression patterns using total RNA extracted from matched samples of formalin-fixed paraffin-embedded (FFPE) cells and snap frozen cells. BMC Biotechnol. 2007, 7: 36-10.1186/1472-6750-7-36.
Borze I, Guled M, Musse S, Raunio A, Elonen E, Saarinen-Pihkala U, Karjalainen-Lindsberg M-L, Lahti L, Knuutila S: MicroRNA microarrays on archive bone marrow core biopsies of leukemias—Method validation. Leuk Res. 2011, 35: 188-195. 10.1016/j.leukres.2010.08.005.
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2-DDCt method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.
Betel D, Wilson M, Gadow A, Merks DS, Sander C: The microRNA.org resource: targets and expression. Nucleic Acids Res. 2008, 36: D149-D153.
John B, Enright AJ, Aravin A, ETuschl T, Sander C, Merks DS: Human microRNA targets. PLoS Biol. 2004, 2: e363-10.1371/journal.pbio.0020363.
Dweep H, Sticht C, Pandey P, Gretz N: miRWalk–database: prediction of possible miRNA binding sites by “walking” the genes of three genomes. J Biomed Inform. 2011, 44: 839-847. 10.1016/j.jbi.2011.05.002.
We thank Dr Elizabeth Algar and Priscilla Siswara for access to flash frozen bone marrow and gene abnormality identification, and Dr. Pei Tian for facilitating sample procurement. Dr. Jovana Maksimovic we thank for HM450 data analysis.
This work was supported by the Victorian Cancer Agency Grant and National Health & Medical Research Council (NHMRC), Australia [Project Grant number 607306] (N.W.); MCRI is supported by the Victorian Government’s Operational Infrastructure Support Program; the National Health & Medical Research Council Dora Lush Postgraduate Scholarship to L.M.; NHMRC Senior Research Fellowship to R.S.; The Leukaemia Foundation of Australia Phillip Desbrow Postdoctoral Fellowship (N.W.) and Project Grant in Aid( N.W. & R.S.); N.W. was also supported by My Room and the Children’s Cancer Centre Foundation.
The authors declare that they have no competing interests.
MH and MP-B carried out cell sorting, extractions and the accumulation of samples. JN assisted in accumulation of samples, extractions, conversions and assay running, as well as preliminary data analysis. FM and NG facilitated sample procurement and conceived the study. ZC participated in statistical and data analysis, study concepts and helped to draft the manuscript. RS and NW conceived the study, acquired the samples, interpreted the data, participated in drafting the manuscript and provided critical revisions for the approved final version. LM designed the study, acquired the samples, organised sample processing, participated in sample processing, undertook all data analysis and statistical evaluations and drafted the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: HM450 methylation heat map of the top 137 significantly differentially methylated probes between paediatric AML (FAB subtype M5) compared to non-leukaemic controls.(PDF 374 KB)
Additional file 2: Paediatric AML DNA methylation results compared to non-leukaemic methylation for the differentially methylated probes identified using HM450.(PDF 175 KB)
Additional file 3: DNA methylation cross-platform validation between HM450 and SEQUENOM MassARRAY® EpiTYPER®.(PDF 431 KB)
Additional file 5: Patient information for all paediatric AML patients used for HM450 DNA methylation interrogation, SEQUENOM MassARRAY® EpiTYPER® validation and gene expression analysis.(PDF 197 KB)
Additional file 6: HM450 probe methylation values for the 13q4 region for paediatric AML (FAB M5) compared to non-leukaemic controls.(PDF 484 KB)
About this article
Cite this article
Morenos, L., Chatterton, Z., Ng, J.L. et al. Hypermethylation and down-regulation of DLEU2 in paediatric acute myeloid leukaemia independent of embedded tumour suppressor miR-15a/16-1. Mol Cancer 13, 123 (2014). https://doi.org/10.1186/1476-4598-13-123
- DNA ethylation