- Research article
- Open Access
Widespread promoter methylation of synaptic plasticity genes in long-term potentiation in the adult brain in vivo
BMC Genomicsvolume 18, Article number: 250 (2017)
DNA methylation is a key modulator of gene expression in mammalian development and cellular differentiation, including neurons. To date, the role of DNA modifications in long-term potentiation (LTP) has not been explored.
To investigate the occurrence of DNA methylation changes in LTP, we undertook the first detailed study to describe the methylation status of all known LTP-associated genes during LTP induction in the dentate gyrus of live rats. Using a methylated DNA immunoprecipitation (MeDIP)-array, together with previously published matched RNA-seq and public histone modification data, we discover widespread changes in methylation status of LTP-genes. We further show that the expression of many LTP-genes is correlated with their methylation status. We show that these correlated genes are enriched for RNA-processing, active histone marks, and specific transcription factors. These data reveal that the synaptic activity-evoked methylation changes correlates with pre-existing activation of the chromatin landscape. Finally, we show that methylation of Brain-derived neurotrophic factor (Bdnf) CpG-islands correlates with isoform switching from transcripts containing exon IV to exon I.
Together, these data provide the first evidence of widespread regulation of methylation status in LTP-associated genes.
Long-term potentiation (LTP), induced by high-frequency stimulation (HFS), is a model for synaptic plasticity recapitulating the neuronal response to learning . Several studies have shown early changes in kinase activity, transcription factors, mRNA, miRNA, and lncRNAs in response to LTP [2–7]. However, the involvement of the epigenome in LTP has largely been overlooked .
Methylation of DNA, which involves the addition of a methyl group to cytosine (5mC, 5-methylcytosine), is an important regulating process in promoters and CpG-islands where it can modulate gene expression (reviewed in ). Although increased DNA methylation often leads to decreased expression, the relationship between methylation and expression is not fully understood .
DNA methylation has been shown to be required for mammalian development , and for neuronal cell differentiation [12, 13]. Interestingly, studies have also shown that DNA methylation plays a role in the adult brain. For example, knockout of de novo DNA methylation proteins, such as Dnmt1 and Dnmt3a, led to learning deficits in mice . Changes to DNA methylation in Bdnf, an important synaptic plasticity gene, have been shown to regulate Bdnf expression , and remodel chromatin after Bdnf stimulation , a stimulus known to evoke synaptic plasticity. DNA methylation is also a feature of contextual learning  and electroconvulsive stimulation  where CpG-islands were identified as having changed methylation status, revealing a potential role for methylation/demethylation in post-synaptic neurons.
Although DNA-methylation has been known to be involved in LTP , to date, no study has attempted to observe the scale of methylation changes in LTP. In a previous study of RNA expression, we observed robust increases in expression of Tet3 , an enzyme involved in DNA demethylation. Tet3 oxidises 5mC to 5hmC, which is the first step in the removal of the methyl group from 5mC [20, 21] and is required in homeostatic plasticity . Whereas all Tet-proteins catalyse the same reaction , an expression analysis showing an increase in only Tet3 suggests it is the primary contributor to the oxidation . Taking these observations together, we set out to investigate the role of methylation in LTP.
Using a rat model, we present the first study exploring changes in DNA methylation in LTP-associated genes in response to HFS. An analysis of methylation state in relation to gene expression and chromatin remodeling suggests that these genes are already primed for transcription by active chromatin. Together, these data suggest an active regulatory role for DNA methylation and histone modification in synaptic plasticity and memory formation.
Long-term potentiation induction through HFS
The long-term potentiation was induced as described previously . Briefly, the left Dentate gyrus (DG) in the hippocampus of urethane-anesthetised rats were subjected to high-frequency stimulation (HFS), using a well-established pattern paradigm, of medial perforant path input to granule cells of the DG. This stimulation resulted in a sharp and lasting increase of the perforant path-evoked field excitatory post-synaptic potential. Matched RNA-seq data was downloaded from ArrayExpress (E-MTAB-3375) and analysed to find 1,110 unique differential expressed genes using a stringent expression threshold described in Methods. While the RNA-seq study utilised the matched contralateral unstimulated DG, this study used baseline test-pulse stimulation as the control group and compared the methylation changes to 3 different time points after HFS. These time points were 30 min, 2 h, and 5 h post-HFS. Only the 6 stimulated samples were matched between the methylation array and the previous RNA-seq study. The expression data from the unstimulated naïve control group from the RNA-seq study was used as a surrogate to investigate the relation between expression and methylation in the base-line test stimuli group.
Design of methylation array of LTP associated gene promoters
To examine the involvement of methylation changes in LTP, we designed a MeDIP-array based on the differentially expressed genes described above. Briefly, probes were designed for the promoter and gene-associated CpG-islands and CpG-shores for each differentially expressed gene, covering a domain ±3 kb around the transcription start site (TSS) (Fig. 1a). The design led to each probe hybridising to a specific region characterised by a CpG-island, CpG-shore, promoter, or a combination of the above (Fig. 1b). The resulting merged domains consisted of 5612 unique regions surrounding the promoter of the LTP-associated genes previously described.
Methylation array reveals early changes in methylation around LTP genes
To investigate changes in methylation status resulting from HFS, we used limma  to determine the differential methylation status of probes (DMP) or collapsed regions (DMR) (see Methods) between HFS stimulated samples (30 min, 2 h, and 5 h) compared to test-pulse stimulation controls.
Few changes (11) in probe intensity were observed after 30 min while 1919 and 1474 probes were differentially methylated (FDR ≤ 0.25) at 2 and 5 h, respectively (Fig. 2a). Interestingly, when investigating differentially methylated regions, 48 regions displayed methylation changes after 30 min while 699 and 448 regions showed changes after 2 and 5 h (Fig. 2b).
To explore the direction of methylation changes, probes and regions were divided into their specific locations revealing an overall mean loss of methylation in CpG-associated regions, while methylation in promoters and shores increased (Additional file 1: Figure S1a, b). Overall, the majority (88% at 30 min, 71% at 2 h, and 75% at 5 h) of DMRs lose part of their methylation with stimulation (Fig. 2c; p-value ≤ 0.05). For example, immediate early responder gene Arl4d , is highly differentially expressed at 30 min, displays a loss of methylation in the promoter region. The differentially methylated regions can be found in Additional file 2: Table S1.
The loss of methylation after HFS is consistent with our previous study  showing upregulation of the majority of genes after HFS, as methylation is mainly known to negatively correlate with gene expression. However, the loss of methylation at 30 min is not explained by the increase of Tet3 expression, which increases only 2 h after HFS (Additional file 3: Figure S2). Moreover, the increase in DNA-methylation could be attributed to the increased expression of Dnmt3a (Additional file 3: Figure S2), which has previously been shown to be required for memory formation .
We further investigated the specificity of the methylation response in LTP by comparing the expression of our studied genes with changes reported in response to a single electroconvulsive stimulus , and contextual fear conditioning . Although our studied genes only made up a subset of genes having differential methylation in their proximity in the other studies, 24% (110/463 genes) of our differentially methylated genes displayed methylation changes in response to other stimuli (Fig. 2e, f).
Taken together, these data suggest that changes to methylation occur in LTP, revealing a high specificity of methylation and demethylation events post-HFS.
Methylation changes correlates with gene expression, and correlated genes are enriched for RNA processing and binding
Whereas site-specific differential methylation reveals changes at a certain moment, this analysis misses trends observed over multiple time points. Furthermore, it is important to investigate the relationship between methylation and expression over time. Previous studies have shown that the loss of methylation correlates with increases in gene expression, although the direction of correlation can be uncertain .
To investigate how methylation changes drive gene expression in LTP, we used the matched RNA-seq (E-MTAB-3375) counts per million reads (cpm), and the resulting logFC changes between each time point, from our previous study . Since our methylation study utilised a baseline test stimuli as control instead of the contralateral DG as in the RNA-seq study, we chose the unstimulated naïve rat control as expression control. Comparison between the logFC between expression and methylation reveals that most genes/regions are negatively correlated (Additional file 4: Figure S3a). In contrast, when correlating the normalised expression with the normalised methylation values, we observed 825 regions that showed no correlation, 1605 regions that showed low correlation, 2417 regions that showed intermediate correlation, and 485 regions that showed high correlation between methylation changes and expression changes (Additional file 4: Figure S3a). When also accounting for p-value < 0.05, we observe only 455 regions showing high correlation with their corresponding gene, the majority (274/455, or 60%) of which are negatively correlated (Fig. 3b and c). Two genes that showed high negative correlation between expression and methylation are Arl4d (Fig. 3d) and Arc (Fig. 3e), both immediate early responder genes (IEGs), which show increased expression after only 30 min post-HFS. In Arc, the promoter associated CpG-island has previously been described to have altered methylation changes in response to spatial behavioural tasks in the hippocampus . Moreover, several other genes show high linear correlation between their methylation and expression (Additional file 5: Figure S4), such as Egr2-4 (Additional file 5: Figure S4a), Etv5, Zswim5 and Ube2g1 (Additional file 5: Figure S4b), and genes with later stage expression such as Stam, Bhlhe40, and Jak2 (Additional file 5: Figure S4c). Highly correlated (absolute(r) ≥ 0.8) genes showed enrichment for RNA-processing (Fig. 3f) and 3′-UTR binding (Fig. 3g) compared to the rest of the genes in the array.
In summary, the majority of observed methylation changes post-HFS correlate with gene expression negatively whereas loss of methylation correlates with increased gene expression.
Expression/methylation-correlated genes are enriched for motifs involved in brain function and development
To further understand the basis for the inconsistency between promoter methylation and gene expression, we examined the difference in the characteristics between LTP-associated genes where promoter methylation and expression were correlated and those that were not. We hypothesised that a decrease in methylation around the promoter could facilitate transcription factor binding leading to increased gene expression.
In support with of this hypothesis, a motif search of expressed transcription factors between correlated and non-correlated gene promoters led to the discovery of enriched (FDR ≤ 0.05) CG-rich motifs 1 kb upstream of TSS in correlated genes corresponding to Tcfap2c, E2F3, Klf4, SP4, and Egr1 (Fig. 4a). We observed a higher fraction of CpG-islands, and a small but statistically significant (μ = 0.53 vs. 0.49 correlated vs. non-correlated, p-value = 1.6e-08) increase of GC-content in our correlated vs. non-correlated genes (Fig. 4b, c), potentially confounding our motif analysis. To correct for GC-content and CpG-islands in this analysis, we only compared motifs present in promoters overlapping a CpG-islands in our two groups (Additional file 6: Figure S5). Although, some of these motifs were no longer enriched after filtering for promoters overlapping a CpG-island, Egr1 and SP4 remained (Additional file 6: Figure S5a). These genes are involved in brain function as drivers in synaptic plasticity and as a regulator of AMPA receptor subunits, respectively [27, 28]. In support of these findings, after filtering the number of CpG-island per transcripts and the CG-frequency of the promoters were almost identical (Additional file 6: Figure S5b), although, a slight enrichment in CG-content in correlated genes still remained (Additional file 6: Figure S5c). To further investigate the bias in GC-content and our results, we investigated the relationship between GC-content in probes (Additional file 7: Figure S6a), regions (Additional file 7: Figure S6b), and the methylated regions correlation with gene expression (Additional file 7: Figure S6c). We observed no clear relationship between the GC-content and correlation between methylation and expression changes. This suggests that these motifs found in this analysis are not relics of the increased CpG-islands, but are important components of LTP, as exemplified by SP4 and Egr1.
Active chromatin marks show enrichment over expression/methylation-correlated genes
To further investigate the differences between the expression/methylation-correlated and non-correlated genes, we wanted to explore the chromatin state surrounding the genes in both groups. The chromatin landscape regulates gene expression through modifications of the histone tails. To address this question, we examined publically available datasets reflecting the normal adult chromatin state of rat brain (see Methods). Genes correlated with their methylation status had a higher enrichment of active histone marks, such as H3K4me1, H3K27ac, H3K4me2, and H4K5ac, around their TSS compared to non-correlated genes (Fig. 5a). Moreover, expression/methylation-correlated genes have enrichment of Matrin3 and Pol II, while they are depleted in 5hmc modifications, and H3K9me3 around the TSS (Fig. 5a). These enrichments are also present over the whole gene body (Fig. 5b), although the highest difference appears around the TSS, which drives the significance in H4K5ac, Matrin 3, and Pol II. Depletion of either H3K9me3 or 5hmc was not found when investigating the whole gene body (Fig. 5b).
To investigate any bias in our histone mark analysis, we compared the expression distribution between expression/methylation-correlated and non-correlated genes. We observe an already increased expression in the naïve rat dentate gyrus of the expression/methylation-correlated genes (Additional file 8: Figure S7) (Control; p-value = 0.014, meancorrelated = 21 RPKM, meannoncorrelated = 14 RPKM). This could potentially explain the enrichment of active chromatin marks in the expression/methylation-correlated group in naïve brain.
This expression difference increases with stimulation and time (Additional file 8: Figure S7) (30 min: meancorrelated = 35 RPKM, meannoncorrelated = 23 RPKM; 2 h: meancorrelated = 44 RPKM, meannoncorrelated = 26 RPKM; 5 h: meancorrelated = 50 RPKM, meannoncorrelated = 29 RPKM) showing that the correlated genes have a sharper expression increase than the non-correlated genes. Taken together, these data suggest that genes with expression correlated with their methylation status are already transcriptionally more active than non-correlated genes, have an enrichment of active chromatin marks, and have a higher proportion of enriched neuro-specific motifs in their promoters.
CpG-island methylation correlates with the isoform switching observed in the Bdnf loci in LTP
Methylation at different promoters of the same gene can influence which isoform is expressed. To examine whether such a mechanism was employed in LTP, we focused on the expression of the Bdnf gene, which plays an important role in synaptic plasticity [16, 29]. Previously, demethylation of Bdnf has been shown to occur in the prefrontal cortex in response to LTP . The existence of multiple Bdnf isoforms is well known , and switches in expressed Bdnf isoforms in fear conditioning is due to their differential methylation status . The rat Bdnf loci consists of three CpG-islands, for which we designed probes along with corresponding CpG-shores (Fig. 6a). In our LTP model, 4 major isoforms of Bdnf are observed, with low expression of Bdnf in the control samples, one isoform lacking expression after 30 min while all isoforms are expressed after 2 h and 5 h post-HFS (Fig. 6b). These isoforms corresponds to exon I, exon II, exon IV, and exon VI . Investigating these isoforms with their corresponding methylation status, we found two isoforms (exon I and exon IV) that show expression patterns consistent with their methylation status, suggesting that methylation may be employed here as a mechanism for isoform selection (Fig. 6c). The other two major isoforms (exon II and exon VI) expressed showed no correlation between methylation and expression, which suggests expression of these isoforms may be regulated by another mechanism.
The regulation of gene expression is a highly complex and tightly controlled process involving DNA methylation, transcription factors, and chromosome organisation.
The role of DNA methylation in LTP has so far been understudied. Here we present the first large-scale methylation study of temporal methylation changes in LTP and correlate the methylation status with matched RNA-seq data for LTP-associated genes. We find that only a few regions show rapid changes in DNA-methylation 30 min after LTP induction, while the majority of changes happen in the intermediate (2 h) to late term (5 h). We also show that genes whose expression correlates with their proximal methylation are enriched for RNA-binding and RNA-processing. Furthermore, we reveal specific transcription factor motifs are enriched in correlated gene promoters, independent of CpG-islands, and we show that these genes are enriched for active chromatin marks in steady state, indicating that these genes already show signs of transcription which are further altered by methylation changes.
The most common changes to promoter methylation is demethylation, supporting our previous observation of enhanced expression for the majority of LTP associated genes . Although methylation levels are generally negatively correlated with expression, this is not always the case . We observed the expression of genes that were both negatively and positively correlated with their methylation level. Although we do not observe any specific relationship between the location of the methylation changes and the correlation, changes may occur in regions not covered by our array. For example, the methylation of gene bodies can influence the expression levels . Moreover, methylation changes in enhancer regions are associated with downstream changes to gene expression . Further studies combining whole-genome bisulphite sequencing with ChIP-seq will be able to clarify the relationship.
Transcription factor binding sites (TFBS) in the promoters are important regulators of gene expression, and DNA methylation changes at TFBS can alter their binding , although it is not considered a general property for transcription factor binding . Here we observe certain expressed transcription factors that are enriched in genes correlated with methylation changes. Specifically, we see enrichment of Egr1, which is highly expressed 30 min post-HFS. Moreover, we observe SP4 and Nfe2I2, two transcription factors involved in expression of AMPA  and NMDA  type glutamate receptors, respectively. These transcription factors are known to have transcriptional control over important receptors involved in LTP. Here we show that their motifs are also enriched in the correlated genes, suggesting that they also control expression of many identified LTP-associated genes. We also observe enrichment of Tcfap2c, a transcription factor forming a complex with Cited2, which regulates Cebpa expression . Cited2 is a gene highly expressed 30 min post-HFS , and, together with Cebpb (differentially expressed after 2 and 5 h post HFS), is a neuronal activity dependent transcription factor . However, the enrichment of Tcfap2c is possibly artifactual, as the enrichment diminished when selection was restricted to genes containing CpG-islands. Other TFs enriched in the promoter of correlated genes are Foxj3, Mafb, Srebf1, and a depletion of Sox10, all which have been linked to either brain development or differentiation [39–42].
Furthermore, we investigated the steady-state chromatin landscape between the correlated and non-correlated genes. Through this analysis, in correlated genes we observed enrichment of active chromatin marks, such as H3K4me1, H3K227ac, H3K4me2, H4K5ac and Matrin3 and PolII, and a slight decrease of inactive chromatin marks near the TSS, as shown by H3K9me3. The enrichment of active chromatin marks was also present over the whole gene body. This suggests that the genes responding to methylation changes are already primed for transcription by active chromatin.
Interestingly, the expression/methylation correlated genes showed a depletion of 5hmC compared to non-correlated genes, suggesting that the correlated genes are more methylated around their TSS than non-correlated genes. This supports our previous observation that Tet3 is responsible for the demethylation of 5mC in response to HFS.
We also provide evidence that CpG-island methylation around gene promoters correlates with the expression of specific Bdnf isoforms from transcripts containing exon IV to transcripts with exon I. However, this did not explain the expression changes of Bdnf isoforms containing exon II or exon VI.
Although we observe methylation changes to be a widespread phenomenon in LTP, the technique used, the limited sample size, and examined regions in this study suggests that the results provide an incomplete representation of the actual methylation changes in LTP. As our MeDIP-array relies on antibodies to capture the methylated regions, we are unable to get an absolute measurement of the methylation changes on a base resolution. Furthermore, the MeDIP-array can only distinguish large changes in methylation, suggesting we might miss more subtle changes. Instead, utilising whole-genome or targeted bisulphite sequencing should provide a more accurate representation, and provide detail of even miniscule difference, of the methylation changes observed in LTP.
As our study only covers the first few hours after HFS, we cannot say if the observed DNA methylation changes persists in the neurons or if the DNA methylation is more transient and goes back to normal later after HFS. Although we believe that the observed changes to DNA methylation in synaptic plasticity genes go back to normal after the required genes have been expressed, longer time points are needed to elucidate for how long these changes persists. Finally, further studies are needed to investigate how the methylation, hydroxymethylation, and chromatin landscape change in LTP in response to HFS.
In conclusion, we show changes to the methylation landscape around synaptic plasticity genes during the first hours after LTP induction. Furthermore, we show that the histone landscape around correlated genes is already primed for transcription in steady-state, suggesting that an active chromatin landscape is required for DNA methylation to occur. This first large-scale study of DNA-methylation in LTP suggests a major role for epigenetic modifications in the activity-dependent control of gene expression underlying long-term synaptic plasticity in the adult mammalian brain.
The animals and electrophysiology were the same as describe in previous study . Briefly, after HFS to the dentate gyrus (DG) animals were decapitated at 30 min, 2 h, and 5 h, and the DG was micro-dissected and snap-frozen on dry ice. For the control group, a baseline test stimuli was applied instead of the HFS. This study used two animals per group. The field excitatory post-synaptic response for the different groups used in this study can be seen in Additional file 9: Figure S8a. After tissue homogenisation, the DNA was extracted using Allprep RNA/DNA isolation kit (Qiagen).
Each differentially expressed gene from previous publication  with modified cut-off of absolute log2FC ≥ 0.9 and a FDR ≤ 0.05 was used in the array design. The FDR was obtained through the Benjamini Hochberg multiple testing corrections of p-values. Additionally, genes not annotated in Ensembl Rn4 69 release but previously characterised in LTP, such as Bdnf and Jarid2 were reanalysed using Refseq and included in the array along with the gene associated CpG-islands and shores. The coverage of each promoter with an associated CpG-island was set to ±1 kb from the TSS. Each CpG-island lacking gene was only covered at 0.5 kb upstream and downstream from the TSS, while the CpG-island coverage was extended with ±1.5 kb in each direction to accommodate for CpG-shores.
The DNA methylation microarray was prepared following the Agilent microarray analysis of methylated DNA immunoprecipitation protocol V2.2. 5-methylcysteine antibody (A-1014-050) was ordered from Epigenetek. The hybridisation of samples for the array was conducted at the Ramaciotti Centre for Genomics at UNSW, where the arrays also were scanned.
The arrays were read and analysed using limma . Briefly, loess normalisation was used to normalise within arrays while quantile normalisation was used between the two arrays. Control probes were removed before analysis. Before quantile normalisation, the probe intensity of each sample was compared, and MDS plots of quantile normalised probes and regions were plotted (Additional file 9: Figure S8b-d). Each probe was given a unique location ID based on their respective region as in Fig. 1b. A contrast matrix was set up for differential methylation analysis between the control and stimulated samples, and probes were called differentially methylated if they had a p-value ≤ 0.05 and a FDR ≤ 0.25. For the region analysis, probes were averaged over their respective region as above in accordance with Fig. 1b. The regions were called differentially methylated regions (DMR) if they had a p-value ≤ 0.05 and a FDR ≤ 0.25.
Comparing LTP differentially methylated genes with other models of neuronal stimulation
Gene lists with closest gene from differentially methylated areas were downloaded from published studies investigating the methylome changes of synchronous activation by electroconvulsive stimulation , and contextual fear condition with and without shock  from their respective supplementary material.
Since the downloaded gene lists contained mice genes, these were converted to rat gene names by first converting the mice gene names into orthologous rat ensembl ids using biomaRt. After that, the rat ensembl ids were converted to corresponding rat name using the ensembl 69 release. All differentially methylated regions were collapsed into their corresponding gene.
Only genes included in our array were investigated for the other stimuli. The resulting intersections were plotted using UpsetR .
Gene expression values (cpm) was imported from previous study . Each gene expression profile was correlated with its respective methylated region using the R-package Hmisc  with Pearson correlation. The effect of the correlation between gene expression and methylation was designated as following: r < 0.1 = no effect, 0.1 < r < 0.30 = small effect, 0.30 < r < 0.70 = intermediate effect, r > 0.70 = large effect. Only genes with a r > 0.70 and a P-value < 0.05 were called correlated in this study. Gene enrichment was done using TopGO on genes with a r ≥ 0.8. These highly correlated genes were compared to non-correlated genes in the array as a background. Both the biological and molecular function enrichment was conducted using the classic algorithm and fisher statistics from the TopGO package. For comparison of LogFC values between expression and methylation, logFC expression values were imported from Maag et al. 2015, and plotted against logFC methylation between each time point and control.
Motif discovery and GC-content
Investigation for motifs present in correlated vs. non-correlated genes was done using PWMEnrich . For the promoter region comparison region 1 kb upstream from the promoter were investigated for motifs. A background object was created using makePWMLognBackground for all Mus musculus motifs for the non-correlated promoters. For enhancers, the background was created the same way except for using enhancers neighbouring non-correlated genes as background. The promoter GC-content was calculated to investigate GC-bias, and the groups were compared using the student t-test. For the motif search, the correlated genes were enriched compared to the non-correlated genes using motifEnrichment with each specific background. Motifs were summarised using the groupReport function, and the top 20 motifs were plotted. Motif names were extracted from the R-package MotifDB. The P-values were adjusted using Benjamini & Hochberg adjustment.
Probe and region GC-bias
All probes and regions were investigated for GC-bias between differentially methylated and non-differentially methylated probes or regions. The t-statistics from limma was plotted against the GC-content for each differential methylation comparison. Furthermore, the GC-content was also compared to the correlation between expression and methylation.
Histone marks and transcription factors analysis over correlated genes
Briefly, the reads from public ChIP-seq control fastq files were aligned to the rat genome using Bowtie2 . The aligned reads were then converted to bam, sorted, and indexed with Samtools . Ngsplot  was used to calculate the average coverage of correlated and non-correlated genes over TSS ± 3 kb or the whole gene body ± 3 kb from TSS and TES. Significance of enrichment was calculated using the Wilcox test for log2 loci enrichment. The p-value for each histone mark was adjusted to account for multiple testing using Benjamini-Hochberg False discovery rate correction. H3K27ac, H2K4me2, H3K4me1, and Matrin 3 was obtained from a growth hormone (GH)- expressing rat pituitary cell line , H3K9me3 from rat hippocampal tissue , PolII from rat neuronal cultures , 5hmc-capture from rat brainstem  and H4K5ac from rat dorsal striatum .
Expression differences between expression/methylation-correlated and non-correlated genes
To investigate the expression differences between the correlated and non-correlated genes, the RNA-seq count data was first converted to RPKM using the gene length from featureCount of the Rattus_norvegicus.RGSC3.4.69.gtf. This allowed for comparison between the two groups. Student’s t-test was used to observe significance between the correlated and non-correlated groups for each time point.
Bdnf isoform switching by methylation
Due to Bdnf being absent from the Ensembl 69 Rn4, we counted reads with Rsem  using Refseq database for Rn4. Using the TPM from the 4 most prominent isoforms, we calculated the percentage of each isoform from the total gene expression. The TPM for each isoform was then correlated with the methylation of their closest CpG-island on the Bdnf loci to observe isoform-specific methylation events.
Differentially methylation probes
Differentially methylation regions, cpm, counts per million
High-frequency stimulation, TSS, transcription start site
Immediate early responder genes
Methylated DNA immunoprecipitation
Transcription factor binding sites
Whitlock JR, Heynen AJ, Shuler MG, Bear MF. Learning induces long-term potentiation in the hippocampus. Science. 2006;313:1093–7.
Malinow R, Madison DV, Tsien RW. Persistent protein kinase activity underlying long-term potentiation. Nature. 1988;335:820–4.
Jones MW, Errington ML, French PJ, Fine A, Bliss TVP, Garel S, et al. A requirement for the immediate early gene Zif268 in the expression of late LTP and long-term memories. Nat Neurosci. 2001;4:289–96.
Wibrand K, Panja D, Tiron A, Ofte ML, Skaftnesmo KO, Lee CS, et al. Differential regulation of mature and precursor microRNA expression by NMDA and metabotropic glutamate receptor activation during LTP in the adult dentate gyrus in vivo. Eur J Neurosci. 2010;31:636–45.
Maag JLV, Panja D, Sporild I, Patil S, Kaczorowski DC, Bramham CR, et al. Dynamic expression of long noncoding RNAs and repeat elements in synaptic plasticity. Front Neurosci. 2015;9:351.
Håvik B, Røkke H, Dagyte G, Stavrum AK, Bramham CR, Steen VM. Synaptic activity-induced global gene expression patterns in the dentate gyrus of adult behaving rats: Induction of immunity-linked genes. Neuroscience. 2007;148:925–36.
Ryan MM, Mason-Parker SE, Tate WP, Abraham WC, Williams JM. Rapidly induced gene networks following induction of long-term potentiation at perforant path synapses in vivo. Hippocampus. 2011;21:541–53.
Sweatt JD. The Emerging Field of Neuroepigenetics. Neuron. 2013;80:624–32.
Jones PA. Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet Nature Publishing Group. 2012;13:484–92.
Wan J, Oliver VF, Wang G, Zhu H, Zack DJ, Merbs SL, et al. Characterization of tissue-specific differential DNA methylation suggests distinct modes of positive and negative gene expression regulation. BMC Genomics BioMed Central. 2015;16:49.
Okano M, Bell DW, Haber DA, Li E. DNA Methyltransferases Dnmt3a and Dnmt3b Are Essential for De Novo Methylation and Mammalian Development. Cell. 1999;99:247–57.
Takizawa T, Nakashima K, Namihira M, Ochiai W, Uemura A, Yanagisawa M, et al. DNA methylation is a critical cell-intrinsic determinant of astrocyte differentiation in the fetal brain. Dev Cell. 2001;1:749–58.
Mohn F, Weber M, Rebhan M, Roloff TC, Richter J, Stadler MB, et al. Lineage-specific polycomb targets and de novo DNA methylation define restriction and potential of neuronal progenitors. Mol Cell. 2008;30:755–66.
Feng J, Zhou Y, Campbell SL, Le T, Li E, Sweatt JD, et al. Dnmt1 and Dnmt3a maintain DNA methylation and regulate synaptic function in adult forebrain neurons. Nat Neurosci. 2010;13:423–30.
Lubin FD, Roth TL, Sweatt JD. Epigenetic Regulation of bdnf Gene Transcription in the Consolidation of Fear Memory. J Neurosci. 2008;28:10576–86.
Martinowich K, Hattori D, Wu H, Fouse S, He F, Hu Y, et al. DNA methylation-related chromatin remodeling in activity-dependent BDNF gene regulation. Science. 2003;302:890–3.
Halder R, Hennion M, Vidal RO, Shomroni O, Rahman R-U, Rajput A, et al. DNA methylation changes in plasticity genes accompany the formation and maintenance of memory. Nat Neurosci Nature Publishing Group. 2015;19:102–10.
Guo JU, Ma DK, Mo H, Ball MP, Jang M-H, Bonaguidi MA, et al. Neuronal activity modifies the DNA methylation landscape in the adult brain. Nat Neurosci Nature Publishing Group. 2011;14:1345–51.
Levenson JM, Roth TL, Lubin FD, Miller CA, Huang IC, Desai P, et al. Evidence that DNA (cytosine-5) methyltransferase regulates synaptic plasticity in the hippocampus. Journal of Biological Chemistry. 2006;281:15763–73.
Iqbal K, Jin S-G, Pfeifer GP, Szabó PE. Reprogramming of the paternal genome upon fertilization involves genome-wide oxidation of 5-methylcytosine. Proc Natl Acad Sci U S A. 2011;108:3642–7.
Ito S, D’Alessio AC, Taranova OV, Hong K, Sowers LC, Zhang Y. Role of Tet proteins in 5mC to 5hmC conversion, ES-cell self-renewal and inner cell mass specification. Nature Nature Publishing Group. 2010;466:1129–33.
Yu H, Su Y, Shin J, Zhong C, Guo JU, Weng Y-L, et al. Tet3 regulates synaptic transmission and homeostatic plasticity via DNA oxidation and repair. Nat Neurosci Nature Publishing Group. 2015;18:836–43.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43:e47–7.
Wibrand K, Messaoudi E, Håvik B, Steenslid V, Løvlie R, Steen VM, et al. Identification of genes co-upregulated with Arcduring BDNF-induced long-term potentiation in adult rat dentate gyrus in vivo. European Journal of Neuroscience. 2006;23:1501–11.
Morris MJ, Adachi M, Na ES, Monteggia LM. Selective role for DNMT3a in learning and memory. Neurobiol Learn Mem. 2014;115:30–7.
Penner MR, Roth TL, Chawla MK, Hoang LT, Roth ED, Lubin FD, et al. Age-related changes in Arc transcription and DNA methylation within the hippocampus. Neurobiol Aging. 2011;32:2198–210.
Priya A, Johar K, Nair B, Wong-Riley MTT. Specificity protein 4 (Sp4) regulates the transcription of AMPA receptor subunit GluA2 (Gria2). Biochimica et Biophysica Acta (BBA) - Molecular Cell Research. 2014;1843:1196–206.
Worley PF, Christy BA, Nakabeppu Y, Bhat RV, Cole AJ, Baraban JM. Constitutive Expression of Zif268 in Neocortex Is Regulated by Synaptic Activity. Proc Natl Acad Sci U S A. 1991;88:5106–10.
Bramham CR, Messaoudi E. BDNF function in adult synaptic plasticity: The synaptic consolidation hypothesis. Prog Neurobiol. 2005;76:99–125.
Sui L, Wang Y, Ju L-H, Chen M. Epigenetic regulation of reelin and brain-derived neurotrophic factor genes in long-term potentiation in rat medial prefrontal cortex. Neurobiology of Learning and Memory. Elsevier Inc. 2012;97:425–40.
Aid T, Kazantseva A, Piirsoo M, Palm K, Timmusk T. Mouse and rat BDNF gene structure and expression revisited. J Neurosci Res. 2007;85:525–35.
Yang X, Han H, De Carvalho DD, Lay FD, Jones PA, Liang G. Gene Body Methylation Can Alter Gene Expression and Is a Therapeutic Target in Cancer. Cancer Cell. 2014;26:577–90.
Heyn H, Vidal E, Ferreira HJ, Vizoso M, Sayols S, Gomez A, et al. Epigenomic analysis detects aberrant super-enhancer DNA methylation in human cancer. Genome Biol. 2016;17. https://genomebiology.biomedcentral.com/about.
Maurano MT, Wang H, John S, Shafer A, Canfield T, Lee K, et al. Role of DNA Methylation in Modulating Transcription Factor Occupancy. Cell Reports [Internet]. 2015;12:1184–95.
Medvedeva YA, Khamis AM, Kulakovskiy IV, Ba-Alawi W, Bhuyan MSI, Kawaji H, et al. Effects of cytosine methylation on transcription factor binding sites. BMC Genomics BioMed Central. 2014;15:1.
Baxter PS, Bell K, Hasel P, Kaindl AM. Synaptic NMDA receptor activity is coupled to the transcriptional control of the glutathione system. Nature. 2015;2015;6:6761.
Xu B, Qu X, Gu S, Doughman Y-Q, Watanabe M, Dunwoodie SL, et al. Cited2 is required for fetal lung maturation. Dev Biol. 2008;317:95–105.
Sun W, Choi S-H, Park SK, Kim SJ, Noh MR, Kim EH, et al. Identification and characterization of novel activity-dependent transcription factors in rat cortical neurons. J Neurochem. 2007;100:269–78.
Landgren H, Carlsson P. FoxJ3, a novel mammalian forkhead gene expressed in neuroectoderm, neural crest, and myotome. Dev Dyn. 2004;231:396–401.
Yu WM, Appler JM, Kim YH, Nishitani AM, Holt JR, Goodrich LV. A Gata3–Mafb transcriptional network directs post-synaptic differentiation in synapses specialized for hearing. Elife. 2013;2:e01341.
Liu J, Magri L, Zhang F, Marsh N, Albrecht S, Huynh J, et al. Chromatin landscape defined by repressive histone methylation during oligodendrocyte differentiation. J Neurochem. 2015;134:96–6.
Chen Y, Bang S, McMullen MF, Kazi H, Talbot K, Ho M-X, et al. Neuronal Activity-Induced Sterol Regulatory Element Binding Protein-1 (SREBP1) is Disrupted in Dysbindin-Null Mice-Potential Link to Cognitive Impairment in Schizophrenia. Mol Neurobiol. 2016. [Epub ahead of print].
Gehlenborg N. UpSetR: A More Scalable Alternative to Venn and Euler Diagrams for Visualizing Intersecting Sets. 2016.
Harrell Jr FE, Dupont WCFC, et al. M. Hmisc: Harrell Miscellaneous. 2015.
Stojnic R, Diez D. PWM enrichment analysis. 2015.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Meth. 2012;9:357–9.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25:2078–9.
Shen L, Shao N, Liu X, Nestler E. ngs.plot: Quick mining and visualization of next-generation sequencing data by integrating genomic databases. BMC Genomics. 2014;15:284.
Skowronska-Krawczyk D, Ma Q, Schwartz M, Scully K, Li W, Liu Z, et al. Required enhancer-matrin-3 network interactions for a homeodomain transcription program. Nature. 2014;514:257–61.
Hunter RG, Murakami G, Dewell S, Seligsohn M, Baker MER, Datson NA, et al. Acute stress and hippocampal histone H3 lysine 9 trimethylation, a retrotransposon silencing response. Proc Natl Acad Sci U S A. 2012;109:17657–62.
Saha RN, Wissink EM, Bailey ER, Zhao M, Fargo DC, Hwang J-Y, et al. Rapid activity-induced transcription of Arc and other IEGs relies on poised RNA polymerase II. Nat Neurosci Nature Publishing Group. 2011;14:848–56.
Szulwach KE, Li X, Li Y, Song C-X, Wu H, Dai Q, et al. 5-hmC–mediated epigenetic dynamics during postnatal neurodevelopment and aging. Nat Neurosci. 2011;14:1607–16.
Cadet JL, Jayanthi S, McCoy MT, Ladenheim B, Saint-Preux F, Lehrmann E, et al. Genome-wide profiling identifies a subset of methamphetamine (METH)-induced genes associated with METH-induced increased H4K5Ac binding in the rat striatum. BMC Genomics BioMed Central. 2013;14:1.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323.
Research Council of Norway (grants 204861 and 199355), Bergen Medical Research Foundation, and the University of Bergen.
Availability of data and materials
The DNA-methylation data generated during the current study are available in the ArrayExpress repository under the accession number E-MTAB-4919.
JLVM, MED, and KW designed the experiments. JLVM designed the array, and conducted the bioinformatic analysis. DCK prepared the samples for the arrays. DP conducted the HFS. TP helped with the statistical analysis. KW extracted the DNA. CRB, KW and MED provided means for the study, and along with JLVM wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
The rat experiments were approved by the Norwegian Committee for Animal Research and were carried out in accordance with the European Community Council Directive of 24 November 1986 (86/609/EEC). The outbred Sprague–Dawley rats used were commercially obtained from Taconic Europe (Eiby, Denmark; model number SD-M).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Region specific log fold changes after stimulation. Violin plots showing the differentially methylated (a) probes and (b) regions for each time point compared to control. The number over each violin represents the mean logFC per probe or region. (PDF 1006 kb)
Differentially methylated regions for each time point. The differentially methylated regions discovered at 30 min, 2 h, and 5 h post-HFS. (XLSX 233 kb)
Expression changes of methylation and de-methylation genes. Log2 fold changes in expression for Dnmt1, Dnmt3, Tet1, and Tet3 between stimulated and non-stimulated dentate gyrus. Data adapted from Maag et al. 2015. Control group is the naïve unstimulated samples from the RNA-seq study. Number on each bar shows the Benjamini Hochberg adjusted FDR values. (PDF 125 kb)
Correlation between expression and methylation changes after HFS. Gene log fold changes of DE genes were imported from Maag et al. 2015. For each time point, all DE genes logFC were plotted against the logFC methylation after HFS vs. control. Each gene has multiple corresponding methylation regions based on figure 1b. Lines were plotted at ± 0.9 logFC expression and ± 0.1 logFC methylation. Each resulting quadrant shows the number of genes present. (PDF 118 kb)
Examples of genes showing linear negative correlation between normalised expression and methylation. Genes showing linear correlation between methylation and expression through all time points from left to right, (a) Egr2, Egr3, Egr4. (b) Etv5, Zswim4, Ube2g1. (c) Stam, Bhlhe40, Jak2. (PDF 575 kb)
Motif search in promoter region after removing all non-CpG-island overlapping promoters. (a) Motifs enriched in significantly methylation/expression correlated 1 kb upstream of the TSS compared to non-correlated genes in genes with CpG-islands present in the promoter region. (b) Number of CpG-islands present in 1 kb upstream from the promoter in both groups compared to number of promoters in each groups. (c) Boxplots to the right show the GC frequency of the significantly and non-significantly correlated genes. μ represents the mean for each group. P-value was calculated using the students’ t-test. (PDF 2830 kb)
Investigation of GC-content in differentially methylated probes and regions reveals no bias between differential methylation and GC-content. GC-content of each (a) probe and each (b) region compared to the t-statistics from limma when comparing to control. (c) GC-content of each region compared to the correlation value when compared to gene expression. (PDF 2187 kb)
Differences in expression levels between methylation/expression correlated and non-correlated genes. RPKM for correlated and non-correlated genes for each time-point. The lighter coloured boxplots show the non-correlated group while the darker coloured boxplots shows the methylation/expression correlated group. Statistics was calculated using student’s t-test. (PDF 134 kb)
Group post-synaptic potential changes in response to HFS, and average array intensity and sample clustering. (a) The differences in the groups’ field excitatory post-synaptic potential after HFS. The data was taken from Maag et al. 2015, and included all 3 samples per group. The control group in this graph is the baseline test stimuli used as control in the present study (b) Distribution of probe intensity per samples before quantilie normalisation. MDS plot of (c) probes and (d) regions after quantile normalisation. (PDF 1131 kb)