- Open Access
DNA methylation alterations in iPSC- and hESC-derived neurons: potential implications for neurological disease modeling
Clinical Epigenetics volume 10, Article number: 13 (2018)
Genetic predisposition and epigenetic alterations are both considered to contribute to sporadic neurodegenerative diseases (NDDs) such as Parkinson’s disease (PD). Since cell reprogramming and the generation of induced pluripotent stem cells (iPSCs) are themselves associated with major epigenetic remodeling, it remains unclear to what extent iPSC-derived neurons lend themselves to model epigenetic disease-associated changes. A key question to be addressed in this context is whether iPSC-derived neurons exhibit epigenetic signatures typically observed in neurons derived from non-reprogrammed human embryonic stem cells (hESCs).
Here, we compare mature neurons derived from hESC and isogenic human iPSC generated from hESC-derived neural stem cells. Genome-wide 450 K-based DNA methylation and HT12v4 gene array expression analyses were complemented by a deep analysis of selected genes known to be involved in NDD. Our studies show that DNA methylation and gene expression patterns of isogenic hESC- and iPSC-derived neurons are markedly preserved on a genome-wide and single gene level.
Overall, iPSC-derived neurons exhibit similar DNA methylation patterns compared to isogenic hESC-derived neurons. Further studies will be required to explore whether the epigenetic patterns observed in iPSC-derived neurons correspond to those detectable in native brain neurons.
The origin and pathophysiology of sporadic neurodegenerative diseases (NDDs) such as idiopathic Parkinson’s disease (iPD) still remain enigmatic. In this regard, epigenetic alterations of genes including α-synuclein (SNCA) could contribute to the pathophysiology of this devastating disorder or the individual susceptibility .
Previously, we have identified significantly reduced DNA methylation levels of SNCA intron 1 in PD brains . In SK-N-SH cells, DNA demethylation of the SNCA intron 1—using the DNA methyltransferase inhibitor 5-Azacytidine—was associated with an increase in SNCA mRNA and α-synuclein protein expression . This is an important finding as elevated α-synuclein protein levels are directly linked to the development of PD [3, 4]. To date, these findings could not be verified in follow-up studies  and/or in native neurons derived from iPD patients compared to controls due to the inaccessibility of living disease-affected cells.
In this context, in vitro differentiated neural cells from patient-derived induced pluripotent stem cells (iPSCs) could provide a tool to study human neurons. Human iPSCs have indeed already been employed to dissect the mechanisms of several mostly inherited NDDs [6,7,8,9]. iPSC-based models of PD, mainly focusing on familial forms of the disease, displayed some interesting and relevant phenotypes . Consequently, patient-specific iPSC-derived neurons could be used to establish in vitro models for iPD allowing epigenetic and gene expression studies. Recently, a comprehensive study comprising epigenetic analyses in PD iPSC models revealed an aberrant epigenome in iPSC-derived dopaminergic neurons (DAn) of iPD and LRRK2 PD patients . These data suggest that cell-specific and presumptive disease-associated epigenetic changes can be retrieved from iPSC-derived patient-specific neurons. Nevertheless, the process of induced reprogramming is based on erasing the existing epigenetic state of the cell of origin [12, 13]. One key question emerging in the context of pluripotent stem cell-based disease modeling is whether iPSC-derived neurons display epigenetic patterns resembling those of neurons derived from non-reprogrammed human embryonic stem cells (hESCs). To address this question, we employed an isogenic stem cell system to compare DNA methylation patterns of hESC- and iPSC-derived mature neurons both generated from the same parental hESC line. Such an isogenic approach should eliminate differences due to individual genetic variation. Specifically, we assessed (i) genome-wide DNA methylation of > 485.000 CpGs in potential regulatory elements such as promoters or enhancers, (ii) gene expression of > 47.000 transcripts, and (iii) DNA methylation and gene expression changes of several NDD-related genes, in particular, SNCA.
Genome-wide and single gene DNA methylation patterns of isogenic hESC- and iPSC-derived neurons are very similar
The current study addressed the question whether iPSC-derived neurons display similar DNA methylation patterns compared to hESC-derived neurons. As a first step, we compared isogenic hESC- and iPSC-derived neurons, which allow the analysis of DNA methylation differences independently of individual genetic variations. To that end, we already differentiated thoroughly the characterized and stable populations of hESC-derived stable neural stem cells (NSCs) from the female hESC line I3 into 6-week-old neurons (Fig. 1). These hESC-derived NSCs, also referred to as long-term self-renewing neuroepithelial stem cells (lt-NES cells) differentiate predominantly into GABAergic neurons with a posterior identity corresponding to the ventral hindbrain area . We reprogrammed these NSCs to iPSC and subsequently differentiated them into three clonal populations of iPSC-derived neurons, which were compared with neurons generated from the parental hESC (Fig. 1). A genome-wide analysis was carried out using Illumina 450 K bead arrays measuring DNA methylation levels at > 450.000 CpG sites covering promoters and putative regulatory domains of all designable RefSeq genes.
A Pearson correlation analysis of gene-linked CpG methylation revealed at least a 96% concordance of iPS-NSC and iPS-Neurons compared to their hESC-derived counterparts (Fig. 2a). The interclonal variance was small; NSC and neurons derived from the three different iPSC clones showed a concordance for DNA methylation levels of at least 97% (Fig. 2a).
Differentially methylated CpGs (DMCGs) in iPSC-derived neurons compared to their hESC-derived counterparts (defined by a methylation difference > 20%) were analyzed in more detail; 8.680 (42.35%) of these DMCGs were hypomethylated, and 11.816 (57.65%) were hypermethylated in iPS- compared to hES-Neurons (Fig. 2b). Transcription start sites were not overrepresented in these DMCGs compared to the overall 450 K array annotation (Fig. 2c upper panels). More hypomethylated CpGs were found in CpG islands (CGI) and more hypermethylated CpGs in non-CGI DNA sequences (Fig. 2c lower panels). We also analyzed the standard deviation (> 0.2) of single DMCG in iPS-NSC (Fig. 2d) and iPS-Neurons (Fig. 2e) to identify an interclonal variation of DMCG. We identified only 3.051 and 3.264 interclonal DMCG sites with a standard deviation of > 0.2 for iPS-NSC and iPS-Neurons, respectively. We found a high correlation of variation between iPSC-derived NSC and neurons (Fig. 2d–f). To identify gene loci showing the most prominent alterations, we carried out an interindividual analysis (hES-Neurons vs. each iPSC-derived neuronal clone) applying a more stringent threshold of > 0.5 methylation difference (Additional file 1 A–C). In line with the previous results, we could identify only a few DMCGs (n = 1010 up to n = 1802) being differentially methylated comparing hES-Neurons to all three iPSC-derived neuronal clones (Additional file 1 A–C). The data indicate that (i) different iPSC neuronal clones exhibit almost the same amount of DMCG when compared to their hESC counterparts and (ii) that the differences are not biased toward certain genomic sites, keeping in mind that the distribution of CpGs on the 450 K array across the genome is not equally distributed per se. Using Fisher’s test, we could, however, identify a highly significant variation for methylation differences on the X chromosome in general (Additional file 1 D).
Consistent gene expression patterns of isogenic hESC- and iPSC-derived neurons
To compare gene expression patterns of hESC- and iPSC-derived neurons, we performed HT12v4 array experiments covering more than 47.000 transcripts. We tested the technical variation of the assay by hybridizing all the neuronal samples in duplicates. We observed high fidelity illustrated by correlation values of greater than 98% between each pair of replicates and a high interclonal correlation of hESC- and iPSC-derived neurons (Additional file 2 A–D). For all downstream comparisons, replicate 1 of each sample was used.
Pearson correlation analysis revealed at least a 96% concordance of gene expression profiles of iPSC- and hESC-derived neurons; interclonal variance was minimal (< 4%, Fig. 3a, Additional file 2 E, F). These data revealed that neurons generated from iPSC are very similar to their hESC-derived counterparts at the gene expression level. Hierarchical cluster analysis revealed that hESC-derived and iPSC-derived neurons cluster together, as do their parental hESC-derived and iPSC-derived NSC. These data indicate that differentiation-associated expression changes outweigh potential changes due to the reprogramming procedure (Fig. 3b). Differentially expressed genes (DEG) in iPS-Neurons were mostly associated with transcription regulation and homeobox genes. Overall, the coefficient of determination for gene expression levels in iPS-NSC and iPS-Neurons amounted to R2 = 0.75 (data not shown).
Taken together, our data indicate that isogenic hESC-derived and iPSC-derived neurons generated via a highly standardized intermediate NSC population  display similar DNA methylation and gene expression patterns. Thus, the reprogramming process does not result in major alterations of the methylation and gene expression patterns in iPSC-derived neurons. It is important to note though that the iPSCs in our experimental setting were generated from hESC-derived NSC, while patient-derived iPSCs are typically generated from fibroblasts or peripheral blood mononuclear cells. To assess how methylation levels observed in our system relate to the data obtained from neurons generated from fibroblast-derived iPSC, we performed a comparative in silico analysis including hESC-derived neural cells from another hESC line (H9) and fibroblast-derived neuronal populations from controls, iPD and LRRK2 PD patients [11, 15]. Hypothesis-free variant principal component analysis (PCA) of the 450 K data revealed that pluripotent stem cells, fibroblast-iPSC-derived neurons, H9 hESC-derived neural cells, and I3 hESC-derived neural cells clustered in distinct tiers (Fig. 4 and Additional file 3). PC2 was most likely associated with cell type, whereas the biggest variance (PC1) could not be attributed to a single factor. The segregation pattern of the different neuronal populations suggests that the individual genetic background (patient-derived vs. H9 ESC-derived vs. I3 ESC-derived) and differences in the neuronal differentiation protocols have a much higher impact on the methylation landscape than the reprogramming process itself. In line with this notion, all isogenic neurons generated from I3 ESC and I3 ESC-NSC-derived iPSC formed a tight cluster, which segregated from fibroblast- and H9 ESC-derived neurons (Fig. 4). We could also identify a weak, but significant correlation regarding PC1 and batch effects (p < 0.0002, R2 = 0.73). From our observations, we conclude that a thorough comparison of neurons generated from ESC-NSC- and fibroblast-derived iPSC would require isogenic human ESC and iPSC. Such comparative studies might be extended to iPSC-derived vs. central nervous system (CNS) tissue-derived neurons isolated from the same donor, which could in the end resolve to what extent the epigenetic landscape of reprogrammed neurons reflects that of primary human neurons.
DNA methylation changes during differentiation of isogenic hESC- and iPSC-derived neurons
To assess the acquisition of methylation patterns during neuronal differentiation in depth, we pooled the data of the hESC- and iPSC-derived NSC and compared them to the pooled data of hESC- and iPSC-derived neurons. We performed a paired t test (p < 0.001) applying a less conservative methylation difference of 10%, which identified 1.314 DMCG only (Fig. 5a, Additional file 4). Only 41 CpGs were hypomethylated, while the vast majority (1.273 CpGs) was hypermethylated in neurons compared to NSC (Fig. 5b). We did not identify an obvious enrichment for DMCG comparing gene sub-regions such as transcription start sites, promoter, and gene body regions in hES- and iPS-Neurons compared to NSC (Additional file 5), but less DMCG in CpG island regions compared to other annotations such as CpG shelfs or shores (Fig. 5c).
DNA methylation levels of SNCA intron 1 increased during neuronal differentiation in isogenic hESC- and iPSC-derived neurons
We next used MiSeq amplicon sequencing of bisulphite-treated genomic DNA (Bi-PROF; ) to validate our data obtained from the 450 K assay and to specifically assess DNA methylation at a single base level. We choose dosage-sensitive genes known to be involved in different neurological disorders, including APP (gene duplication, Alzheimer’s disease, AD), SNCA (gene duplication/triplication, Parkinson’s disease, PD) , and PMP22, which has been associated with Charcot-Marie-Tooth 1A (gene duplication) and several other CNS disorders [18, 19]. Epigenetic modification of such dosage-sensitive genes could represent a missing link between familial and sporadic forms of NDD . Furthermore, we chose MIR886, which was found to be differentially methylated in PD  and GNAS, a gene with a highly complex imprinted expression pattern  (Additional file 6). All aforementioned samples (hES-Neuron and iPS-Neuron clones 1–3) were included in the analysis. Pearson correlation coefficients demonstrated a correlation of 0.92 for individual CpGs detected by 450 K bead array and Bi-PROF, respectively. After a successful validation of single CpGs obtained from the 450 K assay, we analyzed amplicons located at the promoter regions and CpG islands of the candidate genes (Additional file 6). A comparison of hESC- and iPSC-derived neurons revealed that most regulatory regions exhibited unchanged DNA methylation patterns; the DNA methylation levels of APP, GNAS, MIR886, and SNCA promoter and intron 1 were similar comparing hESC- and iPSC-derived neurons (Fig. 6, Additional file 7). Average methylation levels over all CpGs were also comparable between hESC- and iPSC-derived NSC (Fig. 6, Additional file 7). The pattern within the SNCA intron 1 displayed higher mean methylation levels over all analyzed CpGs compared to the canonical promoter region in line with the previous analysis of bulk native brain tissue (Additional file 7) [2, 5]. These data show that DNA methylation levels at promoter regions of the selected candidate genes are similar in ESC- and iPSC-derived neurons.
DNA methylation alterations during neuronal differentiation of NSCs were observed specifically at the SNCA intron 1, but not for the other selected and analyzed NDD-associated genes. HES-NSC generated from I3 hESC and iPS-NSC generated from I3 hES-NSC-derived iPSC exhibited slightly higher, though non-significant (iPSC vs. iPS-NSC p = 0.08) DNA methylation levels of SNCA intron 1 than their parental PSC populations and hESC of the line I6; differentiation of the iPSC-derived NSC line I3 resulted in a further significant increase in DNA methylation (ESC vs. iPS-Neurons p = 0.03, iPSC vs. iPS-Neurons p = 0.02, iPS-NSC vs. iPS-Neurons p = 0.02; Additional file 8). To obtain more insight into DNA methylation changes in the SNCA intron 1 across neuronal differentiation, we analyzed NSC and neurons differentiated for 2, 4, and 6 weeks (Fig. 7a). Mean DNA methylation levels over all analyzed CpG sites of SNCA intron 1 increased during a 6-week period of neuronal differentiation in both hESC-derived (2.3-fold) and iPSC-derived neurons (3.5-fold, p = 0.016; Fig. 7a). While DNA methylation levels increased at each CpG upon neuronal differentiation, CpGs 3, 5, and 20 showed prominent changes (p = 0.001, p = 0.002, and p = 0.008, respectively; Fig. 7a). Remarkably, even at this single CpG level, methylation changes in iPSC- and hESC-derived neurons were very similar and showed a comparable pattern (Fig. 7a). Interestingly, increasing DNA methylation of SNCA intron 1 across neuronal differentiation was not associated with decreased but increased SNCA expression in both hESC- and iPSC-derived neurons (Fig. 7b).
Impact of epigenetic alterations in NDD
There is growing evidence that epigenetic alterations might play a role in NDD [1, 22]. In particular, epigenetic modifications of dosage-sensitive genes such as SNCA, a key player in PD, could be important for the phenotypic features, individual susceptibility, and the variable course of NDD . However, the exact impact of epigenetic modifications in NDD is still unknown. To date, epigenetic studies have been hampered by the inaccessibility of disease-affected cells. Previous studies on bulk brain tissue aimed at exploring whether DNA methylation changes of the SNCA gene could contribute to the development of iPD [2, 5], but the results were not conclusive. iPSC-derived neurons could represent an attractive experimental tool to study epigenetic alterations in PD . However, the reprogramming process itself is based on a major epigenetic “reset” with the resulting iPSC exhibiting epigenetic patterns closely related to those of ESC [24, 25]. This is presumably associated with a loss of both, epigenetic somatic memory (i.e., tissue-specific epigenetic patterns of the donor cells used for iPSC generation) and age-associated epigenetic signatures [26,27,28]. The latter is particularly relevant as NDD typically occur in late age. A global epigenetic rearrangement during the generation of iPSC could largely eliminate epigenetic alterations resulting from prior noxious exposures. IPSC might thus not fully reflect acquired patient-specific epigenetic changes. On the other hand, epigenetic alterations are not only due to environmental impacts but also reflect changes in cis, i.e., changes in methylation patterns as a downstream result of genetic variants .
Isogenic stem cell systems
It is evident that successful and reliable in vitro modeling of epigenetic changes depends on the faithful replication of the epigenetic signatures observed in non-reprogrammed cells. The comparative assessment of epigenetic patterns in reprogrammed versus non-reprogrammed cells is impeded by the interindividual genetic variability. Isogenic systems can be used to resolve this problem. Indeed, while this manuscript was in preparation, several studies started to address the issue of reprogramming-associated epigenetic differences and their biological significance in isogenic stem cell systems [24, 30, 31]. Teichroeb et al. compared the transcriptome of isogenic hESC (male hESC line H9) and iPSC using an array-based analysis (HT-12v3) . They found a striking similarity (99.6%) of the transcriptome of isogenic hESC and iPSC . These results were supported by Mallon et al. who compared array-based DNA methylation (Illumina Meth.27) and gene expression (Agilent 4x44k One Color) patterns of isogenic hESC (male ESC-line H1) and iPSC . They found no major differences in global gene expression patterns and an accordance of DNA methylation patterns of 98–99% . Choi et al. compared isogenic iPSC and hESC derived from two male hESC lines (HUES2 and HUES3) . Employing global expression profiling by RNA sequencing, they found that only 49 genes were differentially expressed in these isogenic cell populations . DNA methylation pattern differences (reduced representation bisulphite sequencing, RRBS) accounted for only approximately 4% .
We were interested in extending these comparative studies from pluripotent stem cells to their neuronal derivatives. To that end, we compared DNA methylation patterns of reprogrammed and non-reprogrammed mature neurons derived from the same hESC line. Employing an isogenic human neural stem cell system, we found that the genome-wide DNA methylations (and gene expression profiles) in hESC- and iPSC-derived neurons are remarkably similar. Furthermore, iPSC-derived cells displayed only minor interclonal variability. These results are in line with several studies on epigenetic and gene expression patterns in isogenic pluripotent stem cell populations [24, 30, 31]. The need for an isogenic system is also underlined by the results of our comparative in silico analysis including hESC-derived neural cells from another hESC line (H9) and fibroblast-derived neuronal populations from controls, iPD, and LRRK2 PD patients [11, 15]. These data indicate that the individual genetic background (patient-derived vs. H9 ESC-derived vs. I3 ESC-derived) and differences in the neuronal differentiation protocols have a much higher impact on the methylation landscape than the reprogramming process itself.
Looking at five NDD-related candidate genes within our isogenic system, we also observed a high similarity between hESC- and iPSC-derived neural cell methylation patterns at the single base level. Furthermore, we detected a high DNA methylation similarity of hESC- and iPSC-derived neurons concerning the imprinted locus GNAS. Interestingly, we found a distinct increase of DNA methylation in SNCA intron 1 upon neuronal differentiation up to 6 weeks in vitro. ESC- and iPSC-derived neurons acquired a highly comparable pattern during neuronal differentiation. However, this differentiation-associated increase in methylation was paralleled by increased expression of SNCA. It remains currently unclear whether DNA methylation changes upon neuronal differentiation of iPSC- and hESC-derived neurons indeed influence SNCA gene expression levels.
Limitations of the current study and future prospects
The current study focused—as a proof of principle—on evaluating whether iPSC-derived neurons from one hESC-line exhibit similar DNA methylation patterns compared to isogenic hESC-derived neurons. A clear limitation of the current study is the lack of a profound DNA methylation analysis comparing pluripotent and differentiated cells, except for the SNCA locus. The current study did also not address the question whether certain CpG methylation patterns were erased or preserved during the reprogramming process which has to be determined using additional isogenic pluripotent, hESC- and iPSC-derived neural cells. Moreover, the complex issue of variation in genomic imprinting and X chromosome inactivation, too, would require more detailed analyses including, e.g., assessment of allele-specific expression, H3K27me accumulation, and XIST RNA levels, which are going beyond the scope of this study. Despite these limitations, our results suggest that the reprogramming process itself has no major impact on methylation patterns in iPSC-derived neurons. Yet, it remains to be determined how in vitro methylation patterns of hESC- and iPSC-derived neurons correlate to those of neurons derived from the native brain tissue. This is particularly relevant in the light of recent findings revealing substantial transcriptional differences between iPSC-derived and primary midbrain dopamine neurons .
iPSC-derived neurons represent a valuable tool for modeling neurological disorders in vitro. While a number of studies have shown that the epigenetic status of iPSC closely mimics that of non-reprogrammed pluripotent stem cells (i.e., hESC), it had remained largely unclear whether neurons differentiated from iPSC adopt methylation patterns similar to those of neurons derived from non-reprogrammed cells. Using an isogenic system, this study demonstrates that ESC- and iPSC-derived neurons are highly similar with respect to DNA methylation and gene expression. While further studies are required to explore how DNA methylation in iPSC-derived neurons correlates to methylation patterns in primary neurons derived from the CNS tissue, our results suggest that the reprogramming process itself does not confound methylation patterns in iPSC-derived neurons. Thus, iPSC-derived neurons should represent a suitable tool for studying the impact of genetic variants on the epigenome of neurons.
HESCs (line I3) were maintained and differentiated to rosette-forming neural stem cells (NSCs) according to standard protocols and differentiated to neurons as described previously (Additional file 9) . iPSCs were obtained and validated according to established protocols. In brief, NSCs were transduced with lentiviral TetON-system in order to generate a stable inducible cell line harboring the reprogramming factors OCT4 and KLF4 (pLVXTP-Tet-On (Clontech), FUW-OCT4, and FUW-KLF4 (Addgene)). Cells were induced to reprogram by the addition of 1 μg/ml doxycycline (Sigma-Aldrich) and cultured on irradiated mouse embryonic fibroblasts in Knockout Dulbecco’s modified Eagle’s medium containing 20% serum replacement, 1% non-essential amino acids, 1 mM L-glutamine, 0.1 mM ß-mercaptoethanol, and 4 ng/ml FGF2. Doxycycline was withdrawn upon colony formation, and doxycycline-independent colonies were mechanically isolated and further propagated. iPSC lines were differentiated to NSC and cultured as described . Terminal differentiation of NSC was performed on Geltrex-coated dishes (Life technologies) in DMEM/F12 and Neurobasal (B27 supplement 1:50, penicillin/streptomycin 1:100, cAMP 300 ng/ml, and 1.6 g/l glucose) mixed at a 1:1 ratio. Neurons were cultured for 2, 4, and 6 weeks, and media were changed every second day (Additional file 9).
Extraction of DNA and RNA
Genomic DNA was isolated using the DNeasy Blood & Tissue Kit (Qiagen) according to the manufacturer’s protocol. RNA was isolated using the RNeasy Mini Kit (Qiagen) according to the manufacturer’s protocol.
cDNA synthesis and semiquantitative PCR
One microliter of total RNA was used for cDNA synthesis with the iScript cDNA Synthesis Kit (Bio-Rad Laboratories) according to the manufacturer’s instructions. A final concentration of 150 ng/μl of cDNA was used for real-time PCR with Power SYBR Green PCR Master Mix (Applied Biosystems) according to the manufacturer’s instructions. Possible contamination with genomic DNA was tested using samples processed without reverse transcriptase. Gene expression levels were normalized to human β-actin using the following primers (final concentration 10 μM): SNCA forward 5′-cgacgacagtgtggtgtaaag-3′ and SNCA reverse 5′-aaatgttggaggagcagtgg-3′.
Bisulphite treatment and single CpG analysis
Bisulphite treatment of genomic DNA was performed with the EZ DNA Methylation-Gold Kit (Zymo Research) according to the manufacturer’s protocol. Amplicons were generated using region-specific primers (final concentration 10 μM each) with the recommended adaptors (Additional file 6). Purified PCR products were pooled in an equimolar ratio and sequenced after cluster formation on a MiSeq instrument benchtop sequencer with the sequencing-by-synthesis technology  according to the manufacturer’s protocol. Runs were set for “Generate FASTQ only” workflow in Illumina Experiment Manager. ESC and iPSC were analyzed as described elsewhere .
Infinium HumanMethylation450 BeadChip
Genome-scale DNA methylation profiles were generated using Illumina’s Infinium HumanMethylation450 Beadchip assay (Illumina). The assay allows the determination of DNA methylation levels at > 450.000 CpG sites covering promoters and putative regulatory domains of all designable RefSeq genes. The Infinium Methylation Assay was performed according to the manufacturer’s instructions. The fluorescently stained chips were imaged using an Illumina HiScan scanner.
Infinium HT12v4 BeadChip
Genome-wide gene expression profiles were generated using Illumina’s Infinium HT12v4 BeadChip assay (Illumina). The assay allows the determination of gene expression of > 47.000 probes covering designable RefSeq genes. The Infinium Expression Assay was performed according to the manufacturer’s instructions. The scanning was performed using an Illumina HiScan. Data from these images were loaded into Illumina’s GenomeStudio software suite (GS) and exported as a GS export file.
Genome-wide methylation raw data were pre-processed using R statistics software suites minfi and RnBeads  to extract the data, subtract the background, and to normalize the data using internal controls present on the chips. Only CpGs with a detection p value < 0.01 in all samples were included (n = 484,932 of 485,577). All samples were analyzed as individual samples (n = 1). For genome-wide expression analysis, raw data was background subtracted and quantile normalized using R and the lumi software suite. Only probes with a detection p value < 0.01 in all samples were included (n = 20,586 of 39,050). All NSCs were analyzed as individual samples (n = 1); neuronal samples were analyzed in duplicates (n = 2). Further statistical analysis was performed with R statistics software or SPSS version 22. p values were based on a (paired) two-sided t test. Associations between variables (expression and DNA methylation) were measured by Pearson’s rank correlation.
For amplicon-based single gene CpG methylation analysis, sequencing reads obtained from the Illumina MiSeq were processed and assigned to the reference sequence and respective sample using specific tags in the universal portion of the adaptors. Methylation level and patterns were assessed using multiple sequence alignment using BiQ Analyzer HT . All samples were analyzed as individual samples (n = 1) excluding all reads with a maximum fraction of unrecognized CpG sites of 10%. Data are shown as mean ± SD.
Probes identified in differential methylation/expression analyses were mapped to corresponding gene IDs according to the Illumina annotation. GO enrichment analysis was performed using the GOrilla tool (http://cbl-gorilla.cs.technion.ac.il/).
Publicly available data
Published and pre-processed 450 K data from Fernandez-Santiago et al.  were downloaded from GEO accession number GSE51921. Non-publicly available data of non-dopaminergic neurons was kindly provided by R. Fernandez-Santiago. Published and pre-processed 450 K data from Kim et al.  were downloaded from GEO accession number GSE38217. For the comparison of ESC-derived and fibroblast-iPSC-derived neurons, we included all samples with a detection p value of < 0.01 and excluded all CpGs having some potential SNP bias (hES-NSC and hES-Neurons n = 4, ESC line H9, NPC line H9 and DAn line H9 n = 2, iPSC iPD and LRRK PD n = 2, iPSC control n = 3, DAn control and LRRK2 PD n = 4, DAn iPD n = 6, non-dopaminergic neurons control, LRRK2 PD and iPD n = 3). Principal component analysis was performed in R using the “prcomp” function. Obtained eigenvalues for principal components 1 and 2 were used for generation of PCA plot. The fibroblast-iPSC-derived non-dopaminergic neurons from the study of Fernández-Santiago et al. were cultured in DMEMF12/Neurobasal medium supplemented with 0.5xB27, 0.5xN2, and 2 mM GlutaMax and penicillin-streptomycin . HESC- and iPSC-derived non-dopaminergic neurons from the isogenic stem cell model were cultured with slightly higher doses of N2 and B27 (1×) and in the presence of cyclic adenosine monophosphate (cAMP, 0.1×). Dopaminergic neurons (DAn) from the study of Fernández-Santiago et al. were generated using lentiviral transduction of LIM homeobox transcription factor alpha (LMX1A), co-cultured with a mouse with mouse PA6 feeding cells overexpressing sonic hedgehog (SHH) and cultured in the presence of fibroblast growth factor 8 (FGF8) . Neurons were cultured for 30 days, and DAn exhibited a dopaminergic marker expression ranging from 44 to 67% (tyrosine hydroxylase (TH)/β-Tubulin III (TUJ1) ratio) . Dopaminergic neural precursor cells from the study of Kim et al. (ESC line H9) were generated culturing hES colonies on irradiated PA6 or MA5 cells stably over-expressing SHH and cultured in insulin/transferrin/selenium/ascorbic acid/basic fibroblast growth factor (bFGF) medium . Terminal differentiation of dopaminergic NSC into DAn was induced in the absence of bFGF but in the presence of brain-derived neurotrophic factor (BDNF), glial cell-derived neurotrophic factor (GDNF), and cAMP . DAn were cultured for 42 days and exhibited a purity of 43% .
Landgrave-Gomez J, Mercado-Gomez O, Guevara-Guzman R. Epigenetic mechanisms in neurological and neurodegenerative diseases. Front Cell Neurosci. 2015;9:1–11.
Jowaed A, Schmitt I, Kaut O, Wüllner U. Methylation regulates alpha-synuclein expression and is decreased in Parkinson’s disease patients’ brains. J Neurosci. 2010;30:6355–9.
Ibáñez P, Bonnet A-M, Débarges B, Lohmann E, Tison F, Pollak P, et al. Causal relation between alpha-synuclein gene duplication and familial Parkinson’s disease. Lancet. 2004;364:1169–71.
Singleton AB, Farrer M, Johnson J, Singleton A, Hague S, Kachergus J, et al. Alpha-Synuclein locus triplication causes Parkinson’s disease. Science. 2003;302:841.
de Boni L, Tierling S, Roeber S, Walter J, Giese A, Kretzschmar H a. Next-generation sequencing reveals regional differences of the α-synuclein methylation state independent of Lewy body disease. NeuroMolecular Med. 2011;13:310–20.
Koch P, Breuer P, Peitz M, Jungverdorben J, Kesavan J, Poppe D, et al. Excitation-induced ataxin-3 aggregation in neurons from patients with Machado-Joseph disease. Nature. 2011;480:543–6.
Cooper O, Seo H, Andrabi S, Guardia-Laguarta C, Graziotto J, Sundberg M, et al. Pharmacological rescue of mitochondrial deficits in iPSC-derived neural cells from patients with familial Parkinson’s disease. Sci Transl Med. 2012;4:141ra90.
Kondo T, Asai M, Tsukita K, Kutoku Y, Ohsawa Y, Sunada Y, et al. Modeling Alzheimer’s disease with iPSCs reveals stress phenotypes associated with intracellular Aβ and differential drug responsiveness. Cell Stem Cell. 2013;12:487–96.
Donnelly CJ, Zhang P-W, Pham JT, Heusler AR, Mistry NA, Vidensky S, et al. RNA toxicity from the ALS/FTD C9ORF72 expansion is mitigated by antisense intervention. Neuron. 2013;80:415–28.
Jacobs BM. Stemming the hype: what can we learn from iPSC models of Parkinson’s disease and how can we learn it? J Parkinsons Dis. 2014;4:15–27.
Fernandez-Santiago R, Carballo-Carbajal I, Castellano G, Torrent R, Richaud Y, Sanchez-Danes A, et al. Aberrant epigenome in iPSC-derived dopaminergic neurons from Parkinson’s disease patients. EMBO Mol Med. 2015;7:1529–46.
Medvedev SP, Pokushalov E a, Zakian SM. Epigenetics of pluripotent cells. Acta Nat. 2012;4:28–46.
Bilic J, Izpisua Belmonte JC. Concise review: induced pluripotent stem cells versus embryonic stem cells: close enough or yet too far apart? Stem Cells. 2012;30:33–41.
Koch P, Opitz T, Steinbeck J a, Ladewig J, Brüstle O. A rosette-type, self-renewing human ES cell-derived neural stem cell with potential for in vitro instruction and synaptic integration. Proc Natl Acad Sci U S A. 2009;106:3225–30.
Kim M, Park Y-K, Kang T-W, Lee S-H, Rhee Y-H, Park J-L, et al. Dynamic changes in DNA methylation and hydroxymethylation when hES cells undergo differentiation toward a neuronal lineage. Hum Mol Genet. 2014;23:657–67.
Gries J, Schumacher D, Arand J, Lutsik P, Markelova MR, Fichtner I, et al. Bi-PROF: bisulfite profiling of target regions using 454 GS FLX Titanium technology. Epigenetics. 2013;8:765–71.
Conrad B, Antonarakis SE. Gene duplication: a drive for phenotypic diversity and cause of human disease. Annu Rev Genomics Hum Genet. 2007;8:17–35.
Dracheva S, Davis KL, Chin B, Woo D a, Schmeidler J, Haroutunian V. Myelin-associated mRNA and protein expression deficits in the anterior cingulate cortex and hippocampus in elderly schizophrenia patients. Neurobiol Dis. 2006;21:531–40.
Ohsawa Y, Murakami T, Miyazaki Y, Shirabe T, Sunada Y. Peripheral myelin protein 22 is expressed in human central nervous system. J Neurol Sci. 2006;247:11–5.
Masliah E, Dumaop W, Galasko D, Desplats P. Distinctive patterns of DNA methylation associated with Parkinson disease: identification of concordant epigenetic changes in brain and peripheral blood leukocytes. Epigenetics. 2013;8:1030–8.
Turan S, Bastepe M. The GNAS complex locus and human diseases associated with loss-of-function mutations or epimutations within this imprinted gene. Horm Res Paediatr. 2013;80:229–41.
Urdinguio RG, Sanchez-Mut JV, Esteller M. Epigenetic mechanisms in neurological diseases: genes, syndromes, and therapies. Lancet Neurol. 2009;8:1056–72.
Torrent R, De Angelis Rigotti F, Dell’Era P, Memo M, Raya A, Consiglio A. Using iPS cells toward the understanding of Parkinson’s disease. J Clin Med. 2015;4:548–66.
Choi J, Lee S, Mallard W, Clement K, Tagliazucchi GM, Lim H, et al. A comparison of genetically matched cell lines reveals the equivalence of human iPSCs and ESCs. Nat Biotechnol. 2015;33:1173–81.
Guenther MG, Frampton GM, Soldner F, Hockemeyer D, Mitalipova M, Jaenisch R, et al. Chromatin structure and gene expression programs of human embryonic and induced pluripotent stem cells. Cell Stem Cell. 2010;7:249–57.
Lo Sardo V, Ferguson W, Erikson GA, Topol EJ, Baldwin KK, Torkamani A. Influence of donor age on induced pluripotent stem cells. Nat Biotechnol. 2016;35:1–8.
Hochedlinger K, Jaenisch R. Induced pluripotency and epigenetic reprogramming. Cold Spring Harb Perspect Biol. 2015;7:a019448.
Li M, Belmonte JCI. Ground rules of the pluripotency gene regulatory network. Nat Rev Genet. 2017;18:180–91.
McDonald KL, Rapkins RW, Olivier J, Zhao L, Nozue K, Lu D, et al. The T genotype of the MGMT C>T (rs16906252) enhancer single-nucleotide polymorphism (SNP) is associated with promoter methylation and longer survival in glioblastoma patients. Eur J Cancer. 2013;49:360–8.
Teichroeb JH, Betts DH, Vaziri H. Suppression of the imprinted gene NNAT and X-chromosome gene activation in isogenic human iPS cells. PLoS One. 2011;6:e23436.
Mallon BS, Hamilton RS, Kozhich O a, Johnson KR, Fann YC, Rao MS, et al. Comparison of the molecular profiles of human embryonic and induced pluripotent stem cells of isogenic origin. Stem Cell Res. 2014;12:376–86.
Xia N, Zhang P, Fang F, Wang Z, Rothstein M, Angulo B, et al. Transcriptional comparison of human induced and primary midbrain dopaminergic neurons. Sci Rep. 2016;6:20270.
Schmitt I, Kaut O, Khazneh H, de Boni L, Ahmad A, Berg D, et al. L-dopa increases α-synuclein DNA methylation in Parkinson’s disease patients in vivo and in vitro. Mov Disord. 2015;13:1794–801.
Assenov Y, Müller F, Lutsik P, Walter J, Lengauer T, Bock C. Comprehensive analysis of DNA methylation data with RnBeads. Nat Methods. 2014;11:1138–40.
Lutsik P, Feuerbach L, Arand J, Lengauer T, Walter J, Bock C. BiQ Analyzer HT: locus-specific analysis of DNA methylation by high-throughput bisulfite sequencing. Nucleic Acids Res. 2011;39:W551–6.
We thank Jasmin Kirch, Katharina Faber, and Christina Lo Porto for their assistance in MiSeq analyis. We thank Dr. R. Fernández-Santiago for the provision of the DNA methylation datasets on the non-dopaminergic neurons used in his study .
The study was supported by the BMBF project 01-KU1216F (DEEP, GG); the BONFOR program of the University Clinic of Bonn (LdB); the German Center for Neurodegenerative Diseases (Deutsches Zentrum für neurodegenerative Erkrankungen, DZNE); the German Federal Ministry of Education and Research (BMBF grant 01GNO813; CH, PK), within the framework of the e:Med research funding concept (BMBF grant # 01ZX1314A; OB); and by the BMBF/ANR through the EpiPD (Epigenomics of Parkinson’s disease) project, under the auspices of the bilateral Epigenomics of Common and Age-related Diseases Programme (grant no. 01KU1403B; UW).
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its additional files.
Ethics approval and consent to participate
Approval for the derivation of pluripotent cell lines for disease and therapeutic research. Ethikkommission, Medizinische Fakultät Bonn, Universität Bonn, serial no. 275/08.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
A–C Interindividual analysis of DMCG comparing hES-Neurons to each iPSC-derived neuronal clone. D X chromosomal-based DMCG analysis using Fisher’s test. Methylation deltas were calculated as indicated and CpGs above applying a moderate (0.2) or high (0.5) threshold. (PDF 437 kb)
A–D High correlation of technical replicates (hES-Neuron and iPS-Neuron clones 1–3) on the HT12v4 gene expression arrays. AU, arbitrary unit; Rep, replicate. E Highly similar expression profile comparing hES-NSC and iPS-NSC. AU, arbitrary unit. F Gene expression correlation profiles comparing hES-Neurons and iPS-Neurons. AU, arbitrary unit. (PDF 202 kb)
Comparison of mean Pearson DNA methylation correlation coefficients. Comparisons were made including all ESC- and iPSC-derived NSC and neurons from the human isogenic stem cell system (isogenic hESC- and iPSC-derived NSC and neurons), all samples from the study of Fernández-Santiago et al.  except fibroblasts, and all samples from the study of Kim et al. . iPSC, induced pluripotent stem cell; LRRK2, leucin-rich repeat kinase 2; PD, Parkinson’s disease; iPD, idiopathic PD; DAn, dopaminergic neuron; n, number; NSC, neural stem cell; hES(C), human embryonic stem (cell); I3, ESC-line I3; H9, ESC-line H9. (PDF 49 kb)
Table of the top 20 differentially methylated CpGs (DMCG) in hESC- and iPSC-derived neurons (grouped together) compared to hESC- and iPSC-derived NSC (grouped together). Data was based on a paired t test (p < 0.001) applying a minimum methylation difference of 10%. Chr., chromosome. (PDF 47 kb)
Mean DNA methylation levels of different genomic regions based on the genome-wide 450 K array data [AU ± SD]. AU, arbitrary unit; SD, standard deviation; hES, human embryonic stem cell; NSC, neural stem cell; iPS, induced pluripotent stem cell. (PDF 39 kb)
Table of analyzed single genes and primers. (PDF 40 kb)
Mean DNA methylation levels (%) of single genes (APP, GNAS, MIR886, PMP22, SNCA promoter, SNCA intron 1) of hESC-derived and iPSC-derived NSC and neurons. Analysis is based on individual samples. Data is shown as mean ± SD. (PDF 40 kb)
DNA methylation levels of SNCA intron 1 in pluripotent stem cells (ESC line I3 and I6 (n = 1 each), iPSC line I3 (n = 3), ES- and iPS-NSC (line I3, n = 1 and n = 3, respectively), and ES- and iPS-Neurons (line I3, n = 1 and n = 3, respectively). ESC versus iPS-Neurons p = 0.03, iPSC vs. iPS-Neurons p = 0.02, and iPS-NSC vs. iPS-Neurons p = 0.02. (PDF 41 kb)
A–C Immunocytochemical analysis of iPSC generated from hESC-derived NSC staining positive for the pluripotency-associated markers Oct4, Tra-1-60, and Tra-1-81. D–F Upon differentiation, these cells give rise to all three germ layers and express appropriate markers for endoderm (D; AFP), mesoderm (E, SMA), and ectoderm (F; TUBB3). G, H iPSC-derived NSC expressing the early neuroectodermal markers nestin, Dach1, ZO1, and PLZF. I, J 6-week-old iPSC-derived neurons staining for the neuronal marker MAP2 (I) and the neurotransmitter GABA (J). Nuclei are counterstained with DAPI. Scale bars: A–F, 75 μm; G–J, 10 μm. (PDF 3601 kb)