Skip to main content

Hypermethylation and down-regulation of DLEU2 in paediatric acute myeloid leukaemia independent of embedded tumour suppressor miR-15a/16-1

Abstract

Background

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.

Methods

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.

Results

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.

Conclusions

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.

Background

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[1]. 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)[2] and more recently World Health Organisation (WHO)[3] 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[6] 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[6]. 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[4]. Despite such emerging findings in an adult context, the utility of individual DNA methylation disruptions in paediatric AML has yet to be fully evaluated[6].

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[12]. Human miRNA are initially transcribed (pri-miRNA), and processed by several complexes to form a 70 bp hairpin-loop (pre-miRNA)[13]. 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[16]) 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[21]. 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[22]. 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[22]). Additionally, miR-193a targeting KIT, and miR-14b targeting CREB have been identified in adult investigations as specifically controlled by DNA methylation (reviewed in[22]). 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, 2426]. 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[27]. AML subtype M5 also has a high proportion of cytogenetically normal (CN-AML) patients[27], and those with complex karyotypes[28]. A combination of these factors add to the overall unfavourable outcome of paediatric M5 diagnosis[29]. 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[31]. Henceforth this will be referred to as the DLEU2 promoter.

Figure 1
figure 1

Regional interrogation of the DLEU2 gene including significantly differentially methylated probes identified by HM450 analysis. Bottom : Distribution of HM450 methylation probes across the DLEU2 region of 13q4. Genes located in this region include: DLEU1, DLEU2, miR15a/16-1 microRNA cluster, TRIM13, KCNRG and miR-3613. Top : HM450 probes identified as significantly differentially methylated between paediatric AML (FAB subtype M5) and non-leukaemic specimens have been plotted against genomic location. The leukaemic group (n = 16) refers to diagnostic bone marrow from paediatric patients. Non-leukaemic group (n = 11) consists of CD sorted cell populations (CD19+, CD33+, CD34+, CD45+) and patient remission specimens. Differential methylation is concentrated in the promoter region of the DLEU2/Alt1 long transcript, which also falls into the body region of DLEU1 within 3 CpG islands (Chr 13: 50, 690,000-50,708,000). Mean methylation β values and 95% CI are shown. Data for leukaemic samples are red and non-leukaemic samples in green. Methylation values range from 0 (0%, no detected methylation) to 1.0 (100% fully methylated). Significantly differentially methylated regions between leukaemic and non-leukaemic samples to p < 0.001 are highlighted with **.

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[32]. 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).

Figure 2
figure 2

Expression analysis of DLEU2 and miR-15/16 in paediatric AML. Here we interrogate the gene, miRNA and precursor miRNA expression in paediatric AML compared to non-leukaemic. This interrogation includes DLEU2 and embedded miR-15a/16-1 on chromosome 13q4. The leukaemic group refers to diagnostic bone marrow from paediatric patients. The non-leukaemic group consists of CD sorted cell populations (CD19+, CD33+, CD34+, CD45+) and patient remission specimens, and is represented by the dashed line at Y = 1. Fold Change (FC) is plotted using normalized data and the 2-ΔΔCt method ± SD, and shows the fold change calculated from the means of each group. A. Gene expression including DLEU1, DLEU2 and TRIM13 in paediatric AML (n = 10) compared to non-leukaemic (n = 13) expression. DLEU2 shows a significant down-regulation in AML compared to non-leukaemic expression (0.07 FC, p = 0.014 represented by **), however there is no significant change in expression for TRIM13 or DLEU1. B. Mature microRNA expression, including primary precursor transcript (PRI) and alternate miRNA expression (*), from the miR-15a/16-1 miRNA cluster embedded within DLEU2 for paediatric AML (n = 28, including the 10 used in Figure 2A) compared to non-leukaemic specimens (n = 30, including the 13 used in Figure 2A). Here the miR-15a/16-1 PRI transcript is 3.03-fold higher in expression compared to non-leukaemic expression. Additionally, miR-16-1* (2.52 FC), miR-15a* (2.24 FC) and miR-15a (1.5 FC) also show increases in expression in paediatric AML. No significant change in expression was observed for miR-16 (0.94 FC) in paediatric AML compared to non-leukaemic expression.

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[34]. 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[32]. 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[35]. 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, 3642]. 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).

Figure 3
figure 3

Interrogation of paediatric AML by clinically defined cytogenetic and FAB subtypes alongside DLEU2 promoter methylation. Paediatric AML diagnostic cytogenetic and subtyping analyses were specifically investigated, including those with known gene abnormalities and CN-AML cases. We investigate here the DNA Methylation of DLEU2/Alt1 promoter region to 95% CI. DNA methylation values range from 0.0 (0% no detected methylation) to 1.0 (100% fully methylated). The leukaemic group refers to diagnostic bone marrow from paediatric patients. Non-leukaemic group consists of CD sorted cell populations (CD19+, CD33+, CD34+, CD45+) and patient remission specimens. A. DNA methylation for DLEU2 HM450 probes (cg12883980, cg20529344, cg5394800) according to cytogenetic type, showing heterogeneous DNA methylation outcomes. Of note, a subset of patients with observable chromosome 11 and trisomy 8 abnormalities have a reduced DNA methylation compared to all other AML abnormalities. B. A subset of DLEU2 DNA methylation results for paediatric AML chromosome 11 and trisomy 8 abnormality patients from Figure 3A fall into a range akin to non-leukaemic specimens. Grouping these patients, regardless of clinical subtyping (M5a, M5b or M1/M2/M4), forms an additional sub-group, denoted here as ‘t(11)/+8’. The DNA methylation values obtained for probe cg12883980 are depicted, as an indication of the methylation at the DLEU2 promoter. t(11)/+8 are not differentially methylated compared to non-leukaemic specimens (t(11)/+8 sub-group Mean Methylation (MM) = 0.38, non-leukaemic MM = 0.32. p = 0.22). M5a, M5b and M1/M2/M4 AMLs (with the removal of t(11)/+8 patients) have become increasingly hypermethylated compared to non-leukaemic methylation (M5a MM = 0.81, M5b MM = 0.78, M1/M2/M4 MM = 0.87. p < 0.0001), and significantly different from t(11)/+8 (M5a/b p < 0.05; M1/M2/M4 p < 0.0001).

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[29]. 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[49]. 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).

Figure 4
figure 4

Gene, primary, miRNA and miRNA* expression for paediatric AML defined through clinical classification and DLEU2 Methylation subtyping. Mature microRNA expression, primary precursor transcript (PRI) and alternate miRNA isoform expression (*) from the miR-15a/16-1 miRNA cluster embedded within DLEU2 for paediatric AML patients (n = 26: 12 M5a, 5 M5b, 4 M1/M2/M4, 5 t(11)/+8 sub-group) all compared to non-leukaemic specimens (n = 30). Leukaemic groups refer to diagnostic bone marrow from paediatric patients. Non-leukaemic group consists of CD sorted cell populations (CD19+, CD33+, CD34+, CD45+)and patient remission specimens. Linear Fold Change (FC) is plotted using normalized data and the 2-ΔΔCt method ± SD, and shows the fold change calculated from the means of each group. DLEU2 gene expression is down-regulated in all subtypes. A. Fold change in expression comparing non-leukaemic to t(11)/+8 subtype. No significant differences in RNA expression are observed between non-leukaemic specimens and this subgroup. B. Fold change in expression comparing patients from subtype M5a to M1/M2/M4. DLEU2 is down-regulated in all subtypes (as previously described in Figure 2), and primary precursor for miR-15a also appears down-regulated (non-significant). M1/M2/M4 groupings do not show any mature miRNA expression changes (defined as >2-fold difference from non-leukaemic). Subtype M5a shows a 2.1-fold (±0.8 SD) increase in miR-15a* and a 2.81-fold (±0.4 SD) increase in miR-16-1*. C. FAB subtype M5b shows up-regulation of miR-15a PRI compared to non-leukaemic expression (20.29-fold (±0.6 SD) p < 0.001), with miR-15a* (4.87-fold ±1 SD) and miR-16-1* (9.86-fold ±1 SD) also up-regulated. M5b additionally shows a significant up-regulation of miR-15a PRI compared to t(11)/+8 sub-group samples (p < 0.05), and also from M5a and M1/M2/M4 sub-groups (p < 0.001).

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.

Conclusions

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.

Methods

Samples

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[51], CCRF-CEM[52]), paediatric AML (Kasumi-1[53], THP-1[54], MV-4-11 and AML-193[55]) 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[23]. 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[56] and analysed using LIMMA[57]. 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:

Δβ = mean β value Leukaemic - mean β value Non - Leukaemic

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[58].

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[34]. 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[34]. 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[61] 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[64]. 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[58].

Statistical methods

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.

Abbreviations

AML:

Acute myeloid leukaemia

CN-AML:

Cytogenetically normal acute myeloid leukaemia

DLEU2:

Deleted in lymphocytic leukaemia 2

DMP:

Differentially methylated probe

FAB:

French-American-British leukaemia classification method

HM450:

Illumina humanmethylation450 bead chip array

miR:

microRNA/miRNA

Pri-miRNA:

Primary microRNA

TSS:

Transcriptional start site.

References

  1. 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.

    Article  CAS  PubMed  Google Scholar 

  2. 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.

    Article  CAS  PubMed  Google Scholar 

  3. 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.

    Article  CAS  PubMed  Google Scholar 

  4. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. 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.

    Article  PubMed Central  PubMed  Google Scholar 

  6. 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.

    Article  PubMed Central  PubMed  Google Scholar 

  7. Gardiner-Garden M, Frommer M: CpG islands in vertebrate genomes. J Mol Biol. 1987, 196: 261-282. 10.1016/0022-2836(87)90689-9.

    Article  CAS  PubMed  Google Scholar 

  8. 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.

    Article  CAS  PubMed  Google Scholar 

  9. Bird AP: DNA methylation patterns and epigenetic memory. Genes Dev. 2002, 16: 6-21. 10.1101/gad.947102.

    Article  CAS  PubMed  Google Scholar 

  10. Bissels U, Bosio A, Wagner W: MicroRNAs are shaping the hematopoietic landscape. Haematologica. 2012, 97: 160-167. 10.3324/haematol.2011.051730.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Calin GA, Croce CM: MicroRNA signatures in human cancers. Nat Rev Cancer. 2006, 6: 857-866. 10.1038/nrc1997.

    Article  CAS  PubMed  Google Scholar 

  12. 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.

    Article  CAS  PubMed  Google Scholar 

  13. Calin GA: MicroRNAs and cancer: what we know and what we still have to learn. Genome Med. 2009, 1: 78-10.1186/gm78.

    Article  PubMed Central  PubMed  Google Scholar 

  14. 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.

    Article  CAS  PubMed  Google Scholar 

  15. Zhao G, Yu D, Weiss MJ: MicroRNAs in erythropoiesis. Curr Opin Hematol. 2010, 17: 155-162.

    CAS  PubMed  Google Scholar 

  16. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  18. 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.

    Article  CAS  PubMed  Google Scholar 

  19. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. 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.

    Article  CAS  PubMed  Google Scholar 

  21. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. 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.

    Article  CAS  PubMed  Google Scholar 

  23. 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.

    Article  CAS  PubMed  Google Scholar 

  24. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. 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.

    Article  PubMed  Google Scholar 

  26. 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.

    Article  CAS  PubMed  Google Scholar 

  27. 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.

    Article  PubMed  Google Scholar 

  28. 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.

    Article  CAS  PubMed  Google Scholar 

  29. 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.

    Article  CAS  PubMed  Google Scholar 

  30. 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.

    Article  CAS  PubMed  Google Scholar 

  31. 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.

    Article  CAS  PubMed  Google Scholar 

  32. 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.

    Article  CAS  PubMed  Google Scholar 

  33. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. 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.

    Article  PubMed Central  PubMed  Google Scholar 

  36. 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.

    Article  CAS  PubMed  Google Scholar 

  37. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. 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.

    Article  CAS  PubMed  Google Scholar 

  39. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. 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.

    Article  PubMed Central  PubMed  Google Scholar 

  41. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  42. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  43. 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.

    Article  CAS  PubMed  Google Scholar 

  44. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  45. 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.

    Article  CAS  PubMed  Google Scholar 

  46. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  47. 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.

    Article  CAS  PubMed  Google Scholar 

  48. 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.

    Article  CAS  PubMed  Google Scholar 

  49. 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.

    Article  CAS  PubMed  Google Scholar 

  50. 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.

    Article  CAS  PubMed  Google Scholar 

  51. 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.

    Article  CAS  PubMed  Google Scholar 

  52. Kaplan J, Mastrangelo R, Peterson WDJ: Childhood lymphoblastic lymphoma, a cancer of Thymus-derived lymphocytes. Cancer Res. 1974, 34: 521-525.

    CAS  PubMed  Google Scholar 

  53. 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.

    CAS  PubMed  Google Scholar 

  54. 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.

    Article  CAS  PubMed  Google Scholar 

  55. 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.

    CAS  PubMed  Google Scholar 

  56. Du P, Kibbe WA, Lin SM: lumi: a pipeline for processing Illumina microarray. Bioinformatics. 2008, 24: 1547-1548. 10.1093/bioinformatics/btn224.

    Article  CAS  PubMed  Google Scholar 

  57. Smyth GK: limma: Linear Models for Microarray Data. 2005, New York, USA: Springer

    Google Scholar 

  58. 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.

    Article  PubMed Central  PubMed  Google Scholar 

  59. 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.

    Article  PubMed Central  PubMed  Google Scholar 

  60. 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.

    Article  CAS  PubMed  Google Scholar 

  61. 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.

    Article  CAS  PubMed  Google Scholar 

  62. Betel D, Wilson M, Gadow A, Merks DS, Sander C: The microRNA.org resource: targets and expression. Nucleic Acids Res. 2008, 36: D149-D153.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  63. 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.

    Article  PubMed Central  PubMed  Google Scholar 

  64. 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.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

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.

Financial support

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.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Leah Morenos.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

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

12943_2013_1333_MOESM1_ESM.pdf

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)

12943_2013_1333_MOESM2_ESM.pdf

Additional file 2: Paediatric AML DNA methylation results compared to non-leukaemic methylation for the differentially methylated probes identified using HM450.(PDF 175 KB)

12943_2013_1333_MOESM3_ESM.pdf

Additional file 3: DNA methylation cross-platform validation between HM450 and SEQUENOM MassARRAY® EpiTYPER®.(PDF 431 KB)

Additional file 4: CpG locations captured by SEQUENOM MassARRAY® EpiTYPER®, including the HM450 probes.(PDF 512 KB)

12943_2013_1333_MOESM5_ESM.pdf

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)

12943_2013_1333_MOESM6_ESM.pdf

Additional file 6: HM450 probe methylation values for the 13q4 region for paediatric AML (FAB M5) compared to non-leukaemic controls.(PDF 484 KB)

Additional file 7: Correlation between DLEU2/Alt1 promoter DNA methylation and DLEU2 gene expression.(PDF 447 KB)

12943_2013_1333_MOESM8_ESM.pdf

Additional file 8: HM450 probe methylation and expression dynamics of paediatric AML compared to non-leukaemic controls for the alternate miR-16 cluster, miR-15b/16-2 embedded in SMC4 on chromosome 3.(PDF 599 KB)

Additional file 9: Comparison of individual patient miRNA miR-15a/16-1 cluster expression.(PDF 517 KB)

Additional file 10: Top Gene Ontology terms using GOrilla gene ontology identification tool.(PDF 308 KB)

12943_2013_1333_MOESM11_ESM.pdf

Additional file 11: DLEU2 promoter DNA methylation divided into paediatric AML patient clinical diagnostic observations.(PDF 513 KB)

12943_2013_1333_MOESM12_ESM.pdf

Additional file 12: Survival and relapse curves of paediatric AML, including a newly defined sub-group based on DLEU2 promoter DNA methylation.(PDF 540 KB)

Additional file 13: DLEU2 DNA methylation comparison of the three ‘non-leukaemic’ sample types.(PDF 457 KB)

12943_2013_1333_MOESM14_ESM.pdf

Additional file 14: Probes used in SEQUENOM MassARRAY® EpiTYPER® validation analysis of DLEU2/Alt1 HM450 promoter probes located on Chromosome 13.(PDF 283 KB)

12943_2013_1333_MOESM15_ESM.pdf

Additional file 15: miR-15a/16-1 cluster TaqMan® microRNA and primary precursor expression assays utilized, including Normalization references.(PDF 295 KB)

Additional file 16: TaqMan® gene expression assays utilized.(PDF 174 KB)

12943_2013_1333_MOESM17_ESM.xlsx

Additional file 17: The top 500 potential target genes identified for each differentially expressed miR-15a/16-1 cluster microRNA member in paediatric AML.(XLSX 26 KB)

Authors’ original submitted files for images

Rights and permissions

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 credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1476-4598-13-123

Keywords