CREBBP is a target of epigenetic, but not genetic, modification in juvenile myelomonocytic leukemia
- Silvia Fluhr1, 2Email author,
- Melanie Boerries3, 4, 5,
- Hauke Busch3, 4, 5,
- Aikaterini Symeonidi3,
- Tania Witte6,
- Daniel B Lipka6,
- Oliver Mücke6,
- Peter Nöllke1,
- Christopher Felix Krombholz1,
- Charlotte M Niemeyer1, 5,
- Christoph Plass†5, 6 and
- Christian Flotho†1, 5
© Fluhr et al. 2016
Received: 16 February 2016
Accepted: 27 April 2016
Published: 5 May 2016
Juvenile myelomonocytic leukemia (JMML) is a myeloproliferative neoplasm of childhood whose clinical heterogeneity is only poorly represented by gene sequence alterations. It was previously shown that aberrant DNA methylation of distinct target genes defines a more aggressive variant of JMML, but only few significant targets are known so far. To get a broader picture of disturbed CpG methylation patterns in JMML, we carried out a methylation screen of 34 candidate genes in 45 patients using quantitative mass spectrometry.
Five of 34 candidate genes analyzed showed recurrent hypermethylation in JMML. cAMP-responsive element-binding protein-binding protein (CREBBP) was the most frequent target of epigenetic modification (77 % of cases). However, no pathogenic mutations of CREBBP were identified in a genetic analysis of 64 patients. CREBBP hypermethylation correlated with clinical parameters known to predict poor outcome.
This study supports the relevance of epigenetic aberrations in JMML pathophysiology. Our data confirm that DNA hypermethylation in JMML is highly target-specific and associated with higher-risk features. These findings encourage the development of prognostic markers based on epigenetic alterations, which will be helpful in the difficult clinical management of this heterogeneous disease.
KeywordsJuvenile myelomonocytic leukemia DNA methylation Epigenetics CREBBP
Findings and conclusions
Risk assessment in juvenile myelomonocytic leukemia (JMML) is important for clinical decision-making as this aggressive myeloproliferative neoplasm is generally unresponsive to conventional chemotherapy and most (but not all) patients need early allotransplantation [1, 2]. Several clinical parameters have been well-established to predict poor outcome , but the picture is less clear at the molecular level. Specifically, if and how the set of genetic mutations in leukemic cells contributes to clinical risk remains a matter of ongoing debate [2–4]. Recent research on epigenetic dysregulation in JMML has provided a growing body of evidence that high-risk cases of JMML are characterized by specific CpG island hypermethylation which integrates the clinical and genetic risk factors identified so far [5–7]. To provide a broader picture of epimutations in JMML, we carried out a DNA methylation analysis of 34 candidate genes in leukemic cells from 45 children with JMML.
Aberrant methylation of CREBBP, MPO, SLC12A8, HIC1, and TCF4 in JMML
We defined a gene as hypermethylated in a JMML sample if the average CpG methylation in the region of analysis exceeded three standard deviations above the mean observed in granulocytes of 11 healthy control subjects. Using this definition, 34 of 44 JMML cases (77 %) showed hypermethylation of the cAMP-responsive element-binding protein-binding protein (CREBBP) gene (the assay was uninformative in one case), 20 of 45 (44 %) were hypermethylated in the MPO gene, and 19 of 45 (42 %) carried a hypermethylated SLC12A8 gene (Fig. 1b). Most CpG units analyzed in the HIC1 and TCF4 genes had normal methylation in JMML, but each region contained a single CpG unit near its border with significant hypermethylation in JMML (Fig. 1c). This observation may hint at aberrant methylation in neighboring regions not covered by the assay. A group-wise comparison of methylation levels using Mann-Whitney test with Bonferroni correction showed for each gene a highly significant difference between the “hypermethylated JMML” group and healthy controls but no difference between “normal methylation JMML” and controls (Fig. 1b, c). The DNA methylation level of the remaining 29 genes was similar between JMML samples and healthy granulocytes or the difference was so small that no biological significance would be expected (Additional file 1: Table S1).
The analysis demonstrates that CpG hypermethylation in JMML is highly locus-specific. Although a significant proportion of JMML cases are characterized by aberrant DNA methylation, the hypermethylation events are confined to a limited array of genetic regions. This observation fits well with three previous studies where 11 of 15 candidate regions did not show any aberrant methylation in 127 cases analyzed , or where hypermethylation was restricted to only one of two CpG islands of the RASA4 locus  or one of three CpG islands of the AKAP12 locus . Contrary to the two-hit hypothesis formulated above, none of the nine 7q22 genes exhibited frequent hypermethylation in JMML, whether monosomy 7 was present or not.
High-resolution DNA methylation analysis of the CREBBP CpG island in JMML
Somatic CREBBP mutations are uncommon in JMML
We next studied if the CREBBP locus was not only subject to epigenetic modification in JMML, but also to genetic mutation. We analyzed leukemic granulocyte DNA from 64 children with JMML (17 of which overlapped with the 44-patient series used for epigenetic analysis) for CREBBP exon sequence variants using the Agilent SureSelect v4.0 capture technique and deep-sequencing, resulting in an average CREBBP exon depth of 82 reads. Eighty percent of CREBBP exons were covered by at least 50 reads. No insertions or deletions in the CREBBP gene were found in any JMML sample. After removing common non-pathogenic single-nucleotide polymorphisms documented in public databases such as dbSNP, we identified 10 non-synonymous sequence variants in 11 cases of JMML (Additional file 1: Table S2). All except one (c.A493G) were listed in at least one of four exome variation resources (Exome Variant Server, ExAC browser, 1000 Genomes, dbSNP) with minor allele frequencies of less than 0.1 %. Three variants (c.A5933G, c.A2728G, c.G2941A) corresponded to alterations associated with Rubinstein-Taybi syndrome [15, 16]. However, indicative clinical features (facial dysmorphism, thumb or toe malformation, intellectual disability) were documented in none of the 11 children. Importantly, we did not detect any of the CREBBP histone acetyltransferase domain mutations described as pathogenic somatic events in lymphoblastic leukemia or lymphoma [17–20]. DNA from paired non-hematopoietic tissue was available in 5 of 11 JMML cases with non-synonymous CREBBP sequence variants, enabling us to test the germline status. In all cases, the variants were present in non-leukemic tissue. The pathogenic significance of these germline CREBBP variants in JMML patients without syndromic features remains unclear. Taken together, the data indicate that the coding sequence of CREBBP is not a relevant target of somatic alteration in JMML, in agreement with the absence of CREBBP mutations from myeloproliferative neoplasms in adults .
CREBBP methylation is associated with established hematologic and clinical parameters indicating an aggressive JMML phenotype
It was repeatedly shown that DNA hypermethylation at specific target genes in JMML is associated with clinical parameters indicating poor prognosis [5–7, 10]. To test if this relationship was also true for CREBBP methylation, we focused on CpG sites #8 to #10 because these sites were unmethylated in healthy individuals but had high methylation variation in JMML. The level of methylation at these sites was significantly associated with older age and elevated fetal hemoglobin at diagnosis (Fig. 2d and Additional file 1: Table S3), two main predictors of reduced survival [1, 2]. In addition, methylation at CpG sites #8 to #10 was associated with the category of Ras pathway mutation (Additional file 1: Table S3), being considerably higher in JMML cases with PTPN11 mutation or neurofibromatosis type 1 (NF1) than in JMML with KRAS, NRAS, or CBL mutation. The relevance of the mutation category to disease course and outcome of JMML remains controversial [2–4, 22], but higher CREBBP methylation in JMML with PTPN11 mutation or NF1 is consistent with the prevailing opinion that these subgroups are more aggressive. Variable correlation with age, high fetal hemoglobin, mutation, or karyotype but no other parameters was also found for MPO, SLC12A8, and HIC1 CpG 3359. The small group size and heterogeneous treatment precluded meaningful outcome analyses. However, we noted that five events of death occurred in the CREBBP hypermethylation group (5/34 patients, 15 %) compared to no such event in the group with normal CREBBP methylation (0/10 patients); the median follow-up time was 6.0 years for all 44 patients. In summary, the data strengthen the concept that higher-risk cases of JMML carry a distinct DNA methylation phenotype. Future studies are needed to elucidate the functional role of coordinated DNA hypermethylation at specific loci in the JMML-initiating cell.
We used clinical material from German JMML patients registered in the European Working Group of MDS in Childhood (EWOG-MDS) studies “98” and “2006” after obtaining informed consent and ethics approval (University of Freiburg, reference number EK247/05). Patients with Noonan syndrome (PTPN11 or KRAS germline mutation) were excluded. Clinical characteristics are given in Additional file 1: Table S3.
High-resolution quantitative DNA methylation analysis
Leukemic granulocytes from 45 patients (bone marrow, N = 36; peripheral blood, N = 9) were enriched by Ficoll centrifugation (Biochrom). Genomic DNA was isolated using the Gentra Puregene kit (Qiagen) and bisulfite converted using the EZ DNA methylation kit (Zymo Research). The MassARRAY EpiTYPER assay (Agena Bioscience) was used to assess CpG island methylation at 34 target regions (Fig. 1a and Additional file 1: Table S1). Primer sequences are listed in Additional file 1: Table S4.
CREBBP sequence analysis
Library construction and exon region capture of fragmented genomic DNA isolated from leukemic granulocytes of 64 patients were done using TruSeq (Illumina) and SureSelect v4.0 (Agilent). 100-base pair paired-end sequencing was performed on HiSeq2000 instruments (Illumina). The reads were aligned to the hg19 reference genome. CREBBP single-nucleotide polymorphisms and indels were identified according to standard recommendations (Genome Analysis Toolkit, Broad Institute).
juvenile myelomonocytic leukemia
cAMP-responsive element-binding protein-binding protein
neurofibromatosis type 1
We are grateful to A. Meier for technical assistance. Grant support: German Research Foundation (CRC 992-C05 to CF, SPP1463 FL345/4-2 to SF and CF, CRC 850 to AS).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Locatelli F, Nöllke P, Zecca M, Korthof E, Lanino E, Peters C, et al. Hematopoietic stem cell transplantation (HSCT) in children with juvenile myelomonocytic leukemia (JMML): results of the EWOG-MDS/EBMT trial. Blood. 2005;105(1):410–9.View ArticlePubMedGoogle Scholar
- Locatelli F, Niemeyer CM. How I treat juvenile myelomonocytic leukemia. Blood. 2015;125(7):1083–90.View ArticlePubMedGoogle Scholar
- Flotho C, Kratz CP, Bergsträsser E, Hasle H, Starý J, Trebo M, et al. Genotype-phenotype correlation in cases of juvenile myelomonocytic leukemia with clonal RAS mutations. Blood. 2008;111(2):966–7.View ArticlePubMedGoogle Scholar
- Niemeyer CM, Strahm B, Dworzak M, de Moerloose B, Hasle H, Stary J, et al. JMML revisited: role and outcome of hematopoietic stem cell transplantation in subtypes of juvenile myelomonocytic leukemia. Blood. 2012;120(21):955.Google Scholar
- Olk-Batz C, Poetsch AR, Nöllke P, Claus R, Zucknick M, Sandrock I, et al. Aberrant DNA methylation characterizes juvenile myelomonocytic leukemia (JMML) with poor outcome. Blood. 2011;117(18):4871–80.View ArticlePubMedGoogle Scholar
- Poetsch AR, Lipka DB, Witte T, Claus R, Nöllke P, Zucknick M, et al. RASA4 undergoes DNA hypermethylation in resistant juvenile myelomonocytic leukemia. Epigenetics. 2014;9(9):1252–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Sakaguchi H, Muramatsu H, Okuno Y, Makishima H, Xu Y, Furukawa-Hibi Y, et al. Aberrant DNA methylation is associated with a poor outcome in juvenile myelomonocytic leukemia. PLoS One. 2015;10(12):e0145394.View ArticlePubMedPubMed CentralGoogle Scholar
- Claus R, Pfeifer D, Almstedt M, Zucknick M, Hackanson B, Plass C, et al. Decitabine induces very early in vivo DNA methylation changes in blasts from patients with acute myeloid leukemia. Leuk Res. 2013;37(2):190–6.View ArticlePubMedGoogle Scholar
- Kratz CP, Emerling BM, Donovan S, Laig-Webster M, Taylor BR, Thompson P, et al. Candidate gene isolation and comparative analysis of a commonly deleted segment of 7q22 implicated in myeloid malignancies. Genomics. 2001;77(3):171–80.View ArticlePubMedGoogle Scholar
- Wilhelm T, Lipka DB, Witte T, Wierzbinska JA, Fluhr S, Helf M, et al. Epigenetic silencing of AKAP12 in juvenile myelomonocytic leukemia. Epigenetics. 2016; in press.
- Kung AL, Rebel VI, Bronson RT, Ch’ng LE, Sieff CA, Livingston DM, et al. Gene dose-dependent control of hematopoiesis and hematologic tumor suppression by CBP. Genes Dev. 2000;14(3):272–7.PubMedPubMed CentralGoogle Scholar
- Zimmer SN, Lemieux ME, Karia BP, Day C, Zhou T, Zhou Q, et al. Mice heterozygous for CREB binding protein are hypersensitive to gamma-radiation and invariably develop myelodysplastic/myeloproliferative neoplasm. Exp Hematol. 2012;40(4):295–306.View ArticlePubMedPubMed CentralGoogle Scholar
- Sanchez-Molina S, Estaras C, Oliva JL, Akizu N, Asensio-Juan E, Rojas JM, et al. Regulation of CBP and Tip60 coordinates histone acetylation at local and global levels during Ras-induced transformation. Carcinogenesis. 2014;35(10):2194–202.View ArticlePubMedGoogle Scholar
- Tsunoda T, Takagi T. Estimating transcription factor bindability on DNA. Bioinformatics. 1999;15(7–8):622–30.View ArticlePubMedGoogle Scholar
- Coupry I, Roudaut C, Stef M, Delrue MA, Marche M, Burgelin I, et al. Molecular analysis of the CBP gene in 60 patients with Rubinstein-Taybi syndrome. J Med Genet. 2002;39(6):415–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Bartsch O, Kress W, Kempf O, Lechno S, Haaf T, Zechner U. Inheritance and variable expression in Rubinstein-Taybi syndrome. Am J Med Genet A. 2010;152A(9):2254–61.View ArticlePubMedGoogle Scholar
- Mullighan CG, Zhang J, Kasper LH, Lerach S, Payne-Turner D, Phillips LA, et al. CREBBP mutations in relapsed acute lymphoblastic leukaemia. Nature. 2011;471(7337):235–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Pasqualucci L, Dominguez-Sola D, Chiarenza A, Fabbri G, Grunn A, Trifonov V, et al. Inactivating mutations of acetyltransferase genes in B-cell lymphoma. Nature. 2011;471(7337):189–95.View ArticlePubMedPubMed CentralGoogle Scholar
- Huether R, Dong L, Chen X, Wu G, Parker M, Wei L, et al. The landscape of somatic mutations in epigenetic regulators across 1,000 paediatric cancer genomes. Nat Commun. 2014;5:3630.View ArticlePubMedPubMed CentralGoogle Scholar
- Nielsen JS, Sedgwick C, Shahid A, Zong Z, Brumme ZL, Yu S, et al. Toward personalized lymphoma immunotherapy: identification of common driver mutations recognized by patient CD8+ T cells. Clin Cancer Res. 2015. [Epub ahead of print]
- Andersen CL, Hasselbalch H, Gronbaek K. Lack of somatic mutations in the catalytic domains of CREBBP and EP300 genes implies a role for histone deacetylase inhibition in myeloproliferative neoplasms. Leuk Res. 2012;36(4):485–7.View ArticlePubMedGoogle Scholar
- Yoshida N, Yagasaki H, Xu Y, Matsuda K, Yoshimi A, Takahashi Y, et al. Correlation of clinical features with the mutational status of GM-CSF signaling pathway-related genes in juvenile myelomonocytic leukemia. Pediatr Res. 2009;65(3):334–40.View ArticlePubMedGoogle Scholar