CBFB-MYH11 hypomethylation signature and PBX3 differential methylation revealed by targeted bisulfite sequencing in patients with acute myeloid leukemia
© Hajkova et al.; licensee BioMed Central Ltd. 2014
Received: 25 March 2014
Accepted: 6 September 2014
Published: 30 September 2014
Studying DNA methylation changes in the context of structural rearrangements and point mutations as well as gene expression changes enables the identification of genes that are important for disease onset and progression in different subtypes of acute myeloid leukemia (AML) patients. The aim of this study was to identify differentially methylated genes with potential impact on AML pathogenesis based on the correlation of methylation and expression data.
The primary method of studying DNA methylation changes was targeted bisulfite sequencing capturing approximately 84 megabases (Mb) of the genome in 14 diagnostic AML patients and a healthy donors' CD34+ pool. Subsequently, selected DNA methylation changes were confirmed by 454 bisulfite pyrosequencing in a larger cohort of samples. Furthermore, we addressed gene expression by microarray profiling and correlated methylation of regions adjacent to transcription start sites with expression of corresponding genes.
Here, we report a novel hypomethylation pattern, specific to CBFB-MYH11 fusion resulting from inv(16) rearrangement that is associated with genes previously described as upregulated in inv(16) AML. We assume that this hypomethylation and corresponding overexpresion occurs in the genes whose function is important in inv(16) leukemogenesis. Further, by comparing all targeted methylation and microarray expression data, PBX3 differential methylation was found to correlate with its gene expression. PBX3 has been recently shown to be a key interaction partner of HOX genes during leukemogenesis and we revealed higher incidence of relapses in PBX3-overexpressing patients.
We discovered new genomic regions with aberrant DNA methylation that are associated with expression of genes involved in leukemogenesis. Our results demonstrate the potential of the targeted approach for DNA methylation studies to reveal new regulatory regions.
KeywordsAcute myeloid leukemia CBFB-MYH11 DNA methylation Targeted bisulfite sequencing PBX3
Changes in DNA methylation patterns are a known hallmark of acute myeloid leukemia (AML) and underlie AML pathogenesis . DNA methylation in patients with AML has been studied extensively and may reflect either specific molecular abnormalities or characterize a group of patients without an apparent molecular aberration. Specific translocations such as PML-RARA, AML1-ETO (RUNX1-RUNX1T1), MLL translocations or CBFB-MYH11 fusion, as well as CEBPA, NPM1, IDH1/IDH2, DNMT3A, TET2 and RUNX1 mutations have been described to display distinct methylation signatures -. These epigenetic profiles are usually accompanied by specific gene expression features. Studying genes that are epigenetically deregulated in different groups of patients may contribute to a more detailed understanding of pathways involved in the leukemic transformation. Importantly, the effect of DNA methylation changes is greatly dependent on the location of differentially methylated regions (DMRs) . New approaches using next-generation sequencing enable studying of DMRs scattered throughout the genome and targeted bisulfite sequencing offers a reasonably balanced ratio between cost and informativeness (number of CpGs covered) . The link between gene expression and DNA methylation data is needed to find pathologically relevant DNA methylation changes, especially because many (or even the majority of) DMRs reflect the tissue of origin and not leukemia (cancer) specific changes .
In this study, 84 megabases (Mb) of 14 AML genomes and one CD34+ pool of cells from healthy donors were captured for DNA methylation and gene expression profiling. The aim was to identify differentially methylated genes with potential impact on AML pathogenesis based on the correlation of methylation and expression data.
AML patients with CBFB-MYH11 fusion (CBFB - Core-binding factor, beta subunit; MYH11 - Myosin, heavy chain 11) resulting from inv(16) rearrangement clustered together in a hierarchal DNA methylation and expression analysis. The majority of differentially methylated regions unique for CBFB-MYH11 patients were hypomethylated and genes assigned to such regions were previously described as overexpressed in inv(16) AML .
PBX3 (pre-B-cell leukemia homeobox 3), recently demonstrated as an important cofactor of HOXA9 in leukemogenesis , was validated as a gene whose gene expression levels correlated with DNA methylation of its putative regulatory region across AML subtypes. The importance of PBX3 is underlined by the fact that PBX3-overexpressing patients relapse more frequently. In summary, we discovered new genomic regions affected by aberrant DNA methylation that are associated with expression of genes implicated in leukemogenesis.
Inv(16) methylation and expression cluster
Characteristics of AML patients
CBFB-MYH11, FLT3- ITD
de novo AML
de novo AML
NPM1, DNMT3A, FLT3- ITD
AML with multilineage dysplasia
AML with multilineage dysplasia
NPM1, DNMT3A, FLT3- ITD
de novo AML
de novo AML
de novo AML
NPM1, FLT3- ITD
de novo AML
de novo AML
de novo AML
AML after MDS
AML after MDS
pool of 4 CD34+ healthy control cells
Clustering of gene expression data for 8 884 genes with detection P value ≤ 0.05 in all 14 AML patients and 4 CD34+ healthy control cells confirmed the inv(16) cluster (Figure 1B). Consistent with the DNA methylation data, no other characteristics formed clusters. The selection of these 14 AML patients was done randomly throughout various clinical subgroups (de novo, secondary after MDS, secondary after breast cancer, AML with dysplasia) with the aim to find common DNA methylation changes among diverse AML subtypes. Therefore, we cannot make any general conclusions with regard to the existence/non-existence of specific DNA methylation/gene expression profiles other than those demontrated in inv(16) AML, which were sufficiently verified in a validation cohort (see below).
Hypomethylation signature of inv(16) AML patients
We extracted genomic regions being uniquely differentially methylated in CBFB-MYH11 patients. There was a clear tendency towards hypomethylation with 125 out of 182 regions (69%) displaying lower DNA methylation levels compared to the healthy donors' CD34+ pool of cells. All CBFB-MYH11 DMRs were uploaded to GREAT  and enrichment for genes previously described as upregulated in inv(16) AML patients was observed (reported in ref. , ID: VALK_AML_CLUSTER_9; Additional file 2: Table S1). This enrichment set comprised 10 genomic regions assigned to 6 genes - MN1, SPARC, ST18, DHRS3, FAM171A and BAHCC1 (see Additional file 3: Table S2). MN1, SPARC, ST18, FAM171A and DHRS3 were chosen for hypomethylation validation by 454 bisulfite pyrosequencing. As we were primarily interested in hypomethylation associated with overexpression, BAHCC1 was excluded from the methylation/expression validation, because its expression levels were undetectable in AML as well as in healthy donors' samples according to microarray expression data. For MN1 and SPARC, two regions per gene were studied.
As we used either PB or BM as a starting material, we investigated whether there is a concordance between DNA methylation results from PB and BM in 10 AML patients with both materials at diagnosis available. We found high correlation between PB and BM results for all of the studied regions (R = 0.96).
In summary, ST18 is the only overexpressed gene in inv(16) compared to other AML subtypes and healthy donors, irrespectively of the RQ-PCR system used.
Correlation of methylation data with expression
For each transcription start site (TSS), a window was defined from 5 kb up- to 1 kb downstream of that TSS. Targeted regions overlapping this window were associated with the gene of the given TSS. GENCODE v14 gene annotation was used for assigning DNA methylation of target regions to corresponding genes. Expression and methylation data were then correlated using Spearman rank tests. The small sample size did not allow for filtering based on P values and thus an ad hoc measure was employed, requiring strong anti-correlation (R ≤-0.7), and a change of methylation between control and at least one AML sample 2-fold or greater (with either the control or at least one of the AMLs having methylation ratio 0.3 or greater). These strict filtering parameters resulted in a list of 163 genomic regions assigned to 130 unique genes (see Additional file 5: Table S3). Among these genes, we chose PBX3 and GFI1 to validate the association between DNA methylation and gene expression changes observed in this cohort of 14 AML patients.
PBX3 and GFI1 differential methylation is involved in their expression regulation
The DMRs associated with PBX3 are located downstream (-160 to -1451 bp) of annotated CpG island (for exact location see Additional file 5: Table S3). We focused on a DMR encompassing a TAF1 binding site (chr9:128 510 974-128 511 259 according to GRCh37/hg19). TAF1 encodes the largest subunit of TFIID and this subunit binds to core promoter sequences encompassing TSSs.
Regarding GFI1, we examined a 409 bp long region within the GFI1 promoter (chr1:92 952 229-92 952 637). Average DNA methylation levels of healthy controls (n = 14) fluctuated around 3.86%. Very similar to that, DNA methylation levels of AML within normal gene expression range were 4.54%. Only a minority of AML samples (3/64) showed decreased GFI1 expression levels. However, all 3 GFI1-downregulated patients displayed higher methylation levels - 22.8%, 28.9% and 19%, respectively. On the contrary, upregulation of GFI1 was not connected with any significant changes of DNA methylation levels, probably due to the fact that normal DNA methylation levels were already low, therefore any further decrease would not have a functional role.
PBX3 expression levels and their impact on prognosis of AML patients
Because PBX3 is one of the four genes, whose common expression signature was described as having an impact on overall survival (OS) in AML patients , we decided to evaluate its expression levels in terms of OS and relapse-free survival (RFS). Only patients receiving standard curative therapy and those who did not die during the first induction were included in this analysis. Altogether 40 AML patients were assessed, 21/40 had low and 19/40 had high expression levels. Low and high expression levels were defined as a change in expression of at least one order of magnitude as well as 2-fold from healthy donors' average PBX3 expression.
We also performed multivariate analysis for both OS and RFS. Firstly, we tested the following parameters by univariate analysis: age, white blood count (WBC), complete remission after induction therapy (CR after induction), good/intermediate/poor prognosis according to cytogenetics, PBX3 expression and FLT3-ITD status. For OS, only CR after induction was statistically significant in both univariate and multivariate analyses (P < 0.001). For RFS, parameters significant in univariate testing (WBC, CR after induction, FLT3-ITD status and PBX3 expression) were evaluated by Cox regression. Only CR after induction and PBX3 expression retained statistical significance (P = 0.002 and P = 0.028, respectively) in multivariate testing.
Furthermore, we evaluated whether there is an association between PBX3 expression levels and presence of MLL-PTD or belonging to prognostically adverse AML subgroup. No such correlation was found, on the other hand there was no PBX3-overexpressing patient among AML patients within the prognostically favourable subgroup (0/6 cases from prognostically favourable versus 26/46 cases from prognostically intermediate or adverse subgroup overexpress PBX3, 0% versus 57% respectively, P = 0.02). Leukemic transformation mediated by MLL-fusion proteins has been suggested to be dependent on the presence of PBX3 expression . In our cohort, all three MLL-translocated patients, two with MLLT3-MLL (MLL-AF9) and one with MLL-MLLT1 (MLL-ENL), had upregulated PBX3 expression, but the low number of MLL translocations limited statistical testing of this association.
Here, we report a CBFB-MYH11, i.e. inv(16), specific hypomethylation that may play a role in upregulation of some previously described inv(16) overexpressed genes . We analysed DNA methylation and expression data of MN1, SPARC, ST18 and DHRS3 in 55 AML patients and 10 healthy controls. Lower methylation levels of these genes in inv(16) patients versus other AML and healthy donors were confirmed. When measured by TaqMan gene expression assays, inv(16)-specific overexpression of MN1, SPARC, ST18 and DHRS3 was found with respect to non-inv(16) AML, but only in ST18 in comparison with healthy donors. This was inconsistent with the Illumina microarray expression data as well as previously published data . Therefore, we re-measured the results using SybrGreen RQ-PCR and we obtained different values. In this RQ-PCR experiment, changes in expression levels (between inv(16) AML and healthy donors) were also significant for MN1 and SPARC. We excluded both the role of PCR nonspecificity and DNA contamination. Interestingly, MN1 and SPARC primers for SybrGreen RQ-PCR, and MN1 and SPARC probes on the expression microarray, are localized within the same exons, while MN1 and SPARC TaqMan probes have different, exon-exon localizations. We cannot claim that this is the only reason for the different results, this issue definitively needs deeper examination in future studies. It is of interest that the results obtained by SybrGreen RQ-PCR are in agreement with the publicly available data (GSE34823) of the study of Bletiere et al. .
Average gene expression levels of MN1, SPARC and ST18 extracted from the above-mentioned dataset are higher in inv(16) AML in comparison with healthy donors' bone marrow (4-times for MN1 and SPARC, 7-times for ST18), and DHRS3 expression is basically the same in both groups.
Further, we extracted the data from The Cancer Genome Atlas (TCGA, https://tcga-data.nci.nih.gov/tcga/) and they also confirmed a link between hypomethylation of MN1 and DHRS3 regulatory regions and their overexpression in inv(16) AML when compared with AML samples with normal karyotype (healthy controls data were not available). Regulatory regions corresponding to remaining genes lacked information of their methylation status in TCGA due to the absence of appropriate CpGs in HumanMethylation450 BeadChip used at TCGA study . This supports the profitability of studying DNA methylation using targeted bisulfite sequencing, which provides more complex coverage than microarray based techniques.
The hypomethylation pattern that we discovered in inv(16) AML patients is remarkable also with respect to the very recently published data of Mandoli and colleagues . For the first time, their study revealed the involvement of CBFB-MYH11 not only in repression but as well in transcriptional activation. The direct involvement of CBFB-MYH11 in overexpression of MN1, ST18 and SPARC is supported with 2-fold downregulation of these genes upon CBFB-MYH11 knockdown as reported in their work. However, none of the 1874 high-confidence CBFB-MYH11 binding sites  overlaps with any of the hypomethylated regions reported here. DHRS3 was among the genes upregulated upon CBFB-MYH11 knockdown, which is in agreement with its disputable upregulation in inv(16) AML. There were great differences in localization of hypomethylated regions with respect to TSSs of individual genes. With regard to MN1 and SPARC, the hypomethylation was located not far from their TSSs (for location see Additional file 3: Table S2), which makes the assumption of their role in the expression of these genes more straightforward. Moreover, MN1 and SPARC potential regulatory regions overlap regions of active chromatin (enhancers) in mobilized CD34+ cells as observed in the EpiGenome Browser (http://epigenomegateway.wustl.edu/browser) suggesting a role of these regions in transcription regulation. On the contrary, differentially hypomethylated sites assigned by GREAT to ST18 and DHRS3 are placed much farther from their TSSs, specifically approximately 275 kilobases (kb) downstream for ST18 and 277 kb upstream for DHRS3. ST18 and DHRS3 assigned regulatory regions are placed within chromatin marked with low transcription activity and enhancer, respectively (in mobilized CD34+ cells, data from EpiGenome Browser).
MN1 expression levels have been shown to stratify prognosis of cytogenetically normal (CN) AML patients and its overexpression is connected with a poor outcome of CN-AML patients ,. Nevertheless, inv(16) AML patients are generally associated with a good prognosis , in spite of their frequent MN1 overexpression. Functional studies have proved that overexpression of MN1 cooperates with inv(16) in developing AML in vivo and that neither inv(16) or MN1 alone are capable of promoting leukemia . According to our results it seems that hypomethylation is present uniquely in inv(16) AML patients in both MN1 assigned regions (none of the other AML subtypes or healthy controls displayed MN1 hypomethylation in our cohort). So it supports the theory that upregulation of some genes that might involve MN1 is crucial for inv(16) leukemogenesis and hypomethylation may be therefore needed to ensure stable overexpression of critical genes. Apparently the mechanism of MN1 upregulation is different in non-inv(16) AML patients or potentially hypomethylation of other regulatory areas located elsewhere may be involved.
Correlation between targeted DNA methylation and microarray expression data of 14 AML patients and a healthy controls' CD34+ pool revealed PBX3 differential methylation and gene expression. PBX3 (pre-B-cell leukemia homeobox 3) is part of the three amino acid loop extension (TALE) family of transcription factors, which include the products of the Pbx and Meis genes and are capable of heterodimerization with the Hox proteins . Recently, PBX3 was reported to have a synergistic effect with HOXA9 in transforming normal hematopoietic progenitor cells in vitro as well as in vivo . Moreover, PBX3 is one of the four genes (HOXA6, HOXA9, PBX3 and MEIS1), whose common expression signature was shown to influence overall survival in CN-AML . All evidence points to an important role of PBX3 in leukemogenesis. This is the first report uncovering DNA methylation as a plausible regulator of PBX3 expression. We found a strong negative correlation between levels of PBX3 methylation and expression in 8 healthy donors' samples and 30 AML patients at diagnosis (P < 0.0001 and P = 0.002 for upregulation and downregulation, respectively). Localization of PBX3 differential methylation overlaps TAF1 binding site according to ENCODE ChIP-Seq data from UCSC genome browser. TAFs (TBP-associated factors) create a stable complex with TBP (TATA-binding protein) and RNAPII to form a preinitiation complex, so we may assume that DNA methylation status of TAF1 binding site can directly influence the accessibility of DNA for transcription enzymes. The probability of transcription initiation is possibly dependent on whether the DNA methylation is low with high expression rates or DNA methylation is high with decreased expression or finally intermediate DNA methylation corresponding to in-between expression levels. As PBX3 has a CpG island (CGI) overlapping its TSS, we also looked at its methylation status. Based on targeted methylation data, there was no methylation present either in AML or healthy controls. Therefore, methylation status of downstream located control element rather than CGI methylation is most likely crucial for PBX3 expression. Further we focused on potential prognostic significance of PBX3 expression in terms of overall survival (OS) and incidence of relapse. High PBX3 expression levels were not related to different OS compared to AML patients with low PBX3 expression, however relapse rates were significantly higher in PBX3-overexpressing patients by both univariate and multivariate testing. This suggests more aggressive phenotype/course of disease of these patients, which is not reflected in the OS probably due to the early and effective treatment of relapses - often followed by bone marrow transplantation. We also showed that PBX3 overexpression did not occur in AML patients within cytogenetically favourable subgroup.
We validated the methylation/expression correlation stated in the Additional file 5: Table S3 also for GFI1. Moreover, the observed correlations are further supported by the presence of genes, for which the role of DNA methylation is already published such as MPO, CEBPα, DAPK1, IRF8 and PRDX8,.
We found a new hypomethylation signature specific for inv(16) AML patients that may be responsible for overexpression of some genes that are crucial for inv(16) pathogenesis. MN1 gene is likely to be a key gene involved in the pathogenesis of inv(16) AML and hypomethylation in the regulatory region near its TSS in inv(16) AML patients was confirmed, even on the basis of publicly available data from TCGA. Furthermore, we explored new regulatory region for PBX3 and association of its methylation with PBX3 expression. Therefore, targeted bisulfite sequencing represents a convenient approach in terms of genome coverage and informativeness with a great potential to reveal new regulatory regions of genes involved in leukemic transformation.
Materials and methods
For targeted bisulfite sequencing, 14 AML patients at diagnosis (see Table 1) and pooled CD34+ cells from 4 healthy donors were sequenced. Genes selected based on targeted bisulfite sequencing results were examined using 454 bisulfite pyrosequencing in a larger cohort of AML patients (their characteristics are given in Additional file 7: Table S4). Informed consent was obtained from all patients and healthy blood donors enrolled in the study. The study was approved by the IHBT Institutional Ethics Committee according to the Helsinki Declaration.
Mononuclear cells (MNC) from peripheral blood (PB) or bone marrow (BM) of the AML patients at diagnosis were separated by Ficoll gradient centrifugation (Histopaque, Sigma-Aldrich, Steinheim, Germany). CD34+ cells were harvested from buffy coats of healthy blood donors using MicroBead kits (Miltenyi Biotec GmbH, Bergish Gladbach, Germany). The CD34+ pool was created by mixing of 4 individual healthy blood donors' separated cells (all of them men aged 42 to 58 years old, median age 45.5). DNA and RNA were extracted using AllPrep DNA/RNA Mini Kit (Qiagen, Hilden, Germany). Bisulfite conversion was performed from 1 μg of DNA by EpiTect Bisulfite Kit (Qiagen) and eluted into 40 μl of EB buffer. cDNA was prepared using M-MLV RT (Moloney Murine Leukemia Virus Reverse Transcriptase, Promega, Madison, WI, USA).
Targeted bisulfite sequencing
Preparation of targeted bisulfite libraries started with 3 μg of genomic DNA and was carried out using SureSelectXT Human Methyl-Seq kit (Agilent, Agilent Technologies, Santa Clara, CA, USA) according to the manufacturer's instructions. Libraries were multiplexed into 4 pools and each pool was sequenced on 2 HiSeq-2000 (Illumina, San Diego, CA, USA) lanes using 105 bp paired-end sequencing reads with average coverage of 83 - ranging from 46 to 131.
4 bisulfite pyrosequencing
Bisulfite-treated (BS) DNA was subjected to 2-round PCR. In the 1st round of PCR, loci-specific primers with M13 universal tails were used to amplify regions of interest. Subsequently, primers specific to M13 universal tails now tailed with 454 -specific sequencing primers and a unique barcode sequence (MID) were applied to the 2nd PCR. Loci-specific primers were designed with Methyl Primer Express v1.0 software (Applied Biosystems Inc. Foster City, CA, USA; see Additional file 8: Table S5 for primer sequences). HotStarTaq DNA polymerase (Qiagen) and manufacturer's recommended PCR reaction conditions were used for amplification. 2 μl of BS DNA was added to the 1st PCR and 1 μl of 100× diluted 1st round PCR product was subjected to the 2nd PCR. PCR cycling conditions were as follows: 1st round PCR - initial denaturation (15 min at 95°C), followed by 35 cycles of denaturation (30s at 94°C), annealing (30s at Ta °C, Ta - annealing temperature, see Additional file 8: Table S5) and extension (1 min at 72°C), final extension (10 min at 72°C); 2nd round of PCR - initial denaturation (15 min at 95°C), followed by 26 cycles of denaturation (30s at 94°C), annealing (30s at 60°C) and extension (1 min at 72°C), final extension (5 min at 72°C). All amplicons after 2nd round of PCR (up to 288 for one run) were purified using Agencourt AMPureXP magnetic beads (Beckman Coulter, Fullerton, CA, USA) and Biomek® FXP Laboratory Automation Workstation (Beckman Coulter). Precise concentration of amplicons were determined using Quant-iT™ PicoGreen dsDNA Assay Kit (Life Technologies, Carlsbad, CA, USA) and amplicons were equimolarly pooled to obtain amplicon library with 109 fragments/μl concentration. Subsequent procedures were carried out according to 454 amplicon sequencing manuals (454 Life Sciences, Roche Applied Science, Branford, CT, USA) on the GS Junior sequencer (Roche). An average overall coverage of 225 reads was observed (ranging from 53 to 659 for individual amplicons).
mRNA microarray profiling
Gene expression profiles from 14 AML patients and 4 CD34+ cells of healthy controls were generated by HumanHT-12 v4 Expression BeadChip Kit (Illumina). The chip scanning was done with a BeadStation 500 instrument (Illumina).
Quantitative real-time PCR (RQ-PCR)
The expression levels of selected genes were assessed with TaqMan Gene Expression Assays (Life Technologies) - see Additional file 9: Table S6 for individual Assay IDs. GAPDH was utilized as a house-keeping gene. Amplification was carried out with TaqMan Universal Master Mix II (Life Technologies) and recommended cycling conditions. SybrGreen RQ-PCR was performed using QuantiTect SYBR Green PCR Kit (Qiagen) and pre-designed KiCqStart® SYBR® Green primers (Sigma-Aldrich). Each sample was run in duplicates on a StepOne instrument (Life Technologies).
The presence of internal tandem duplication (ITD) in the juxtamembrane (JM) and tyrosine kinase 1 (TKD1) domains (exons 12-14) of FLT3 gene and the presence of CBFB-MYH11 fusion transcript at diagnosis was detected as described previously . Further, we examined mutations in NPM1, CEBPA and DNMT3A and intragenic MLL abnormalities such as partial tandem duplications (MLL-PTD) by direct sequencing ,.
For cytogenetic analyses and fluorescence in situ hybridization (FISH) the samples of bone marrow were cultivated for 24 hrs in medium RPMI 1640 with 10% of fetal calf serum. Twenty G-banded Wright-Giemsa stained mitoses, if available, were evaluated. The karyotypes were described following ISCN 2013 nomenclature. For precise identification of chromosomal aberrations, FISH with locus specific DNA probes (Vysis, Downers Grove, IL, USA) and multicolor FISH with color kit probes and ISIS computer analysis (both from Metasysteme, Altlusheim, Germany) were used.
Data processing and statistics
Data from targeted bisulfite sequencing were processed and evaluated using freely available programs: (i) FastQC  (quality control of reads), (ii) Trimmomatic  (removal of adapters/primer-dimers and bases with low-quality scores), (iii) Bismark  (methylation-aware alignment of reads to the reference genome and computation of methylation ratios). Differentially methylated target regions were assigned to genes using the on-line annotation tool GREAT . Quantile normalization and subtraction of background was applied to the raw microarray expression data in BeadStudio Data Analysis Software (Illumina). Raw data from 454 pyrosequencing were processed using a filter template to relax the stringency of the original valley filter (kindly provided by Dr. Esteban Czwan, Roche). This step was necessary due to the lower complexity of bisulfite treated DNA containing long stretches of homopolymers. The filter template is available on-line (Additional file 10) and its usage is described in Additional file 11. Data from 454 pyrosequencing were aligned to a reference in GS Amplicon Variant Analyzer (AVA) (Roche) software and DNA methylation levels were assessed using the web-based software BISMA .
Kaplan-Meier curves and two-sided log-rank test were used to estimate the overall survival and to compare differences between survival curves. The relations between qualitative parameters were compared in contingency tables using Fisher's exact test. For analyses of quantitative data, medians were detected and non-parametric two-tailed Mann-Whitney tests were performed. All these tests were conducted at a level of significance of 0.05 using GraphPad Prism4 software (GraphPad Software, San Diego, CA, USA). Cox regression analysis. was performed applying the SPSS statistical software (SPSS Inc., Chicago, IL, USA).
HH performed targeted and 454 bisulfite analysis and gene expression assays, analysed and interpreted data and wrote the manuscript; MHYF did the bioinformatics and statistical analyses and wrote the manuscript; CH interpreted data and wrote the manuscript; JS selected patients for 454 bisulfite sequencing, provided and interpreted clinical data and reviewed the manuscript; CS chose patients for targeted bisulfite analysis, provided and interpreted clinical data and reviewed the manuscript; JM performed molecular genetics analysis and reviewed the manuscript; ZK and MDM performed and evaluated the gene expression microarrays and reviewed the manuscript; AK, MV and OF performed molecular genetics analysis and reviewed the manuscript, KM provided cytogenetic data and reviewed the manuscript; PC provided clinical data and reviewed the article; VB supervised targeted bisulfite analysis and wrote the manuscript. All authors critically reviewed the manuscript and read and approved the final version.
We thank Dinko Pavlinic for help with targeted bisulfite library preparation and Esteban Czwan from Roche for help with raw 454 bisulfite data acquisition. This study was supported by the project (Ministry of Health, Czech Republic) for conceptual development of research organization (00023736, IHBT), by the grant ERDF OPPK CZ.2.16/3.1.00/24001 and was a part of the COST Action BM0801 (EuGESMA).
- Sonnet M, Claus R, Becker N, Zucknick M, Petersen J, Lipka DB, Oakes CC, Andrulis M, Lier A, Milsom MD, Witte T, Gu L, Kim-Wanner S, Schirmacher P, Wulfert M, Gattermann N, Luebbert M, Rosenbauer F, Rehli M, Bullinger L, Weichenhan D, Plass C: Early aberrant DNA methylation events in a mouse model of acute myeloid leukemia. Genome Med. 2014, 6: 34-10.1186/gm551.PubMed CentralView ArticlePubMedGoogle Scholar
- 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 A: DNA methylation signatures identify biologically distinct subtypes in acute myeloid leukemia. Cancer Cell. 2010, 17 (1): 13-27. 10.1016/j.ccr.2009.11.020.PubMed CentralView ArticlePubMedGoogle Scholar
- Genomic and epigenomic landscapes of adult de novo acute myeloid leukemia. N Engl J Med. 2013, 368 (22): 2059-2074. 10.1056/NEJMoa1301689.Google Scholar
- Akalin A, Garrett-Bakelman FE, Kormaksson M, Busuttil J, Zhang L, Khrebtukova I, Milne TA, Huang Y, Biswas D, Hess JL, Allis CD, Roeder RG, Valk PJM, Lowenberg B, Delwel R, Fernandez HF, Paietta E, Tallman MS, Schroth GP, Mason CE, Melnick A, Figueroa ME: Base-pair resolution DNA methylation sequencing reveals profoundly divergent epigenetic landscapes in acute myeloid leukemia. PLoS Genet. 2012, 8 (6): e1002781-10.1371/journal.pgen.1002781.PubMed CentralView ArticlePubMedGoogle Scholar
- Jones PA: Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet. 2012, 13 (7): 484-492. 10.1038/nrg3230.View ArticlePubMedGoogle Scholar
- Ivanov M, Kals M, Kacevska M, Metspalu A, Ingelman-Sundberg M, Milani L: In-solution hybrid capture of bisulfite-converted DNA for targeted bisulfite sequencing of 174 ADME genes. Nucleic Acids Res. 2013, 41 (6): e72-10.1093/nar/gks1467.PubMed CentralView ArticlePubMedGoogle Scholar
- Sproul D, Kitchen RR, Nestor CE, Dixon JM, Sims AH, Harrison DJ, Ramsahoye BH, Meehan RR: Tissue of origin determines cancer-associated CpG island promoter hypermethylation patterns. Genome Biol. 2012, 13 (10): R84-10.1186/gb-2012-13-10-r84.PubMed CentralView ArticlePubMedGoogle Scholar
- Valk P, Verhaak R, Beijen M, Erpelinck C, van Doorn-Khosrovani S, Boer J, Beverloo H, Moorhouse M, van der Spek P, Lowenberg B, Delwel R: Prognostically useful gene-expression profiles in acute myeloid leukemia. N Engl J Med. 2004, 350 (16): 1617-1628. 10.1056/NEJMoa040465.View ArticlePubMedGoogle Scholar
- Li Z, Zhang Z, Li Y, Arnovitz S, Chen P, Huang H, Jiang X, Hong G, Kunjamma RB, Ren H, He C, Wang C, Elkahloun AG, Valk PJM, Doehner K, Neilly MB, Bullinger L, Delwel R, Lowenberg B, Liu PP, Morgan R, Rowley JD, Yuan C, Chen J: PBX3 is an important cofactor of HOXA9 in leukemogenesis. Blood. 2013, 121 (8): 1422-1431. 10.1182/blood-2012-07-442004.PubMed CentralView ArticlePubMedGoogle Scholar
- Suzuki R, Shimodaira H: Pvclust: an R package for assessing the uncertainty in hierarchical clustering. Bioinformatics. 2006, 22 (12): 1540-1542. 10.1093/bioinformatics/btl117.View ArticlePubMedGoogle Scholar
- McLean CY, Bristor D, Hiller M, Clarke SL, Schaar BT, Lowe CB, Wenger AM, Bejerano G: GREAT improves functional interpretation of cis-regulatory regions. Nat Biotechnol. 2010, 28 (5): 495-U155. 10.1038/nbt.1630.View ArticlePubMedGoogle Scholar
- Dickson GJ, Liberante FG, Kettyle LM, O'Hagan KA, Finnegan DPJ, Bullinger L, Geerts D, McMullin MF, Lappin TRJ, Mills KI, Thompson A: HOXA/PBX3 knockdown impairs growth and sensitizes cytogenetically normal acute myeloid leukemia cells to chemotherapy. Haematologica. 2013, 98 (8): 1216-1225. 10.3324/haematol.2012.079012.PubMed CentralView ArticlePubMedGoogle Scholar
- de la Bletiere DR, Blanchet O, Cornillet-Lefebvre P, Coutolleau A, Baranger L, Genevieve F, Luquet I, Hunault-Berger M, Beucher A, Schmidt-Tanguy A, Zandecki M, Delneste Y, Ifrah N, Guardiola P: Routine use of microarray-based gene expression profiling to identify patients with low cytogenetic risk acute myeloid leukemia: accurate results can be obtained even with suboptimal samples. BMC Med Genet. 2012, 5: 6-Google Scholar
- Mandoli A, Singh AA, Jansen PWTC, Wierenga ATJ, Riahi H, Franci G, Prange K, Saeed S, Vellenga E, Vermeulen M, Stunnenberg HG, Martens JHA: CBFB-MYH11/RUNX1 together with a compendium of hematopoietic regulators, chromatin modifiers and basal transcription factors occupies self-renewal genes in inv(16) acute myeloid leukemia. Leukemia. 2014, 28 (4): 770-778. 10.1038/leu.2013.257.View ArticlePubMedGoogle Scholar
- Heuser M, Beutel G, Krauter J, Doehner K, von Neuhoff N, Schlegelberger B, Ganser A: High meningioma 1 (MN1) expression as a predictor for poor outcome in acute myeloid leukemia with normal cytogenetics. Blood. 2006, 108 (12): 3898-3905. 10.1182/blood-2006-04-014845.View ArticlePubMedGoogle Scholar
- Langer C, Marcucci G, Holland KB, Radmacher MD, Maharry K, Paschka P, Whitman SP, Mrozek K, Baldus CD, Vij R, Powell BL, Carroll AJ, Kolitz JE, Caligiuri MA, Larson RA, Bloomfield CD: Prognostic importance of MN1 transcript levels, and biologic insights from MN1-associated gene and MicroRNA expression signatures in cytogenetically normal acute myeloid leukemia: a cancer and leukemia group B study. J Clin Oncol. 2009, 27 (19): 3198-3204. 10.1200/JCO.2008.20.6110.PubMed CentralView ArticlePubMedGoogle Scholar
- Grimwade D, Walker H, Oliver F, Wheatley K, Harrison C, Harrison G, Rees J, Hann I, Stevens R, Burnett A, Goldstone A: The importance of diagnostic cytogenetics on outcome in AML: analysis of 1,612 patients entered into the MRC AML 10 trial. Blood. 1998, 92 (7): 2322-2333.PubMedGoogle Scholar
- Mrozek K, Bloomfield CD: Clinical significance of the most common chromosome translocations in adult acute myeloid leukemia. JNCI Monogr. 2008, 2008 (39): 52-57. 10.1093/jncimonographs/lgn003.View ArticleGoogle Scholar
- Carella C, Bonten J, Sirma S, Kranenburg TA, Terranova S, Klein-Geltink R, Shurtleff S, Downing JR, Zwarthoff EC, Liu PP, Grosveld GC: MN1 overexpression is an important step in the development of inv(16) AML. Leukemia. 2007, 21 (8): 1679-1690. 10.1038/sj.leu.2404778.View ArticlePubMedGoogle Scholar
- Chang C, Brocchieri L, Shen W, Largman C, Cleary M: Pbx modulation of Hox homeodomain amino-terminal arms establishes different DNA-binding specificities across the Hox locus. Mol Cell Biol. 1996, 16 (4): 1734-1745.PubMed CentralView ArticlePubMedGoogle Scholar
- Itonaga H, Imanishi D, Wong Y, Sato S, Ando K, Sawayama Y, Sasaki D, Tsuruda K, Hasegawa H, Imaizumi Y, Taguchi J, Tsushima H, Yoshida S, Fukushima T, Hata T, Moriuchi Y, Yanagihara K, Miyazaki Y: Expression of myeloperoxidase in acute myeloid leukemia blasts mirrors the distinct DNA methylation pattern involving the downregulation of DNA methyltransferase DNMT3B. Leukemia. 2014, 28 (7): 1459-1466. 10.1038/leu.2014.15.View ArticlePubMedGoogle Scholar
- Schneider M, Szaumkessel M, Richter J, Ammerpohl O, Hansmann ML, Küppers R, Siebert R, Giefing M: The PRDX2 gene is transcriptionally silenced and de novo methylated in Hodgkin and Reed-Sternberg cells of classical Hodgkin lymphoma. Blood. 2014, 123 (23): 3672-3674. 10.1182/blood-2014-02-553263.View ArticlePubMedGoogle Scholar
- Jiang D, Hong Q, Shen Y, Xu Y, Zhu H, Li Y, Xu C, Ouyang G, Duan S: The diagnostic value of DNA methylation in leukemia: a systematic review and meta-analysis. PLoS One. 2014, 9 (5): e96822-10.1371/journal.pone.0096822.PubMed CentralView ArticlePubMedGoogle Scholar
- Markova J, Trnkova Z, Michkova P, Maaloufova J, Stary J, Cetkovsky P, Schwarz J: Monitoring of minimal residual disease in patients with core binding factor acute myeloid leukemia and the impact of C-KIT, FLT3, and JAK2 mutations on clinical outcome. Leuk Lymphoma. 2009, 50 (9): 1448-1460. 10.1080/10428190903085951.View ArticlePubMedGoogle Scholar
- Pitiot AS, Santamaria I, Garcia-Suarez O, Centeno I, Astudillo A, Rayon C, Balbin M: A new type of NPM1 gene mutation in AML leading to a C-terminal truncated protein. Leukemia. 2007, 21 (7): 1564-1566. 10.1038/sj.leu.2404679.View ArticlePubMedGoogle Scholar
- Fuchs O, Provamikova D, Kocova M, Kostecka A, Cvekova P, Neuwirtova R, Kobylka P, Cermak J, Brozinova J, Schwarz J, Markova J, Salaj P, Klamova H, Maaloufova J, Lemez P, Novakova L, Benesova K: CEBPA polymorphisms and mutations in patients with acute myeloid leukemia, myelodysplastic syndrome, multiple myeloma and non-Hodgkin's lymphoma. Blood Cells Mol Dis. 2008, 40 (3): 401-405. 10.1016/j.bcmd.2007.11.005.View ArticlePubMedGoogle Scholar
- Markova J, Michkova P, Burckova K, Brezinova J, Michalova K, Dohnalova A, Maaloufova JS, Soukup P, Vitek A, Cetkovsky P, Schwarz J: Prognostic impact of DNMT3A mutations in patients with intermediate cytogenetic risk profile acute myeloid leukemia. Eur J Haematol. 2012, 88 (2): 128-135. 10.1111/j.1600-0609.2011.01716.x.View ArticlePubMedGoogle Scholar
- Caligiuri M, Strout M, Schichman S, Mrozek K, Arthur D, Herzig G, Baer M, Schiffer C, Heinonen K, Knuutila S, Nousiainen T, Ruutu T, Block A, Schulman P, PedersenBjergaard J, Croce C, Bloomfield C: Partial tandem duplication of ALL1 as a recurrent molecular defect in acute myeloid leukemia with trisomy 11. Cancer Res. 1996, 56 (6): 1418-1425.PubMedGoogle Scholar
- Shiah H, Kuo Y, Tang J, Huang S, Yao M, Tsay W, Chen Y, Wang C, Shen M, Lin D, Lin K, Tien H: Clinical and biological implications of partial tandem duplication of the MLL gene in acute myeloid leukemia without chromosomal abnormalities at 11q23. Leukemia. 2002, 16 (2): 196-202. 10.1038/sj.leu.2402352.View ArticlePubMedGoogle Scholar
- Krueger F, Kreck B, Franke A, Andrews SR: DNA methylome analysis using short bisulfite sequencing data. Nat Methods. 2012, 9 (2): 145-151. 10.1038/nmeth.1828.View ArticlePubMedGoogle Scholar
- Lohse M, Bolger AM, Nagel A, Fernie AR, Lunn JE, Stitt M, Usadel B: RobiNA: a user-friendly, integrated software solution for RNA-Seq-based transcriptomics. Nucleic Acids Res. 2012, 40 (W1): W622-W627. 10.1093/nar/gks540.PubMed CentralView ArticlePubMedGoogle Scholar
- Krueger F, Andrews SR: Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics. 2011, 27 (11): 1571-1572. 10.1093/bioinformatics/btr167.PubMed CentralView ArticlePubMedGoogle Scholar
- Rohde C, Zhang Y, Reinhardt R, Jeltsch A: BISMA - Fast and accurate bisulfite sequencing data analysis of individual clones from unique and repetitive sequences. BMC Bioinformatics. 2010, 11: 230-10.1186/1471-2105-11-230.PubMed CentralView ArticlePubMedGoogle Scholar
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/4.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.