- Open Access
Epigenetic control of epilepsy target genes contributes to a cellular memory of epileptogenesis in cultured rat hippocampal neurons
Acta Neuropathologica Communications volume 5, Article number: 79 (2017)
Hypersynchronous neuronal excitation manifests clinically as seizure (ictogenesis), and may recur spontaneously and repetitively after a variable latency period (epileptogenesis). Despite tremendous research efforts to describe molecular pathways and signatures of epileptogenesis, molecular pathomechanisms leading to chronic epilepsy remain to be clarified. We hypothesized that epigenetic modifications may form the basis for a cellular memory of epileptogenesis, and used a primary neuronal cell culture model of the rat hippocampus to study the translation of massive neuronal excitation into persisting changes of epigenetic signatures and pro-epileptogenic target gene expression. Increased spontaneous activation of cultured neurons was detected 3 and 7 days after stimulation with 10 μM glutamate when compared to sham-treated time-matched controls using calcium-imaging in vitro. Chromatin-immunoprecipitation experiments revealed short-term (3 h, 7 h, and 24 h) and long-term (3 d and 2 weeks) changes in histone modifications, which were directly linked to decreased expression of two selected epilepsy target genes, e.g. excitatory glutamate receptor genes Gria2 and Grin2a. Increased promoter methylation observed 4 weeks after glutamate stimulation at respective genes suggested long-term repression of Gria2 and Grin2a genes. Inhibition of glutamatergic activation or blocking the propagation of action potentials in cultured neurons rescued altered gene expression and regulatory epigenetic modifications. Our data support the concept of a cellular memory of epileptogenesis and persisting epigenetic modifications of epilepsy target genes, which are able to turn normal into pro-epileptic neurons and circuits.
Temporal lobe epilepsy (TLE) is the most frequent focal epilepsy in humans and often associated with an initial precipitating injury, such as brain trauma, inflammation or prolonged febrile seizures, followed by a clinically silent latency period before onset of chronic recurrent seizures [1, 2]. Despite tremendous research efforts to describe molecular pathways and signatures of epileptogenesis, pathophysiological mechanisms leading to chronic epilepsy remain to be clarified. Gene expression profiling studies identified many genes to be differentially expressed in chronic TLE. Comprehensive gene pathway analysis suggested inflammation, stress, synaptic transmission and plasticity to play a role [3,4,5,6,7]. More recent data suggested epigenetic mechanisms, including DNA methylation, histone modifications, chromatin remodeling and non-coding RNAs as key upstream mechanisms deregulating gene expression, thereby promoting epileptogenesis [8,9,10,11,12,13,14]. These studies identified locus-specific as well as genome-wide alterations in DNA methylation patterns in different animal models or human surgical brain samples, compatible with a compromised gene regulation machinery underlying epileptogenesis. Furthermore, hyper- and hypoacetylation of histone H4 have been reported for several gene promoters following an initial precipitating injury, e.g. status epilepticus, and seem to be independent of the analyzed model [15,16,17,18]. Histone modifications are, therefore, of particular interest and capable to promote the epileptogenic process. Here we asked if neuronal hyperexcitation alters the epigenetic machinery of hippocampal neurons towards the previously described pro-epileptogenic cellular signature.
We intended to induce rhythmic hyperexcitation in cultured hippocampal neurons with 10 μM glutamate, to be documented by live-cell calcium imaging [19, 20]. Chromatin immunoprecipitation was performed at various time intervals to study epigenetic histone modifications, i.e. H4 acetylation as well as H3K4, H3K9 and H3K27 trimethylation. Well-characterized epilepsy genes were then investigated as potential targets of a pro-epileptogenic cellular signature in our simplistic cell culture model. Blockage of glutamatergic signaling by D,L-AP5 and NBQX and of the propagation of action potentials by TTX was performed to provide evidence for the principal role of neuronal excitation as trigger of the epigenetic machinery. This experimental strategy was designed to help answering the question if synchronized neuronal hyperexcitation is capable of inducing long-lasting epigenetic signatures and facilitating a cellular memory of epileptogenesis (CME).
Materials and methods
Animals and tissue preparation
Adult Wistar rats were obtained from Charles River (Sulzfeld, Germany), bred and maintained at the local animal center in breeding cages under controlled environmental conditions (12 h light/dark cycle, 20–23 °C, 50% relative humidity, drinking and feeding ad libitum). Newborn or up to two-day-old male and female offspring were used for the in vitro model. All animal experiments have been approved by the local animal care and use committee (TS-1/13) and were in accordance with the European Communities Council Directive and German Animal Welfare Act (54–2532.1-23/09, Directive 2010/63/EU).
Preparation of cell suspensions and dispersed hippocampal cell culture
Cell suspensions were prepared from rat hippocampus as described previously . Coverslips (for immunofluorescence staining and calcium imaging) and culture dishes (for any other application) were coated with poly-D-lysine hydrobromide (Sigma-Aldrich, Taufkirchen, Germany). Cell suspension was preplated onto an uncoated flask and incubated at 37 °C in 5% CO2 for 1 h. During this time glial cells settled down and adhered to the bottom of the flask, while neurons remained in the supernatant. Supernatant was collected thereafter and centrifuged at 800 rpm for 8 min at room temperature. Cell pellets were resuspended and cultured in serum-free Neurobasal-A medium supplemented with 2% B27, 0.5 mM GlutaMAX and 1% penicillin-streptomycin (all Life Technologies, Darmstadt, Germany). Cells were plated on poly-D-lysine coated dishes or coverslips at a density of 2.5 × 105 onto 3.5 cm2. Cells were maintained at 37 °C in a fully humidified incubator containing 5% CO2. After 24 h Cytosine β-D-arabinofuranoside hydrochloride (AraC; Sigma-Aldrich) was added to inhibit proliferation of remaining glial cells. Neurons were maintained in dispersed culture with the original media up to 40 days in vitro (DIV).
Glutamate treatment of cell cultures was performed as described elsewhere [22, 23]. At 12 DIV, culture media was replaced by a physiological treatment solution (145 mM NaCl, 2.5 mM KCl, 10 mM HEPES [pH 7.4], including 10 mM glucose, 2 mM CaCl2, 1 mm MgCl2, and 2 μM glycine for control cultures, adding 10 μM glutamate for stimulation of the glutamate group). Neuronal cultures were exposed to treatment solution for 10 min, washed with treatment solution 3 times and replaced by original culture media again until the end of the experiment. 1 μM TTX (Sigma-Aldrich, Taufkirchen, Germany) was added during glutamate treatment to inhibit action potential discharges through interference with voltage-gated sodium channels. 10 μM 2,3-dihydroxy-6-nitro-7-sulfamoyl-benzo-quinoxaline-2,3-dione (NBQX, Signal-Aldrich) and 50 μM D-amino-5-phosphonovaleric acid (D,L-AP5, Sigma-Aldrich) were added to block excitatory synapses (i.e. AMPA/kainate and NMDA receptors) and prevent recurrent glutamatergic activity. After glutamate treatment neurons were kept in culture up to 4 weeks for downstream applications.
Cell viability assay
Neuronal viability was determined 2 and 4 weeks after glutamate treatment with propidium iodide (PI) and fluorescein diacetate (FDA) (both Thermo Scientific, Dreieich, Germany). FDA diffuses through the membrane of viable cells where it is metabolized into the green fluorescent fluorescein. PI cannot pass the viable cell membrane but intercalates with DNA in necrotic cells. Neurons were incubated with 500 nM PI and 160 nM FDA for 5 min in Neurobasal A medium and washed 2 times with Neurobasal A medium. Semi-quantitative measurements of PI/FDA stained cells were performed in triplicate using a microcomputer imaging system (ColorView II CCD camera, Cell^F imaging software, Olympus, Tokyo, Japan) equipped to an Olympus XI70 microscope. PI and FDA positive cell bodies were tagged on the computer screen and manually counted in four independent visual fields at 20× objective magnification. Cell viability was determined three independent times in glutamate-treated cultures and time matched sham controls.
Immunofluorescence staining was performed 2 and 4 weeks after glutamate treatment. Neurons cultured on coverslips were fixed in 4% formaldehyde. Fixed cells were pre-treated with blocking solution containing tris-buffered saline (TBS; 0.05 M Tris-Cl pH 7.6) with 1% bovine serum albumin (BSA, Biochrom, Berlin, Germany), 2% fish skin gelatin (Sigma-Aldrich) and Triton X-100 (0.1%, Sigma-Aldrich, Steinheim, Germany) for 2 h. Primary antibodies were diluted in blocking solution and incubated overnight at 4 °C. GFAP (mouse monoclonal anti-GFAP, clone GA5, mAb3670, Cell Signaling, Danvers, USA) was used as glial marker at a dilution of 1:300 and MAP2C (chicken polyclonal anti-MAP2, ab5392, abcam, Cambridge, UK) as neuronal marker at a dilution of 1:2000. Secondary fluorophore-conjugated antibodies were diluted in blocking solution and incubated for 4 h at room temperature. Secondary antibodies were labelled with FITC (1:200, rabbit anti-chicken, ab6749, abcam) and Cy3 (1:100, goat anti-mouse, Dianova). Nuclei were counterstained with DAPI (Sigma-Aldrich) at a dilution of 1:1000 for 15 min. Semi-quantitative measurements of neuronal and glial cell numbers were performed in triplicate using a microcomputer imaging system (ColorView II CCD camera, Cell^F imaging software,) equipped to an Olympus XI70 inverted fluorescence microscope (Olympus, Tokyo, Japan). Immunohistochemically stained neuronal and glial cell bodies were tagged on the computer screen and manually counted at 20× objective magnification in four independent visual fields.
Calcium imaging experiments were performed in triplicate during acute glutamate stimulation and two time points after glutamate treatment (3 d and 7 d). Coverslips were incubated with 5 μM Fluo-4 AM (Thermo Fisher) in recording solution (144 mM NaCl, 2.5 mM KCl, 2.5 mM CaCl2, 2.5 mM MgCl2, 10 mM Glucose, 10 mM HEPES, pH 7.4, 320 mOsm) for 30 min at 37 °C. After coverslips were washed twice with PBS, they were placed into a perfusion chamber filled with 500 μl recording solution. A Nikon TI-Eclipse inverted fluorescence microscope (Nikon, Tokio, Japan) with 10×, 0.45 NA objective and Perfect Focus System was used for all imaging experiments. Fluo-4 was excited by a Nikon Intensilight C-HGFI lamp (Neutral Density Filter 16) through excitation filters of 455–485 nm and a 495 nm dichroic long-pass mirror. The emitted light was passed through an emission bandpass filter of 500–545 nm (Semrock, Rochester, NY) and was projected onto a cooled EM-CCD camera (iXonEM DU-885, Andor). Imaging was controlled and recorded by NIS Elements (Nikon).
Images were recorded at 100 ms intervals and at an exposure time of 70 ms. Image stacks were converted into tagged image file format. All image analysis was performed using custom-written routines in MATLAB (The MathWorks, Natick, MA). Regions of interest (ROIs) were automatically detected using the background-determined feature point detection as described in Sbalzarini et al.  (intensity parameter w = 0.1–0.8%, approximate cell radius r = 5.0 pixel). For each ROI the relative fluorescence trace was calculated . Spike estimation on these traces was performed defining a spike as an event of more than 3–6 standard deviations in amplitude within 30 time frames.
Gene expression analysis
RNA was extracted from primary neuronal cultures at 5 different time points after glutamate treatment (3 h, 7 h, 24 h, 3 d, and 14 d) using TRIzol Reagent (Life technologies; according to the manufacturer’s instructions) followed by DNase treatment (Life technologies). First-strand cDNA synthesis was performed using the SuperScript II Reverse Transcriptase Kit (Life technologies; according to the manufacturer’s instructions). Quantitative real-time PCR was performed using the 7500 Fast Real-Time PCR System (Life technologies) with Power SYBR Green PCR Master Mix (Life technologies) as fluorescent dye according to the manufacturer’s protocol. Primers used for amplification of cDNA are specified in Table 1. GAPDH quantification was used as internal control for normalization. Fold difference of mRNA levels were calculated using the ΔΔCt method. No-template controls for each primer were included on every plate and melt curve analysis was performed to exclude unspecific amplification. All PCR reactions were performed in triplicate and repeated three independent times.
Chromatin immunoprecipitation assay
Chromatin immunoprecipitation (ChIP) was performed according to a modified X-ChIP protocol by abcam. Cultured neurons were cross-linked in 0.8% formaldehyde for 10 min at room temperature at 5 different time points after glutamate treatment (3 h, 7 h, 24 h, 3 d, 14 d). Crosslinking reaction was stopped by adding glycine to a final concentration of 0.125 M and centrifuged at 1000 rpm for 8 min at 4 °C. Pellet was washed twice in cold PBS and snap frozen in liquid nitrogen. Cells were lysed in SDS-lysis buffer (50 mM Tris-HCl [pH 8.0], 10 mM EDTA dihydrate, 1% SDS) for 1 h. Lysate was sheared using a Diagenode Bioruptor on high power for 4 × 5 cycles (30 s on-off). Equal amounts of chromatin lysate for glutamate and control group (1500–3000 ng) were diluted ten times with RIPA Buffer (50 mM Tris-HCl [pH 8:0], 150 mM NaCl, 2 mM EDTA dihydrate [pH 8.0], 1% NP-40, 0.5% Sodium deoxycholate, 0.1% SDS) and pre-cleared with Dynabeads Protein G magnetic beads (Life technologies) for 1 h. Immunoprecipitation was performed overnight with the corresponding primary IgG antibody along with the Dynabeads Protein G magnetic beads. Immune complexes were washed three times with ChIP wash buffer (0.1% SDS, 1% Triton X-100, 2 mM EDTA dihydrate [pH 8.0], 150 mM NaCl, 20 mM Tris-HCl [pH 8.0]), once with ChIP final wash buffer (0.1% SDS, 1% Triton X-100, 2 mM EDTA dihydrate [pH 8.0], 500 mM NaCl, 20 mM Tris-HCl [pH 8.0]) and were then eluted from the beads with 120 μl elution buffer (1% SDS, 0,1 M NaHCO3, pH 9.0). Protein-DNA cross-links were reverted by overnight incubation at 65 °C and DNA was extracted using the QIAquick PCR Purification Kit (Qiagen) according to the manufacturer’s protocol. Immunoprecipitated DNA was preamplified using the SsoAdvanced PreAmp Supermix (Biorad, Ismaning, Germany) according to the manufacturer’s instructions. DNA quantification was done by quantitative real-time PCR using the same protocol as mentioned above. The following antibodies were used: 2 μg H3K4me3 (07–473, Millipore, Molsheim, France), 2 μg H4ac (06–866, Millipore), 4 μg H3K27me3 (07–449, Millipore), 6 μg H3K9me3 (ab8898, abcam) and IgG (Millipore). Primers used for preamplification and qPCR of immunoprecipitated DNA are shown in Table 2. Hist1H4B was used as control primer for normalization for H4ac, H3K4me3, MyoD for H3K27me3, and Zfp12 for H3K9me3. All steps were performed at 4 °C up to and including the immunoprecipitation step. PBS, SDS-lysis buffer and RIPA buffer were supplemented with protease inhibitors (EDTA-free protease inhibitor cocktail tablets; Roche) and deacetylase inhibitors (Active Motif, La Hulpe, Belgium) immediately before application. Non-specific IgG antibody was included as a negative control to exclude unspecific interaction. All ChIP experiments were repeated three independent times for each antibody.
Genomic DNA was extracted from primary neuronal culture 4 weeks after glutamate treatment using the QIAamp DNA Micro Kit (Qiagen, Hilden, Germany), followed by bisulfite conversion of 1 μg of genomic DNA using the EpiTect Bisulfite Kit (Qiagen) according to the manufacturer’s protocol. Regions of interest were amplified using the TaKaRa EpiTaq HS Kit (TaKaRa Clontech, Otsu, Japan) and cloned using the TOPO TA Cloning Kit (Life technologies) doing blue/white screening. White colonies were selected and grown in LB medium overnight. The plasmid was purified using the GeneJET Plasmid Miniprep Kit (Thermo Scientific, Darmstadt, Germany) and clones were sequenced by Sanger sequencing (GATC Biotech, Cologne, Germany). Sequences were quality controlled and aligned using the CLC sequence viewer (CLC bio) and Quantitation tool for Methylation Analysis software (RIKEN). Primers used to amplify bisulfite converted DNA are summarized in Table 3.
Statistical analysis was performed using GraphPad Prism6 (GraphPad Software, San Diego, CA). Graphs were generated using MATLAB (The MathWorks Inc., Natick) or GraphPad. Error bars and specified values represent standard deviation (SD). Depending on number of groups analyzed, sample size and respective distribution, statistical comparisons were made with Fisher’s exact test (site-specific DNA methylation), one-way ANOVA and Dunnett’s correction for multiple comparisons (gene expression, calcium imaging and ChIP data), or Mann-Whitney test (IF/cell counts, DNA methylation), respectively. In all figures, asterisks indicate a level of significance of p < 0.05.
Transient exposure to glutamate has been shown to induce spontaneous ictal-like activity in hippocampal neuronal cultures. In the present study we used this simplistic neuronal cell culture model to analyze how massive excitation of neuronal membranes modifies the epigenetic machinery and translates into persisting gene expression changes. Rat hippocampal neurons were cultured 12 DIV and then exposed to 10 μM glutamate for 10 min. Cultures were kept up to 4 weeks under standard conditions for further use of specific downstream applications (see below and summary in Fig. 1).
Viability of dissociated rat hippocampal neurons
Using PI/FDA stainings for assaying cell viability, 77% (+/−1%) of cells survived 4 weeks after transient 10 μM glutamate exposure, compared to 78% (+/− 2%) in time-matched sham-treated controls (Mann-Whitney test, p > 0.05; Fig. 2a and a’). These cells established a complex neuronal network. Double immunofluorescence with the neuronal marker MAP2C and glial marker GFAP identified 89% (+/− 1%) of neurons in the cultured cell population compared to 92% (+/− 1%) in time-matched sham-treated controls 4 weeks after glutamatergic stimulation (Mann-Whitney test, p > 0.05; Fig. 2b and b’).
Glutamate induced alterations in neuronal excitability, development of epileptogenic networks
Next we asked if glutamatergic excitation changed the intrinsic network activity of cultured hippocampal neurons. Using live cell calcium (Ca2+) imaging with Fluo4-AM dye, we monitored glutamate-evoked recurrent spontaneous Ca2+ transients as surrogate marker for spiking activity in neuronal networks 3d and 7d after glutamate stimulation. Twenty minute traces of spontaneous Ca2+ activity, were simultaneously captured from neuronal populations of up to 200 cells at single cell resolution. Isolated episodes of low frequency and low amplitude Ca2+ oscillations could be observed in dissociated hippocampal neurons. Following glutamate exposure a significant increase in intracellular Ca2+ uptake was identified (Fig. 3a). This could be reverted to baseline, but not completely blocked, by either NBQX/AP-5 or TTX application (Fig. 3a and b). Over time we detected a significant increase in amplitude and frequency of spikes (from 3d on), synchronization of neuronal firing over large networks of neurons as well as a shift from single spikes towards bursts of spikes following glutamate stimulation (from 7d on) compared to sham-treated time-matched controls (Fig. 3c). Cells that were co-treated with either NBQX/AP5 or TTX and glutamate showed little synchronization and no burst activity similar to sham controls (data not shown). Taken together, our data is compatible with a simplistic pathogenic model of epileptogenesis.
Glutamate receptor subunit expression in epileptogenesis
Altered expression or functional mutations affecting glutamate receptor subunit composition have been previously linked with epilepsy. We studied epigenetic gene expression of ionotropic α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptor subunit Gria2 and N-methyl-D-aspartate (NMDA) receptor subunit Grin2a at 5 different time points after glutamate stimulation (i.e. 3 h, 7 h, 24 h, 3 d, and 2 w). Using qPCR analyses, we identified a consistent reduction in Gria2 and Grin2a gene expression at all experimental time points. Gria2 mRNA levels were significantly decreased after 3 h and declined further until 2 weeks after glutamate treatment (Fig. 4a; one-way ANOVA and Dunnett’s post-hoc test p < 0.0001). Likewise, Grin2a mRNA levels decreased consecutively from 3 h to 2 weeks after glutamate stimulation (Fig. 5a; p < 0.0001). These indings support previous data describing similar changes in subunit composition of glutamate receptors in experimental and human epilepsy (Table 4).
Epigenetic control of glutamatergic mechanisms in epileptogenesis
We then asked whether altered Gria2 and Grin2a expression can be linked to suppressive epigenetic histone modifications [15, 16, 26]. We observed a significant and permanent decrease in H4ac levels within the Gria2 promoter up to 2 weeks post injury, both up- and downstream the transcriptional start site (TSS) including the 5′ untranslated region (5’UTR) and part of exon 1 of the Gria2 gene (Fig. 4b and c, left panel, locus Gria2 A, one-way ANOVA and Dunnett’s post-hoc test p < 0.0001; locus Gria2 C, p < 0.01). NMDA receptor subunit Grin2a also showed permanent decrease in promoter associated acetylation levels of histone H4 at a region corresponding to the most upstream of three alternative TSSs (Fig. 5b and c, left panel, locus Grin2a A, p < 0.0001).
No changes in H3K4 trimethylation (H3K4me3; p > 0.05), another gene activating histone modification, were found at any given locus of the Gria2 and Grin2a promoter following glutamate injury compared to sham controls (Figs. 4c and 5c, second panel on left, p > 0.05).
In contrast, repressive H3K9 trimethylation (H3K9me3) increased in a time-dependent manner peaking at 3 h and 24 h after glutamate treatment providing an inhibitory signal at the promoter structures of Gria2 and Grin2a (Figs. 4c and 5c, right panel, p < 0.05). H3K27me3 also increased significantly at Gria2 and Grin2a promoters (Fig. 4c and 5c, second panel on right, p < 0.01), thereby adding, with a slight delay over H3K9me3, another repressive signal at these gene structures. Both signals were transient and returned to control levels at 3 d after glutamate exposure. As long-term gene silencing may be mediated through DNA methylation we performed bisulfite cloning and sequencing of the Gria2 and Grin2a promoters including TSSs, 5′ UTRs and parts of exon 1 at 4 weeks following glutamate stimulation and in time-matched sham-treated controls. We identified increased DNA methylation at both CpG and non-CpG sites within promoter regions of both Gria2 (Fig. 4d) and Grin2a (Fig. 5d).
Taken together our data indicated that glutamate induced massive neuronal excitation induced changes in histone modifications and DNA methylation at the promoter regions of ionotropic glutamate receptor subunits compatible with a collaborative epigenetic regulation of gene expression.
Inhibition of excitatory glutamatergic signaling rescues epigenetic modifications of epilepsy candidate genes
To test whether epigenetic alterations were induced by neuronal hyperactivity in our cell culture we performed inhibition studies either interfering with excitatory synapses or propagation of action potentials (Fig. 6). Neuronal cultures were simultaneously treated with glutamate and NBQX/AP5 or TTX, respectively. Calcium imaging performed 7 d after treatment showed little synchronization, no burst activity (data not shown) and significantly decreased spike frequency and amplitude, similar to sham controls (Fig. 6a). Cell viability assays and quantification of neuronal and glial cells in culture identified cell survival rates similar to time-matched sham-treated control cultures two weeks following treatment (data not shown). TTX and NBQX/AP5 prevented glutamate induced gene repression of Gria2 and Grin2a (Fig. 6b). As an indicator of epigenetic changes H4ac levels were determined 3 h and 3 d after treatment with TTX or NBQX/AP5. Glutamate induced decrease of H4ac at Gria2 and Grin2a could be abolished by TTX and by NBQX/AP5 (Fig. 6c, with no difference to sham-treated controls. Taken together, the results demonstrated a rescue of aberrant glutamate receptor subunit gene expression and regulatory epigenetic marks after treatment with TTX or NBQX/AP5.
Understanding the sequence and timing of molecular epileptogenesis is mandatory for targeted therapies and disease-modifying intervention. The hippocampus is key to ictogenesis in TLE and likely encodes structural, molecular and functional signaling pathways promoting chronic epilepsy. We asked whether and how massive hyperexcitation of neurons translates into aberrant and persisting pro-epileptogenic gene expression and upstream epigenetic modifications, thereby contributing to a cellular memory of epileptogenesis (CME). The divergent cellular composition of CNS tissue with glia, neurons and mesenchyme makes it challenging to unravel such complex issue . A “simplistic” neuronal cell culture model of rat hippocampal neurons was specifically used to address this question.
Following transient glutamatergic stimulation  we studied cellular and molecular changes in principal hippocampal neurons at different time points up to 4 weeks after glutamatergic stimulation. We recorded rhythmic neuronal activation at 7 days after glutamatergic stimulation, and identified complex epigenetic alterations leading to decreased expression of excitatory glutamate receptor genes Gria2 and Grin2a. Inhibition of ionotropic glutamatergic signaling and propagation of action potentials with NBQX/AP5 and TTX, respectively, during glutamate stimulation rescued aberrant gene expression and epigenetic modifications in cultured neurons. Moreover, the time-dependent development of epileptiform neuronal activation was blocked by this treatment. Gria2 and Grin2a are well recognized candidate genes of epileptogenesis [28, 29]. Previous studies identified transcriptional regulation of Gria2 by epigenetic mechanisms (Table 4; ). Grin2a expression was previously linked to HDAC2 activity and H4K12acetylation in animal models of Alzheimer’s disease . Moreover, aberrant DNA methylation at the GRIN2A locus was described in patients with major depression . Our model revealed fast decrease in both Gria2 and Grin2a gene expression following glutamate induced neuronal hyperactivity. Downregulation of glutamate receptor subunits was initiated within 3 h following glutamate exposure and remained stable thereafter. No downregulation of Gria2 and Grin2a was observed upon inhibition of glutamatergic excitation or propagation of action potentials with NBQX/AP5 or TTX, respectively. Our data were in line with previous findings showing that glutamate receptor subunit composition can be adjusted to neuronal activity within minutes or hours . Long-term downregulation of Gria2 and Grin2a are suggested to contribute to lasting changes in AMPA and NMDA receptor properties and downstream pro-epileptogenic events including neuronal death and functional network reorganization [33,34,35]. Changes in synaptic strength all over the network may thus have formed the basis for the development of rhythmic neuronal network hyperactivity as shown in our model and monitored by calcium imaging. Future studies are required to gain further insight into signaling over and synchronization of the neuronal network in the present model. Temporal resolution of our imaging experiment was too low to reliably dissect spike timing and temporal recruitment of cells into the active neuronal network. Moreover extended imaging of glutamate-treated cultures at earlier (e.g. 24 h and 48 h) and later time-points (e.g. 14 d and 4 weeks) will be needed to understand the physiological properties of the cells in relation to observed molecular changes.
In our particular interest to understand regulatory events of altered gene expression in epileptogenesis, we observed long-lasting decrease in histone 4 acethylation (H4ac) at Gria2 and Grin2a promoters. Histone acetylation promotes a chromatin structure permissive to gene expression . There are 26 sites of acetylation on a nucleosome and histone acetylation is dynamically regulated by histone acetyltransferases (HATs) and the antagonistic effects of histone deacetylases (HDACs) to adjust the level of transcription. Both Gria2 and Grin2a were previously shown to be sensitive to H4 deacetylation by HDACs .
Beside histone deacetylation we identified a rapid but transient increase in inhibitory H3K9me3 over the first 3 h following the glutamatergic excitation. Although H3K9me3 is mostly recognized for its role in stabilizing heterochromatin , recent studies provided evidence that H3K9me3 contributes also to short-term silencing of actively transcribed genes [38, 39]. There was also a delayed and transient increase in H3K27me3, another transcriptional repressing histone modification, occurring 7 h hours to one day after glutamate stimulation. Fast but transient changes in histone methylation indicate that they mediate transcriptional repression, but are not used for long-term gene silencing. Therefore, we analyzed DNA methylation patterns 4 weeks after glutamate stimulation. DNA methylation is commonly established at a much slower rate than histone modifications and rather contribute to stabilizing silencing marks than to initiate transcriptional repression . Machnes et al. previously showed that kainic acid exposure to mouse hippocampal slices as well as kainic acid induced status epilepticus in rats caused a mild increase in DNA methylation at the 5′ regulatory region of the Gria2 gene . Promoter methylation correlated with suppression of Gria2 mRNA expression. Also, in patients with major depressive disorder altered DNA methylation had been reported for three distinct CpG sites along the GRIN2A promoter . Here we identified increased CpG and non-CpG methylation in the Gria2 promoter after glutamate injury and seizure-like events in our model suggesting a role in long-term down-regulation of the gene.
Time dependent epigenetic patterns may correlate with disease progression
The present study allowed us to analyze long-term dynamics of epigenetic mechanisms from the time point of glutamate-induced massive neuronal hyperactivity through a latency period until the onset and presumed maintenance of recurrent hypersynchronous activity (Fig. 7). We detected both, transient gene specific and long-lasting epigenetic modifications. Fast changes detected in histone acetylation corresponded well to early changes in gene expression at both gene promoters. This is in line with other studies showing that histone acetylation turnover mediates fast responses to environmental signals. We determined fast and transient H3K9me3 at the promoter and first exon of Gria2 and Grin2a accompanying transcriptional repression, which suggested that H3K9me3 is involved in the initiation of transcriptional repression, but not maintenance. Alterations in H3K27me3 emerged later within 1 day after glutamate stimulation and were transient. These findings also confirm previous hypotheses that histone modifications are not stable enough to mediate long-term transcriptional changes. DNA methylation patterns are in general regarded less dynamic than histone modifications, but long-lasting. This could be confirmed in our cell culture model by increased CpG methylation at both the Grin2a and Gria2 promoters.
Our data linked the pathogenic model of gene repression in epileptogenesis with the regulatory epigenetic machinery of the cell. Persisting downregulation of both target genes was initiated through loss of H4 acetylation, accompanied by transient increase in H3K9me3, and subsequently stabilized by establishing H3K27me3. Finally, gene-specific repression is locked by DNA promoter methylation, leading to a cellular memory of epileptogenesis (Fig. 7). Our model may be further used to decipher the signaling pathway from an excited neuronal membrane to the cell nucleus, addressing key regulatory events of the epigenetic machinery for other epilepsy target genes of interest, asking for communalities as shared pathomechanism of epileptogenesis, and also looking into epigenetic mechanisms of previously described gene activation (e.g. RE1 Silencing Transcription Factor/Rest , neurotransmitter receptors mGluR1 and 4 , or Ca2+ channels ). Different mechanisms including Ca2+ signaling (Fig. 6e; [44, 45]) and metabolic regulation of transcription factors have been already proposed [46, 47] and may be the missing link between neuronal hyperactivity, epigenetic gene regulation and CME.
Blumcke I, Coras R, Miyata H, Ozkara C (2012) Defining clinico-neuropathological subtypes of mesial temporal lobe epilepsy with hippocampal sclerosis. Brain Pathol 22:402–411
Pitkanen A, Lukasiuk K (2011) Mechanisms of epileptogenesis and potential treatment targets. Lancet Neurol 10:173–186
Becker AJ, Chen J, Zien A, Sochivko D, Normann S, Schramm J, Elger CE, Wiestler OD, Blumcke I (2003) Correlated stage- and subfield-associated hippocampal gene expression patterns in experimental and human temporal lobe epilepsy. Eur J Neurosci 18:2792–2802
Elliott RC, Miles MF, Lowenstein DH (2003) Overlapping microarray profiles of dentate gyrus gene expression during development- and epilepsy-associated neurogenesis and axon outgrowth. J Neurosci 23:2218–2227
Gorter JA, van Vliet EA, Aronica E, Breit T, Rauwerda H, Lopes da Silva FH, Wadman WJ (2006) Potential new antiepileptogenic targets indicated by microarray analysis in a rat model for temporal lobe epilepsy. J Neurosci 26:11083–11110
Hendriksen H, Datson NA, Ghijsen WE, van Vliet EA, da Silva FH, Gorter JA, Vreugdenhil E (2001) Altered hippocampal gene expression prior to the onset of spontaneous seizures in the rat post-status epilepticus model. Eur J Neurosci 14:1475–1484
Lukasiuk K, Dabrowski M, Adach A, Pitkanen A (2006) Epileptogenesis-related genes revisited. Prog Brain Res 158:223–241
Belhedi N, Perroud N, Karege F, Vessaz M, Malafosse A, Salzmann A (2014) Increased CPA6 promoter methylation in focal epilepsy and in febrile seizures. Epilepsy Res 108:144–148
Debski KJ, Pitkanen A, Puhakka N, Bot AM, Khurana I, Harikrishnan KN, Ziemann M, Kaspi A, El-Osta A, al LK (2016) Etiology matters - genomic DNA methylation patterns in three rat models of acquired epilepsy. Sci Rep 6:25668
Kobow K, Jeske I, Hildebrandt M, Hauke J, Hahnen E, Buslei R, Buchfelder M, Weigel D, Stefan H, al KB (2009) Increased reelin promoter methylation is associated with granule cell dispersion in human temporal lobe epilepsy. J Neuropathol Exp Neurol 68:356–364
Martinowich K, Hattori D, Wu H, Fouse S, He F, Hu Y, Fan G, Sun YE (2003) DNA methylation-related chromatin remodeling in activity-dependent BDNF gene regulation. Science 302:890–893
Miller-Delaney SF, Bryan K, Das S, McKiernan RC, Bray IM, Reynolds JP, Gwinn R, Stallings RL, Henshall DC (2015) Differential DNA methylation profiles of coding and non-coding genes define hippocampal sclerosis in human temporal lobe epilepsy. Brain 138:616–631
Ryley Parrish R, Albertson AJ, Buckingham SC, Hablitz JJ, Mascia KL, Davis Haselden W, Lubin FD (2013) Status epilepticus triggers early and late alterations in brain-derived neurotrophic factor and NMDA glutamate receptor Grin2b DNA methylation levels in the hippocampus. Neuroscience 248C:602–619
Williams-Karnesky RL, Sandau US, Lusardi TA, Lytle NK, Farrell JM, Pritchard EM, Kaplan DL, Boison D (2013) Epigenetic changes induced by adenosine augmentation therapy prevent epileptogenesis. J Clin Invest 123:3552–3563
Huang Y, Doherty JJ, Dingledine R (2002) Altered histone acetylation at glutamate receptor 2 and brain-derived neurotrophic factor genes is an early event triggered by status epilepticus. J Neurosci 22:8422–8428
Jia YH, Zhu X, Li SY, Ni JH, Jia HT (2006) Kainate exposure suppresses activation of GluR2 subunit promoter in primary cultured cerebral cortical neurons through induction of RE1-silencing transcription factor. Neurosci Lett 403:103–108
Sng JC, Taniura H, Yoneda Y (2006) Histone modifications in kainate-induced status epilepticus. Eur J Neurosci 23:1269–1282
Tsankova NM, Kumar A, Nestler EJ (2004) Histone modifications at gene promoter regions in rat hippocampus after acute and chronic electroconvulsive seizures. J Neurosci 24:5603–5610
Cohen E, Ivenshitz M, Amor-Baroukh V, Greenberger V, Segal M (2008) Determinants of spontaneous activity in networks of cultured hippocampus. Brain Res 1235:21–30
Grienberger C, Konnerth A (2012) Imaging Calcium in Neurons. Neuron 73:862–885
Kiese K, Jablonski J, Boison D, Kobow K (2016) Dynamic regulation of the adenosine kinase gene during early postnatal brain development and maturation. Front Mol Neurosci 9:99
Srinivas KV, Jain R, Saurav S, Sikdar SK (2007) Small-world network topology of hippocampal neuronal network is lost, in an in vitro glutamate injury model of epilepsy. Eur J Neurosci 25:3276–3286
Sun DA, Sombati S, DeLorenzo RJ (2001) Glutamate injury-induced epileptogenesis in hippocampal neurons: an in vitro model of stroke-induced "epilepsy". Stroke 32:2344–2350
Sbalzarini IF, Koumoutsakos P (2005) Feature point tracking and trajectory analysis for video imaging in cell biology. J Struct Biol 151:182–195
Jia H, Rochefort NL, Chen X, Konnerth A (2011) In vivo two-photon imaging of sensory-evoked dendritic calcium signals in cortical neurons. Nat Protoc 6:28–35
Machnes ZM, Huang TC, Chang PK, Gill R, Reist N, Dezsi G, Ozturk E, Charron F, O'Brien TJ, al JNC (2013) DNA methylation mediates persistent epileptiform activity in vitro and in vivo. PLoS One 8:e76299
Cai X, Evrony GD, Lehmann HS, Elhosary PC, Mehta BK, Poduri A, Walsh CA (2014) Single-cell, genome-wide sequencing identifies clonal somatic copy-number variation in the human brain. Cell Rep 8:1280–1289
Dingledine R (2012) Glutamatergic mechanisms related to epilepsy: ionotropic receptors. In: Noebels JLAM, Rogawski MA, Olsen RW, Delgado-Escueta AV (eds) Jasper's basic mechanisms of the epilepsies 4th edn. Oxford university press, City, pp, pp 122–131
Lemke JR, Lal D, Reinthaler EM, Steiner I, Nothnagel M, Alber M, Geider K, Laube B, Schwake M, al FK (2013) Mutations in GRIN2A cause idiopathic focal epilepsy with rolandic spikes. Nat Genet 45:1067–1072
Graff J, Rei D, Guan JS, Wang WY, Seo J, Hennig KM, Nieland TJ, Fass DM, Kao PF, Kahn M et al (2012) An epigenetic blockade of cognitive functions in the neurodegenerating brain. Nature 483:222–226
Kaut O, Schmitt I, Hofmann A, Hoffmann P, Schlaepfer TE, Wullner U, Hurlemann R (2015) Aberrant NMDA receptor DNA methylation detected by epigenome-wide analysis of hippocampus and prefrontal cortex in major depression. Eur Arch Psychiatry Clin Neurosci 265:331–341
Bellone C, Nicoll RA (2007) Rapid bidirectional switching of synaptic NMDA receptors. Neuron 55:779–785
Casillas-Espinosa PM, Powell KL, O’Brien TJ (2012) Regulators of synaptic transmission: roles in the pathogenesis and treatment of epilepsy. Epilepsia 53:41–58
Liu SJ, Zukin RS (2007) Ca2+−permeable AMPA receptors in synaptic plasticity and neuronal death. Trends Neurosci 30:126–134
Tanaka H, Grooms SY, Bennett MV, Zukin RS (2000) The AMPAR subunit GluR2: still front and center-stage. Brain Res 886:190–207
Graff J, Tsai LH (2013) Histone acetylation: molecular mnemonics on the chromatin. Nat Rev Neurosci 14:97–111
Stewart MD, Li J, Wong J (2005) Relationship between histone H3 lysine 9 methylation, transcription repression, and heterochromatin protein 1 recruitment. Mol Cell Biol 25:2525–2538
Lu Y, Chu A, Turker MS, Glazer PM (2011) Hypoxia-induced epigenetic regulation and silencing of the BRCA1 promoter. Mol Cell Biol 31:3339–3350
Martin C, Zhang Y (2005) The diverse functions of histone lysine methylation. Nat Rev Mol Cell Biol 6:838–849
Bird A (2002) DNA methylation patterns and epigenetic memory. Genes Dev 16:6–21
Spencer EM, Chandler KE, Haddley K, Howard MR, Hughes D, Belyaev ND, Coulson JM, Stewart JP, Buckley NJ, al KA (2006) Regulation and role of REST and REST4 variants in modulation of gene expression in in vivo and in vitro in epilepsy models. Neurobiol Dis 24:41–52
Pitsch J, Schoch S, Gueler N, Flor PJ, van der Putten H, Becker AJ (2007) Functional role of mGluR1 and mGluR4 in pilocarpine-induced temporal lobe epilepsy. Neurobiol Dis 26:623–633
Becker AJ, Pitsch J, Sochivko D, Opitz T, Staniek M, Chen CC, Campbell KP, Schoch S, Yaari Y, Beck H (2008) Transcriptional upregulation of Cav3.2 mediates epileptogenesis in the pilocarpine model of epilepsy. J Neurosci 28:13341–13353
Hardingham GE, Bading H (2002) Coupling of extrasynaptic NMDA receptors to a CREB shut-off pathway is developmentally regulated. Biochim Biophys Acta 1600:148–153
Kokubo M, Nishio M, Ribar TJ, Anderson KA, West AE, Means AR (2009) BDNF-mediated cerebellar granule cell development is impaired in mice null for CaMKK2 or CaMKIV. J Neurosci 29:8901–8913
Ivanova D, Dirks A, Fejtova A (2016) Bassoon and piccolo regulate ubiquitination and link presynaptic molecular dynamics with activity-regulated gene expression. J Physiol 594:5441–5448
Ivanova D, Dirks A, Montenegro-Venegas C, Schone C, Altrock WD, Marini C, Frischknecht R, Schanze D, Zenker M, al GED (2015) Synaptic activity controls localization and function of CtBP1 via binding to bassoon and piccolo. EMBO J 34:1056–1077
Walczak A, Szczepankiewicz AA, Ruszczycki B, Magalska A, Zamlynska K, Dzwonek J, Wilczek E, Zybura-Broda K, Rylski M, al MM (2013) Novel higher-order epigenetic regulation of the Bdnf gene upon seizures. J Neurosci 33:2507–2511
Crosio C, Heitz E, Allis CD, Borrelli E, Sassone-Corsi P (2003) Chromatin remodeling and neuronal response: multiple signaling pathways induce specific histone H3 modifications and early gene expression in hippocampal neurons. J Cell Sci 116:4905–4914
Kobow K, Kaspi A, Harikrishnan KN, Kiese K, Ziemann M, Khurana I, Fritzsche I, Hauke J, Hahnen E, et al CR (2013) Deep sequencing reveals increased DNA methylation in chronic rat epilepsy. Acta Neuropathol 126:741–756
Miller-Delaney SF, Das S, Sano T, Jimenez-Mateos EM, Bryan K, Buckley PG, Stallings RL, Henshall DC (2012) Differential DNA methylation patterns define status epilepticus and epileptic tolerance. J Neurosci 32:1577–1588
We kindly thank B. Rings for her expert technical assistance.
Our work was supported by the German Research Council (DFG Bl 421/3–1) as part of the European Science Foundation EUROCORES Programme EuroEPINOMICS, and by the European Union’s Seventh Framework Program (DESIRE project, grant agreement #602531). JH received an MD thesis scholarship through the Interdisciplinary Center for Clinical Research (IZKF) of the Faculty of Medicine of the Friedrich-Alexander-University Erlangen-Nürnberg.
Availability of data and materials
The datasets used and analyzed in the current study are available from the corresponding author on reasonable request.
All animal experiments have been approved by the local animal care and use committee (TS-1/13) and were in accordance with the European Communities Council Directive and German Animal Welfare Act (54–2532.1-23/09, Directive 2010/63/EU).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.