Skip to main content

Tracking of epigenetic changes during hematopoietic differentiation of induced pluripotent stem cells

Abstract

Background

Differentiation of induced pluripotent stem cells (iPSCs) toward hematopoietic progenitor cells (HPCs) raises high hopes for disease modeling, drug screening, and cellular therapy. Various differentiation protocols have been established to generate iPSC-derived HPCs (iHPCs) that resemble their primary counterparts in morphology and immunophenotype, whereas a systematic epigenetic comparison was yet elusive.

Results

In this study, we compared genome-wide DNA methylation (DNAm) patterns of iHPCs with various different hematopoietic subsets. After 20 days of in vitro differentiation, cells revealed typical hematopoietic morphology, CD45 expression, and colony-forming unit (CFU) potential. DNAm changes were particularly observed in genes that are associated with hematopoietic differentiation. On the other hand, the epigenetic profiles of iHPCs remained overall distinct from natural HPCs. Furthermore, we analyzed if additional co-culture for 2 weeks with syngenic primary mesenchymal stromal cells (MSCs) or iPSC-derived MSCs (iMSCs) further supports epigenetic maturation toward the hematopoietic lineage. Proliferation of iHPCs and maintenance of CFU potential was enhanced upon co-culture. However, DNAm profiles support the notion that additional culture expansion with stromal support did not increase epigenetic maturation of iHPCs toward natural HPCs.

Conclusion

Differentiation of iPSCs toward the hematopoietic lineage remains epigenetically incomplete. These results substantiate the need to elaborate advanced differentiation regimen while DNAm profiles provide a suitable measure to track this process.

Background

Several protocols have been described for directed differentiation of iPSCs toward hematopoietic progenitor cells (iHPCs) [1,2,3]. These differentiation protocols frequently utilize a step-wise exposure of cells to specific cytokine cocktails that mimic multiple steps of mesoderm commitment and hematopoietic commitment [3,4,5]. iHPCs have been shown to express hematopoietic markers and to have a multipotent differentiation potential toward various hematopoietic lineages [6,7,8,9]. However, the efficiency of differentiation in terms of cell counts of resulting hematopoietic cells is rather low. More importantly, so far, iHPCs were not able to show long-term engraftment and thus iHPCs cannot be used in regenerative medicine [1, 5]. A better understanding of the molecular makeup of iHPCs is therefore important to further adjust differentiation procedures to ultimately achieve fully functional cells.

DNA methylation (DNAm) changes at CG dinucleotides (CpG sites) resemble an important epigenetic modification, which is relevant for normal hematopoietic differentiation into specific lineages [10, 11]. Conversely, DNAm patterns can be utilized to characterize the cellular composition of blood [12, 13]. Recently, Nishizawa et al. compared DNAm patterns among iPSC-lines that revealed varying differentiation potential toward the hematopoietic lineage depending on IGF2 levels in pluripotent cells [14]. This study also described DNAm profiles of iPSCs and iPSC-derived HPCs. However, a systematic comparison of DNAm profiles of iHPCs with normal hematopoietic counterparts was not yet performed.

It is well known that primary HPCs are regulated by tight interaction with their niche, which can be mimicked by co-culture with suitable feeder layers in vitro [15]. Mesenchymal stromal cells (MSCs) were shown to support culture expansion of primary HPCs [16,17,18], and it has even been tested in clinical trials if co-culture of HPCs with MSCs increases hematopoietic engraftment in allogeneic transplantation [19]. Hence, it might be anticipated that co-culture of iHPCs with MSCs might further support culture expansion and hematopoietic commitment of iHPCs. MSC-like cells can also be derived from iPSCs [20,21,22], and we have previously described that such iMSCs support culture expansion of primary HPCs [23]. Thus, iPSCs might be differentiated in parallel toward either the hematopoietic or mesenchymal stromal cell lineage for subsequent co-culture in a syngenic setting. This might provide new perspectives to further increase differentiation efficiency.

In this study, we compared DNAm patterns of iHPCs with various hematopoietic cell types. We demonstrate that in particular patterns of the monocytic lineage are recapitulated, while DNAm patterns remained overall different from primary hematopoietic cells. There was no evidence that additional co-culture with syngenic MSCs or iMSCs supported further epigenetic maturation of iHPCs.

Results

Differentiation of iPSCs toward the hematopoietic lineage is reflected by corresponding DNA methylation changes

Directed differentiation of iPSCs into iHPCs was induced by embryoid body (EB) formation and subsequent stage-specific application of specific cytokine cocktails, as indicated in Fig. 1a [4, 24]. These experiments were done with three iPSC lines that were reprogrammed from MSCs of different donors. After 10 days, cells with typical hematopoietic morphology started to bud off from the endothelial layer and these cells were harvested from the supernatant after 20 days (Additional file 1: Fig. S1a). Cytospins revealed predominantly a rather monocytic morphology (Additional file 1: Fig. S1b). Immunophenotype analysis demonstrated high expression of CD31 and CD45, which is in line with hematopoietic cell fate decisions (Additional file 1: Fig. S1c, d). However, high expression of CD34, which is indicative for hematopoietic progenitor cells, was only observed on a subpopulation at day 6 and was downregulated at day 20. Similar findings were described by other authors before [4], and the results are overall in line with differentiation toward the hematopoietic lineage.

Fig. 1
figure 1

DNA methylation profiles of iPSC-derived hematopoietic progenitor cells. a Schematic representation of the differentiation protocol of iPSCs toward hematopoietic progenitor cells. b Principal component analysis of DNA methylation profiles of our iHPCs (donor 1 and 2) and those of Nishizawa et al. [14], primary CD34+ cells from cord blood (CB) or bone marrow (BM), whole blood, and pluripotent cells. c Scatter plot depicts CpG sites that are either hypermethylated (961 CpG sites, red) or hypomethylated (6075 CpG sites, blue; delta mean β value > 0.5 or < − 0.5) in iHPCs as compared to iPSC (GSE37066). CpG sites associated with promoter regions are highlighted in bold. d Gene ontology analysis of genes with differentially methylated CpG sites in the promoter region. Enrichment of specific categories was calculated by the one-sided Fisher’s exact P value

We have then analyzed DNAm profiles of two iPSC clones after 20 days of differentiation using the Illumina Infinium MethylationEPIC BeadChip. Principal component analysis demonstrated that iHPCs clustered closely together with iPSC-derived hematopoietic progenitor cells of Nishizawa et al. [14] (Fig. 1b). These authors used a differentiation protocol with a different cytokine composition and without hypoxic conditions. Thus, the epigenetic state of iHPCs appears to be independent of the differentiation regimen used. Notably, the iHPCs were overall clearly separated from primary hematopoietic cells, indicating that hematopoietic differentiation was epigenetically incomplete.

To better understand the epigenetic changes that are acquired during differentiation of iPSCs into iHPCs, we filtered for CpGs with more than 50% change in DNAm levels: 961 CpGs were hypermethylated and 6075 CpGs were hypomethylated in iHPCs as compared to iPSCs (Fig. 1c). These DNAm changes were enriched in the gene body and intergenic regions, whereas they were less present in the promoter regions that may have more impact on gene expression (Additional file 2: Fig. S2a). Two hundred twenty hypermethylated and 1493 hypomethylated CpGs were associated with promoter regions (Additional file 3: Tab. S1). These differentially methylated CpG sites show similar methylation patterns in iHPCs as compared to cord blood-derived CD34+ cells (Additional file 4: Fig. S3a). Gene ontology analysis of the associated genes revealed enrichment of hypomethylation in categories of the immune system, blood cell activation, and hemostasis, whereas hypermethylation was rather related to development and cell fate (Fig. 1d).

Subsequently, we analyzed differential DNAm between iHPCs and primary CD34+ cells from umbilical cord blood (50% difference in DNAm levels): 2692 CpGs were higher methylated and 2307 CpGs were less methylated in iHPCs (Fig. 2a). Hypermethylation in iHPCs was enriched at intergenic regions, CpG islands, and their N shore region, whereas hypomethylation was enriched in 5′UTR, gene bodies, and N shore regions (Additional file 2: Fig. S2b). Six hundred fifty-nine hypermethylated and 587 hypomethylated CpGs were associated with promoter regions (Additional file 5: Tab. S2). Gene ontology analysis of the associated genes showed enrichment of hypermethylation in categories of various cellular processes, such as cell division or membrane potential, whereas hypomethylation was related to the immune system, hematopoiesis, development, and metabolic processes (Fig. 2b). These differentially regulated CpG sites do not recapitulate the DNAm patterns in MSCs but they were overall related to DNAm in iPSCs, indicating that these epigenetic modifications of hematopoietic development are not acquired in iHPCs during differentiation (Additional file 4: Fig. S3b). We have then exemplarily analyzed DNAm patterns within the genes CD34 and CD45. Overall, iHPCs acquired a similar DNAm pattern as observed in primary CD34+ progenitor cells, particularly within the CD45 gene (Fig. 2c). Thus, the DNAm profiles clearly recapitulated features of normal hematopoietic differentiation, albeit the overall epigenetic makeup was yet incomplete.

Fig. 2
figure 2

Comparison of DNA methylation profiles of iPSC-derived HPCs versus primary CD34+ cells. a Scatter plot of CpG sites that are either higher methylated (2692 CpG sites, red) or less methylated in iHPCs (2307 CpG sites, blue; delta mean β value > 0.5 or < − 0.5) compared to primary CD34+ cells. CpG sites associated with promoter regions are highlighted in bold (659 CpG sites and 587 CpG sites, respectively). b Gene ontology analysis of genes with differentially methylated CpG sites in the promoter region. Enrichment of specific categories was calculated by the one-sided Fisher’s exact P value. c DNA methylation patterns within the genes CD45 and CD34. Mean DNAm levels are depicted for iPSCs (GSE37066), iHPC d20, and primary CD34+ cells. TSS1500: 1500 bp upstream of transcription start site; UTR: untranslated region. For all analyses, iHPCs were compared to primary CD34+ cells isolated from human cord blood (GSE40799)

Epigenetic changes point toward monocytic differentiation

We subsequently analyzed whether lineage-specific DNAm patterns are recapitulated in iHPCs. Integrative PCA analysis of DNAm profiles of iHPCs and various hematopoietic cell types demonstrated that the epigenetic profiles of iHPCs were rather distinct from all other cell types (Fig. 3a). We then used a deconvolution algorithm that has been developed to estimate the percentages of cell types in blood based on DNAm profiles [13]. The results indicated that the majority of our iHPCs and of the iPSC-derived progenitor cells generated by Nishizawa et al. [14] were classified as monocytes (Fig. 3b). However, this tool has been developed for primary cells and only includes few different cell types. Therefore, we independently filtered for specific sets of CpGs that are characteristically hypomethylated in granulocytes, CD4+ T cells, CD8+ T cells, B cells, NK cells, monocytes, or lymphocytes, as described before [12]. Clustering according to DNAm levels in these cell type-specific CpGs supported the notion that iHPCs recapitulate DNAm patterns of the monocytic lineage (Fig. 3c).

Fig. 3
figure 3

Comparison of DNA methylation profiles of iPSC-derived HPCs with different hematopoietic cell types. a Principal component analysis demonstrates that DNAm profiles of iHPCs are distinct to other cell types. b A deconvolution algorithm of Houseman et al. [13] was used to estimate the composition of different hematopoietic cell types based on the DNAm profiles of our iHPCs (donor 1 and 2) and those of Nishizawa et al. [14]. c Heatmap of DNAm levels at CpGs that are specifically hypomethylated in specific hematopoietic cell types. iHPCs (d20) are most closely related to monocytes (CD34+ cells isolated from human cord blood: GSE40799; B cells, peripheral blood mononuclear cells (PBMC), CD4+ T cells, CD8+ T cells, NK cells, whole blood, granulocytes, monocytes: all GSE35069; and iPSCs: GSE37066). Gray areas indicate that no data is available for this specific CpG site. CB: cord blood

Mesenchymal stromal cells support proliferation of iPSC-derived hematopoietic progenitor cells

Subsequently, we followed the hypothesis that co-culture with MSCs might provide a cellular microenvironment that stimulates further differentiation toward the hematopoietic lineage. To this end, we used MSCs and iPSC-derived MSCs (iMSCs) that were generated from the same donors used for the reprogramming of iPSCs and differentiation into iHPCs (Fig. 4a). iMSCs revealed typical MSC-like morphology, immunophenotype, and three-lineage differentiation potential (Additional file 6: Fig. S4). Upon culture expansion for an additional 14 days on MSC and iMSC stroma or tissue culture plastic (TCP) as control, the vast majority of cells were non-adherent round cells with typical hematopoietic morphology (Fig. 4b). Proliferation of iHPCs was faster under co-culture conditions (Fig. 4c), which is in line with the accelerated proliferation of primary HPCs on a suitable feeder layer [16]. Furthermore, colony-forming unit (CFU) potential was maintained under co-culture conditions with MSCs and iMSCs, but not on TCP (Fig. 4d). However, compared to the iHPCs at day 20 without additional culture expansion, the proportion of cells with CFU forming potential was significantly decreased upon culture expansion. Flow cytometric analysis demonstrated that hematopoietic markers except for cKIT were downregulated by additional culture expansion (Fig. 4e,f). Cellular morphology in cytospins indicated multilineage hematopoietic differentiation, albeit most iHPCs were relatively large with a foam-like appearance of the cytoplasm that might be related to a macrophage phenotype [24] (Fig. 4g). Thus, the stromal support of MSCs or iMSCs does not seem to improve maintenance of hematopoietic progenitors and their multilineage potential as compared to iHPCs without further culture expansion.

Fig. 4
figure 4

Co-culture with MSCs or iPSC-derived MSCs supports growth of iPSC-derived HPCs. a Schematic representation of co-culture experiments of iHPCs with syngenic MSCs and iMSCs (donors 1–3). b Phase contrast images of iHPCs upon expansion for 2 weeks on tissue culture plastic (TCP), MSCs, and iMSCs. Scale bar 100 μm. c Proliferation of iHPCs on TCP, MSCs, or iMSCs was estimated by residual CFSE stain after 5 days. The mean fluorescence intensities (MFI) were normalized to corresponding measurements at d0. Data represent the mean of three independent experiments ± SD (ANOVA was not significant). d Colony-forming unit (CFU) frequency of iHPCs at day 20 and after 14 days of co-culture expansion on TCP, MSCs, or iMSCs. Data represent mean of three independent experiments ± SD; ANOVA: **P < 0.01; ***P < 0.001. e Flow cytometry analysis of iHPCs on day 14 of co-culture as depicted in a. f Quantitative analysis of flow cytometric measurements. GlyA: glycophorin A. g Cytospins of iHPCs at day 34. Morphologies relate to hematopoietic progenitor cells (1), macrophages (2), monocytes (3), and neutrophils (4). Scale bar 20 μm

Additional culture expansion does not result in epigenetic maturation

To estimate the impact of stromal support on the epigenetic makeup of iHPCs, we analyzed DNAm profiles of the expanded cell preparations. Albeit co-cultures resembled a mixture of iHPCs with MSCs or iMSCs, we did not separate the cell types by cell sorting, as the results would be largely affected by sorting strategies and the majority of cells seemed to be iHPC-derived. In fact, PCA demonstrated that the additional culture expansion for 2 weeks with or without stromal support resulted in DNAm profiles that remained similar to those of iHPCs at day 20 (Fig. 5). In tendency, co-cultured iHPCs clustered closer to iMSCs than the non-co-cultured iHPCs due to the residual feeder cells. Either way, the results demonstrated that additional culture expansion of iHPCs either with or without stromal support did not support further epigenetic maturation toward a DNAm profile of the hematopoietic lineage.

Fig. 5
figure 5

Principal component analysis of DNA methylation profiles of iPSC-derived HPCs with and without co-culture. The DNA methylation profiles of iHPCs (d20; donor 1 and 2) and iPSC-derived HPCs from Nishizawa et al. [14] are closely related to those of iHPCs that were further cultured for 2 weeks with or without stromal support. Additional stromal support (iHPC on MSCs, iMSCs, and TCP, respectively) did not enhance further commitment toward the hematopoietic lineage

Discussion

Generation of hematopoietic progenitor cells from iPSCs opens new perspectives for patient-specific disease modeling and drug testing. Furthermore, it may eventually overcome limitations in culture expansion of primary HPCs for cellular therapy. However, the results of this study support the notion that with current differentiation regimen, iHPCs remain epigenetically distinct from their in vivo counterparts.

Optimization of differentiation conditions for iHPCs remains a challenge [3]. Most protocols seem to be biased toward myeloid differentiation, which is in line with results shown here demonstrating that DNAm profiles are related to the monocytic lineage. The immunophenotype of our iHPCs is overall similar to previous studies, albeit, for example, Sontag et al. [24] and Dorn et al. [25] described a more pronounced upregulation of CD43. More recently, several groups suggested alternative protocols that utilize overexpression of specific transcription factors or activation of signal cascades that may better support definitive hematopoiesis resulting in improved engraftment capabilities of hematopoietic progenitors [5, 26, 27]. In the future, it will be interesting to investigate if these approaches give rise to progeny that resemble their primary counterparts more closely on epigenetic levels.

DNA methylation has been demonstrated to be an important regulative mechanism for cell-fate decisions in hematopoietic differentiation [11, 28]. Here, we show that DNAm profiles of our iHPCs are closely related to those by Nishizawa and colleagues [14], although these authors did not use the same differentiation conditions. It will be important to gain further understanding of how epigenetic profiles of iPSC-derived hematopoietic cells are affected by various differentiation regimens. The discrepancy in the epigenetic makeup of iHPCs versus in vivo hematopoietic stem and progenitor cells may at least be partially attributed to other epigenetic signatures that are not directly related to lineage-decisions, but rather to attributes such as tissue-specific and age-associated DNAm patterns [20, 29].

It has been shown that early steps in normal hematopoietic development are reflected by continuous gain or loss of DNAm during differentiation and that these changes are overall inversely correlated with gene expression at key loci [28, 30]. This is in line with our observation that demethylation in iHPCs is particularly enriched in promoter regions of genes that are relevant for hematopoietic differentiation. On the other hand, association of DNAm with gene expression is frequently rather low [28, 31, 32]. It may be too simplistic to consider only DNAm changes as functionally relevant that directly impact on gene expression changes. DNAm patterns are also closely related to other epigenetic marks, such as histone modifications and higher order chromatin conformation. These aspects should be taken into consideration for future studies.

Mesenchymal stromal cells represent an important component of the hematopoietic stem cell niche [33]. They secrete various chemokines and facilitate cell-cell interactions that were shown to support primary HPCs also in vitro [16, 34]. In this study, we demonstrate that MSCs and iMSCs support proliferation and maintenance of CFU potential of iHPCs. However, our results do not indicate that this additional co-culture regimen increases CFU frequency or maturation of a hematopoietic immunophenotype. Furthermore, the additional culture expansion, either with or without stromal support, does not result in epigenetic maturation of iHPCs to their in vivo counterparts. It is conceivable that co-culture needs to be implemented at earlier steps of the differentiation procedure or that it needs to be combined with other culture media.

Conclusion

Differentiation of iPSC toward the hematopoietic lineage is associated with DNAm profiles that point toward normal hematopoietic differentiation, but the epigenetic patterns remain incomplete. It will be important to investigate whether epigenetic assimilation can be improved by alternative regimen for hematopoietic differentiation, and for this purpose, DNAm profiles provide an important measure for evaluation of differentiation procedures.

Methods

Cell culture

Mesenchymal stromal cells were isolated from bone marrow (caput femoris) of patients undergoing orthopedic surgery and culture-expanded in standard medium consisting of Dulbecco’s Modified Eagle Medium (DMEM; 1 g/l glucose; PAA, Pasching, Austria), 1% L-glutamine (PAA), 1% penicillin/streptomycin (PAA), and 10% pooled human platelet lysate (hPL) [35]. The medium was supplemented with 0.1% heparin (5000 IU/ml; Ratiopharm, Ulm, Germany) to prevent coagulation [36].

Induced pluripotent stem cells were generated from three MSC preparations with episomal plasmids [37] and thoroughly characterized as described before [38, 39]. iPSCs were cultured on tissue culture plastic (TCP) coated with vitronectin (0.5 mg/cm2) in StemMACS iPS-Brew XF (all Miltenyi Biotec, Bergisch Gladbach, Germany). Pluripotency was validated by in vitro differentiation and Epi-Pluri-Score (Cygenia GmbH, Aachen, Germany) [40].

Generation of iPSC-derived MSCs (iMSCs) was performed by switching culture conditions to standard hPL-medium and further passaging on 0.1% gelatin-coated plates [20]. Three-lineage differentiation potential of MSCs and iMSCs was validated as described before [41, 42]. Cell images were taken on a digital EVOS FL Auto microscope (Thermo Fisher Scientific, Carlsbad, California, USA).

Hematopoietic progenitor differentiation

Differentiation of iPSCs toward hematopoietic progenitor cells was performed as described before with iPSCs from three different donors [24]. In brief, iPSC colonies were harvested and aggregated into embryoid bodies (EBs), which were cultured in basal differentiation medium consisting of StemPro34 supplemented with 100 U/ml penicillin, 100 μg/ml streptomycin, 2 mM L-glutamine (all Thermo Fisher Scientific), 0.4 mM monothioglycerol, and 50 μg/ml L-ascorbic acid (all Sigma Aldrich, St. Louis, Missouri, USA). The following growth factors were added as indicated in Fig. 1a: 8 ng/ml bone morphogenic protein 4 (BMP4), 10 ng/ml human basic fibroblast growth factor (bFGF), 10 ng/ml vascular endothelial growth factor (VEGF), 100 ng/ml interleukin 6 (IL-6), 100 ng/ml stem cell factor (SCF), 25 ng/ml insulin-like growth factor 1 (IGF1), 30 ng/ml interleukin 3 (IL-3), 20 ng/ml thrombopoietin (TPO), and 10 ng/ml fms-related tyrosine kinase 3 ligand (FLT3L). BMP4, IL-3, IL-6, and SCF were purchased from Miltenyi Biotec, and all other cytokines from Peprotech (Hamburg, Germany). Between day 0 and day 7, cells were cultured in a 5% O2 and 5% CO2 atmosphere. On day 6, cells were moved to 0.1% gelatin-coated TCP, and on day 8, cultures were transferred to normoxia and 5% CO2 atmosphere. On day 20, iHPCs were harvested for analysis or additional co-culture by collecting the supernatant, washing the well twice with PBS and centrifugation of the cells.

Immunophenotypic analysis

Flow cytometric analysis of hematopoietic surface markers was analyzed on a FACS Canto II (BD Biosciences, Franklin Lakes, New Jersey, USA) with the following antibodies: CD3-PerCP-Cy5.5 (clone UCHT1), CD19-FITC (clone HIB19), CD31-PE (clone WM59), CD34-APC (clone 581), CD38-PE (clone HIT2), CD43-FITC (clone 1G10), CD14-PE (clone MOP9), CD66b-PE (clone G10F5) (all BD Bioscience), CD45-APC-Cy7 (clone 2D1), CD117-PE-Cy7 (clone 104D2), CD11c-PE-Cy7 (clone 3.9), HLA-DR-FITC (clone LN3), CD235a-Pacific Blue (clone 6A7M) (all eBioscience, San Diego, California, USA), CD33-APC (clone AC104.3E3; Miltenyi Biotec). EBs were dissociated with Accutase (Stemcell Technologies, Vancouver, British Columbia, Canada) for 10–15 min prior to staining with antibodies. Single cells were incubated with 1% human IgG solution (Privigen, CSL Behring, King of Prussia, Pennsylvania, USA) for 30 min at 4 °C to block unspecific binding.

Immunophenotypic analysis of MSCs and iMSCs was performed with CD14-APC (clone M5E2), CD29-PE (clone MAR4), CD31-PE (clone WM59), CD34-APC (clone 581), CD45-APC (clone HI30), CD73-PE (clone AD2), CD90-APC (clone 5E10) (all BD Biosciences), and CD105-FITC (clone MEM-226; ImmunoTools, Friesoythe, Germany). Data was analyzed using the FlowJo software (Tree Star, Ashland, Oregon, USA).

Cytospin

Cellular morphology of iHPCs was analyzed by cytospins with Diff Quik staining (Medion Grifols Diagnostics, Düdingen, Switzerland) and analyzed with a Leica DMRX microscope (Leica Microsystems, Wetzlar, Germany).

Colony-forming unit assay

After 20 days of differentiation or after additional co-culture with MSCs or iMSCs for 2 weeks (day 34 of differentiation), 1000 iHPCs were seeded in 1 ml of methylcellulose-based medium (HSC-CFU lite with EPO; Miltenyi Biotec) in 24-well plates. Colonies were quantified 14 days later.

Proliferation analysis

Cells were labeled with Carboxyfluorescein N-succinimidyl ester (CFSE; Sigma-Aldrich) and the residual dye was analyzed after 5 days by flow cytometry, as described before [16].

Co-culture conditions

Mesenchymal stromal cells, iPSCs, and iMSCs were generated from the same three donors. After culture expansion, MSCs and iMSCs were irradiated (30 Gy) and used as feeder layers at a density of 50,000 cells/cm2 on 0.1% gelatin-coated dishes. After differentiation, 100,000 iHPCs were seeded onto MSCs, iMSCs, or TCP in StemSpan SFEM (Stemcell Technologies, Cologne, Germany) supplemented with 10 ng/ml fibroblast growth factor 1 (FGF-1), 10 ng/ml SCF, 10 ng/ml TPO (all purchased from Peprotech), and 1% penicillin/streptomycin [15].

DNA methylation analysis

Genomic DNA of iHPCs from two donors was harvested with the NucleoSpin Tissue kit (Macherey-Nagel, Düren, Germany), bisulfite-converted with the EZ DNA Methylation™ Kit (Zymo Research, Irvine, California, USA), and analyzed on the Infinium MethylationEPIC BeadChip (Illumina, San Diego, California, USA). Hybridization and initial data processing with the GenomeStudio (v2011.1) Methylation Module (v1.9.0) was performed at Life&Brain (Bonn, Germany). Raw data has been deposited at Gene Expression Omnibus (GEO) under the accession number GSE119079. In addition, we used previously published DNAm profiles of iHPCs (GSE60811) [14]; MSCs, embryonic stem cells (ESCs), and iPSCs (all GSE37066) [43]; iMSCs (GSE95058) [38]; cord blood-derived and bone marrow-derived CD34+ cells (GSE40799 [44], GSE60923 [14], GSE66459 [45], GSE49618 [46], GSE58477 [47]); whole blood (GSE40279 [48], GSE35069 [49]); and different hematopoietic subsets (GSE35069 [49], GSE45459 [50]).

For further analysis, we focused on CpGs that were represented on both the Infinium HumanMethylation450K BeadChip (in total 485,577 CpGs on the platform) and the Infinium MethylationEPIC BeadChip (in total 867,926 CpGs) with an overlap of 454,181 CpG sites on both platforms. CpGs on X and Y chromosomes were excluded resulting in 443,542 CpGs for further analysis. More than 99% of these filtered CpGs were detected in each iHPC sample. Principal component analysis (PCA) was performed with R; heatmaps and hierarchical clustering (HCL) were generated with the MultiExperiment Viewer (MeV; version 4.9.0). Differentially methylated CpGs were selected by the mean difference in DNAm as indicated in the text. For enrichment analysis, CpG sites were classified according to gene regions and CpG islands as described before [43, 51]. Gene ontology (GO) analysis was performed on genes with differentially methylated CpGs in the promoter region (located in TSS1500, TSS200, and 5′UTR) using the GoMiner software [52]. Categories comprising more than 1000 genes were not considered and similar categories are only listed once.

Statistical analysis

Proliferation, immunophenotypic, and CFU experiments were performed with three independent biological replicates and results are presented as mean ± standard deviation (SD). Statistical significance was estimated by unpaired one-way ANOVA in GraphPad Prism 6. Enrichment of gene ontology categories was calculated by the one-sided Fisher’s exact P value using all genes represented on the array as a reference (within the 443,542 CpGs that were filtered for this analysis as described above). For enrichment analyses in gene regions or CpG islands, statistical significance was estimated by hypergeometric distribution in Microsoft Excel.

Abbreviations

CB:

Cord blood

CFU:

Colony-forming unit

DNAm:

DNA methylation

ESC:

Embryonic stem cell

HPC:

Hematopoietic progenitor cell

iHPC:

Induced pluripotent stem cell-derived hematopoietic progenitor cell

iMSC:

Induced pluripotent stem cell-derived mesenchymal stromal cell

iPSC:

Induced pluripotent stem cell

MSC:

Mesenchymal stromal cell

PCA:

Principal component analysis

TCP:

Tissue culture plastic

TSS:

Transcription start site

UTR:

Untranslated region

References

  1. Ditadi A, Sturgeon CM, Keller G. A view of human haematopoietic development from the Petri dish. Nat Rev Mol Cell Biol. 2017;18(1):56–67.

    Article  CAS  Google Scholar 

  2. Slukvin II. Hematopoietic specification from human pluripotent stem cells: current advances and challenges toward de novo generation of hematopoietic stem cells. Blood. 2013;122(25):4035–46.

    Article  CAS  Google Scholar 

  3. Ackermann M, Liebhaber S, Klusmann JH, Lachmann N. Lost in translation: pluripotent stem cell-derived hematopoiesis. EMBO Mol Med. 2015;7(11):1388–402.

    Article  CAS  Google Scholar 

  4. Kennedy M, Awong G, Sturgeon CM, Ditadi A, LaMotte-Mohs R, Zuniga-Pflucker JC, et al. T lymphocyte potential marks the emergence of definitive hematopoietic progenitors in human pluripotent stem cell differentiation cultures. Cell Rep. 2012;2(6):1722–35.

    Article  CAS  Google Scholar 

  5. Sugimura R, Jha DK, Han A, Soria-Valles C, da Rocha EL, Lu YF, et al. Haematopoietic stem and progenitor cells from human pluripotent stem cells. Nature. 2017;545(7655):432–8.

    Article  CAS  Google Scholar 

  6. Hansen M, Varga E, Aarts C, Wust T, Kuijpers T, von Lindern M, et al. Efficient production of erythroid, megakaryocytic and myeloid cells, using single cell-derived iPSC colony differentiation. Stem Cell Res. 2018;29:232–44.

    Article  CAS  Google Scholar 

  7. Börger AK, Eicke D, Wolf C, Gras C, Aufderbeck S, Schulze K, et al. Generation of HLA-universal iPSCs-derived megakaryocytes and platelets for survival under refractoriness conditions. Mol Med. 2016;22:274–85.

    Article  Google Scholar 

  8. Zhang H, Xue C, Shah R, Bermingham K, Hinkle CC, Li W, et al. Functional analysis and transcriptomic profiling of iPSC-derived macrophages and their application in modeling Mendelian disease. Circ Res. 2015;117(1):17–28.

    Article  CAS  Google Scholar 

  9. Lachmann N, Ackermann M, Frenzel E, Liebhaber S, Brennig S, Happle C, et al. Large-scale hematopoietic differentiation of human induced pluripotent stem cells provides granulocytes or macrophages for cell replacement therapies. Stem Cell Reports. 2015;4(2):282–96.

    Article  CAS  Google Scholar 

  10. Meissner A, Mikkelsen TS, Gu H, Wernig M, Hanna J, Sivachenko A, et al. Genome-scale DNA methylation maps of pluripotent and differentiated cells. Nature. 2008;454(7205):766–70.

    Article  CAS  Google Scholar 

  11. Broske AM, Vockentanz L, Kharazi S, Huska MR, Mancini E, Scheller M, et al. DNA methylation protects hematopoietic stem cell multipotency from myeloerythroid restriction. Nat Genet. 2009;41(11):1207–15.

    Article  Google Scholar 

  12. Frobel J, Bozic T, Lenz M, Uciechowski P, Han Y, Herwartz R, et al. Leukocyte counts based on DNA methylation at individual cytosines. Clin Chem. 2018;64(3):566–75.

    Article  CAS  Google Scholar 

  13. Houseman EA, Accomando WP, Koestler DC, Christensen BC, Marsit CJ, Nelson HH, et al. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics. 2012;13:86.

    Article  Google Scholar 

  14. Nishizawa M, Chonabayashi K, Nomura M, Tanaka A, Nakamura M, Inagaki A, et al. Epigenetic variation between human induced pluripotent stem cell lines is an indicator of differentiation capacity. Cell Stem Cell. 2016;19(3):341–54.

    Article  CAS  Google Scholar 

  15. Walenda T, Bokermann G, Ventura Ferreira M, Pieroth D, Hieronymus T, Neuss S, et al. Synergistic effects of growth factors and mesenchymal stromal cells for expansion of hematopoietic stem and progenitor cells. Exp Hematol. 2011;39(6):617–28.

    Article  CAS  Google Scholar 

  16. Walenda T, Bork S, Horn P, Wein F, Saffrich R, Diehlmann A, et al. Co-culture with mesenchymal stromal cells increases proliferation and maintenance of haematopoietic progenitor cells. J Cell Mol Med. 2010;14(1–2):337–50.

    Article  CAS  Google Scholar 

  17. Robinson SN, Ng J, Niu T, Yang H, McMannis JD, Karandish S, et al. Superior ex vivo cord blood expansion following co-culture with bone marrow-derived mesenchymal stem cells. Bone Marrow Transplant. 2006;37(4):359–66.

    Article  CAS  Google Scholar 

  18. Fei XM, Wu YJ, Chang Z, Miao KR, Tang YH, Zhou XY, et al. Co-culture of cord blood CD34(+) cells with human BM mesenchymal stromal cells enhances short-term engraftment of cord blood cells in NOD/SCID mice. Cytotherapy. 2007;9(4):338–47.

    Article  CAS  Google Scholar 

  19. Lee SH, Lee MW, Yoo KH, Kim DS, Son MH, Sung KW, et al. Co-transplantation of third-party umbilical cord blood-derived MSCs promotes engraftment in children undergoing unrelated umbilical cord blood transplantation. Bone Marrow Transplant. 2013;48(8):1040–5.

    Article  CAS  Google Scholar 

  20. Frobel J, Hemeda H, Lenz M, Abagnale G, Joussen S, Denecke B, et al. Epigenetic rejuvenation of mesenchymal stromal cells derived from induced pluripotent stem cells. Stem Cell Reports. 2014;3(3):414–22.

    Article  CAS  Google Scholar 

  21. Sheyn D, Ben-David S, Shapiro G, De Mel S, Bez M, Ornelas L, et al. Human induced pluripotent stem cells differentiate into functional mesenchymal stem cells and repair bone defects. Stem Cells Transl Med. 2016;5(11):1447–60.

    Article  CAS  Google Scholar 

  22. Diederichs S, Tuan RS. Functional comparison of human-induced pluripotent stem cell-derived mesenchymal cells and bone marrow-derived mesenchymal stromal cells from the same donor. Stem Cells Dev. 2014;23(14):1594–610.

    Article  Google Scholar 

  23. Vasko T, Frobel J, Lubberich R, Goecke TW, Wagner W. iPSC-derived mesenchymal stromal cells are less supportive than primary MSCs for co-culture of hematopoietic progenitor cells. J Hematol Oncol. 2016;9:43.

    Article  Google Scholar 

  24. Sontag S, Forster M, Qin J, Wanek P, Mitzka S, Schuler HM, et al. Modelling IRF8 deficient human hematopoiesis and dendritic cell development with engineered iPS cells. Stem Cells. 2017;35(4):898–908.

    Article  CAS  Google Scholar 

  25. Dorn I, Klich K, Arauzo-Bravo MJ, Radstaak M, Santourlidis S, Ghanjati F, et al. Erythroid differentiation of human induced pluripotent stem cells is independent of donor cell type of origin. Haematologica. 2015;100(1):32–41.

    Article  CAS  Google Scholar 

  26. Tan YT, Ye L, Xie F, Beyer AI, Muench MO, Wang J, et al. Respecifying human iPSC-derived blood cells into highly engraftable hematopoietic stem and progenitor cells with a single factor. Proc Natl Acad Sci U S A. 2018;115(9):2180–5.

    Article  CAS  Google Scholar 

  27. Leung A, Zulick E, Skvir N, Vanuytsel K, Morrison TA, Naing ZH, et al. Notch and aryl hydrocarbon receptor signaling impact definitive hematopoiesis from human pluripotent stem cells. Stem Cells. 2018;36:1004–19.

    Article  CAS  Google Scholar 

  28. Farlik M, Halbritter F, Muller F, Choudry FA, Ebert P, Klughammer J, et al. DNA methylation dynamics of human hematopoietic stem cell differentiation. Cell Stem Cell. 2016;19(6):808–22.

    Article  CAS  Google Scholar 

  29. Weidner CI, Lin Q, Koch CM, Eisele L, Beier F, Ziegler P, et al. Aging of blood can be tracked by DNA methylation changes at just three CpG sites. Genome Biol. 2014;15(2):R24.

    Article  Google Scholar 

  30. Cabezas-Wallscheid N, Klimmeck D, Hansson J, Lipka DB, Reyes A, Wang Q, et al. Identification of regulatory networks in HSCs and their immediate progeny via integrated proteome, transcriptome, and DNA methylome analysis. Cell Stem Cell. 2014;15(4):507–22.

    Article  CAS  Google Scholar 

  31. Kulis M, Heath S, Bibikova M, Queiros AC, Navarro A, Clot G, et al. Epigenomic analysis detects widespread gene-body DNA hypomethylation in chronic lymphocytic leukemia. Nat Genet. 2012;44(11):1236–42.

    Article  CAS  Google Scholar 

  32. Wagner JR, Busche S, Ge B, Kwan T, Pastinen T, Blanchette M. The relationship between DNA methylation, genetic and expression inter-individual variation in untransformed human fibroblasts. Genome Biol. 2014;15(2):R37.

    Article  Google Scholar 

  33. Mendez-Ferrer S, Michurina TV, Ferraro F, Mazloom AR, Macarthur BD, Lira SA, et al. Mesenchymal and haematopoietic stem cells form a unique bone marrow niche. Nature. 2010;466(7308):829–34.

    Article  CAS  Google Scholar 

  34. Wagner W, Roderburg C, Wein F, Diehlmann A, Frankhauser M, Schubert R, et al. Molecular and secretory profiles of human mesenchymal stromal cells and their abilities to maintain primitive hematopoietic progenitors. Stem Cells. 2007;10:2638–57.

    Article  Google Scholar 

  35. Horn P, Bokermann G, Cholewa D, Bork S, Walenda T, Koch C, et al. Impact of individual platelet lysates on isolation and growth of human mesenchymal stromal cells. Cytotherapy. 2010;12(7):888–98.

    Article  CAS  Google Scholar 

  36. Hemeda H, Kalz J, Walenda G, Lohmann M, Wagner W. Heparin concentration is critical for cell culture with human platelet lysate. Cytotherapy. 2013;15(9):1174–81.

    Article  CAS  Google Scholar 

  37. Okita K, Matsumura Y, Sato Y, Okada A, Morizane A, Okamoto S, et al. A more efficient method to generate integration-free human iPS cells. Nat Methods. 2011;8(5):409–12.

    Article  CAS  Google Scholar 

  38. Goetzke R, Franzen J, Ostrowska A, Vogt M, Blaeser A, Klein G, et al. Does soft really matter? Differentiation of induced pluripotent stem cells into mesenchymal stromal cells is not influenced by soft hydrogels. Biomaterials. 2018;156:147–58.

    Article  CAS  Google Scholar 

  39. Willmann CA, Hemeda H, Pieper LA, Lenz M, Qin J, Joussen S, et al. To clone or not to clone? Induced pluripotent stem cells can be generated in bulk culture. PLoS One. 2013;8(5):e65324.

    Article  CAS  Google Scholar 

  40. Lenz M, Goetzke R, Schenk A, Schubert C, Veeck J, Hemeda H, et al. Epigenetic biomarker to support classification into pluripotent and non-pluripotent cells. Sci Rep. 2015;5:8973.

    Article  CAS  Google Scholar 

  41. Koch CM, Suschek CV, Lin Q, Bork S, Goergens M, Joussen S, et al. Specific age-associated DNA methylation changes in human dermal fibroblasts. PLoS One. 2011;6(2):e16679.

    Article  CAS  Google Scholar 

  42. Wagner W, Horn P, Castoldi M, Diehlmann A, Bork S, Saffrich R, et al. Replicative senescence of mesenchymal stem cells: a continuous and organized process. PLoS One. 2008;3(5):e2213.

    Article  Google Scholar 

  43. Koch CM, Reck K, Shao K, Lin Q, Joussen S, Ziegler P, et al. Pluripotent stem cells escape from senescence-associated DNA methylation changes. Genome Res. 2013;23(2):248–59.

    Article  CAS  Google Scholar 

  44. Weidner CI, Walenda T, Lin Q, Wolfler MM, Denecke B, Costa IG, et al. Hematopoietic stem and progenitor cells acquire distinct DNA-hypermethylation during in vitro culture. Sci Rep. 2013;3:3372.

    Article  Google Scholar 

  45. Fernando F, Keijser R, Henneman P, Van der Kevie-Kersemaekers AM, Mannens MM, van der Post JA, et al. The idiopathic preterm delivery methylation profile in umbilical cord blood DNA. BMC Genomics. 2015;16:736.

    Article  Google Scholar 

  46. Cancer Genome Atlas Research N, Ley TJ, Miller C, Ding L, Raphael BJ, Mungall AJ, et al. Genomic and epigenomic landscapes of adult de novo acute myeloid leukemia. N Engl J Med. 2013;368(22):2059–74.

    Article  Google Scholar 

  47. Qu Y, Lennartsson A, Gaidzik VI, Deneberg S, Karimi M, Bengtzen S, et al. Differential methylation in CN-AML preferentially targets non-CGI regions and is dictated by DNMT3A mutational status and associated with predominant hypomethylation of HOX genes. Epigenetics. 2014;9(8):1108–19.

    Article  Google Scholar 

  48. Hannum G, Guinney J, Zhao L, Zhang L, Hughes G, Sadda S, et al. Genome-wide methylation profiles reveal quantitative views of human aging rates. Mol Cell. 2013;49(2):359–67.

    Article  CAS  Google Scholar 

  49. Reinius LE, Acevedo N, Joerink M, Pershagen G, Dahlen SE, Greco D, et al. Differential DNA methylation in purified human blood cells: implications for cell lineage and studies on disease susceptibility. PLoS One. 2012;7(7):e41361.

    Article  CAS  Google Scholar 

  50. Lee ST, Xiao Y, Muench MO, Xiao J, Fomin ME, Wiencke JK, et al. A global DNA methylation and gene expression analysis of early human B-cell development reveals a demethylation signature and transcription factor network. Nucleic Acids Res. 2012;40(22):11339–51.

    Article  CAS  Google Scholar 

  51. Sandoval J, Heyn H, Moran S, Serra-Musach J, Pujana MA, Bibikova M, et al. Validation of a DNA methylation microarray for 450,000 CpG sites in the human genome. Epigenetics. 2011;6(6):692–702.

    Article  CAS  Google Scholar 

  52. Zeeberg BR, Qin H, Narasimhan S, Sunshine M, Cao H, Kane DW, et al. High-Throughput GoMiner, an 'industrial-strength' integrative gene ontology tool for interpretation of multiple-microarray experiments, with application to studies of common variable immune deficiency (CVID). BMC Bioinformatics. 2005;6:168.

    Article  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the Else Kröner-Fresenius-Stiftung (2014_A193), by the Deutsche Forschungsgemeinschaft (DFG; WA1706/8-1 and WA1706/11-1), and by the Interdisciplinary Center for Clinical Research (IZKF; O1-3 to WW and O1-4 to MZ) within the faculty of Medicine at the RWTH Aachen University.

Availability of data and materials

Raw data of DNAm profiles have been deposited at Gene Expression Omnibus (GEO) under the reference ID GSE119079.

Author information

Authors and Affiliations

Authors

Contributions

OC analyzed the data and compared epigenetic profiles. JFro and SR designed the study, conducted the experiments, and analyzed the data. JFra performed the bioinformatic analysis. SS performed the flow cytometric analysis and supported the hematopoietic differentiation. RG generated and characterized the iMSCs. MS supported the hematopoietic differentiation. MZ and WW conceived the study, participated in the data analysis, and coordinated the research. OC and WW wrote the manuscript, and all authors read, edited, and approved the final manuscript.

Corresponding author

Correspondence to Wolfgang Wagner.

Ethics declarations

Ethics approval and consent to participate

Mesenchymal stromal cells were isolated from bone marrow (caput femoris) of patients undergoing orthopedic surgery after written consent and the study was approved by the Ethic Committee of the Use of Human Subjects at the University of Aachen (permit number: EK128/09).

Consent for publication

Not applicable.

Competing interests

WW is cofounder of Cygenia GmbH that can provide service for various epigenetic signatures (www.cygenia.com). Apart from this, the authors have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Figure S1. Differentiation of iPSCs toward hematopoietic progenitor cells. (a) Phase contrast images in the course of hematopoietic differentiation cultures (on days 0, 6, 10, and 20). Scale bar = 500 μm and scale bar in inlet = 100 μm. (b) Exemplary cytospin analysis on day 22 reveals prevailing monocytic morphology. Scale bar = 10 μm. (c) Flow cytometry analysis of hematopoietic differentiation cultures on days 6 and 20. Blots are representative for three independent experiments. (d) Frequencies of CD31, CD34, CD43, CD45, and cKIT populations on day 6 and 20 of differentiation. Data represent the mean of three independent experiments ± SD. (PDF 376 kb)

Additional file 2:

Figure S2. Enrichment of hyper- and hypomethylated CpG sites in iPSC-derived HPCs. Changes in DNA methylation (delta mean β value > 0.5 or < − 0.5) in comparison (a) between iHPCs d20 and iPSCs (GSE37066) or (b) between iHPCs d20 and primary CD34+ cells from human cord blood (GSE40799). The differentially methylated CpGs were classified according to gene regions and in relation to CpG islands. Hypergeometric distribution: *P < 0.05; #P < 0.01; $P < 0.001; §P < 0.0001; ♦P < 0.00001. (PDF 344 kb)

Additional file 3:

Table S1. Differentially methylated CpGs in iPSC-derived HPCs versus iPSCs. Promoter-associated CpG sites that are either hypermethylated (220 CpG sites) or hypomethylated (1493 CpG sites; delta mean β value > 0.5 or < − 0.5) in iHPCs d20 compared to iPSCs (GSE37066) with related genes, gene groups, association to CpG islands, and mean β values of the cell types. (XLSX 119 kb)

Additional file 4:

Figure S3. Comparison of differentially methylated CpG sites across different cell types. Heatmap of DNAm levels at promoter-associated CpG sites that are either at least 50% hypo- or hypermethylated in (a) iHPCs versus iPSCs (corresponding to Fig. 1c) or in (b) iHPCs versus cord blood-derived CD34+ cells (corresponding to Fig. 2a). DNAm levels are compared between MSCs, iPSCs, iHPCs d20, and cord blood-derived CD34+ cells. The heatmaps were sorted by the mean DNAm levels in MSCs. (PDF 126 kb)

Additional file 5:

Table S2. Differentially methylated CpGs in iPSC-derived HPCs versus CD34+ cells. Promoter-associated CpG sites that are either hypermethylated (659 CpG sites) or hypomethylated (587 CpG sites; delta mean β value > 0.5 or < − 0.5) in iHPCs compared to human cord blood-derived CD34+ cells (GSE40799) with related genes, gene groups, association to CpG islands, and mean β values of the cell types. (XLSX 91 kb)

Additional file 6:

Figure S4. Differentiation of iPSCs toward MSCs. (a) Phase contrast images of iPSCs and in the course of differentiation toward iPSC-derived MSCs on day 5, 10, 20, and 30. Scale bar = 100 μm. (b) Flow cytometric analysis of iMSCs, MSCs, and iPSCs. Data is representative of three independent experiments. Autofluorescence is indicated in white. (c) iMSCs can be differentiated into adipocytes (BODIPY staining of fat droplets), osteocytes (Alizarin Red staining) and chondrocytes (Alcian Blue/PAS staining). (PDF 342 kb)

Rights and permissions

Open Access This 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Cypris, O., Frobel, J., Rai, S. et al. Tracking of epigenetic changes during hematopoietic differentiation of induced pluripotent stem cells. Clin Epigenet 11, 19 (2019). https://doi.org/10.1186/s13148-019-0617-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13148-019-0617-1

Keywords