Dissecting the role of aberrant DNA methylation in human leukemia

Chronic Myeloid Leukemia (CML) is a myeloproliferative disorder characterized by the genetic translocation t(9;22)(q34;q11.2) encoding for the BCR-ABL fusion oncogene. However, many molecular mechanisms of the disease progression still remain poorly understood. A growing body of evidence suggests that epigenetic abnormalities are involved in tyrosine kinase resistance in CML, leading to leukemic clone escape and disease propagation. Here we show that, by applying cellular reprogramming to primary CML cells, aberrant DNA methylation contributes to the disease evolution. Importantly, using a BCR-ABL inducible murine model, we demonstrate that a single oncogenic lesion triggers DNA methylation changes which in turn act as a precipitating event in leukemia progression.


Introduction
. Recent evidences suggest that epigenetic alterations such as DNA methylation can contribute to leukemia pathogenesis 4,5,6 . In particular, abnormal epigenetic regulation of specific genes might play an important role in both leukemia development and therapeutic responsiveness. CML is a particular type of leukemia characterized by the presence of the BCR-ABL oncogene. BCR-ABL is a constitutive active tyrosine kinase mediating deregulation of several pathways involved in proliferation and differentiation 7 . Although the expression of BCR-ABL is considered the primarily feature associated with CML onset, other molecular mechanisms contributing to CML development remain to be elucidated. Both hypermethylation of specific genomic loci and genome-wide hypomethylation have been found to correlate with cancer development. Specifically, hypermethylation of tumor suppressor genes has been found to play a crucial role in carcinogenesis by affecting normal cell growth 8 . Further, aberrant DNA methylation has been linked to the onset of leukemic clones resistant to tyrosine kinase inhibitors and deemed responsible for CML propagation and progression 9 . Among the genes found to be hypermethylated in CML and other lymphoid malignancies, and correlating with a poor outcome, are HOXA4 and HOXA5, both acting as differentiating genes during hematopoiesis 10 . Interestingly, the observation that a member of CCAAT/enhancer-binding protein (CEBP) family of transcription factors: DNAdamage-inducible transcript 3 (DDIT3 gene, hypermethylated in CML 11 suggests that aberrant epigenetic regulation of the CEBPs genes may also play a role in leukemia onset 12 .
Taken together, these findings indicate that aberrant DNA methylation might contribute to leukemia development by synergizing with acquired genetic lesions. However, the lack of suitable animal models or other tools to directly address this particular question has constrained a thorough analysis.
Human induced pluripotent stem (iPS) cells can be generated from various cell types, including cancer cells, by overexpressing a defined set of transcription factors including OCT4, SOX2, KLF4 and c-MYC 13,14,15,16,17 . Previous studies have shown that cellular reprogramming might be able to erase tissue-specific DNA methylation and to re-establish an ES-like DNA methylation state 18 . Here, by means of cellular reprogramming followed by in vitro re-differentiation we tested the functional relevance of DNA aberrant methylome in CML development. Reprogramming of CML cells into an iPS-like state was able to erase the cancer-specific DNA methylation signature and to identify a cell population no longer effective in producing CML when subsequently transplanted into immunocompromised mice. Finally, using an inducible BCR-ABL transgenic mouse 19 , we demonstrate in vivo that a single genetic aberration perturbing DNA methylation profile acts as a crucial precipitating event in CML disease development.
positive CML patient. Leukemia cells were transduced as previously reported 13,20 . Two weeks after infection, colonies with typical human ES-like morphology were picked and expanded on mouse embryonic fibroblast feeder layers resulting in stable ES-like cell lines: LiPS1-K562 and LiPS2-K562, both derived from the K562 cell line, and CML-LiPS1 and CML-LiPS2, derived from primary CML cells. Additionally, we included the previously characterized KBM7 cell line and its reprogrammed counterpart in our analysis 15 . Remarkably, reprogrammed primary CML cells still preserved the BCR-ABL oncogene (Supplementary Figure 1C and 15 ). A comprehensive SNP array analysis confirmed that LiPS clones derived from K562 and KBM7 cell lines retained the same genetic alterations as the parental leukemia cells (Supplementary Figure 2 and Supplementary Data 1-2) ruling out the possibility that a essentially normal subclone or contaminating cell was selected during reprogramming.
Having established several LiPS cell lines, we proceeded to test whether cellular reprogramming was sufficient to reset DNA methylation of the parental leukemic cells. Genomic DNA methylation profiles of K562, KBM7, primary CML cells and of the respective LiPS clones were assessed by Reduced Representation Bisulfite Sequencing (RRBS), which has been shown to provide high sensitivity and specificity for detecting cancer-specific changes in DNA methylation not only in CpG islands but also throughout genes and in repetitive regions 21,22 . Compared to human ES cells, CD34 + -derived iPS cells (CD34 + -iPS) and CD34 + -cells, K562 and KBM7 cell lines exhibited widespread hypermethylation throughout the genome, including CpG islands, genes, and promoters ( Figure 1A) as well as across families of repetitive elements. Primary CML cells also demonstrated significant hypermethylation in CpG islands, gene, and promoter regions, although to a lesser extent than in the cell lines. Compared to pluripotent cells, primary CML cells demonstrated hypermethylation across families of repetitive elements, similarly to CD34 + cells ( Figure 1A). DNA methylation levels were substantially decreased throughout the genome in the LiPS cells as compared to their parental leukemic cell types, demonstrating the ability of cellular reprogramming to reset aberrant methylome ( Figure 1A). Overall, the LiPS clones acquired DNA methylation patterns that were similar but not identical to those of ES cell lines and CD34+iPS, supporting the notion of cell-of-origin specific 'epigenetic memory' in many types of induced pluripotent cells 23,24 . Indeed, through genome-scale DNA methylation analysis, we identified LiPS-specific DNA methylation changes of relevant pluripotent and hematopoietic-associated genes. Pluripotent genes such as SALL4 were reactivated in LiPS and exhibited decreased DNA methylation levels, particularly at the promoter regions consistent with the well-established role of DNA demethylation as a key event in cellular reprogramming ( Figure 1B and Supplementary Figure 3) 18 . Epigenetic reprogramming affected not only pluripotency-associated genes but also other genomic loci, i.e., hematopoietic transcription factors or tumor suppressor genes highly methylated across leukemia cells and responding to cellular reprogramming with different extent of LiPS demethylation (Fig. 1B).

Characterization of Differentially Methylated Regions in CML
In order to characterize the changes in DNA methylation that are induced by the reprogramming of CML cells to pluripotency, we used the RRBS data to bioinformatically identify differentially methylated regions (DMRs) between LiPS clones and the cells they were derived from (see Materials and Methods). Among the 500 strongest differentially methylated regions that are associated with gene promoters, the vast majority was hypermethylated in CML cells as compared to the derived LiPS cells (Figure 2A). This observation is consistent with the increased levels of promoter methylation that are common among cancer cells. To better understand the biological functions of the genes that are characterized by differential promoter methylation, we performed gene ontology enrichment analyses. The 500 gene promoters hypermethylated in CML cells as compared to LiPS cells were significantly enriched for functional categories associated with development, differentiation and signaling ( Figure 2B and Supplementary Table 1), suggesting that certain aspects of these stem cell-associated processes may be silenced by DNA methylation in CML cells. In contrast, the 500 gene promoters strongly hypermethylated in LiPS as compared to CML cells were enriched for hematopoietic regulatory gene categories such as lymphocyte activation and immune response ( Figure 2C and Supplementary Table 1), the latter recently shown to be attenuated in iPS cells [25][26] .
Beyond this gene promoter centric analysis of the potential biological functions of differential DNA methylation between CML cells and LiPS cells, we also investigated the co-localization of DNA methylation changes with other epigenetic and regulatory mechanisms throughout the genome. To that end, we repeated the differential DNA methylation analysis of 5kb tiling regions spanning the whole genome, and we used the EpiExplorer software 27 to compare the 500 most differentially methylated regions with the background of all tiling regions. This analysis revealed that DMRs are enriched for epigenetic marks associated with gene regulatory sequences including promoters and enhancers. Further, consistent with previous studies demonstrating a strong association between Polycomb binding in pluripotent cells and hypermethylation in cancer 28,29 , we observed that the CML hypermethylated regions were highly enriched for Polycomb binding sites ( Figure 2D).
Overall, these observations confirm that CML reprogrammed into LiPS cells undergo epigenome reprogramming that erases characteristic patterns of cancer-associated DNA methylation while establishing characteristic DNA methylation patterns of pluripotent cells.

Reprogramming restores myeloid differentiation of CML cells
To assess whether epigenetic reprogramming observed in LiPS cells coincides with phenotypic changes in their oncogenic potential, we applied in vitro and in vivo strategies using both primary human CML cells and the CML model cell line K562, chosen for its high efficiency of engraftment in xenotransplantation approaches. By co-culture with OP9 stroma cells 30 (Supplementary Figure 4A), we demonstrated that CML-LiPS derived CD45 + /CD34 + cells reacquire myeloid and erytroid differentiation capabilities giving rise to CD15 + /CD14 + /GlyA + cells, despite the presence of BCR-ABL oncogene (Supplementary Figure 4A-D). With a similar OP9 co-culture we proved that K562-LiPS derived CD45 + cells exhibit a well-differentiated monocyte/macrophage phenotype in contrast with the leukemia counterpart, locked into an undifferentiated cellular state ( Figure 3A and Supplementary Figure 4E). Indeed, while the parental leukemia and LiPS cells showed low expression of the hematopoietic transcription factor CEBPA, the purified CD45 + LiPS derived cells demonstrated reactivation of this critical differentiation-associated gene locus and increase of its levels of expression (Supplementary Figure 4F-G)). Remarkably, when K562-LiPS derived CD45 + cells and the parental leukemia cells (20,000 per mouse) were injected into immunocompromised NOD.Cg-Prkdc scid Il2 rgtm1Wjl /SzJ (NSG) mice, all mice (=3) transplanted with leukemic cells developed the disease and died within 30 days after the injection. In contrast, seven out of eight mice transplanted with LiPS-derived CD45 + cells, regardless the successful engraftment, did not show any sign associated with a leukemic phenotype, i.e. spleen infiltration as observed in case of K562, during the entire 80 days of monitoring period ( Figure 3B-C).
Finally, to investigate the potential contribution of DNA methylation in CML pathogenesis, we applied 5-Azacytidine (AZA) treatment to primary CML cells. Upon treatment with AZA, we observed a skew in differentiation towards the granulocyte and monocyte/ macrophage lineages, as shown by the presence of CD15 + and CD11b + cells, as compared to untreated cells ( Figure 3D-E) suggesting that DNA methylation abnormalities cooperate in disrupting hematopoietic differentiation during CML development and progression.

BCR-ABL induction reprograms DNA methylation pattern
Our observations prompted us to investigate the functional relationship between the presence or absence of BCR-ABL oncogene and DNA methylation response. To address this question, we took advantage of a transgenic murine model in which BCR-ABL expression can be conditionally induced and reversibly repressed 19 . Indeed previous studies, both in cell culture and in vivo, have shown that overexpression of specific oncogenes such as MLL-AF9 and mutant IDH1 are capable of inducing epigenetic reprogramming comparable to that observed in human disease 31, 32, 33 . In our model, BCR-ABL is under the control of the 3′ enhancer of the murine stem cell leukemia (SCL) gene, thus allowing BCR-ABL expression mainly in the hematopoietic stem sell (HSC) compartment and leading to development of a chronic myeloproliferative disorder similar to human CML 19 . Induction of BCR-ABL expression by withdrawal of tetracycline from the drinking water leads to significant increase of Gr-1 and Mac-1 positive myeloid population (leukemic mice). This phenotype can be reversed when tetracycline is added to the water (rescued mice), thus confirming a direct correlation between BCR-ABL presence and the expansion of that population ( Figure 4A and 19 ). Using this mouse model, we analyzed DNA methylation patterns in the HSCs population harvested from 3 distinct groups of mice: i) the non-induced BCR-ABL mice (control); ii) the induced BCR-ABL mice (leukemic) and iii) the repressed BCR-ABL mice (rescued). Leukemic mice showed moderate increase of CpG islands DNA methylation levels, indicating that a single oncogenic fusion protein can trigger characteristic perturbations of DNA methylation in a highly defined mouse model ( Figure 4B). Upon BCR-ABL activation, the leukemic mice acquired aberrant DNA methylation at several gene loci, one example being the hematopoietic transcription factor Hoxb1 (Figure 4C-D and Supplementary Table 2). Remarkably, repression of BCR-ABL led to widespread reversion of the BCR-ABL-induced changes in DNA methylation, supporting the hypothesis that BCR-ABL triggers DNA methylation changes. This observation could be explained by transcription factors and epigenetic regulator proteins directly downstream of BCR-ABL influencing DNA methylation, but it could also be a reflection of a broader change in epigenetic cell states when BCR-ABL transforms cells.
In summary, these data demonstrate the functional effect of a single genetic lesion in determining the aberrant methylome state, suggesting that epigenetic alterations act as a secondary precipitating event during leukemia development.

5-Azacytidine reduces BCR-ABL cells oncogenic potential
To test whether inhibition of BCR-ABL-induced aberrant DNA methylation could restore normal hematopoiesis during leukemogenesis, six leukemic mice per group were treated for 4 weeks with AZA (2 mg/kg) or with the tyrosine kinase inhibitor Imatinib (200 mg/kg) (IM, first line treatment for CML), respectively. Mice treated with AZA exhibited no expansion of the Mac1 + Gr1 + myeloid population as compared to the untreated counterpart, closely mimicking the results observed for the IM-treated group used as a positive control. Furthermore, both AZA and IM treated mice displayed normal lymphopoiesis, absent in the untreated controls (Supplementary Figure 5A). These data indicate that despite their different mechanisms of action, both AZA, at the epigenetic level by correcting aberrant DNA methylation, and IM at the genetic level by blocking the BCR-ABL oncogene, are capable of re-establishing normal hematopoiesis.
To further evaluate the contribution of aberrant DNA methylation to leukemia development in vivo, 44 irradiated B6.SJL-Ptprc a Pepc b /BoyJ (pep boys) congenic recipient mice were transplanted with murine BCR-ABL leukemic bone marrow cells. Eight to ten mice were randomly assigned to five groups: control (not leukemic); leukemic untreated, treated with AZA (2 mg/kg); or IM (200 mg/kg) or combination of the two drugs IM and AZA. The treatment was administered for 4 weeks (Supplementary Figure 5B) and mice were followed up for a period of up to 150 days. The untreated mice developed leukemia (as expected) with a median survival of 41 days. Strikingly, the AZA treatment provided longer survival (median survival 100 days; p= 0.0008) than the IM treatment (median survival: 56 days; p= 0.03) as compared to the untreated mice. Unfortunately, the combination of AZA and IM, at their standard dosages, resulted in hematological and gastrointestinal toxicity for the mice that did not allow us to make any conclusive observation on the potential benefits of using combination of the drugs to correct the leukemic phenotype. Future studies will be necessary to assess optimal dosage range and scheduling mode to: 1) reduce the toxicity and 2) understand the potential synergic effect of the two drugs ( Figure 5B). Moreover, the AZA treatment efficiently targeted the stem cell compartment (LSK population from the CD45.2+ donor cells), restoring expression and promoting demethylation of loci differentially methylated between control and sick mice (Supplementary Data 2), such as the Growth differentiation factor 3 (Gdf3) (Supplementary Figure 6), a promising target of the BCR-ABL oncogene.
Taken together, these data indicate that DNA methylation aberrations play a critical role in CML development, supporting the relevance of DNA methylation modulating drugs as therapeutic options.

Discussion
This study analyzed the interconnection and functional relevance of genetic and epigenetic alterations in CML. By deriving genetically matched but epigenetically distinct LiPS cell lines from CML cell lines and human primary CML cells, and by using an inducible mouse model of CML, we show that: 1) genetic and epigenetic alterations are both required to maintain the leukemic potential, and 2) that the BCR-ABL fusion protein is able to trigger DNA methylation changes that can act as a secondary event and contribute to leukemia formation. Using nuclear reprogramming as an "epigenetic tool" we were able to erase aberrant DNA methylation from leukemia cells, thereby delaying the onset of the malignancy. Among the leukemia-specific DMRs, less than half were associated with promoters, whereas the vast majority corresponded to specific histone marks such as enhancer-associated H3K4me1 and promoter-associated H3K4me3. Interestingly, genes hypermethylated in LiPS compared to CML cells belonged to functional categories of lymphocyte activation and immune response. It has been shown that iPS cells are defective in their antiviral response as a result of downregulation of key genes involved in immune response 25,26 . Similarly, attenuated antiviral response reported in human embryonic stem (ES) cells is restored by differentiation, suggesting that the impaired antiviral response of iPS cells likely stemmed from the pluripotent state and might promote cell pluripotency 25,26,34 .
Our data show that the mice transplanted with CD45 + cells derived from LiPS do not develop leukemia within the 80 days monitoring time ( Figure 3B). In particular, we demonstrated that leukemia progression was susceptible not only to the genetic inhibition of the driving oncogene but also to the inhibition of the secondary DNA methylation changes, which provides strong evidence for the cooperative role of both genetic and epigenetic alterations in leukemia.
In a related study focusing on glioblastoma, Stricker and colleagues 35 recently demonstrated that reprogramming of glioblastoma cells was also able to erase aberrant DNA methylation. But in contrast to our observations for CML, these induced epigenetic changes were insufficient to compromise oncogenic potential, and neural progenitors derived from the reprogrammed glioblastoma remained highly malignant. This important difference between glioblastoma and CML supports the relevance of aberrant DNA methylation specifically in CML, consistent with the much lower rates of genetic alterations in leukemia compared to glioblastoma 36 and with the history of successful use of DNA methylation inhibitors in leukemias. This conclusion is further supported by comparing our results in CML with a recent study focusing on cells derived from Ewing sarcoma patients, which is a solid tumor driven by an oncogenic fusion protein EWS-FLI1. Upon knock-down, Tomazou et al. observed widespread changes in histone acetylation and histone methylation, whereas DNA methylation patterns were largely unaffected 37 . Consistent with this observation, AZA had a much weaker effect in EWS-FLI1 expressing cells than HDAC inhibitors, whereas the functional impact of AZA on BCR-ABL expressing cells in our mouse model was striking and even exceeded those observed for Imatinib.
Overall, these results highlight the specific relevance of aberrant DNA methylation in CML, which are driven by BCR-ABL expression and are at the same time essential for maintaining a strong leukemic phenotype in a mouse model. This study provides a rationale for broadly exploring the use of demethylating agents as adjuvant treatments in leukemia, even in subtypes with recurrent and well-characterized genetic abnormalities such as BCR-ABL.

Cell cultures and LiPS cells generation
Human K562 leukemia cells (ATCC ® CCL-243™) and U937 cells (ATCC ® CRL-1593.2™) were cultured following the manufacturer's instructions. KBM7 cells were a kind gift from the T. R. Brummelkamp laboratory. They were cultured in DMEM 10% FBS. iPS derived from CD34 + cells were a gift from George Daley laboratory and were maintained in human ES cell media (Invitrogen) in the presence of FGF (Peprotech) and 20% of Replacement Serum (Invitrogen) as described 38 . Primary CML total bone marrow mononuclear cells were obtained after a written informed consent in accordance with the Declaration of Helsinki in agreement with the approved protocol from the review board of the S.Orsola-Malpighi Hospital, Bologna, Italy. Primary samples were cultured in QBSF-60 (Quality Biological), 20% FBS, SCF 50µg per ml, TPO 20µg perml, Flt3 10µg per ml, IL3 10µg per ml. Retroviral infections were performed as described 13,20 . Briefly, cells were transduced with a pool of the four factors Oct4, Sox2, Klf4, and c-Myc in a ratio of 1:1:1:1 for a 72-hour period. At day four post transduction, regular medium was switched to human ES cell media (Invitrogen) in the presence of FGF (Peprotech) and 20% of Replacement Serum (Invitrogen). After 15-16 days, human iPS-like colonies were morphologically identified, picked, and mechanically disaggregated on irradiated mouse embryonic fibroblasts (Globalstem). CML primary cells were treated with 1µM of 5-Azacytidine for 72 hours to allow monocytic differentiation.

In vitro differentiation
Embryoid body (EB) formation was induced by culturing human LiPS in clumps in the absence of FGF in bacterial plates to avoid attachment to the bottom of the plates. To allow spontaneous endoderm formation, after 3-4 days, embryoid bodies were transferred onto 0.2% gelatin-coated glass six-well plates and cultured in differentiation medium (DMEM supplemented with 20% fetal bovine serum, 2 mM l-glutamine, 0.1 mM 2-mercaptoethanol, nonessential amino acids, and penicillin-streptomycin, all from Invitrogen) for 2 to 3 weeks. The medium was changed every other day. For mesoderm/cardiomyocyte differentiation, LiPS cells were maintained on gelatin-coated plates in differentiated medium supplemented with 100 µM ascorbic acid (Sigma).
For ectoderm differentiation, EBs were cultured on laminin (Stemgent) coated six-well plates. Briefly, after 4 days in EB medium as a floating culture, the cells were plated on laminin (Stemgent) coated plates in N2B27 medium (Invitrogen), Sonic Hedgehog (0.1 µg/ µl), and FGF8 (100 ng/ml) (both from Peprotech), and maintained for 3 to 5 weeks in the absence of FGF2. The medium was changed every other day.

Immunofluorescence analysis
Cells were grown in six-well plates and fixed with 4% paraformaldehyde. The following antibodies were used: OCT4 (

Teratoma Formation and Analysis
Between 2 and 10 million iPS cells were mechanically scraped, pelleted, resuspended in 100 µl of hES media, and injected subcutaneously and intramuscularly into NOD.Cg-Prkdc scid Il2rg tm1Wjl /SzJ (NSG). Teratomas were harvested after 9 weeks, fixed overnight with 10% formaldehyde, embedded in paraffin, sectioned, HE stained, and analyzed using a BX40 microscope (Olympus). All animal work and protocols were approved by the Institutional Animal Care and Use Committee of Beth Israel Deaconess Medical Center (Protocol #006-2014) and were performed in accordance with the guidelines set forth in the Guide for the Care and Use of Laboratory Animals of the National Research Council of the National Institutes of Health.

Methylcellulose Colony Assay
Methylcellulose colony formation assays were performed with Methocult H4434 (Stem Cell Technologies) according to the manufacturer's instructions. CD34 + /CD45 + cells sorted from Bone Marrow or from in vitro differentiated CML-LiPS were suspended in culture medium containing methylcellulose and plated onto 35-mm culture dishes at a density of 10 3 (BMderived) or 10 4 (CML-LiPS derived) cells/dish, in triplicates. Colonies were counted after two weeks and processed for differentiation markers by FACS after three weeks.

SNP analysis
The analysis was performed using Affymetrix 6.0 SNP arrays platform and then computationally using the Partek Genomics Suite. Briefly, normalization using global median was first performed for array intensities cross all chips. Copy numbers and allele ratios were then estimated by comparison with a common reference to identify the copy number variations (CNVs) in the parental cell lines as well as the additional variations in copy numbers caused by reprogramming. Genomic segmentation was performed to smooth the data with the following parameters: 10 as the minimum of data points (SNPs), 0.001 as the cutoff for p value, 0.2 as the noise-to-signal cutoff.

Reduced Representation Bisulfite Sequencing
RRBS was performed as it follows. Genomic DNA was isolated from K562, K562-LiPS1, K562-LiPS2, KBM7, KBM7-LiPS1, KBM7-LiPS2, primary CML, CML-LiPS1 and CML-LiPS2, CD34 + iPS and digested with MspI (NEB), a methylation-insensitive enzyme that cuts C'CGG sites independent of the DNA methylation status of the center CpG. Digested DNA was size-selected on a 4% NuSieve 3:1 Agarose gel (Lonza). For each sample, two slices containing DNA fragments of 40-120 bp and 120-220 bp, respectively, were excised from the unstained preparative portion of the gel. These two size fractions were kept apart throughout the procedure including the final sequencing. Pre-annealed Illumina adaptors containing 5'-methyl-cytosine instead of cytosine were ligated to ends of size-selected MspI fragments. Adapter-ligated fragments were bisulfite-treated using the Epitect kit (Qiagen). The products were PCR amplified, size selected, and sequenced on the Illumina GA IIx at a reading length of 36 bp. Sequencing reads were mapped to the human genome assembly NCBI37/hg19 and mm9 for human and mouse data respectively using RRBSMAP 39 allowing 2 mismatches. Methylation levels were called using the BisSNP software 40 . RRBS data for human H1 ES line passage 37 and 38 and CD34 + cells were obtained from http:// www.epigenomeatlas.org. Further processing was done using RnBeads (http://rnbeads.mpiinf.mpg.de/) 41 . In this analysis, promoter methylation was defined as the mean across all scores for CpGs within a window of −1500bp to +500bp of the TSS of Ensembl-annotated genes. Differentially methylated sites were identified using a ranking-based approach taking into account statistical testing as well as absolute and relative differences in methylation levels between two groups of samples. In brief, we defined a combined ranking score for a single CpG as the maximum, i.e. worst, rank among the following three criteria: a) absolute difference in group mean methylation levels, b) absolute log ratios in group mean methylation levels and c) p-values originating from hierarchical linear models 42 . For genomic regions, such as CpG islands, genes, promoters, etc. the combined ranking score was defined analogously based on a) mean differences in mean methylation levels across all CpGs in a particular region b) mean ratios of mean methylation levels across CpGs c) combining the CpG specific p-values using an extension of Fisher's method 43 . Full RnBeads analysis reports, including tables characterizing differential methylation, can be accessed at http://reprogramming-leukemia.computational-epigenetics.org/.

Bisulfite Sequencing
Methylation profiling of CEBPA, gene loci was performed using 1µg of genomic DNA. For the analysis DNA was bisulfite converted using the EpiTect kit (Qiagen). The primers and PCR conditions for bisulfite sequencing are: F-OCT4: TAG TTG GGA TGT GTA GAG  TTT GAGA, R-OCT4: TAA ACC AAA ACA ATC CTT CTA CTC C; F-CEBPA: TAG  TTT YGT TAG TTT GGG GGG TTT

In vitro hematopoietic differentiation
Embryoid bodies were produced by scraping of confluent LiPS cells and cultured in suspension in embryoid body medium (90% DMEM, 10% FBS) for 96 hours. Embryoid bodies were then placed over a feeder layer of confluent OP9 stromal cells and allowed to attach. The medium used for the first 48 hours of differentiation was 50% embryoid body medium and 50% hematopoietic differentiation medium. The hematopoietic differentiation medium was prepared with QBSF-60 serum-free medium (Quality Biological Inc) supplemented with cytokines BMP4 (10ng/ml), VEGF (10ng/ml), SCF (25ng/ml), FGF (10ng/ml), TPO (20ng/ml), and Flt ligand (10ng/ml) (all from Peprotech). After 48h, cells were cultured with hematopoietic differentiation medium, with medium changes every 48h until the end of the differentiation protocol. At day 26, OP9 and embryoid bodies were collected by trypsinization (0.05% trypsin), washed and labeled for the cytofluorimetric analysis.

qRT-PCR
Total RNA was isolated from cells using TRIzol reagent (Invitrogen) and RNeasy Plus Mini Kit (Qiagen). Genomic DNA was eliminated with RNase free DNase 1 (Roche), and total RNA was then subjected to another round of purification by phenol solution with pH4.3 ± 0.2 (Sigma-Aldrich). For pluripotent gene expression analysis, cDNA was generated by SuperScript II reverse transcriptase (Invitrogen) with random hexamers (Invitrogen) from 1 µg of total RNA, and real-time PCR was performed with GoTaq qPCR Master Mix (Promega). Expression of CEBPA was quantified using Taqman one-step RT-PCR master mix (Applied Biosystems) relative to β-Actin. Primers used for CEBPA qRT-PCR were

Transplantation of immunocompromised and B6.SJL-Ptprc a Pepc b /BoyJ
Sorted LiPS-derived CD45 + cells and K562 cells were transferred by intra-orbital injection into sub-lethally irradiated (150 cGy) 6-week old recipient NSG mice (20,000 cells per mice). Transplanted mice were sacrificed for isolation of bone marrow and spleen. Transplantation was performed in quadruplicate for LiPS 1 clone and K562 and in triplicate for LiPS3 clone. 1×10 6 whole bone marrow cells obtained from SCL-tTA/BCR-ABL transgenic mice were transplanted by intra-orbital injection into 44-six weeks old CD45.1 B6.SJL-Ptprca Pepcb/ BoyJ congenic mice lethally irradiated. Tetracycline was maintained in the drinking water for three weeks after transplantation. BCR-ABL expression was induced by withdrawal of tetracycline and two weeks after BCR-ABL induction, treatment was started. Mice were randomly assigned to five groups: Control (Bcr-Abl un-induced); untreated (Bcr-Abl induced); IM (Bcr-Abl induced, treated with Imatinib); AZA (Bcr-Abl induced, treated with 5-Azacytidine); IM+AZA (Bcr-Abl induced, treated with Imatinib and 5-Azacytidine). The IM group received Imatinib at the dosage of 200mg per kg by gavauge, every other day. The AZA group received 5-Azacytidine by intraperitoneal injection at the dosage of 2mg per kg every other day. The IM-AZA group received a combination of Imatinib, 200mg per kg by gavage, and 5-Azacytidine by intraperitoneal injection at the dosage of 2mg per kg every other day. Mice were followed up for 150 days.

Statistical analysis
Survival curves were analyzed using the Kaplan-Meier method. Differences in the survival curves were evaluated with the log-rank test. Values of P ≤ 0.05 were considered statistically significant.
Statistical analyses of qRT-PCR were performed with paired two-tailed Student's t-test.

Figure 1. Reprogramming erases aberrant DNA methylation patterns in LiPS cells
A) DNA methylation profile of K562 and KBM7 cell lines, primary patient-derived CML cells, their corresponding LiPS clones, as well as human ES cells, CD34 + -iPS and primary CD34 + cells from healthy donors. Each column of the heat map represents the average over two biological replicates. The heat map visualizes DNA methylation levels, with red color indicating low levels of DNA methylation and blue color indicating high levels. White color indicates DNA regions not covered by RRBS analysis. For each genomic region type, shown are the 2% of regions that scored highest according to a ranking score in differential analysis Amabile et al. Page 16 comparing the group of K562, KBM7, and primary CML samples to the group including all derived iPS cell lines (see also the Methods section). B) DNA methylation profile of SALL4 gene locus, hypermethylated in all the leukemia cells and demethylated after reprogramming (upper panel). DNA methylation profiles of HOXA5 and BRCA1 gene loci are shown in the middle and lower panels, respectively. Each row represents a sample, while each column represents a CpG dinucleotide. Blue color indicates high DNA methylation levels, and red color indicates low methylation levels. CpG islands are annotated with green bars. Mean promoter and gene methylation levels are indicated at the right.

. BCR-ABL activation induces aberrant DNA methylation pattern in murine HSCs
A) FACS analysis of BCR-ABL transgenic murine bone marrow mononuclear cells in which the BCR-ABL oncogene is not expressed ("control"), is expressed ("leukemic") or repressed after previous expression ("rescued"). Leukemic mice showed moderate expansion of the HSCs population cKit + Sca1 + and a strong expansion of the Mac1 + Gr1 + myeloid population compared to control and rescued animals. B) DNA methylation analysis of CpG islands in the Lin − cKit + Sca1 + population of BCR-ABL transgenic mice with different status of oncogene expression. Zoomed-in region shows example of CpG islands differentially methylated between leukemic, control and rescue mice. Each row represents the average of two biological replicates. Each column represents a complete CpG island. Blue: high methylation, Red: low methylation. White color indicates DNA regions not covered by RRBS analysis. C) Heatmap visualizing genes and promoters hypermethylated in leukemic mice versus control and rescued ones. D) High-resolution view of DNA methylation of the Hoxb1 gene locus, which was hypermethylated in leukemic mice as compared to the control mice. Each row represents the average of two biological replicates; each column represents a single CpG dinucleotide.

Author Manuscript
Author Manuscript Author Manuscript Author Manuscript