Stable DNMT3L overexpression in SH-SY5Y neurons recreates a facet of the genome-wide Down syndrome DNA methylation signature
Epigenetics & Chromatin volume 14, Article number: 13 (2021)
Down syndrome (DS) is characterized by a genome-wide profile of differential DNA methylation that is skewed towards hypermethylation in most tissues, including brain, and includes pan-tissue differential methylation. The molecular mechanisms involve the overexpression of genes related to DNA methylation on chromosome 21. Here, we stably overexpressed the chromosome 21 gene DNA methyltransferase 3L (DNMT3L) in the human SH-SY5Y neuroblastoma cell line and assayed DNA methylation at over 26 million CpGs by whole genome bisulfite sequencing (WGBS) at three different developmental phases (undifferentiated, differentiating, and differentiated).
DNMT3L overexpression resulted in global CpG and CpG island hypermethylation as well as thousands of differentially methylated regions (DMRs). The DNMT3L DMRs were skewed towards hypermethylation and mapped to genes involved in neurodevelopment, cellular signaling, and gene regulation. Consensus DNMT3L DMRs showed that cell lines clustered by genotype and then differentiation phase, demonstrating sets of common genes affected across neuronal differentiation. The hypermethylated DNMT3L DMRs from all pairwise comparisons were enriched for regions of bivalent chromatin marked by H3K4me3 as well as differentially methylated sites from previous DS studies of diverse tissues. In contrast, the hypomethylated DNMT3L DMRs from all pairwise comparisons displayed a tissue-specific profile enriched for regions of heterochromatin marked by H3K9me3 during embryonic development.
Taken together, these results support a mechanism whereby regions of bivalent chromatin that lose H3K4me3 during neuronal differentiation are targeted by excess DNMT3L and become hypermethylated. Overall, these findings demonstrate that DNMT3L overexpression during neurodevelopment recreates a facet of the genome-wide DS DNA methylation signature by targeting known genes and gene clusters that display pan-tissue differential methylation in DS.
DS is the leading genetic cause of intellectual disability and results from trisomy 21 . However, genes outside of chromosome 21 are also altered in DS and differences in gene expression and DNA methylation are observed across the entire genome. Most DS tissues exhibit differentially methylated sites that tend to be hypermethylated when compared to typically developing controls . Furthermore, some genes show pan-tissue differential methylation, which is also skewed heavily towards hypermethylation. Mechanistically, there are a number of genes located on chromosome 21 that belong to pathways related to DNA methylation and have the potential to result in the hypermethylation pattern observed in DS . Functional experimentation into the cause of hypermethylation has demonstrated a key role of the chromosome 21 encoded DNA methyltransferase DNMT3L at select genes . While DNMT3L is catalytically inactive, it is a regulatory factor that binds to and stimulates the de novo methyltransferases DNMT3A and DNMT3B [4,5,6].
Structurally, DNMT3L and DNMT3A form elongated heterotetramers through their C-terminal domains (DNMT3L–DNMT3A–DNMT3A–DNMT3L), and this complex multimerizes on DNA to form nucleoprotein filaments that spread DNA methylation over larger regions, such as CpG islands [7,8,9]. Members of the de novo methyltransferase family (DNMT3A,B,L) contain an ATRX–DNMT3–DNMT3L (ADD) domain that binds to the unmodified histone H3 tail (H3K4me0), to localize DNA methylation to previously unmethylated regions, and this binding is inhibited by methylation at lysine 4 of H3 (H3K4me3) [10,11,12].
In mouse embryonic stem cells (mESCs), DNMT3L has been reported to be either a positive or negative regulator of DNA methylation depending on genomic context ; however, another report found it to only function as a positive regulator . This difference appears to be due to the mESCs cell lines used as well as the methodologies used to assay DNA methylation, since the studies examined different regions of the genome at different levels of resolution. Experiments using transient DNMT3L overexpression in mESCs and somatic cells have demonstrated a role for DNMT3L in assembling a repressive chromatin modifying complex that silences retroviral sequences through de novo DNA methylation and methylation-independent mechanisms . Finally, DNMT3L is also required for the establishment of parental imprinting patterns in developing gametes [16, 17].
While DNMT3L is a prime candidate for the genome-wide DS DNA hypermethylation profiles, the genome-wide effect of DNMT3L overexpression has not yet been characterized. Previously, we profiled post-mortem DS brain and compared to matched controls using WGBS. We observed a genome-wide wide impact on DNA methylation profiles, where ~ 75% of the differentially methylated regions (DMRs) were hypermethylated and mapped to genes involved in one-carbon metabolism, membrane transport, and neurotransmission. Given these observations, we sought to characterize the effect of DNMT3L overexpression on different phases of neuronal development by assaying differentiating SH-SY5Y cells.
SH-SY5Y cells are a human neuronal cell line that have been subcloned three times from the SK-N-SH cell line, which was derived from metastatic cells of a 4-year-old female with neuroblastoma. Unlike many other cancer lines, the SH-SY5Y cell line has a stable karyotype with a dominant ploidy of 2, which includes trisomy 1q and 7, gains on 2p and 17q, and losses on 14q and 22q [18,19,20]. Similar to other cancer cell lines, SH-SY5Y cells have a relatively hypomethylated genome when compared to neurons from human post-mortem brain samples [21, 22]. SH-SY5Y cells display a mixed morphology that is primarily neuroblast-like with short processes; however, there is also a smaller proportion of epithelial-like cells, which likely reflects the cell lines multipotency and an origin from the neural crest. SH-SY5Y cells can be consistently differentiated from a neuroblast-like state into a homogenous population of mature dopaminergic neurons in 18 days . The differentiation protocol is a three-phased reprogramming approach, which involves gradual serum deprivation and the addition of retinoic acid during the first two phases and then a transition to an extracellular matrix where the cells are given a neuronal media that lacks serum. The neuron-like cells become dependent on neurotrophic factors, such as BDNF, for their survival and the epithelial-like cells, which require serum, do not survive. Notably, during this differentiation process, serial splitting via a brief incubation with low strength trypsin allows for the selection of the less adherent neuron-like cells (for the next phase or harvest) and ultimately results in a morphologically pure culture of mature dopaminergic neurons in the final phase. Thus, SH-SY5Y cells serve as an ideal model system to investigate the methylome-wide effects of DNMT3L overexpression on undifferentiated, differentiating, and differentiated human neurons.
Here, we interrogated the contribution of DNMT3L, a chromosome 21 encoded gene, in DS-associated genome-wide changes in DNA methylation. Stable overexpression of DNMT3L in human SH-SY5Y cells provided a model system to test the effects of this known methylation regulator on genome-wide DNA methylation across neuronal differentiation. We found that DNMT3L overexpression recreated a facet of the genome-wide DNA profile observed in DS tissues, including brain. Interestingly, chromatin state analyses revealed that the DNMT3L-induced hypermethylated DMRs were enriched within bivalent domains, whereas the hypomethylated DMRs were enriched within heterochromatin. These results indicate that an increase in DNMT3L copy number during neuronal differentiation leads to chromatin-specific effects on DNA methylation.
Stable DNMT3L Overexpression in SH-SY5Y Cells
SH-SY5Y cells were transfected with either myc-tagged human DNMT3L  and enhanced GFP (eGFP) via a piggyBac transposase  construct or a matched control construct without DNMT3L (Fig. 1a). eGFP-positive cells were selected by flow cytometry and then cultured again to enrich for a cell line with stable DNMT3L overexpression. For the experimental investigations, three developmental phases were assayed: 1) undifferentiated cells (growth), 2) the first phase of differentiation (phase 1), and 3) the final phase of differentiation (phase 3) that represents the pure and mature neuronal culture (Fig. 1b). When compared to the growth phase, phase 1 neurons displayed a relatively elongated cell body with an increased number of longer neurites with neuronal growth cones that were starting to form synapses with their neighboring cells. Notably, phase 1 neurons are the most commonly assayed form of differentiated SH-SY5Y cells in the literature. Phase 2 is a relatively brief intermediate phase for the final selection that occurs in phase 3 and, thus, was not examined. After being plated as a monolayer, phase 3 neurons migrated to form large spheres of highly interconnected neurons, which were similar to neurospheres; however, they are were attached to an extracellular matrix and produced a three-dimensional network with a larger-scale organization that consists of long tracts of bundled processes. Notably, the eGFP-only control cells displayed visibly higher fluorescence levels than those with DNMT3L + eGFP, a result that was observed in all cell lines. Endogenous DNMT3L expression was confirmed by Western blot in all transfected cell lines and compared to a non-transfected cell line (WT). As expected, DNMT3L expression was much higher in the transgenic DNMT3L cell lines and only low basal levels were detected in the WT and eGFP-only control cell lines (Fig. 1c).
Stable DNMT3L overexpression results in global and CpG island hypermethylation
WGBS libraries were constructed using DNA isolated from all cell lines and sequenced to > 10 × coverage, resulting in > 26 million assayed CpGs. There was a significant (p = 0.016) increase in global CpG methylation (Fig. 2a) and a significant (p = 0.004) increase in global CpG island methylation in cells with DNMT3L overexpression (3L) when compared to vector-only (eGFP) controls. Cell developmental phase, herein referred to as “phase”, also had a significant (p = 0.004) effect on global CpG methylation and a significant (p = 0.030) effect on global CpG island methylation. In the growth phase, DNMT3L cells showed 69.70% global methylation and 34.11% CpG island methylation, while the eGFP cells showed 69.43 and 33.82%, respectively. In phase 1 of differentiation, DNMT3L cells showed 70.38% global methylation and 34.39% CpG island methylation, while eGFP cells showed 69.93 and 33.89%, respectively. In phase 3 of differentiation, the DNMT3L cells showed 69.68% global methylation and 33.95% CpG island methylation, while eGFP cells showed 69.33 and 33.52%, respectively. Notably, both global CpG methylation and global CpG island methylation levels were the highest in phase 1 and the lowest in phase 3. Principal component analysis (PCA) revealed that the cells clustered together by phase for both the 20 Kb window and single CpG approaches (Fig. 2b), but clustered by genotype for the CpG island approach (Fig. 2c).
DNMT3L DMRs map to genes related to neurodevelopment, cellular signaling, and gene regulation
To examine the specific genetic loci altered by DNMT3L overexpression during neuronal differentiation, differentially methylated regions (DMRs) were identified using pairwise genotype comparisons of all three phases. The replicate cell lines were utilized as biological replicates and batch was directly adjusted for in the analyses. The significant (permutation p < 0.05) DMRs from each pairwise phase comparison were skewed towards hypermethylation (Additional File 1: Supplementary Table 1). In the growth phase, 62% of the 6746 DMRs were hypermethylated, as well as 71% of the 5954 DMRs in phase 1, and 68% of the 6389 DMRs in phase 3 (Fig. 3). DMRs were mapped to genes and gene ontology (GO) terms were slimmed to identify the least dispensable significant (p < 0.05) terms. Across all phases, when examining the effect of DNMT3L overexpression, these slimmed GO terms represented different neurodevelopmental processes, including synapse, signaling, and biological adhesion. However, there were also phase-specific GO terms, including the glucuronidation hierarchy of GO terms in the phase 1 comparison. While represented by xenobiotic glucuronidation, the glucuronidation hierarchy appears to represent retinoic acid, which is added to the media in phase 1, since the glucuronidation of retinoic acid is involved in neurodevelopment .
To further test the functional relevance of DNA methylation changes resulting from DNMT3L overexpression, we tested the DMRs for enrichments within genomic regions annotated as CpG islands, CpG shores, CpG shelves, or open sea, as well as for enrichments within promoter and gene body annotations. Both the hypermethylated and hypomethylated DMRs from all pairwise phase comparisons were significantly (q < 0.05) enriched within CpG islands (Fig. 4a) and significantly de-enriched within intergenic regions (Fig. 4b). Furthermore, the hypermethylated DMRs from all pairwise phase comparisons were significantly (q < 0.05) enriched within promoters and regions of the gene body; while, the hypomethylated DMRs from all pairwise phase comparisons were significantly (q < 0.05) enriched within introns.
To understand the similarities and differences between the impact of DNMT3L overexpression on each developmental phase, DMRs from all pairwise comparisons were compared for overlap by both genomic coordinate and gene association. Overlaps of the DMRs by gene symbol included 962 in common to all three phases (Fig. 5a), even though only 12 overlapped by genomic coordinate in all three phases. The DMRs from all phase comparisons were then merged by sequence overlap to produce a consensus DMR profile, and hierarchal clustering analysis revealed that the cells clustered primarily by genotype and then developmental phase, and not by cell line of origin (Fig. 5b). The consensus DMR profile highlighted not only the overall skew towards DMR hypermethylation, but also clusters of DMRs that were specific to differentiation phase. A meta p-value of analysis of the GO terms from all pairwise phase comparisons revealed that the least dispensable significant (p < 0.05) terms were largely similar to the individual comparisons and primarily represented neurodevelopment (Fig. 5c).
The hypermethylated DMRs are enriched for bivalent chromatin and DS cross-tissue regions
To further investigate the functional relevance between the hypermethylated and hypomethylated DNMT3L DMRs, enrichment analyses within chromatin state maps from the reference epigenomes were performed. The most significant (q < 0.05) chromHMM chromatin state enrichments overall were for the hypermethylated DMRs in regions of bivalent chromatin (Fig. 6a) marked by H3K4me3 (Fig. 6b) in stem cells. This effect was most pronounced in the differentiated cells. In contrast, the hypomethylated DMRs were most significantly (q < 0.05) enriched within heterochromatic regions (Fig. 6a) marked by H3K9me3 (Fig. 6b) in ESCs.
To test the hypothesis that the DNMT3L DMRs represent a facet of the differentially methylated genes observed in DS, cross-tissue and pan-tissue comparisons with differentially methylated DS sites across diverse tissues were performed [26,27,28,29,30,31,32,33,34,35]. These results revealed a significant (q < 0.05) enrichment of the DNMT3L DMRs within DS differentially methylated sites identified from multiple tissues (Fig. 7a). This enrichment was predominantly for the hypermethylated DNMT3L DMRs from all pairwise phase comparisons. The only significant (q < 0.05) enrichment for the hypomethylated DNMT3L DMRs was from the growth comparison and within fetal frontal cortex CpGs. To further understand the DS differentially methylated genes most impacted by DNMT3L, genes with DNMT3L-associated differential methylation in all three phases of neural differentiation were filtered to include those also differentially methylated in at least 5 DS methylome datasets, and also investigated for overlap with embryonic bivalent chromatin marks from a previous study  (Fig. 7b). In this analysis, the common pan-tissue genes mapping to the DNMT3L-specific DMRs were generally from perinatal or neuronal DS datasets and most also mapped to bivalent chromatin. Notably, in the DNMT3L overexpression cells lines, one of the pan-tissue DMR sites identified maps to the promoter of TBX1 (Fig. 7c). TBX1 is a transcription factor involved in embryonic development and a key gene in 22q11.2 deletion syndrome (DiGeorge syndrome). DiGeorge syndrome shares phenotypic similarities with DS, which include congenital facial and cardiac abnormalities as well as intellectual disability . Additionally, TBX1 is regulated by retinoic acid, which is a key component of the SH-SY5Y differentiation protocol utilized .
Also present amongst the DNMT3L differentially methylated genes were subsets of the 25 previously known pan-tissue and multi-tissue DS genes , which were all hypermethylated. This observation is in line with the previous observation that 24 out of the 25 pan-tissue and multi-tissue DS genes are hypermethylated. The growth comparison contained DMRs mapping to ZNF837, RYR1, BCL9L, VPS37B, and RUNX1. The phase 1 comparison contained a DMR mapping to TEX14. The phase 3 comparison contained DMRs mapping to BCL9L, RFPL2, RYR1, ADAMTS10, MZF1, and RUNX1.
The DNMT3L DMRs also mapped to two clusters of genes previously known to be altered in DS and important in neurodevelopment, specifically the protocadherin and HOX genes clusters . The clustered protocadherin locus is located on chromosome 5 and involved in the establishment of single-cell neuronal identity . Notably, the protocadherin gamma cluster contains one of the previously known 25 pan-tissue DS genes (PCDHGA2) , and differential methylation within the cluster has been observed in diverse DS tissues, which include fetal brain, adult brain, neonatal blood, T-lymphocytes, buccal epithelial cells, and placenta [29,30,31, 33,34,35]. The protocadherin gamma cluster contained two hypermethylated DNMT3L DMRs and one hypomethylated DMR in the phase 1 comparison. The protocadherin beta cluster contained a hypermethylated DNMT3L DMR in the growth comparison and a different hypermethylated DMR in the phase 3 comparison. The clustered HOX genes, which are involved in early developmental patterning, contained DNMT3L DMRs in all phase comparisons. The HOXA gene cluster also contains one of the previously known 25 pan-tissue genes (HOXA4) , and large-scale hypermethylation of the cluster was apparent in DNMT3L-overexpressing cells at all phases of differentiation (Fig. 7d). Differential methylation of the HOX clusters has been observed in a number of tissues, which includes adult brain, fetal brain, neonatal blood, whole-blood, T-lymphocytes, placenta, buccal epithelial cells, and neural induced pluripotent stem cell (iPSC) derivatives [26, 27, 30,31,32,33,34,35]. DNMT3L DMRs mapping to HOXB7, HOXC4, HOXC5, HOXD10, and HOXD12 were present in the growth comparison, while those mapping to HOXD-AS2, HOTAIR (“HOX transcript antisense RNA” for the HOXC cluster), HOXC5, HOXC6, HOXC12, and HOXC13, and HOXD12 were present in the phase 1 comparison, and those mapping to HOTTIP (“HOXA distal transcript antisense RNA”), HOXC-AS2, HOXC-AS3, HOXC12, and HOTTAIR were present in the phase 3 comparison. The presence of DNMT3L DMRs in the HOX gene clusters is also notable given that retinoic acid is known to regulate their expression and chromatin state in mESCs .
While the DNMT3L DMRs overlapped with many of the 106 imprinted genes in the human genome, these results were not significantly higher than expected in enrichment analyses. The growth DMRs mapped to 18 imprinted genes: SNRPN, ZNF597, ZFAT, NTM, LIN28B, GRB10, MIMT1, KCNQ1, ADTRP, ANO1, OSBPL5, FAM50B, ATP10A, SNORD116, MIR296, GLIS3, PPP1R9A, and KCNK9. The phase 1 DMRs mapped to 19 imprinted genes: ATP10A, MKRN3, PHLDA2, RB1, MESTIT1, IGF2, DGCR6L, ZFAT, NTM, MIR296, MAGI2, ZIM2, SNORD116, KCNK9, ZDBF2, NAA60, KCNQ1, GNAS, and KCNQ1OT1. The phase 3 DMRs mapped to 22 imprinted genes: OSBPL5, TP73, MKRN3, WT1, KCNK9, NAA60, PPP1R9A, ATP10A, SNRPN, PLAGL1, GRB10, SNORD116, SLC22A3, MEG8, MEG3, ZIM2, MIMT1, GNAS, UBE3A, NDN, IGF2, and ADTRP. This lack of enrichment is consistent with the literature, as although DNMT3L is required for the establishment of imprinting, it is not required for its maintenance .
Finally, we examined the effect of lot-to-lot variation on DNMT3L overexpression in cell culture. This involved generating an additional set of replicate SH-SY5Y cell lines from a more recent lot of cells than the one used for the first two sets of transfection experiments. While this additional set of cell lines clustered differently from the other two and showed a stronger hypermethylation profile in DNMT3L-overexpressing cells compared to control cells, a replication analysis that included all three transfection experiments still reproduced the main the findings (Additional File 1: Supplementary Table 2 and Additional File 2: Supplementary Note and Supplementary Figures S1–S7). Together, these results demonstrate the reproducibility of the main findings of this study, which is that DNMT3L overexpression results in hypermethylation of genes differentially methylated in DS and regions of bivalent chromatin.
The results of this study confirmed previous observations and provided novel findings that are relevant to understanding the role of DNMT3L in establishing DNA methylation profiles during development and in DS. DNMT3L overexpression in undifferentiated, differentiating, and differentiated SH-SY5Y neurons resulted in a global increase in CpG and CpG island methylation. The DMRs mapped to genes involved in neurodevelopment, cellular signaling and gene regulation. Both DNMT3L genotype and cell differentiation phase were inter-related in the consensus DMR clustering, demonstrating that different regions of the same genes are affected during development. Furthermore, these results demonstrate that DNMT3L overexpression results in differential methylation of cross-tissue and pan-tissue sites within genes identified in 15 previous DS datasets from diverse tissues, with a majority of the pan-tissue genes mapping to bivalent chromatin. However, our results also highlight a directional dichotomy in the DNA methylation changes associated with DNMT3L overexpression. The hypermethylated DNMT3L DMRs showed a cross-tissue DS profile, in contrast to the hypomethylated DMRs that showed a distinct profile related to early neurodevelopment. There were also distinct differences in the chromatin state enrichments for the hypermethylated and hypomethylated DMRs.
The hypermethylated DMRs were enriched within regions of bivalent chromatin marked by H3K4me3. Bivalent chromatin consists of regions marked with both activating H3K4me3 and repressive H3K27me3 and serves to regulate key developmental genes in ESCs by repressing them during pluripotency and poising them for rapid activation upon removal of H3K27me3 during differentiation [40, 41]. Notably, while bivalent promoters in ESCs are unmethylated or hypomethylated, upon differentiation some bivalent promoters lose H3K4me3 and gain DNA methylation [14, 42, 43]. Disruption of bivalent chromatin has been previously implicated in DS [29, 35, 44]. Since DNMT3L recognizes unmethylated H3K4 , a possible mechanism is that regions of bivalent chromatin lose H3K4me3 as they differentiate and become hypermethylated by the excess DNMT3L that recognizes the H3K4me0. This hypothesis is consistent with the observation that phases 1 and 3 DNMT3L DMRs show much stronger enrichments within regions marked by H3K4me3 when compared to the growth phase. This proposed mechanism is also consistent with the growth phase still showing enrichment within regions marked by H3K4me3, as the cells are already committed to the neural crest lineage. Additionally, while bivalent chromatin is generally hypomethylated in normal cells, in cancer cells, bivalent chromatin becomes hypermethylated and is associated with increased developmental gene expression (45, 46). In neural progenitors, many promoters with the H3K27me3 modification gain DNA methylation during differentiation .
In comparison to the hypermethylated DMRs, the hypomethylated DNMT3L DMRs in our study were enriched within regions of heterochromatin marked by H3K9me3. Notably, the observation of DNMT3L causing hypomethylation has also been previously established. Knockdown of DNMT3L in mESCs resulted in 14,107 regions showing a decrease in methylation and 5,724 showing a gain in methylation . This ratio is comparable to our observations and the regions identified significantly (p = 0.0002) overlap with the genomic coordinates for the consensus DNMT3L DMRs. Mechanistically, DNMT3L is known to either release or sequester DNMT3A from heterochromatin, which has been hypothesized to enable preferential targeting of euchromatin . Taken together, we hypothesize that in our study, the release or sequestering of DNMT3A from heterochromatin by excess DNMT3L resulted in the preferential targeting of bivalent chromatin within regions that recently lost H3K4me3 as the cells differentiated.
While DNMT3L overexpression during neuronal differentiation recapitulates an important facet of the epigenetic signature of DS, it is not the only mechanism. The other methylation differences observed in DS brain are likely related to age, medication exposure, and an increase in copy number of the other genes encoding proteins related to DNA methylation and one-carbon metabolism that are also located on chromosome 21 . For example, we observed hypermethylation of the chromosome 21 RUNX1 locus in both the growth phase and phase 3 cells. Previously, we observed hypermethylation of RUNX1 as well as genome-wide hypomethylation of RUNX1-binding sites in neonatal DS when compared to non-DS control blood . However, we did not observe hypomethylation of RUNX1-binding sites in DNMT3L-overexpressing cells, likely because RUNX1 was diploid in SH-SY5Y cells. The results of our study also demonstrate reproducibility in gene ontologies and chromatin states impacted by DNMT3L overexpression, despite variation inherent to individual cell lines. Overall, our findings suggest that DNMT3L overexpression during neuronal differentiation recreates a facet of the aberrant Down syndrome DNA methylation signature by targeting specific chromatin states that regulate genes important for neurodevelopment.
Plasmids and cloning
DNMT3L was cloned into an expression plasmid by Gibson assembly. First, human wild-type DNMT3L cDNA (IMAGE clone 1541874) was PCR amplified from pD3LMyc . PCR amplification primers were design with the complementary overhangs necessary for Gibson assembly using NEBuilder (New England Biolabs). The sequence for the forward primer was: TGTCTCATCATTTTGGCAAAATGGAGCAGAAGCTGATCTCAGAGGAGGAC. The sequence of the reverse primer was: TCACCGCATGTTAGCAGACTTCCTCTGCCCTCGCCACCTCCGCTGCCGCCTAAAGAGGAAGTGAGTTCTGTTGAAAAATACTTG. The reverse primer was designed to include a sequence coding for a self-cleaving P2A peptide inframe of DNMT3L. A PiggyBac-compatible expression plasmid, PbCAG-eGFP, was a gift from Joseph Loturco (Addgene plasmid # 40973) . PbCAG-eGFP was linearized with EcoRI. Both the linearized plasmid and the DNMT3L PCR amplicon were run on an agarose gel and purified in a silica column (Qiagen). Purified fragments were mixed at 1:3 molar ratio and assembled with NEBuilder HiFi DNA Assembly Mater Mix (New England Biolabs) according to the manufacturer’s instructions. Assembled products were then introduced into NEB Turbo E. coli competent bacteria by chemical transformation, plated on agar plates supplemented with ampicillin, and grown at 37 °C according to the manufacturer’s instructions. Colonies were then grown in LB and the plasmid was purified in spin silica columns according to manufacturer’s instructions. Resulting pb-DNMT3L-P2A-eGFP plasmids were confirmed by Sanger sequencing before proceeding to cell culture transfection experiments.
Cell culture and transfection
Low passage SH-SY5Y cells in the undifferentiated growth phase were transfected with the plasmids via Lipofectamine 3000 (Invitrogen) according to the manufacturer’s instructions. After transfection, the cells were grown to a previously established timepoint when transient expression was shown to dissipate, so that only cells with stable integration could be sorted for GFP via flow cytometry. After flow cytometry, the population of > 95% GFP-positive cells were re-cultured, expanded for ~ 5 passages, and then differentiated according to Shipley et al.  except that the newer B-27 plus neuronal cell culture system (Gibco) was used and the retinoic acid was dissolved in DMSO. DNA was isolated from flash frozen cell pellets using the QIAamp DNA Micro kit (Qiagen) according to the manufacturer’s instructions.
Protein was prepared by homogenizing cells in RIPA buffer (Alfa Aesar) and quantitated by Bicinchoninic Acid (Thermo Fisher) following the manufacturers' instructions. Ten micrograms of protein per sample were loaded on Bis–Tris gel (Thermo Fisher Scientific) and transferred to a PVDF membrane (Bio-Rad) according to the manufacturers' instructions. Immunoblotting was visualized on a LiCor Odyssey instrument. Antibodies for immunoblotting were: DNMT3L (Santa Cruz Biotechnology, Catalog # sc-393603) and GAPDH (Advanced ImmunoChemicals Inc., Catalog # 2-RGM2). Molecular weights were determined using known weight markers and reported weights.
WGBS library preparation was performed using the post-bisulfite adaptor tagging (PBAT) method with the terminal deoxyribonucleotidyl transferase-assisted adenylate connector-mediated single-stranded-DNA ligation technique [49, 50] via the Accel-NGS Methyl-Seq DNA Library Kit (Swift Biosciences) with the Methyl-Seq Combinatorial Dual Indexing Kit (Swift Biosciences) according to the manufacturer’s instructions. The library pool was sequenced across 2 lanes on an Illumina NovaSeq 6000 S4 flow cell for 150 bp paired end reads to generate ~ 150 million unique read-pairs (~ 10X coverage) of the genome per sample.
Trimming of adapters and methylation bias, screening for contaminating genomes, alignments to hg38, deduplication, calculation of coverage and insert size metrics, extraction of CpG methylation values, generation of genome-wide cytosine reports (CpG count matrix), and examination of quality control metrics were performed using CpG_Me (https://github.com/ben-laufer/CpG_Me) [35, 51,52,53]. DMR calling as well as most downstream analyses were performed using DMRichR, which utilizes the dmrseq and bsseq algorithms (https://github.com/ben-laufer/DMRichR) [35, 54, 55]. Linear mixed-effects models of the average methylation level were utilized for global and CpG island methylation level analyses, where genotype (DNMT3L or eGFP) and cellular differentiation phase were direct effects, while cell line of origin was a random effect. Transcription factor motif enrichment testing using HOMER , DS cross-tissue enrichment testing via GAT , as well as enrichment testing for histone post-translational modifications (5 marks, 127 epigenomes) and chromatin states (15-state chromHMM model) [58, 59] through LOLA , were all performed as previously described .
Availability of Data and Materials
The datasets generated and/or analyzed during the current study are available in the NCBI Gene Expression Omnibus (GEO) repository, through accession number: GSE168276.
DNA methyltransferase 3L
Differentially methylated regions
Whole genome bisulfite sequencing
Mouse embryonic stem cells
Enhanced green fluorescent protein
Principal component analysis
Induced pluripotent stem cell
Presson AP, Partyka G, Jensen KM, Devine OJ, Rasmussen SA, McCabe LL, McCabe ERB. Current estimate of down syndrome population prevalence in the United States. J Pediatr. 2013;163:1163–8.
Do C, Xing Z, Yu YE, Tycko B. Trans-acting epigenetic effects of chromosomal aneuploidies: lessons from Down syndrome and mouse models. Epigenomics. 2017;9:189–207.
Lu J, Mccarter M, Lian G, Esposito G, Capoccia E, Delli-Bovi LC, Hecht J, Sheen V. Global hypermethylation in fetal cortex of Down syndrome due to DNMT3L overexpression. Hum Mol Genet. 2016;25:1714–27.
Chedin F, Lieber MR, Hsieh C-L. The DNA methyltransferase-like protein DNMT3L stimulates de novo methylation by Dnmt3a. Proc Natl Acad Sci. 2002;99:16916–21.
Chen Z-X, Mann JR, Hsieh C-L, Riggs AD, Chédin F. Physical and functional interactions between the human DNMT3L protein and members of the de novo methyltransferase family. J Cell Biochem. 2005;95:902–17.
Jeltsch A, Jurkowska RZ. Allosteric control of mammalian DNA methyltransferases—a new regulatory paradigm. Nucleic Acids Res. 2016;44:8556–75.
Jurkowska RZ, Anspach N, Urbanke C, Jia D, Reinhardt R, Nellen W, Cheng X, Jeltsch A. Formation of nucleoprotein filaments by mammalian DNA methyltransferase Dnmt3a in complex with regulator Dnmt3L. Nucleic Acids Res. 2008;36:6656–63.
Rajavelu A, Jurkowska RZ, Fritz J, Jeltsch A. Function and disruption of DNA Methyltransferase 3a cooperative DNA binding and nucleoprotein filament formation. Nucleic Acids Res. 2012;40:569–80.
Stepper P, Kungulovski G, Jurkowska RZ, Chandra T, Krueger F, Reinhardt R, Reik W, Jeltsch A, Jurkowski TP. Efficient targeted DNA methylation with chimeric dCas9-Dnmt3a-Dnmt3L methyltransferase. Nucleic Acids Res. 2017;45:1703–13.
Ooi SKT, Qiu C, Bernstein E, Li K, Jia D, Yang Z, Erdjument-Bromage H, Tempst P, Lin S-P, Allis CD, et al. DNMT3L connects unmethylated lysine 4 of histone H3 to de novo methylation of DNA. Nature. 2007;448:714–7.
Otani J, Nankumo T, Arita K, Inamoto S, Ariyoshi M, Shirakawa M. Structural basis for recognition of H3K4 methylation status by the DNA methyltransferase 3A ATRX-DNMT3-DNMT3L domain. EMBO Rep. 2009;10:1235–41.
Zhang Y, Jurkowska R, Soeroes S, Rajavelu A, Dhayalan A, Bock I, Rathert P, Brandt O, Reinhardt R, Fischle W, et al. Chromatin methylation activity of Dnmt3a and Dnmt3a/3L is guided by interaction of the ADD domain with the histone H3 tail. Nucleic Acids Res. 2010;38:4246–53.
Neri F, Krepelova A, Incarnato D, Maldotti M, Parlato C, Galvagni F, Matarese F, Stunnenberg HG, Oliviero S. Dnmt3L antagonizes DNA methylation at bivalent promoters and favors DNA methylation at gene bodies in ESCs. Cell. 2013;155:121.
Veland N, Lu Y, Hardikar S, Gaddis S, Zeng Y, Liu B, Estecio MR, Takata Y, Lin K, Tomida MW, et al. DNMT3L facilitates DNA methylation partly by maintaining DNMT3A stability in mouse embryonic stem cells. Nucleic Acids Res. 2019;47:152–67.
Kao T-H, Liao H-F, Wolf D, Tai K-Y, Chuang C-Y, Lee H-S, Kuo H-C, Hata K, Zhang X, Cheng X, et al. Ectopic DNMT3L triggers assembly of a repressive complex for retroviral silencing in somatic cells. J Virol. 2014;88:10680–95.
Bourc’his, D., Xu, G. L., Lin, C. S., Bollman, B. and Bestor, T. H. . Dnmt3L and the establishment of maternal genomic imprints. Science. 2001;294:2536–9.
Hata K, Okano M, Lei H, Li E. Dnmt3L cooperates with the Dnmt3 family of de novo DNA methyltransferases to establish maternal imprints in mice. Development. 2002;129:1983–93.
Do JH, Kim IS, Park T-K, Choi D-K. Genome-wide examination of chromosomal aberrations in neuroblastoma SH-SY5Y cells by array-based comparative genomic hybridization. Mol Cells. 2007;24:105–12.
Yusuf M, Leung K, Morris KJ, Volpi EV. Comprehensive cytogenomic profile of the in vitro neuronal model SH-SY5Y. Neurogenetics. 2013;14:63–70.
Krishna A, Biryukov M, Trefois C, Antony PMA, Hussong R, Lin J, Heinäniemi M, Glusman G, Köglsberger S, Boyd O, et al. Systems genomics evaluation of the SH-SY5Y neuroblastoma cell line as a model for Parkinson’s disease. BMC Genomics. 2014;15:1154.
Schroeder DI, Lott P, Korf I, LaSalle JM. Large-scale methylation domains mark a functional subset of neuronally expressed genes. Genome Res. 2011;21:1583–91.
Dunaway K, Goorha S, Matelski L, Urraca N, Lein PJ, Korf I, Reiter LT, LaSalle JM. Dental pulp stem cells model early life and imprinted DNA methylation patterns. Stem Cells. 2017;35:981–8.
Shipley MM, Mangold CA, Szpara ML. Differentiation of the SH-SY5Y human neuroblastoma cell line. J Vis Exp. 2016. https://doi.org/10.3791/53193.
Chen F, LoTurco J. A method for stable transgenesis of radial glia lineage in rat neocortex by piggyBac mediated transposition. J Neurosci Methods. 2012;207:172–80.
Reay WR, Cairns MJ. The role of the retinoids in schizophrenia: genomic and clinical perspectives. Mol Psychiatry. 2020;25:706–18.
Henneman P, Bouman A, Mul A, Knegt L, Van Der Kevie-Kersemaekers A-MM, Zwaveling-Soonawala N, Meijers-Heijboer HEJJ, Paul van Trotsenburg AS, Mannens MM, van Trotsenburg ASP, et al. Widespread domain-like perturbations of DNA methylation in whole blood of Down syndrome neonates. PLoS ONE. 2018;13:e0194938.
Bacalini MG, Gentilini D, Boattini A, Giampieri E, Pirazzini C, Giuliani C, Fontanesi E, Scurti M, Remondini D, Capri M, et al. Identification of a DNA methylation signature in blood cells from persons with down syndrome. Aging (Albany, NY). 2015;7:82–96.
Cejas RB, Wang J, Hageman-Blair R, Liu S, Blanco JG. Comparative genome-wide DNA methylation analysis in myocardial tissue from donors with and without Down syndrome. Gene. 2021;764:145099.
Mendioroz M, Do C, Jiang X, Liu C, Darbary HK, Lang CF, Lin J, Thomas A, Abu-Amero S, Stanier P, et al. Trans effects of chromosome aneuploidies on DNA methylation patterns in human Down syndrome and mouse models. Genome Biol. 2015;16:263.
Jones MJ, Farré P, McEwen LM, MacIsaac JL, Watt K, Neumann SM, Emberly E, Cynader MS, Virji-Babul N, Kobor MS. Distinct DNA methylation patterns of cognitive impairment and trisomy 21 in down syndrome. BMC Med Genomics. 2013;6:58.
Jin S, Lee YK, Lim YC, Zheng Z, Lin XM, Ng DPY, Holbrook JD, Law HY, Kwek KYC, Yeo GSH, et al. Global DNA hypermethylation in down syndrome placenta. PLoS Genet. 2013;9:e1003515.
Laan L, Klar J, Sobol M, Hoeber J, Shahsavani M, Kele M, Fatima A, Zakaria M, Annerén G, Falk A, et al. DNA methylation changes in down syndrome derived neural iPSCs uncover co-dysregulation of ZNF and HOX3 families of transcription factors. Clin Epigenetics. 2020;12:9.
El Hajj N, Dittrich M, Böck J, Kraus TFJ, Nanda I, Müller T, Seidmann L, Tralau T, Galetzka D, Schneider E, et al. Epigenetic dysregulation in the developing Down syndrome cortex. Epigenetics. 2016;11:563–78.
Laufer BI, Hwang H, Vogel Ciernia A, Mordaunt CE, LaSalle JM. Whole genome bisulfite sequencing of Down syndrome brain reveals regional DNA hypermethylation and novel disorder insights. Epigenetics. 2019;14:672–84.
Laufer BI, Hwang H, Jianu JM, Mordaunt CE, Korf IF, Hertz-Picciotto I, LaSalle JM. Low-pass whole genome bisulfite sequencing of neonatal dried blood spots identifies a role for RUNX1 in down syndrome DNA methylation profiles. Hum Mol Genet. 2020;11(8):563.
Court F, Arnaud P. An annotated list of bivalent chromatin regions in human ES cells: a new tool for cancer epigenetic research. Oncotarget. 2017;8:4110–24.
Du Q, de la Morena MT, van Oers NSC. (2020) The genetics and epigenetics of 22q11.2 deletion syndrome. Front Genet. 2020;10:1365.
Chen WV, Maniatis T. Clustered protocadherins. Development. 2013;140:3297–302.
De Kumar B, Parrish ME, Slaughter BD, Unruh JR, Gogol M, Seidel C, Paulson A, Li H, Gaudenz K, Peak A, et al. Analysis of dynamic changes in retinoid-induced transcription and epigenetic profiles of murine Hox clusters in ES cells. Genome Res. 2015;25:1229–43.
Bernstein BE, Mikkelsen TS, Xie X, Kamal M, Huebert DJ, Cuff J, Fry B, Meissner A, Wernig M, Plath K, et al. A bivalent chromatin structure marks key developmental genes in embryonic stem cells. Cell. 2006;125:315–26.
Vastenhouw NL, Schier AF. Bivalent histone modifications in early embryogenesis. Curr Opin Cell Biol. 2012;24:374–86.
Meissner A, Mikkelsen TS, Gu H, Wernig M, Hanna J, Sivachenko A, Zhang X, Bernstein BE, Nusbaum C, Jaffe DB, et al. Genome-scale DNA methylation maps of pluripotent and differentiated cells. Nature. 2008;454:766–70.
Galonska C, Charlton J, Mattei AL, Donaghey J, Clement K, Gu H, Mohammad AW, Stamenova EK, Cacchiarelli D, Klages S, et al. Genome-wide tracking of dCas9-methyltransferase footprints. Nat Commun. 2018. https://doi.org/10.1038/s41467-017-02708-5.
Lane AA, Chapuy B, Lin CY, Tivey T, Li H, Townsend EC, Van Bodegom D, Day TA, Wu SC, Liu H, et al. Triplication of a 21q22 region contributes to B cell transformation through HMGN1 overexpression and loss of histone H3 Lys27 trimethylation. Nat Genet. 2014;46:618–23.
Bernhart SH, Kretzmer H, Holdt LM, Jühling F, Ammerpohl O, Bergmann AK, Northoff BH, Doose G, Siebert R, Stadler PF, et al. Changes of bivalent chromatin coincide with increased expression of developmental genes in cancer. Sci Rep. 2016;6:1–18.
Dunican DS, Mjoseng HK, Duthie L, Flyamer IM, Bickmore WA, Meehan RR. Bivalent promoter hypermethylation in cancer is linked to the H327me3/H3K4me3 ratio in embryonic stem cells. BMC Biol. 2020;18:25.
Mohn F, Weber M, Rebhan M, Roloff TC, Richter J, Stadler MB, Bibel M, Schübeler D. Lineage-specific polycomb targets and de novo DNA methylation define restriction and potential of neuronal progenitors. Mol Cell. 2008;30:755–66.
Jurkowska RZ, Rajavelu A, Anspach N, Urbanke C, Jankevicius G, Ragozin S, Nellen W, Jeltsch A. Oligomerization and binding of the Dnmt3a DNA methyltransferase to parallel DNA molecules: Heterochromatic localization and role of Dnmt3L. J Biol Chem. 2011;286:24200–7.
Miura F, Enomoto Y, Dairiki R, Ito T. Amplification-free whole-genome bisulfite sequencing by post-bisulfite adaptor tagging. Nucleic Acids Res. 2012;40:e136.
Miura F, Shibata Y, Miura M, Sangatsuda Y, Hisano O, Araki H, Ito T. Highly efficient single-stranded DNA ligation technique improves low-input whole-genome bisulfite sequencing by post-bisulfite adaptor tagging. Nucleic Acids Res. 2019;47:e85.
Wingett SW, Andrews S. FastQ Screen: A tool for multi-genome mapping and quality control. F1000Research. 2018;7:1338.
Krueger F, Andrews SR. Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics. 2011;27:1571–2.
Ewels P, Magnusson M, Lundin S, Käller M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics. 2016;32:3047–8.
Korthauer K, Chakraborty S, Benjamini Y, Irizarry RA. Detection and accurate false discovery rate control of differentially methylated regions from whole genome bisulfite sequencing. Biostatistics. 2018;20:367–83.
Hansen KD, Langmead B, Irizarry RA. BSmooth: from whole genome bisulfite sequencing reads to differentially methylated regions. Genome Biol. 2012;13:R83.
Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, Cheng JX, Murre C, Singh H, Glass CK. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol Cell. 2010;38:576–89.
Heger A, Webber C, Goodson M, Ponting CP, Lunter G. GAT: a simulation framework for testing the association of genomic intervals. Bioinformatics. 2013;29:2046–8.
Ernst J, Kellis M. ChromHMM: Automating chromatin-state discovery and characterization. ChromHMM: Automating chromatin-state discovery and characterization. Nat Methods. 2012;2012(9):215–6.
Roadmap Epigenomics Consortium, Kundaje A, Meuleman W, Ernst J, Bilenky M, Yen A, Heravi-Moussavi A, Kheradpour P, Zhang Z, Wang J et al. Integrative analysis of111 reference human epigenomes. Nature 2015; 518: 317–329.
Sheffield NC, Bock C. LOLA: enrichment analysis for genomic region sets and regulatory elements in R and Bioconductor. Bioinformatics. 2015;32:587–9.
This work was supported by a National Institutes of Health (NIH) grant [3R01ES015359-10S1] to JML and I. Hertz-Picciotto, a Canadian Institutes of Health Research (CIHR) postdoctoral fellowship [MFE-146824] to BIL, a CIHR Banting postdoctoral fellowship [BPF-162684] to BIL, a Jane Coffin Childs Memorial Fund for Medical Research fellowship to JAG, a NIH T32 Fellowship to JAG [5T32CA108459], and the UC Davis Intellectual and Developmental Disabilities Research Center (IDDRC) [P50HD103526]. The flow cytometry was carried out by the UC Davis Flow Cytometry Shared Resource Laboratory and was supported by the UC Davis Comprehensive Cancer Center Support Grant (CCSG) awarded by the National Cancer Institute [NCI P30CA093373]. The library preparation and sequencing was carried out by the DNA Technologies and Expression Analysis Cores at the UC Davis Genome Center and was supported by a NIH Shared Instrumentation Grant [1S10OD010786-01].
Ethics approval and consent to participate
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.
: Table 1. Testable background regions and significant DMRs for DNMT3L overexpression in A) growth phase, B) phase 1, and C) phase 3. Table 2. Testable background regions and significant DMRs for DNMT3L overexpression in all 3 cell line replicates for A) growth phase, B) phase 1, and C) phase 3 (XLSX 78248 KB)
: Supplementary note. Figures S1–S7. Additional replicate cell line set analysis (PDF 10503 KB)
About this article
Cite this article
Laufer, B.I., Gomez, J.A., Jianu, J.M. et al. Stable DNMT3L overexpression in SH-SY5Y neurons recreates a facet of the genome-wide Down syndrome DNA methylation signature. Epigenetics & Chromatin 14, 13 (2021). https://doi.org/10.1186/s13072-021-00387-7
- Down syndrome
- Whole genome bisulfite sequencing
- Bivalent chromatin
- PiggyBac transgenesis