Selective demethylation of two CpG sites causes postnatal activation of the Dao gene and consequent removal of d-serine within the mouse cerebellum

Background Programmed epigenetic modifications occurring at early postnatal brain developmental stages may have a long-lasting impact on brain function and complex behavior throughout life. Notably, it is now emerging that several genes that undergo perinatal changes in DNA methylation are associated with neuropsychiatric disorders. In this context, we envisaged that epigenetic modifications during the perinatal period may potentially drive essential changes in the genes regulating brain levels of critical neuromodulators such as d-serine and d-aspartate. Dysfunction of this fine regulation may contribute to the genesis of schizophrenia or other mental disorders, in which altered levels of d-amino acids are found. We recently demonstrated that Ddo, the d-aspartate degradation gene, is actively demethylated to ultimately reduce d-aspartate levels. However, the role of epigenetics as a mechanism driving the regulation of appropriate d-ser levels during brain development has been poorly investigated to date. Methods We performed comprehensive ultradeep DNA methylation and hydroxymethylation profiling along with mRNA expression and HPLC-based d-amino acids level analyses of genes controlling the mammalian brain levels of d-serine and d-aspartate. DNA methylation changes occurring in specific cerebellar cell types were also investigated. We conducted high coverage targeted bisulfite sequencing by next-generation sequencing and single-molecule bioinformatic analysis. Results We report consistent spatiotemporal modifications occurring at the Dao gene during neonatal development in a specific brain region (the cerebellum) and within specific cell types (astrocytes) for the first time. Dynamic demethylation at two specific CpG sites located just downstream of the transcription start site was sufficient to strongly activate the Dao gene, ultimately promoting the complete physiological degradation of cerebellar d-serine a few days after mouse birth. High amount of 5′-hydroxymethylcytosine, exclusively detected at relevant CpG sites, strongly evoked the occurrence of an active demethylation process. Conclusion The present investigation demonstrates that robust and selective demethylation of two CpG sites is associated with postnatal activation of the Dao gene and consequent removal of d-serine within the mouse cerebellum. A single-molecule methylation approach applied at the Dao locus promises to identify different cell-type compositions and functions in different brain areas and developmental stages. Electronic supplementary material The online version of this article (10.1186/s13148-019-0732-z) contains supplementary material, which is available to authorized users.


Background
Embryonic and postnatal brain development is a continuous process that involves well-orchestrated changes, including age-and cerebral region-specific gene expression programs. The perinatal period is the main temporal window in which the timely occurrence of programmed molecular changes is critical for correct neural functioning throughout the remainder of life. Epigenetic mechanisms favoring brain development and neuronal plasticity may underlie these fundamental events [1][2][3][4][5]. Accordingly, neuroepigenetics is a growing research field since DNA methylation and chromatin modifications influence gene programs, neuronal plasticity, synaptogenesis and complex behaviors, and epigenetic marks are altered in several neuropsychiatric conditions [6][7][8][9][10][11][12][13][14]. Notably, from recent epigenome-wide studies, it emerged that several genes that undergo perinatal changes in DNA methylation are associated with neuropsychiatric disorders [15][16][17]. However, the mammalian brain is a very heterogeneous structure composed of functionally different areas and a variety of cell populations [18], which may exhibit unique DNA methylation profiles, whose diversity has been poorly investigated to date. In this regard, epigenetic modifications during the perinatal period may potentially drive critical changes in the genes regulating brain levels of D-amino acids [19,20]. In particular, free D-serine (D-Ser) and D-aspartate (D-Asp) have been identified as critical regulators of mammalian glutamatergic neurotransmission by acting as a co-agonist and agonist of N-methyl-D-aspartate receptors (NMDARs), respectively. These atypical amino acids show a unique regional and temporal pattern of emergence in the mammalian brain [21][22][23]. On the day of birth, a substantial amount of D-Ser is present in the mouse cortex, striatum, and cerebellum and then increases selectively in the cortex and striatum, whereas the cerebellar concentration of this D-amino acid declines drastically to trace levels around day 18 after birth [24]. Similarly, D-Asp content significantly decreases in the whole mouse brain a few days after birth, and such degradation is ultimately associated with an age-dependent increase in D-aspartate oxidase (Ddo) expression [19,20,25,26]. Although the specific molecular mechanisms involved are still unclear, the precise spatiotemporal regulation of the genes involved in D-Ser synthesis (serine racemase, Srr) or degradation (D-amino acid oxidase, Dao) and D-Asp degradation (Ddo) is thought to be critical for regulating the correct regional concentrations of these D-amino acids during pre-and postnatal brain development [27]. Accordingly, the expression and activity of the enzymes encoded by these genes have been shown to vary across different brain areas and different perinatal times in humans and rodents [24,26,28,29]. Since D-Asp and D-Ser are believed to play a direct role in regulating NMDAR-related synaptic plasticity, morphology, and function [22,30,31], the expression program of the genes controlling the metabolism of these D-amino acids must be rigorously orchestrated in mammalian brain development, and epigenetic mechanisms are the major candidate directors.
In this work, we provide a comprehensive analysis of mRNA expression and DNA methylation changes occurring at genes regulating the concentrations of D-Ser and D-Asp in three mouse brain areas during postnatal development. Here, we report striking spatiotemporal epigenetic modifications at the Dao gene during neonatal cerebellar development occurring selectively in the astrocytes for the first time. Notably, ultradeep methylation analysis allowed us to track the dynamic evolution of Dao epialleles (specific methyl-CpG arrangements) during ontogenesis in different brain areas and in specific cell types.

Results
Occurrence of free D-Ser and D-Asp in the hippocampus, cortex, and cerebellum during postnatal life We hypothesized that dynamic epigenetic modifications may drive changes in gene expression programs ultimately leading to area-specific perinatal shifts in D-amino acid concentrations. As a first step, we wished to validate our study system in terms of changes in D-Ser and D-Asp levels in different brain areas during mouse brain postnatal development. To this end, we quantified the ratio of the amounts of the D-and total (D +L) forms of both Ser and Asp by HPLC analysis at different postnatal time points (P1-P15-P30-P60) and in different brain regions, including the cortex (CX), hippocampus (HIPP), and cerebellum (CB) (Fig. 1). We found that the D-Ser/ total serine ratio significantly increased from P1 to P15 in both the hippocampus and cortex and that a significantly high ratio persisted throughout the following postnatal phases (one-way ANOVA, hippocampus: F(3, 9) = 166.7, p < 0.0001; cortex: F(3, 7) = 64.26, p < 0.0001) (Fig. 1a, b). In striking contrast to what occurs in the hippocampus and cortex, in the cerebellum, we observed a significant progressive age-dependent reduction in the D-Ser/total serine ratio during the same postnatal development stages (one-way ANOVA, F(3, 8) = 20.75, p = 0.0002; Fig. 1c). At the same postnatal time points in the same brain structures, we also investigated the ratio between D-Asp and total Asp levels by HPLC analysis (Fig. 1d-f). Consistent with previous work [20], the D-Asp/total aspartate ratio dramatically decreased after P1 in the hippocampus, cortex and cerebellum and remained very low in all the following postnatal stages (one-way ANOVA, hippocampus: F(3, 9) = 53.24, p < 0.0001; cortex: F(3, 7) = 28.12, p = 0.0003; cerebellum: F(3,9) = 476.3, p < 0.000; Fig. 1d-f).
Dynamic demethylation of two specific relevant CpG sites at the Dao promoter determines programmed Dao gene activation and a decrease in D-Ser levels in the cerebellum Here, we investigated the DNA methylation and mRNA expression of the Dao gene in the same brain regions and postnatal stages for which HPLC detection was performed. Methylation analysis was performed by targeted bisulfite sequencing with high coverage (approximately 100,000 reads per sample). First, we analyzed the average methylation level at the Dao promoter. We first focused on a region of 389 nucleotides of the Dao promoter (region PR1, Fig. 2a) that contains four CpG sites (+ 7; + 101; + 217; + 334), all of which are located downstream of the transcriptional start site (TSS). The average methylation levels at P1, P15, P30, and P60 in all analyzed brain regions were compared (Fig. 2b). No significant differences were found in HIPP and CX during development, and these structures displayed a high (approximately 75%) nearly stable degree of methylation during the different stages. In contrast, remarkable differences in the average DNA methylation were observed in CB. We found a significantly (one-way ANOVA; p ≤ 0.001) higher degree of methylation at postnatal day 1 (71.5% ± 0.03) compared to P15 (43.5% ± 0.02), P30 (44.8% ± 0.004), and P60 (53.1% ± 0.003). No significant differences were observed in the comparison between P15, P30, and P60. Additionally, we analyzed mRNA expression levels in HIPP, CX, and CB during brain development (Fig. 2c). Our data indicated very low Dao mRNA levels in both HIPP and CX, which showed no significant changes during  The numbers of the CpG sites refer to the putative transcriptional start site (TSS), indicated with + 1. White box indicates the putative upstream regulatory region; black box indicates the first exon; and gray box indicates the first intron. Black arrows at the top of the map specify the position of the primers used for bisulfite amplification. PR1 denotes the promoter region downstream of the TSS analyzed in all brain areas. PR2 indicates the promoter region upstream of the TSS explored in the cerebellum area. The Dao sequence was retrieved by Ensembl with the following accession number: ENSMUSG00000042096. b Average Dao methylation (%) is shown for HIPP, CX, and CB at the P1, P15, P30, and P60 developmental stages (n = 3 mice/stage). Comparisons between developmental stages in each area were performed using one-way ANOVA followed by Tukey's multiple comparison post-hoc test. *p ≤ 0.001. c mRNA expression levels are reported for each analyzed brain area during brain development. Dao mRNA expression was normalized to the mean of two housekeeping genes and is expressed as 2 -ΔCt values, * p ≤ 0.001 (one-way ANOVA followed by Tukey's multiple comparison post-hoc test). d Graphs represent the average methylation at single CpG sites in promoter region 1 (PR1), downstream of the TSS, in HIPP, CX, and CB during ontogenesis. e Methylation levels at CpG sites located in promoter region 2 (PR2, see a) are shown. The methylation trend of the seven CpG sites is indicated at each developmental stage. f A comprehensive view of single CpG methylation levels upstream and downstream of the Dao promoter region in the cerebellum during the transition from P1 to P15. Nine CpG sites are indicated, and their position is referred to as the TSS. CpGs with average methylation < 30% are shown in green, while CpG sites with average methylation > 50% are indicated in red ontogenesis. Conversely, a significant (one-way ANOVA; p ≤ 0.001) increase in Dao mRNA expression was found in the cerebellum during brain development (Fig. 2c). We then evaluated the methylation at a single-CpG resolution during ontogenesis in HIPP, CX, and CB (Fig. 2d). In the CB area in particular, a striking significant (one-way ANOVA; p ≤ 0.001) decrease in the degree of methylation was observed specifically at the CpG + 7 and + 101 sites from P1 to P15 (from 52.8 ± 0.06 to 8.2 ± 0.02%), and this low level of methylation persisted at P30 and P60. In contrast, the average methylation at CpG + 217 and CpG + 334 was elevated and remained rather constant over time in all brain areas except at the CpG + 217 site in CB, which showed a transient decrease at P15, switching back to high levels (above 80%) of methylation at P30 and P60. To validate the above results and determine whether other CpG sites located upstream of the TSS were subjected to similar dynamic demethylation in CB, we extended the bisulfite analysis to the PR2 region (− 234/+ 137 nts, Fig. 2a). This region encompasses the TSS and includes both the abovementioned + 7 and + 101 sites and five additional CpG sites located upstream (− 208, − 177, − 162, − 129, and − 102; Fig. 2a) in a single amplicon. The results confirmed the decrease in methylation at the + 7 and + 101 sites and showed that the methylation drift was essentially limited to these sites just downstream of the TSS (Fig. 2e, f). Thus, the very large selective reduction of DNA methylation observed in the cerebellum at the + 7 and + 101 sites (Fig. 2d, e) was functionally consistent with the robust activation of the Dao gene and ultimately with the dramatic reduction of the D-Ser concentration detected in the mouse cerebellum at P30 (Fig. 1c).
Single-CpG resolution analysis of 5-mCpG and 5-hmCpG dynamics at the Dao regulatory region argues for the involvement of active DNA demethylation DNA methylation and demethylation dynamics depend on the competitive opposite activities of DNA methyltransferases (DNMTs) and Ten-eleven translocation (TET) enzymes [32]. 5-hmC represents an intermediate product of TET-mediated active DNA demethylation [33,34]. However, eventual presence of 5-hydoxymethylcytosine (5-hmC) could not be detected by bisulfite analysis. Therefore, we evaluated the levels of 5-hmC at the Dao gene in CB at P1 and P15 by oxidative bisulfite sequencing [35]. Results, shown in Fig. 3, revealed the presence of high amount of hydroxymethylation at CpG + 7 and CpG + 101. Consistently, these sites underwent substantial demethylation at later stages (Fig. 2d). Although CpG + 217 exhibited high levels of 5-hmC, at this specific site 5-mC levels increased at P15 while lower levels of 5-hmC were detected at this stage. These findings suggest that during the brain development, the activity of TET enzymes prevailed over DNMTs activity at the + 7 and + 101 CpGs while the opposite dynamics occurred at the + 217 CpG. Strikingly, at the + 334 CpG, 5methylcytosine was very stable and did not undergo either hydroxymethylation or demethylation over time. Overall, our data strongly support the existence of an active demethylation process at selected sites and poor dynamics of 5-hmC at surrounding stably methylated CpG site (+ 334).

Dao demethylation in the cerebellum at postnatal day 15 mostly occurs in astrocytes
Previous immunohistochemical analyses demonstrated that during postnatal stages, Dao activation occurs prevalently in cerebellar astrocytes and to a lesser extent in neurons [36,37]. However, during prenatal and neonatal stages, the mouse brain undergoes radical changes in the cell type composition, and these modifications are distinct and characteristic for each brain area [5,18]. To verify whether the Dao gene undergoes demethylation in cerebellar astrocytes, we performed cell separation experiments from cerebellum tissue at the P1 and P15 time points. Then, we analyzed Dao promoter (PR1 region; Fig. 2a) methylation in purified astrocytes, neurons and a fraction containing microglia, oligodendrocytes, and endothelial cells (MOE fraction) (Fig. 4a). The results showed a decrease in the degree of methylation between days P1 and P15 mainly at the + 7 and + 101 CpGs and in astrocytes compared to neurons and other cell types (Fig. 4b). Because at these postnatal stages astrocytes are likely active proliferating cells [38], these data do not exclude passive demethylation. However, oxidative bisulfite sequencing analysis performed in whole cerebellum tissue (Fig. 3) makes us confident that a process of active demethylation is mainly responsible for the methylation decrease occurring at the + 7 and + 101 CpGs in cerebellar astrocytes.
Epiallele distribution analysis in HIPP, CX, and CB during brain development reveals area-and stage-specific epialleles at the Dao promoter We recently provided proof of concept that studying DNA methylation at specific loci at single-molecule resolution (epiallele distribution analysis) allows one the tracking of the spatiotemporal evolution of cell-to-cell methylation differences in a given cell population [19,39,40]. By using this method, requiring a large number of sequence reads per sample, it is possible to quantitatively and qualitatively analyze the arrangements of all methylated CpG at a specific locus at single-molecule level. Each specific arrangement is defined here as an "epiallele." We applied epiallele analysis to investigate whether epialleles at the Dao gene locus evolve in specific selected patterns depending on the brain developmental stage. After ultradeep bisulfite sequencing, an average of 100,000 sequences/sample were obtained and analyzed via the ampliMethProfiler pipeline (https://sourceforge.net/projects/amplimethprofiler) [41], which generated a BIOM format table containing the counts of Dao epialleles for all samples. The BIOM file was subjected to a rarefaction procedure to avoid bias due to different numbers of reads per sample and was used to perform the further analyses. Considering the Dao region including four CpG sites (PR1, Fig. 2a), the number of potential epialleles types was 16 (2 4 ). All the possible epialleles (16) were found to be present in each brain area and developmental stage but in different proportions. For each brain region, we performed principal component analysis (PCA) to verify the presence of patterns associated with various developmental stages (Fig. 5a). Strikingly, although the average levels of methylation in HIPP and CX were found to be almost identical during development   (Fig. 5a). In both HIPP and CX, a highly heterogeneous epiallelic composition was observed on postnatal day 1, which converged to relatively homogeneous epiallele types during brain development (Fig. 5a). For the cerebellum, the same analysis gave rise to more distinct grouping of each developmental stage, especially at P30 and P60, which was clearly influenced by the quantitative variation at these stages. Together, these data revealed non-stochastic, wellorchestrated, DNA methylation remodeling that occurs in the cells of different brain areas during development. We then performed PCA by integrating all brain areas and developmental time points (Fig. 5b). We confirmed that the The overall PCA analysis of brain region differences based on the epiallelic composition is shown, including all developmental stages. On the left, the score plot displays the epigenetic correlation of Dao epialleles between HIPP (green), CX (blue), and CB (red). On the right, the loading plot of the PCA displays the impact of each of 16 Dao epialleles on the brain area distribution evolution of the epiallelic distribution in the cerebellum was distinct compared to that in HIPP and CX (Fig. 5b). Interestingly, in very early stages (P1), the cerebellum shared epiallelic patterns with HIPP and CX, whereas at later stages, the epiallelic profiles evolved in a distinct manner (Fig. 5b). None of the above information could be retrieved from average standard methylation analyses alone, demonstrating the high potential of epiallelic distribution analysis.
DNA methylation and mRNA expression analyses of the Srr gene in HIPP, CX, and CB at different developmental times We next performed DNA methylation analysis and mRNA expression analysis of the Srr gene, encoding the D-Ser-synthesizing enzyme [42][43][44]. The Srr promoter region includes a very dense CpG island. We analyzed a region of 350 bp in length including 30 CpG sites located just downstream of the TSS (Fig. 6a). Very low levels of CpG methylation were detected in all analyzed brain areas and developmental stages, with average levels between 0.14% and 0.32%, (Fig. 6b, c). This finding was expected given the presence of a dense CpG island in the Srr promoter. No significant differences were found between different areas or during postnatal brain development. We then performed an mRNA expression analysis of Srr (Fig. 6d). In both HIPP and CB, we found significant (one-way ANOVA; p ≤ 0.001) differences in mRNA levels between the P1 and P30 stages. In CX, a significant (one-way ANOVA; p ≤ 0.001) increase in Srr mRNA levels was found in the comparison between P1 and the other developmental times. We conclude that the observed increase in mRNA levels is not controlled by a change in the CpG methylation state, at least in the analyzed region. Due to the low degree of methylation detected at each CpG site, epiallele analysis was considered uninformative and was therefore not performed. Interestingly, the increase in Srr mRNA levels (Fig. 6d) and stable low Dao expression (Fig. 2c) correlated well with the increase in D-Ser levels measured by HPLC in CX (and transiently in HIPP) over time (Fig. 1a, b). More strikingly, in CB, the increase in Srr levels (Fig. 6d) over time was counterbalanced by a simultaneous robust increase in Dao expression (Fig. 2c), resulting in a net decrease in D-Ser levels (Fig. 1c).
DNA methylation and mRNA expression of the Ddo gene in different brain areas at different developmental stages We previously reported that demethylation of Ddo in whole brain leads to increased Ddo mRNA expression and downregulation of cerebral D-Asp levels [19,20]. Here, we further investigated whether regional Ddo demethylation and changes in the D-Asp concentration occur in specific brain areas and postnatal stages. We analyzed the DNA methylation state and mRNA expression of the Ddo gene in HIPP, CX, and CB at P1, P15, P30, and P60. A region of 405 bp upstream of the TSS, including six CpG sites (− 363, − 330, − 318, − 242, − 175, − 125), was analyzed (Fig. 7a). In CX, neither the average Ddo methylation level nor single-CpG methylation within the Ddo promoter changed at any of the examined developmental stages (Fig. 7b, c). However, a slight increase in mRNA expression was found during development (Fig. 7d). An evident and significant decrease (oneway ANOVA; p = 0.001) in average methylation was found in HIPP during the transition from P1 to P60 (Fig. 7b). Interestingly, all analyzed CpG sites underwent demethylation from P1 to P15, and the methylation level then remained low over time (Fig. 7c). Consistent with the DNA methylation data, Ddo mRNA levels significantly (one-way ANOVA; p ≤ 0.001) increased from P1 to P60 (Fig. 7d, Additional file 2: Figure S2). In CB, a tendency of the average Ddo methylation to decrease was observed from P1 to P15, which then slightly increased at P30 and P60 (Fig. 7b). Accordingly, a gradual significant (one-way ANOVA; p ≤ 0.001) increase in mRNA expression was found over time in the cerebellum area ( Fig. 7d; Additional file 2: Figure S2). Interestingly, when we evaluated the average single-CpG methylation of Ddo, we found that all six CpG sites presented a trend of decreased methylation from the P1 to P15 stages, and the methylation level then remained consistently low at the − 125 and − 175 CpG sites (Fig. 7c). Conversely, all the other distal CpGs underwent a slight increase in methylation (Fig. 7c), possibly explaining the apparent tendency of higher average methylation levels observed from developmental stages P15 to P60 (Fig. 7b). These data highlight a decisive role of the degree of methylation at the − 125 and − 175 CpG sites in the control of Ddo expression during brain development for the first time. Consistent with the observed Ddo mRNA increase in HIPP, CX, and CB during development ( Fig. 7d; Additional file2: Figure S2), a considerable decrease in D-Asp levels was observed in all analyzed brain areas in the temporal window of P1-P15 (Fig. 1d-f). We previously reported a Ddo epiallele distribution analysis of the developing whole brain [19]. Here, we performed Ddo epiallele analysis of specific areas over time. Hierarchical cluster analysis (HCA) showed that developmental times were distinguishable by their epiallelic compositions mainly in cerebellum (Additional file 1: Figure S1).

Discussion
D-Ser is a physiological endogenous co-agonist of NMDA receptors at central excitatory synapses [45]. D-Ser can influence many NMDAR-dependent functions, including brain development [46,47], and behaviors such as cognition and social interaction. Although the role of free D-Asp has been less addressed, it has been recently shown that this D-amino acid influences NMDAR-mediated transmission [48]. In the brain, extracellular D-Asp activates NMDARs inducing considerable L-glutamate release in the prefrontal cortex of freely moving mice through the presynaptic activation of NMDA, AMPA/kainate, and mGlu5 receptors [49]. D-Aspartate is enriched in the embryonic brain of rodents and humans and its concentration strongly decreases after birth [50]. The proper control of brain D-amino acid levels has a clinical implication in both pathogenesis and therapeutic approach of different clinical conditions. For instance, serine deficiency syndromes include epilepsy, spasticity, and neurocognitive symptoms; these symptoms may partially be ameliorated by serine supplementation [51] In addition, D-serine deregulation is involved in the positive, negative, and cognitive symptoms of schizophrenia [52] and the inhibition of DAO appears to be a viable strategy to increase D-serine level and to have therapeutic potential in schizophrenia [53].
Based on the basic and translational interest of free Damino acids in brain physiology and pathology, we have The analyzed CpG island in the Srr promoter region is shown. White box represents the putative upstream regulatory region; black box indicates the first exon; and gray box indicates the first intron. The analyzed Srr amplicon is located between the black arrows on the top, indicating the forward and reverse primers used for bisulfite amplification. The positions of the analyzed CpGs refer to the TSS (+ 1). The Srr sequence was retrieved by Ensembl with the accession number: ENSMUSG00000001323. b Average methylation (%) at single-CpG level is indicated for each brain region during ontogenesis. c The average methylation (%) of the Srr gene in all brain regions at all analyzed developmental stages is reported. Statistical analyses were performed using one-way ANOVA followed by Tukey's multiple comparison post-hoc test. Alpha was considered significant at ≤ 0.001. d The spatiotemporal distribution of Srr mRNA levels is reported in histograms. Srr mRNA expression is normalized to the mean values for two housekeeping genes and expressed as 2 -ΔCt values. Statistical analyses were performed using one-way ANOVA followed by Tukey's multiple comparison post-hoc test. * p ≤ 0.001 here investigated the spatiotemporal changes in the DNA methylation and expression of genes involved in the modulation of D-Ser and D-Asp, along with contextual changes in the levels of these neurochemicals. Methylation analysis was performed via a single-molecule approach (ultradeep methylation analysis), as a proxy for single-cell methylation analysis, in three different brain regions (HIPP, CX, and CB) at postnatal days P1, P15, P30, and P60. We found programmed selective demethylation at two specific CpG sites in genes involved in D-Ser and D-Asp metabolism during specific developmental stages and in specific brain areas. Dynamic demethylation phenomena resulted in progressive activation of these genes and a dramatic area-specific modulation of D-Ser and D-Asp The average Ddo methylation is presented for HIPP, CX and CB at each developmental stage. Statistical analyses were performed using one-way ANOVA followed by Tukey's multiple comparison posthoc test. Alpha was considered significant at ≤ 0.001. c The average methylation (%) calculated at all seven CpGs is shown for each brain region during ontogenesis. d Ddo mRNA levels are reported for all brain regions and developmental times. Ddo mRNA expression is normalized to the mean values for two housekeeping genes and expressed as 2 -ΔCt values. One-way ANOVA followed by Tukey's multiple comparison post-hoc test was used to assess statistical significance. *p ≤ 0.001 levels during postnatal brain development. Environmentally induced alterations in physiological DNA methylation/demethylation turnover at specific genes in the brain during the perinatal period have been linked to the subsequent onset of neuropsychiatric diseases [54,55]. Pioneering studies demonstrated the crucial role of postnatal programmed demethylation of a single or few CpG sites in the glucocorticoid receptor gene (exon 1 7 promoter). The absence of this demethylation phenomenon, caused by environmental sources, leads to persistent alteration of the HPA axis, resulting in the transmission of the adaptive response across generations [14,56]. It has been proposed that effects on DNA methylation such as those described in the present work, even at single CpG sites, serve as an intermediate mechanism that imprints dynamic environmental experiences upon the fixed genome, resulting in stable alterations in phenotype [56].
In this work, we found that the Dao gene underwent programmed demethylation selectively in CB in a specific time window (P1-P15). The concomitant Dao activation and consequent postnatal D-Ser degradation likely represent physiological events required for correct postnatal cerebellum development. Accordingly, it has been shown that the major impact on D-Ser levels observed in mice lacking the Dao gene is limited to the cerebellum, where D-Ser levels remain high until adult life [57]. Dynamic demethylation of the Dao gene occurred specifically at only two CpG sites (+ 7 and + 101) located just downstream of the TSS, and this phenomenon supports gene activation. The possibility that the astrocyte fraction (assuming that these cells exhibit demethylated + 7 and + 101 CpG sites from P1 onward) may prevail over time in the cerebellum, providing false evidence of demethylation, was considered. Following cell separation, a decreased degree of methylation was still observed in purified cerebellar astrocytes collected at different stages. Therefore, our findings indicate that loss of methylation, rather than simple astrocyte number expansion, was mainly responsible for the observed demethylation. The upstream signaling mechanisms as well as the mechanisms driving specificity of targeted methylation and demethylation events remain unknown. Nevertheless, it is now established that mammalian cells, including neural cells, harbor the complete enzymatic machinery for not only methylating CpGs (and non-CpGs) but also actively demethylating DNA [58,59]. In this work, we provide evidence that hydroxymethylcytosine, an intermediate compound of TET-mediated demethylation, is specifically enriched at the CpG sites undergoing demethylation at later developmental stages. In contrast, poor methylation dynamics was detected at the + 334 CpG, where the 5-methylcytosine was very stable and did not undergo either hydroxymethylation or demethylation over time. Our data strongly argue for the occurrence of an active demethylation process at Dao gene in cerebellum during development.
Through epiallele distribution analysis, we also investigated how specific methylation patterns are generated in different stages and brain areas during development. The evolution of Dao epialleles was distinct between different stages and different brain areas in mice. This phenomenon was often not associated with changes in average methylation levels but was instead likely due to the reconfiguration of methylation patterns among the cells of each area at each developmental time point. Our recent epiallele analysis of the human Dao gene based on methylation classes performed in different areas of postmortem brain [40] showed that the epiallele distribution at Dao is distinct in the cerebellum compared to other areas, indicating the conservative nature of epiallele evolution across species.
To obtain a complete picture of the methylation events occurring at D-Ser-modulating genes, we also performed methylation and expression analysis of the Srr gene in selected mouse brain regions for the first time in this work. Consistent with a recent human study [40], no changes in Srr gene methylation (stably hypomethylated) were found in mouse brain, although substantial regionally and temporally specific expression modifications were detected. On the other hand, increasing Srr mRNA levels and stable low Dao expression correlated well with the progressive increases in D-Ser levels observed in CX and HIPP over time. Interestingly, in the cerebellum, the increase in Srr mRNA levels over time was counterbalanced by a simultaneous robust increase in Dao expression, resulting in a net decrease in D-Ser levels. For Ddo, the D-Asp-degrading enzyme, in a recent study conducted on the whole brain, we found that this gene is strongly activated at early postnatal stages by contextual demethylation of a large (approximately 400 bp) region surrounding the TSS, accompanied by a postnatal decrease in D-Asp levels [19,20]. Here, we have investigated separate brain areas and found that demethylation occurred in HIPP and CB but, surprisingly, not in CX. In CX, the observed slight increase in Ddo mRNA levels can likely be attributed to mechanisms other than DNA methylation changes. Of particular interest, strong Ddo gene activation was prominent in CB, where only two CpG sites (− 125 and − 175) remained demethylated over time, suggesting a critical role of these two specific sites.
Taken together, our data highlight the importance of post-natal changes of DNA methylation and hydroxymethylation at few CpG sites in the proper control of D-Ser and D-Asp levels with potential clinical implications.

Conclusion
Our data strongly encourage the investigation in the near future of whether failure of timely programmed active demethylation at relevant CpG sites may affect the physiological decreases in D-Ser and D-Asp levels thought to be required for correct brain development [22,27,45,[60][61][62][63][64]. Coherently, altered levels of D-Ser and D-Asp have been found in different neuropsychiatric conditions [50,[65][66][67][68]. It will be very important to investigate possible environmental factors affecting correct demethylation at the Dao and Ddo genes and whether altered methylation of the Dao and Ddo genes may be associated with neuropsychiatric conditions.
In conclusion, we hypothesize that eventual dysfunction of postnatal changes of DNA methylation/hydroxymethylation dynamics at few CpG sites, indirectly establishing the proper D-Ser and D-Asp brain levels, may have profound clinical impact in the etio-pathogenesis and treatment of neurodevelopmental disorders.

Animals
All experiments were performed on male animals. C57BL/ 6 J mice were purchased from The Jackson Laboratory. Mice were housed in groups (n = 4 or 5) in standard cages (29 × 17.5 × 12.5 cm) at constant temperature (21-24°C) and maintained on a 12/12 h light/dark cycle, with food and water ad libitum. All research involving animals was performed in accordance with the European directive 86/ 609/EEC governing animal welfare and protection, which is acknowledged by the Italian Legislative Decree no. 26 (March 14, 2014). Animal research protocols were also reviewed and approved by the local animal care committee at University of Naples "Federico II." All efforts were made to minimize the animal's suffering.

Mouse tissue collection
Hippocampus, cortex, and cerebellum were collected from C57BL/6 J mice at different developmental stages, including the following time points: post natal day (P) P1 (n = 3), P15 (n = 3), P30 (n = 3), and P60 (n = 3). Mice were killed, and all brain regions were dissected out within 20 s on an ice-cold surface. All tissue samples were pulverized in liquid nitrogen and stored at 80°C for subsequent processing.

High-performance liquid chromatography analysis
Brain tissue samples were analyzed as previously reported [69] with minor modifications [20,70]. Samples were homogenized in 1:10 (w/v) 0.2 M TCA, sonicated (3 cycles, 10 s each), and centrifuged at 13,000 g for 20 min. The precipitated protein pellets were stored at − 80°C for protein quantification, whereas the supernatants were neutralized with NaOH and subjected to precolumn derivatization with o-phthaldialdehyde (OPA)/N-acetyl-L-cysteine in 50% methanol. Enantiomer derivatives were then resolved on a Simmetry C8 5 m reversed-phase column (Waters, 4.6 × 250 mm), in isocratic conditions (0.1 M sodium acetate buffer, pH 6.2, 1% tetrahydrofuran, 1 mL/min flow rate). A washing step in 0.1 M sodium acetate buffer, 3% tetrahydrofuran, and 47% acetonitrile was performed after every single run. Identification and quantification of D-Ser and Dasp was based on retention times and peak areas, compared with those associated with external standards. The identity of D-Asp peak was confirmed by adding known amount of external standards, and by the selective degradation catalyzed by an active hDDO recombinant enzyme (kind gift by Dr. Homma and Dr. Katane). The samples were added with 10 μg of the enzyme, incubated at 30°C for 4 h and then derivatized. Total protein content of homogenates was determined using the Bradford assay method after resolubilization of the TCA precipitated protein pellets. The detected total levels of the amino acids in homogenates were normalized by the total protein content. The content of free D-Ser and D-Asp was expressed as ratio between D-and total (D + L) amount of both Ser and Asp.
Isolation of neurons, astrocytes, and enriched microglia/ oligodendrocyte/endothelial cells Cerebral cortex and cerebellum were dissected from six to eight C57/BL6 mice at P1 and P15 and dissociated into a single-cell suspension using the Neural Tissue Dissociation Kit Postnatal Neurons (Milteny Biotec), according to the manufacturer's protocol, with minor modifications [71,72]. Briefly, after carefully removing the meninges, brain tissue was weighed and cut into small pieces in ice-cold Hank's buffered salt solution (HBSS, Gibco). The minced tissue were then dissociated enzymatically in a solution of Enzyme Mix 1 for 25 min at 37°C and then mechanically after adding Enzyme Mix 2. Finally, the sample was applied to a cell strainer, and the resulting single-cell suspension was centrifuged at 300×g for 10 min. Mouse astrocytes, neurons, and other nonneuronal cells were sequentially separated from the same brain samples with a magnetic activated cell sorting approach (MACS®). ACSA2-expressing astrocytes were enriched by positive selection using antibodyconjugated magnetic beads (Milteny Biotec). Neuronal cells were enriched by a negative depletion of nonneuronal cells using the (mouse) Neuron Isolation Kit (Milteny Biotec). Briefly, up to 10 7 dissociated cells were resuspended in 80 μL of cold Dulbecco's phosphatebuffered saline containing 0.5% bovine serum albumin (DPBS-BSA buffer), incubated with 10 μL FcR Blocking Buffer for 10 min, and then with 10 μL ACSA-2 MicroBeads for 15 min, always at 4°C. After washing, cells were centrifuged to remove excess beads from the solution. The pellet was suspended with 500 μL of DPBS-BSA buffer and the suspension was applied to the appropriate MACS column fitted in MACS Midi magnetic cell separator. The flow-through containing the unlabeled negative fraction was collected for subsequent cell separation.
Following column removal from the magnetic separator, bound astrocytes were flushed out with 1 mL of DPBS-BSA, centrifuged, and the pellet stored at − 80°C for further analysis. Cell number in the negative fraction was determined and up to 10 7 dissociated cells were resuspended in 80 μL DPBS-BSA. Then, they were first incubated with 20 μL of Non-Neuronal Cell Biotin-Antibody Cocktail for 5 min and then with 20 μL di Anti-Biotin MicroBeads for 10 min, always at 4°C. After resuspension in 500 μL buffer, the solution was eluted on MACS column. The flow-through containing the negative fraction of unlabeled neuronal cells was centrifuged and the pellet stored at − 80°C. Following column removal from the magnetic separator, bound non-neuronal cells (MOE fraction: mainly microglia, oligodendrocytes, and endothelial cells) were eluted in 1 mL buffer, centrifuged and stored. Cell identity and purity of astrocyte and neuronal fractions was confirmed in immunocytochemical experiments (data not shown).

DNA extraction from tissues
DNA was extracted using Dneasy Blood & Tissue Kit (Qiagen, Hilden, Germany), following the manufacturer's instructions. DNA was quality checked using NanoDrop 2000, (Thermo Scientific) and quantified using Qubit 2.0 Fluorometer with the dsDNA broad range assay kit (Invitrogen, Q32850).

RNA extraction
Total RNA was extracted by RNeasy mini kit (QIAGEN) according to the manufacturer's instructions. The integrity of the RNA was assessed by denaturing agarose gel electrophoresis (presence of sharp 28S, 18S, and 5S bands) and spectrophotometry (NanoDrop 2000, Thermo Scientific). Total RNA was purified to eliminate potentially contaminating genomic DNA using recombinant DNase (QIAGEN).

Bisulfite and oxidative bisulfite conversion and amplicon library preparation
Genomic DNA (1 μg) was converted by sodium bisulfite with EZ DNA Methylation Kit (Zymo Research) and eluted in 50 μL H 2 O according to the manufacturer's instruction. Oxidative bisulfite experiments for hydroxymethylcytosine detection at single nucleotide level were performed using TrueMethyl oxBS module (Nugen, Tecan, California USA) following the manufacturer's instruction). Quality and amount (1 μg) of genomic DNA was evaluated by Nanodrop and Qubit instruments. DNA was first purified with magnetic beads and then treated with oxidant solution. Both bisulfite-and ox-bisulfite-treated samples underwent to double amplification strategy to generate an amplicon library. The first PCR step was performed using bisulfite-specific primers. Reactions were performed in 30 μL total volumes: 3 μL 10× reaction buffer, 0.6 μL of 10 mM dNTP mix, 0.9 μL of 5 mM forward and reverse primers, 3.6 μL MgCl2 25 mM, 2-4 μL bisulfite template DNA, 0.25 μL FastStart Taq, and H 2 O up to the final volume. Sequences of gene-specific primers, along with individual PCR conditions, are reported in Table 1. The second step of PCR was used to add multiplexing indices and Illumina sequencing adaptors to first amplicons. Second PCR step was performed in 50 μL total volumes: 5 μL 10× reaction buffer, 1 μL dNTP mix, 5 μL forward and reverse "Nextera XT" primers (Illumina, San Diego, CA), 6 μL 25 mM MgCl 2 , 5 μL of first PCR product, 0.4 μL FastStart Taq, and H 2 O up to the final volume. After both PCR steps, a purification phase was performed using AMPure purification magnetic Beads (Beckman-Coulter, Brea, CA) following the manufacturer's protocol. All amplicons were quantified using Qubit® 2.0 Fluorometer. A library of amplicons derived from bisulfite-and ox-bisulfitetreated genomic DNAs was obtained pooling amplicons at equimolar ratio. Library was diluted to final concentration of 8 pM. Phix control libraries (Illumina) were combined with normalized library [10% (v/v)] to increase diversity of base calling during sequencing. Amplicons' library was subjected to sequencing using V3 reagents kits on Illumina MiSeq system (Illumina). Paired-end sequencing was performed in 281 cycles per read (281 × 2). An average of approximately 100,000 reads/sample were obtained. To estimate the rate of bisulfite conversion, fully unmethylated M13mp18 double-strand DNA (New England BioLabs) was added in representative samples.

Sequence handling and bioinformatics analyses
Paired-end reads were obtained from Illumina Miseq sequencer platform. Using PEAR tool [73], reads were assembled together with a minimum of 40 overlapping residues as threshold. FASTQ assembled reads were then converted in FASTA format using PRINSEQ tool [74]. Sequences derived from bisulfite-and ox-bisulfite-treated DNAs were analyzed with ampliMethProfiler pipeline software (https://sourceforge.net/projects/amplimethprofiler) [36], specifically designed for deep-targeted bisulfite amplicon sequencing. ampliMethProfiler produces quality filtered FASTA files for all samples and directly extracts methylation average and methylation profiles. To determine the average 5-hmC levels, the average values at single-CpGs obtained from sequence analysis of DNAs treated with oxidant plus bisulfite were subtracted to the average methylation values obtained from DNA treated with bisulfite only. As output, the pipeline generates also a summary file with information about the number of reads passing filters, the methylation percentage of each C in CpG sites, and the bisulfite efficiency for each C in non-CpG sites. ampliMeth-Profiler produces a tabular format file (BIOM format) containing the number of methylation profiles (epialleles) for all samples. The BIOM table was normalized for the same number of sequence/sample through a rarefaction procedure using QIIME [75]. Rarefied BIOM table was then processed using JMP to obtain principal component analysis (PCA) and hierarchical cluster analysis (HCA) based on epialleles distribution.

Statistical analyses
Methylation average data are expressed as means ± standard error. Comparisons between groups were performed using one-way ANOVA with α significance level ≤ 0.001. mRNA expression levels are reported as 2 -ΔCt and analyzed by one-way ANOVA. All statistical analyses were performed using GraphPad (GraphPad Prism Software, Inc., La Jolla, CA, USA www.graphpad.com/ guides/prism/7/statistics/index.htm).

Additional files
Additional file 1: Figure S1. Ddo epiallelic distribution analysis in all analyzed brain regions during ontogenesis. Hierarchical cluster based on the epiallelic composition of each developmental stage is presented for HIPP, CX and CB. Developmental stages are indicated with different colors (P1 = red; P15 = blue; P30 = yellow; P60 = green). (JPG 87 kb) Additional file 2: Figure S2. Regression analysis associating the degree of methylation with mRNA expression levels of the Ddo gene. Regression analysis associating DNA methylation and mRNA expression of Ddo gene is shown in each analyzed brain area during development. Ddo mRNA expression is normalized to the mean values for two housekeeping genes and expressed as 2 -ΔCt values. The three mice for each time point (P1 = red; P15 = blue; P30 = yellow; P60 = green) are indicated. Statistical analyses were performed using Pearson correlation. * p ≤ 0.05 (PDF 86 kb)