Next Article in Journal
Arabidopsis SUMO E3 Ligase SIZ1 Interacts with HDA6 and Negatively Regulates HDA6 Function during Flowering
Previous Article in Journal
Targeting Mitochondrial Metabolism as a Strategy to Treat Senescence
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental Epileptogenesis in a Cell Culture Model of Primary Neurons from Rat Brain: A Temporal Multi-Scale Study

1
Department of Neuropathology, Universitätsklinikum Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), 91054 Erlangen, Germany
2
Department of Psychiatry and Psychotherapy, Universitätsklinikum Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), 91054 Erlangen, Germany
3
NGS Core Unit, Institute of Human Genetics, Universitätsklinikum Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), 91054 Erlangen, Germany
4
Department of Epileptology, University Hospital Bonn, 53127 Bonn, Germany
*
Author to whom correspondence should be addressed.
These authors contributed equally.
Cells 2021, 10(11), 3004; https://doi.org/10.3390/cells10113004
Submission received: 10 September 2021 / Revised: 27 October 2021 / Accepted: 2 November 2021 / Published: 3 November 2021
(This article belongs to the Section Cells of the Nervous System)

Abstract

:
Understanding seizure development requires an integrated knowledge of different scales of organization of epileptic networks. We developed a model of “epilepsy-in-a-dish” based on dissociated primary neuronal cells from neonatal rat hippocampus. We demonstrate how a single application of glutamate stimulated neurons to generate spontaneous synchronous spiking activity with further progression into spontaneous seizure-like events after a distinct latency period. By computational analysis, we compared the observed neuronal activity in vitro with intracranial electroencephalography (EEG) data recorded from epilepsy patients and identified strong similarities, including a related sequence of events with defined onset, progression, and termination. Next, a link between the neurophysiological changes with network composition and cellular structure down to molecular changes was established. Temporal development of epileptiform network activity correlated with increased neurite outgrowth and altered branching, increased ratio of glutamatergic over GABAergic synapses, and loss of calbindin-positive interneurons, as well as genome-wide alterations in DNA methylation. Differentially methylated genes were engaged in various cellular activities related to cellular structure, intracellular signaling, and regulation of gene expression. Our data provide evidence that a single short-term excess of glutamate is sufficient to induce a cascade of events covering different scales from molecule- to network-level, all of which jointly contribute to seizure development.

1. Introduction

Epilepsy is not a single disease, but highly heterogeneous with diverse clinical syndromes, associated etiologies (i.e., structural, genetic, infectious, metabolic, immune, unknown), and many different mechanisms contributing to the epileptic phenotype and known frequent comorbidities. Temporal lobe epilepsy (TLE) is the most common epilepsy syndrome in adults. Hippocampal sclerosis (HS) is a frequent structural brain lesion identified in patients with drug-resistant focal epilepsy submitted to surgical treatment [1]. It is histopathologically characterized by segmental neuronal cell loss in the hippocampal formation [2]. Clinical disease presentation in patients with TLE often follows a characteristic pattern. An initial precipitating injury is reported in many of these patients, such as prolonged febrile seizures or status epilepticus (SE), followed by a clinically silent latency period before the onset of spontaneous, recurrent, and often drug-resistant seizures [3]. A prevailing hypothesis to best anticipate this phenomenon is secondary epileptogenesis [4,5]. A precipitating injury triggers a molecular cascade turning a normal brain into an epileptic disease condition [6]. A growing number of experimental studies specifically addressed the structural and molecular basis of specific disease pathomechanisms in epileptogenesis, e.g., the imbalance between excitatory and inhibitory neurotransmission [7], the role of fast-spiking interneurons in seizure initiation and network synchronization [8], aberrant neuronal circuitries and plasticity [9], seizure-induced aberrant neurogenesis [10,11], neuroinflammation [12,13], or energy metabolism [14], amongst many other more. Epigenetic alterations were more recently also proposed as a fundamental pathomechanism in epilepsy. Several studies have identified locus-specific and genome-wide alterations in DNA methylation patterns in experimental epilepsy models [15,16] or human surgical brain samples [17,18]. Aberrant DNA methylation was linked to the compromised execution of gene expression programs and thereby suggested to account for many known structural and functional changes in the epileptic brain [15,16].
While seizures are commonly defined as the clinical manifestation of an abnormal, excessive, hypersynchronous discharge of a population of cortical neurons [19], it remains unclear to date what the minimum cellular and molecular requirements are that induce and promote the complex sequence of events leading to seizure development. Many experimental animal models have been developed in the past to induce and study acute seizures or epileptogenesis [20], e.g., following neonatal hyperthermia [21,22], hypoxia [23,24,25,26,27], traumatic brain injury [28,29,30], tetanus toxin injection [31,32,33], pilocarpine induced SE [34,35], kainic acid induced SE [36], or electrical kindling [37]. However, animal models have a limited ability to recapitulate human pathologies even though they appear to symptomatically mimic human disease. Furthermore, they are highly cost and labor- and time-consuming. Additionally, the complexity of the mammalian brain limits the possibilities of addressing questions of cause and consequence, and it remains challenging to attribute specific disease pathomechanisms over time to individual cell populations. In vitro models, including acute hippocampal slice preparations, organotypic hippocampal slice cultures [38], and iPSC-based brain organoid models [39] may be used to address some of these limitations. Moreover, dissociated primary neuronal cultures and neuronal cell lines are used to study basic electrophysiological properties and screen, e.g., for neurotoxicity of certain compounds [40,41,42].
We recently developed a simple and powerful cell culture model of epileptogenesis using a single transient glutamate-induced neuronal mass discharge to transform a primary dissociated cell culture from neonatal rat hippocampus into a repetitively firing synchronous seizure-like network, i.e., “epilepsy-in-a-dish” [43]. To understand seizure development, an integrated understanding of different scales of organization of epileptic networks across time is needed. However, research domains in epilepsy are commonly studied either alone or only at a single time point. The present multi-scale study aimed to quantitatively investigate network and single-cell activity and structural remodeling processes, including dendritic sprouting with or without new synapse formation. We further used whole-genome bisulfite sequencing (WGBS) to study the regulatory role of DNA methylation on biological pathways implicated in the control of cell structure and function. With this model, we can monitor and interfere with the timeline of molecular, structural, and functional changes during epileptogenesis at the neuronal cell and network level.

2. Materials and Methods

2.1. Animals and Tissue Preparation

Animal and tissue preparation was performed as previously described [43]. In short, adult male and female Wistar rats were obtained from Charles River (Sulzfeld, Germany), bred and maintained at the local animal housing facility in breeding cages under controlled environmental conditions (12 h light/dark cycle, 20–23 °C, 50% relative humidity, drinking and feeding ad libitum). Male and female rat pups (P0–P2) were used. All animal experiments were approved by the local animal care and use committee (TS-1/13) and followed the European Communities Council Directive and German Animal Welfare Act (54–2532.1-23/09, Directive 2010/63/EU).

2.2. Cell Culture

Primary neuronal cell cultures from rat neonatal hippocampus were prepared as described before [43]. Briefly, cell suspensions from dissociated newborn rat hippocampi (P0–P2) were preplated for 1 h at 37 °C, 5% CO2 onto an uncoated flask to allow settlement of glial cells and their adhesion to the flask. Remaining cells from the supernatant were harvested 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 coverslips (Ø 1cm; Greiner Bio-One, Kremsmünster, Austria) at a density of 2.5 × 105 cells/well. Cells were maintained at 37 °C in a fully humidified incubator with 5% CO2. After 24 h, 6 µM Cytosine β-D-arabinofuranoside hydrochloride (AraC; Sigma-Aldrich, Taufkirchen, Germany) was added to inhibit proliferation of remaining glial cells. Neurons were continuously maintained in dispersed culture with their original media supplemented with AraC. After 12 days in vitro, the original culture medium was transiently replaced by physiological treatment solution (10 mM HEPES (pH 7.4), 145 mM NaCl, 10 mM glucose, 1 mM MgCl2, 2.5 M KCl, 2 mM CaCl2, 2 µM glycine) with or without 10 µM glutamate [44,45]. Cultures were then washed three times, resupplied with the original culture medium with AraC, and kept at 37 °C in a fully humidified incubator containing 5% CO2 until further experiments.

2.3. Viability Assay

Cellular viability was determined 7 and 15 days after glutamate treatment with propidium iodide (PI) and fluorescein diacetate (FDA). Cells were incubated with 500 nM PI and 160 nM FDA for 10 min in Neurobasal A medium at 37 °C. Semiquantitative measurements of PI/FDA-stained cells were performed using a microcomputer imaging system (ColorView II CCD camera, Cell^F imaging software) 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 routinely determined in glutamate-treated cultures and time-matched sham controls.

2.4. Immunofluorescence

Immunofluorescence staining for Sholl analysis was performed 1, 3, 5, 7, and 15 days after glutamate or sham treatment, for synapse quantification after 3 and 7, as well as for interneuron detection 7 days following initial treatment. Neurons cultured on coverslips were fixed in 4% paraformaldehyde and blocked in 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, Taufkirchen, Germany) and Triton X-100 (0.1%; Sigma-Aldrich) for 2 h at room temperature. Primary antibodies were diluted in blocking solution and applied overnight at 4 °C. Secondary fluorophore-conjugated antibodies were diluted in blocking solution and applied for 4 h at room temperature. Nuclei were dyed with 4′,6-diamidino-2-phenylindole (DAPI; Sigma-Aldrich, Taufkirchen, Germany, 1:1000) diluted in blocking solution. For Sholl analysis, primary chicken anti-Map2 (Abcam, Cambridge, UK, Cat# ab5392, RRID:AB_2138153, 1:2000) and mouse anti-Gfap (Cell Signaling Technology, Frankfurt a.M., Germany, Cat# 3670, RRID:AB_561049, clone GA5, 1:300) and secondary Alexa-647 donkey anti-mouse (Abcam, Cambridge, UK, Cat# ab150107, RRID:AB_2890037, 1:1000) and Alexa-488 goat anti-chicken (Thermo Fisher Scientific, Bremen, Germany, Cat# A-11039, RRID:AB_2534096, 1:1000) were used. For interneuron quantification, primary rabbit anti-Calbindin (Synaptic Systems, Göttingen, Germany, Cat# 214 002, RRID:AB_2068199) and secondary Alexa-555 goat anti-rabbit (Thermo Fisher Scientific, Bremen, Germany, Cat# A-21428, RRID:AB_2535849, 1:1000) were used. For synapse quantification, primary rabbit anti-Vglut1 (Synaptic Systems, Göttingen, Germany, Cat# 135 308, RRID:AB_2864787, 1:500), chicken anti-Vgat (Synaptic Systems, Göttingen, Germany, Cat# 131 006, RRID:AB_2619820, 1:500), and mouse anti-Bassoon (Enzo Life Sciences, Lörrach, Germany, Cat# SAP7F407, RRID:AB_2313990, 1:300), together with secondary Alexa-488 goat anti-chicken and Alexa-555 goat anti-rabbit (Thermo Fisher Scientific, Bremen, Germany, Cat# A-21428, RRID:AB_2535849, 1:1000), as well as Alexa-647 donkey anti-mouse (Abcam, Cambridge, UK, Cat# ab150107, RRID:AB_2890037, 1:1000) were used.

2.5. Sholl Analysis

Immunofluorescence microscopy for Sholl analysis was performed with the Zeiss MN Observer microscope (Zeiss, Jena, Germany) and ZEN 2012 SP5 software. Singular neurons were tracked in the periphery of cultures, avoiding overlay of different cells and false-positive attribution of dendrites. All images were calibrated to scale for Sholl analysis afterwards. After processing of images, we used the Sholl analysis plug-in v3.7.0 [46] for Fiji software v1.51s [47] to describe dendritic arbors on singular neurons. A total of 60 neurons were analyzed per time point after glutamate treatment with 30 neurons per treatment, with or without glutamate, respectively.

2.6. Synapse Quantification

Synapses were visualized using a Zeiss Observer Z1 LSM 780 microscope (Zeiss, Jena, Germany) equipped with 405/488/555/647 HXP 120 V Fluorescence lamp, 405/488/555/647 LASOS Argon Laser, and 594 Zeiss In Tune Laser, using ZEN 2012 SP5 software. Images were taken with a 1 mm pinhole using a 63× objective covering an area of 134 µm2. All images were acquired under standardized conditions. For synapse quantification, five non-overlapping images per coverslip and six independent preparations per experimental condition were used. Co-localized presynaptic (Bassoon-positive) and postsynaptic (Vglut1- or Vgat-positive) puncta within a mask were counted manually using the cell counter plug-in for Fiji (1.51s “Multi-point” tool provided by Fiji). In a second approach, we used the automated Fiji plug-in Synapse Counter. During quantification, settings for synapse counting were uniformly chosen for all images within an experiment. Results of automatic puncta quantification were always validated by visual inspection and manual counting for a subset of images. All coverslips were coded and processed in a blinded manner. We normalized all puncta counts to the median number of puncta detected on control neurons for a particular coverslip to combine data acquired from different coverslips in different experiments. An average of 50,000 synapses was analyzed per experimental condition and time point.

2.7. Calcium Imaging

Calcium imaging was performed 1, 3, 5, 7, 9, and 15 days after glutamate or sham treatment. Cultures were incubated for 30 min at 37 °C with 5 µM Fluo-8 AM (Abcam, Cambridge, UK,) in 1 µg/µL DMSO. The osmolarity of modified Cohen’s Recording Solution (1x SBM: 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; [48]) was adjusted to the medium’s osmolarity and 5 µM Fluo-8 AM was added. After incubation, coverslips were placed into a perfusion chamber filled with 450 µL recording solution. Nikon Eclipse Ti’s inverted fluorescence microscope (Nikon, Tokyo, Japan) with a Lambda SC shutter was used for the recording. An area of 819.2 µm2 was detected with a Plan Apo 10×, 0.45 NA objective, and perfect focus system in every experiment. Fluo-8 AM was excited by Nikon Intensilight C-HGFI Lamp (Neutral Density Filter 16) through a 455–485 nm excitation and a 495 nm dichroic long-pass mirror. Emitted light was passed through a 488 nm filter Cube GFP (Semrock, Rochester, NY, USA). Detection was performed with a cooled EM-CCD camera (iXonEM DU-885; DU-887, Andor Oxford Instruments, Belfast, Northern Ireland, UK). The culture was recorded with the NIS Elements (Nikon, Tokyo, Japan) software in two different settings: no binning and 200 ms exposure time or 2 × 2 binning and 50 ms exposure time. For all experimental conditions, 5–10 biological replicates from at least three different preparations were analyzed.
Cellular calcium signals were quantified as the average light intensity of manually selected regions of interest (ROI; Fiji 1.51s) covering cell bodies. Per-pixel standard deviation was used to identify neurons (ROI) with significant changes in calcium signal. Signals from 80–150 neurons were analyzed in every recording session. Individual traces were exported and analyzed in a custom-developed software designed by V.G. in LabView (National Instruments, Austin, TX, USA) [49].

2.8. Whole Genome Bisulfite Sequencing

WGBS paid service was performed at the Core Unit NGS (FAU Erlangen, Germany) and included DNA extraction, bisulfite conversion, library preparation, sequencing, and data analysis. Sequencing was performed with Illumina HiSeq2500 (Illumina, San Diego, CA, USA). Raw paired-end FASTQ files (per sample, flowcell, read direction, and lane) were trimmed with Trim Galore v0.6.2 using cutadapt v2.3, autodetecting and removing Illumina adapter sequences. Trimmed FASTQ files were then merged with zcat/pigz into a single FASTQ file pair per sample. For bisulfite-sensitive alignment and differential methylation calling, Bismark v0.22.1 with hisat2 v2.1.0 mapping was used [50]. Reads were aligned to the rn6 reference genome, with Ensembl rnor6.0.95 as an annotation reference. Differential methylation was examined using the Bioconductor package DSS v2.32.0 [51] from the methylation counts as reported by Bismark. Upon pairwise comparison of biological groups, differentially methylated positions (DMPs) were called and exported (FDR < 0.05). Then, differentially methylated regions (DMRs) with a minimum number of 4 consecutive DMPs were called using the bumphunter package [52]. DMRs (p < 0.001) that mapped to enhancers, promotors (2 kb upstream of transcriptional start site/TSS), or gene bodies were used to identify differentially methylated genes (DMGs). The biological function of gene networks was analyzed using Qiagen Ingenuity Pathway Analysis (IPA) software (Qiagen, Hilden, Germany).

2.9. Statistics

Statistical analysis was performed using GraphPad Prism7 (GraphPad Software, San Diego, CA, USA). All points shown in this work represent an average from at least three independent measurements in three independent preparations, including epileptogenic and sham control cultures obtained from the same preparations (i.e., sister cultures). Depending on sample size and distribution, paired or unpaired two-tailed t-test (2 groups), ANOVA (≥3 groups), or respective non-parametric tests and corrections for multiple comparisons were performed. Details are provided in the respective Results section. In all figures, asterisks indicate a level of significance of p < 0.05, and error bars and specified values represent mean and standard error of the mean (SEM), if not stated otherwise.

3. Results

This section may be divided by subheadings. It should provide a concise and precise description of the experimental results, their interpretation, as well as the experimental conclusions that can be drawn.

3.1. “Epilepsy-in-a-Dish”

We previously developed an in vitro model of epileptogenesis based on a self-organized network of interconnected neurons and a few glial cells [43]. Starting from dissociated primary neuronal cells from neonatal rat hippocampus, we showed that a single and transient (10 min) application of physiological levels of 10 µM glutamate stimulated neurons to generate spontaneous synchronous spiking activity progressing into spontaneous seizure-like events (SLE). However, a comprehensive characterization and temporal resolution of network activity dynamics following glutamate injury were missing yet. In the present study, we performed a temporal series of live cell calcium imaging to functionally characterize our model and quantify the dynamics of the potential epileptogenic process. Calcium imaging using Fluo8-AM green fluorescent calcium-binding dye was performed 1, 3, 5, 7, 9, and 15 days after glutamate injury (Figure 1), and the neuronal activity of up to 150 cells per visual field was quantified. Neuronal cultures showed no synchronous activity one day after treatment. After three days, single spikes were detected in individual neurons. At five days after glutamate treatment, synchronous spikes, polyspikes, and a few SLEs were observed. At seven days, 60–70% of all glutamate-treated cultures showed SLEs. Almost all cultures showed synchronous spiking activity compared with time-matched sham-treated controls, which remained silent (Figure 1a,b). SLEs became significantly longer in the second week after glutamate treatment, suggesting further progression of the epileptogenic process (Mann–Whitney test, mean SLE duration1st week = 49 ± 10 s, n = 18, SLE duration2nd week = 111 ± 6 s, n = 11, p = 0.0014; Figure 1c).
Computational analysis and visual comparison of in vitro SLEs with human intracranial stereo-EEG recordings of seizures identified strong similarities and a related sequence of events with defined onset, progression, and termination (Figure 1d). SLEs typically showed three prominent features: a build-up phase with increasing spiking activity, followed by a direct current (DC) shift with on-top oscillatory activity starting from high frequency then slowing down (=chirps in power spectrum), and post-SLE depression (Figure 1b,d).

3.2. Structure Determines Function

We next thought to analyze whether the functional differences between glutamate-treated and control cultures resulted from structural changes of cultured cells or altered cellular assembly of the networks under study. Therefore, we monitored cultured neurons (Map2 immunoreactivity) for up to 15 days after glutamate treatment and quantified neurite length and sprouting using Sholl analysis (Figure 2). With time, glutamate-treated neurons displayed significantly longer processes (paired two-tailed t-test, mean lengthCtrl1d = 241.8 ± 70.2 µm, lengthGlu1d = 263.1 ± 76.9 µm, p = 0.33; lengthCtrl3d = 185.9 ± 59.8 µm, lengthGlu3d = 241.1 ± 58.7 µm, p = 0.0003; lengthCtrl5d = 187.6 ± 54.9 µm, lengthGlu5d = 241.1 ± 68.9 µm, p = 0.0014; lengthCtrl7d = 227.4 ± 37.2 µm, lengthGlu7d = 303.8 ± 74.0 µm, p < 0.0001; lengthCtrl15d = 272.9 ± 68.5 µm, lengthGlu15d = 289.5 ± 101.7 µm, p = 0.62).
When we plotted the number of neurite intersections against distance from soma: neurons following glutamate treatment displayed fewer intersections close to the soma, but more intersections with increasing distance from the soma, as compared with time-matched sham-treated controls (paired two-tailed t-test, mean branchingCtrl1d = 7.0 ± 4.3, Glu1d = 4.9 ± 3.3, p < 0.0001; Ctrl3d = 4.6 ± 2.8, Glu3d = 3.5 ± 2.3, p < 0.0001; Ctrl5d = 4.8 ± 2.9, Glu5d = 5.2 ± 3.6, p = 0.62; Ctrl7d = 7.8 ± 4.7, Glu7d = 5.9 ± 3.8, p < 0.0001; Ctrl15d = 5.9 ± 4.1, Glu15d = 4.0 ± 3.0, p < 0.0001; Figure 2a,b).
We next performed in-depth manual and automated quantification of glutamatergic and GABAergic synapses. We identified a relative loss of active GABAergic synapses stained for vesicular GABA transporter (Vgat) over excitatory synapses stained positive for vesicular glutamate transporter 1 (Vglut1) seven days after glutamate treatment compared with time-matched sham controls (unpaired two-tailed t-test, %VgatCtrl7d = 43.4 ± 3.2, %VgatGlu7d = 24.0 ± 3.1, p < 0.001). This difference was not yet detectable three days after epileptogenesis induction, i.e., following glutamate exposure (unpaired two-tailed t-test, mean %VgatCtrl3d = 35.0 ± 1.2, %VgatGlu3d = 37.9 ± 3.5, p > 0.05). Only synapses co-labeled for Bassoon (Bsn, red), a scaffolding protein involved in organizing the presynaptic cytoskeleton of active axon terminals, and Vgat or Vglut1 were included in the measurements (Figure 3a–c, Figure S1).
To examine whether the loss in GABAergic synapses resulted from a loss in GABAergic interneurons or rather from a change in synaptic plasticity, immunofluorescence staining was performed for neuronal microtubule-associated protein 2 (Map2), glial fibrillary acidic protein (Gfap), and the interneuronal marker protein Calbindin (Calb1). A general decline in Calb1-positive cells relative to overall cell numbers was observed in cultures seven days after glutamate treatment (unpaired two-tailed t-test, %Calb1Ctrl7d = 43.1 ± 1.6, %Calb1Glu7d = 36.7 ± 2.5, p = 0.031; Figure 3d,e). Taken together, the data suggest a time-dependent shift towards increased excitability.

3.3. Epigenetic Regulation of Causative Gene Networks

Activity-driven alterations of gene expression patterns induce structural and functional changes in neurons, critical for long-term adaptations of brain circuits and indispensable for proper brain function [53,54]. To test whether functional and structural changes in our model could be linked with potentially gene-regulatory DNA methylation signatures, we performed whole-genome bisulfite sequencing (WGBS) in primary neuronal cell cultures seven days after glutamate treatment (n = 4) and in time-matched controls (n = 5).
From pairwise comparisons, 538 differentially methylated positions (DMP) were identified with an FDR < 0.05 (Figure 4a). Principal component analysis for dimensionality reduction and hierarchical clustering showed that all samples clustered in their respective treatment groups (control—blue, glutamate–red; Figure 4b). Both increase and decrease in DNA methylation across the genome were observed in the disease model (Figure 4b,c). Next, differentially methylated regions (DMRs, p < 0.001, n = 279) were called, which uniquely mapped to 130 differentially methylated genes (DMGs). We then performed functional pathway analysis using Qiagen’s ingenuity pathway analysis software. There was an enrichment of DMG in the axonal guidance signaling pathway in line with the dominant structural changes that we detected in neuronal networks following glutamate treatment. Next, functional interaction networks of proteins encoded by hyper- (green) and hypomethylated genes (red) were identified and mapped to the cell as a structural and functional unit. Thereby the spatial location of protein nodes was shifted into their corresponding subcellular locations (Figure 4c). DMGs and their encoded proteins were found to be part of a broad range of cellular structures and functional processes including extracellular matrix (ECM) structure; plasma membrane composition, including membrane-bound and soluble signaling factors, ion channels, transporters, and receptors; components of the actin cytoskeleton; intracellular signaling and metabolism; or proteins involved in chromatin structure and function (Table 1). Many of these proteins play a significant role in nervous system development, function, and disease (Figure 4c, pink bordered icons; Table 1). Our data suggest that genomic DNA methylation is altered early on in the epileptogenic process and may play a central role in regulating pro-epileptogenic neuronal gene expression and downstream cellular structure and function throughout the entire disease process.

4. Discussion

Epileptogenesis, i.e., the process that leads to the development of spontaneous recurrent seizures and disease progression, cannot be studied in the human brain in situ. Only in a subset of patients undergoing epilepsy surgery for treatment of their focal drug-resistant epilepsy do we have access to human brain tissue. However, the described molecular, structural, and functional changes in such tissue are representative of a late (chronic) disease stage and do not allow conclusions to be drawn on when these changes occurred. Therefore, experimental models are needed to study the disease course. Here, we describe a simple but powerful model providing insights into the full process of epileptogenesis in an exceptionally short time frame covering all steps from initial precipitating injury to silent latency period, development of spontaneous recurrent seizure-like events (SLEs), and disease progression. We performed an integrated study addressing multiple scales from molecules to cellular and network structure and neurophysiological function. Our data show striking and broad similarities between our cellular model system with human epilepsy across scales including similar seizure pattern and duration, cell structure, and a specific phenotype-associated gene-regulatory epigenetic profile.
Rodent models of chronic epilepsy have been used extensively in the past to study epileptogenesis. Still, they have been proven to be labor- and cost-intensive and challenging to study in sufficient temporal, cellular, and molecular resolution. Acute or organotypic slice cultures are often used for more detailed studies on seizure mechanisms, as they sustain much of the necessary circuitry. While acute preparations are viable for a few hours and valuable for studying acute provoked seizures, organotypic slice cultures are sustainable for multiple weeks. They can be used to study temporal processes, including circuit development, synaptic plasticity, and axonal sprouting. However, interpretation of molecular data from brain tissue with its high cellular heterogeneity remains challenging, as results cannot easily be attributed to a specific cell type. Approaches that help address this issue include cell sorting before molecular analysis, microfluidics and subsequent single-cell sequencing, or culturing select cell populations, all of which have specific technical and biological advances and limitations. Here, we used primary neuronal cells from rat HC that, after maturation in vitro, self-assembled into 3D functional networks (without the support of an artificial 3D matrix) and developed spontaneous recurrent epileptiform activity following a single transient stimulus with physiological levels of glutamate. With a stable 92–96% neuronal phenotype, our cultures could be used for a nearly unbiased study of neuronal features of seizure generation.
Epileptogenesis was monitored following glutamate exposure using live-cell calcium imaging at different time points. After 24 h, no specific neuronal activity was detected except for episodic random unsynchronized spikes. At 3 days after the transient glutamate exposure, synchronous spiking activity was detected and remained prominent at 5, 7, 9, and 15 days. SLEs emerged 5 days after treatment and were maintained in the 7-, 9-, and 15-day cultures. Cell type attribution of detected signals was always confirmed by Map2/Gfap immunofluorescence of all calcium-imaged cultures. We identified specific but extremely rare events of glial calcium activity (i.e., low frequency, high amplitude, slowly propagating along glial processes), which only affected individual cells and never occurred in temporal relation with synchronized spiking or SLE development across the network. SLEs in our cultures showed strong similarities with human intracranial EEG seizures, both in terms of duration, stability for around 1–2 min, and according to seizure pattern. The similarities between experimental and human seizure power spectra were striking despite the differences in resolution. Our results indicate the development of spontaneous recurrent epileptiform activity after a short silent period and progressive worsening with more prolonged SLEs over time, thus recapitulating the presumed general pathogenic time course of epileptogenesis. In our study design, no continuous measurement of epileptiform activity was possible. Preparations of cultures on (high-density) multi-electrode arrays would allow measurements of field potentials to examine network activity at higher temporal resolution, however not at single-cell resolution. Other complementary studies could involve whole-cell patch-clamp approaches, enabling manipulations of individual cells and determining downstream small network effects if several cells over the network are monitored simultaneously.
Structural changes of neuronal cells during epileptogenesis were studied using Map2 immunofluorescence staining and Sholl analysis. The neurite length of neurons in glutamate-treated cultures was significantly increased early on in the epileptogenic process. This notable change in cellular structure persisted over time, together with a higher degree of terminal branching. Altered cell morphology in our model matched structural changes previously described ex vivo in humans [55] and rats [16], where dendritic and axonal sprouting have been identified and suggested to contribute to the proepileptogenic reorganization of hippocampal circuits. In addition, we identified an increase in active glutamatergic over GABAergic synapses at 7 but not 3 days after glutamate treatment compared with time-matched sham controls. Variability in relative Vglut1- and Vgat-positive synapse quantities may have been related to individual preparations. We aimed to control this technical factor by analyzing sister cultures when comparing glutamate-treated to sham controls. The observed shift in synaptic composition towards an increase in active glutamatergic synapses was evident, however, and may contribute to an increase in excitability of cell cultures.
To further investigate DNA methylation changes, a more permanent regulator for gene expression, we performed WGBS of cultures 7 days after initial precipitating injury. Principal component analysis and hierarchical clustering revealed a highly distinct pattern in genomic DNA methylation for glutamate- versus sham-treated neurons. Thereby, our data support extensive epigenomic changes even as short as 7 days after the treatment, confirming previous reports of DNA methylation dynamics in epileptogenesis [15,16,18,43,56,57]. Our data further support the notion of a leading regulatory role in pathologic gene expression in epilepsy. Both acquisition and loss of DNA methylation were observed to target genic and intergenic parts of the genome. These findings were in line with previous reports from Kobow et al. in vivo [15,16]. Next, differentially methylated CpGs were clustered to identify biologically meaningful differentially methylated regions (DMR) and then mapped to gene bodies, promoters, and enhancers. Thereby, 130 genes were targeted by differential methylation in epileptogenesis and encoded proteins with a broad range of functions throughout the entire cell, from ECM and cell membrane components to signaling molecules in the cytoplasm and structural proteins of the cytoskeleton, as well as nuclear factors involved in chromatin structure and function (Table 1). Functional pathway analysis specifically highlighted axon and neurite growth. This was in line with the broad structural changes and increased neurite outgrowth that we detected in our model. While we identified no further enrichment of differentially methylated genes (DMGs) in Kyoto Encyclopedia of Genes and Genomes (KEGG) or Reactome pathways, which is likely due to an annotation bias in these databases, ingenuity pathway analysis revealed several functional protein interaction networks. In-depth analysis of specific candidate genes highlighted various roles in brain development, function, and disease, including neurodegenerative disease, autism spectrum disorders, intellectual disability, developmental brain malformations, and epilepsy. We detected hypermethylation of Potassium Voltage-Gated Channel Subfamily H Member 5 (Kcnh5), an outward-rectifying, non-inactivating potassium channel that regulates neurotransmitter and hormone release. Diseases associated with KCNH5 mutations include early infantile epileptic encephalopathy (EIEE). We also found hypermethylation of collagen type XIV, alpha 1 (Col14a1), an ECM component that is highly expressed in the brain and has been linked to familial adult myoclonic epilepsy (FAME). We further identified hypermethylation of doublecortin (Dcx). Dcx is a microtubule-associated protein involved in neuronal migration. Mutations in human DCX have been identified as being causative in double cortex syndrome and other cortical malformations associated with epilepsy. Hence, although we lack a transcriptomic profile of our cultures, these data support a gene regulatory role of DNA methylation, which could explain many of the structural and functional changes seen in our model of epileptogenesis. The diversity of gene functions potentially affected by DNA methylation changes in our model are in line with the multifactorial nature of epilepsy and associated comorbidities and its link to brain development, function, and aging.
Characteristic genomic methylation signatures for different disease stages have been reported in various species, etiologies, and pathologies [15,16,17,18,57,58]. A comprehensive integration and comparison of genomic DNA methylation data obtained in the present study (WGBS), with published data from different rat epilepsy models (Methyl-Seq), and human focal epilepsy (Methyl-Seq, Illumina arrays) has not yet been performed. Such integrated analysis with the aim of identifying common epileptogenic mechanisms would be challenging with regard to species differences, annotation biases, platform variability, normalization algorithms, cutoff points for selecting regulated genes, analysis of different brain structures and cells, use of variable insults to trigger epileptogenesis, selection of time points for sampling after the insult, and characterization of epilepsy phenotype at the time of sampling. However, comparison of DNA methylation profiles from Methyl-Seq in the epileptic hippocampus obtained from three different rat epilepsy models (i.e., pilocarpine, amygdala stimulation, and traumatic brain injury) distinguished controls from epileptic animals but failed to identify a joint methylation signature in all three different models and few regions common to any two models despite the fact that the study design involved control for the same species, same time point/age of sampling, same tissue, same profiling platform, and the same sequencing facility. Debski et al. concluded that their data provided evidence that genome-wide alteration of DNA methylation signatures is a general pathomechanism associated with epileptogenesis and epilepsy in experimental animal models, but the broad pathophysiological differences between models were reflected in distinct etiology-dependent DNA methylation patterns [15]. Monitoring DNA methylation at 7 days after glutamate injury in our model showed that epigenetic alterations were not confined to single genes, extending previous reports of locus-specific epigenetic gene regulation in this model [43]. Instead, our data further confirmed that genome-wide aberrant DNA methylation is a conserved fundamental gene regulatory mechanism in the pathogenesis of epilepsy. Future studies will be needed to identify specific and reproducible DNA-methylation-based biomarkers of epileptogenesis and the risk of developing seizures after an initial precipitating injury—both in vitro and in vivo.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cells10113004/s1, Figure S1: Z-stack of a confocal image showing co-localization of Vglut1 (blue) and Vgat (green) with Bassoon (red) marking active synaptic terminals

Author Contributions

Conceptualization, K.K. and I.B.; methodology, J.J. and L.H.; calcium imaging analysis, V.G.; whole genome bisulfite sequencing and data analysis, A.B.E. and S.U.; resources, A.F.; writing—original draft preparation, K.K.; writing—review and editing, all authors; visualization, V.G. and K.K.; supervision, K.K.; project administration, K.K.; funding acquisition, I.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by European Union FP7 HEALTH, project “DESIRE”, grant agreement number 602531.

Institutional Review Board Statement

The study was approved by the local animal care and use committee (TS-1/13) following the European Communities Council Directive and German Animal Welfare Act (54–2532.1-23/09, Directive 2010/63/EU).

Informed Consent Statement

Not applicable.

Data Availability Statement

The WGBS data presented in this study will be made openly available in GEO under accession number GSEXXX.

Acknowledgments

We kindly thank Birte Rings for her expert technical assistance. We further thank Laura Tassi from the Claudio Munari Epilepsy Surgery Center of the Niguarda Hospital in Milan, Italy, for providing the intracerebral data on human seizures.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Blumcke, I.; Spreafico, R.; Haaker, G.; Coras, R.; Kobow, K.; Bien, C.G.; Pfäfflin, M.; Elger, C.; Widman, G.; Schramm, J.; et al. Histopathological Findings in Brain Tissue Obtained during Epilepsy Surgery. N. Engl. J. Med. 2017, 377, 1648–1656. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Blumcke, I.; Thom, M.; Aronica, E.; Armstrong, D.D.; Bartolomei, F.; Bernasconi, A.; Bernasconi, N.; Bien, C.G.; Cendes, F.; Coras, R.; et al. International consensus classification of hippocampal sclerosis in temporal lobe epilepsy: A Task Force report from the ILAE Commission on Diagnostic Methods. Epilepsia 2013, 54, 1315–1329. [Google Scholar] [CrossRef] [PubMed]
  3. Sloviter, R.S. Hippocampal epileptogenesis in animal models of mesial temporal lobe epilepsy with hippocampal sclerosis: The importance of the "latent period" and other concepts. Epilepsia 2008, 49 (Suppl. 9), 85–92. [Google Scholar] [CrossRef] [PubMed]
  4. Morrell, F. Secondary epileptogenesis in man. Arch. Neurol. 1985, 42, 318–335. [Google Scholar] [CrossRef] [PubMed]
  5. McNamara, J.O. Emerging insights into the genesis of epilepsy. Nature 1999, 399, A15–A22. [Google Scholar] [CrossRef] [PubMed]
  6. Pitkänen, A.; Lukasiuk, K. Molecular and cellular basis of epileptogenesis in symptomatic epilepsy. Epilepsy Behav. 2009, 14, 16–25. [Google Scholar] [CrossRef] [PubMed]
  7. Buckingham, S.C.; Campbell, S.L.; Haas, B.R.; Montana, V.; Robel, S.; Ogunrinu, T.; Sontheimer, H. Glutamate release by primary brain tumors induces epileptic activity. Nat. Med. 2011, 17, 1269–1274. [Google Scholar] [CrossRef]
  8. Magloire, V.; Mercier, M.S.; Kullmann, D.M.; Pavlov, I. GABAergic Interneurons in Seizures: Investigating Causality With Optogenetics. Neuroscientist 2019, 25, 344–358. [Google Scholar] [CrossRef] [Green Version]
  9. Goldberg, E.M.; Coulter, D.A. Mechanisms of epileptogenesis: A convergence on neural circuit dysfunction. Nat. Rev. Neurosci. 2013, 14, 337–349. [Google Scholar] [CrossRef] [Green Version]
  10. Blumcke, I.; Schewe, J.C.; Normann, S.; Brustle, O.; Schramm, J.; Elger, C.E.; Wiestler, O.D. Increase of nestin-immunoreactive neural precursor cells in the dentate gyrus of pediatric patients with early-onset temporal lobe epilepsy. Hippocampus 2001, 11, 311–321. [Google Scholar] [CrossRef]
  11. Jessberger, S.; Zhao, C.; Toni, N.; Clemenson, G.D., Jr.; Li, Y.; Gage, F.H. Seizure-associated, aberrant neurogenesis in adult rats characterized with retrovirus-mediated cell labeling. J. Neurosci. 2007, 27, 9400–9407. [Google Scholar] [CrossRef]
  12. Kim, S.Y.; Buckwalter, M.; Soreq, H.; Vezzani, A.; Kaufer, D. Blood-brain barrier dysfunction-induced inflammatory signaling in brain pathology and epileptogenesis. Epilepsia 2012, 53 (Suppl. 6), 37–44. [Google Scholar] [CrossRef] [Green Version]
  13. Vezzani, A.; French, J.; Bartfai, T.; Baram, T.Z. The role of inflammation in epilepsy. Nat. Rev. Neurol. 2011, 7, 31–40. [Google Scholar] [CrossRef] [Green Version]
  14. Patel, M. A Metabolic Paradigm for Epilepsy. Epilepsy Curr. 2018, 18, 318–322. [Google Scholar] [CrossRef] [Green Version]
  15. Debski, K.J.; Pitkanen, A.; Puhakka, N.; Bot, A.M.; Khurana, I.; Harikrishnan, K.N.; Ziemann, M.; Kaspi, A.; El-Osta, A.; Lukasiuk, K.; et al. Etiology matters—Genomic DNA Methylation Patterns in Three Rat Models of Acquired Epilepsy. Sci. Rep. 2016, 6, 25668. [Google Scholar] [CrossRef]
  16. Kobow, K.; Kaspi, A.; Harikrishnan, K.N.; Kiese, K.; Ziemann, M.; Khurana, I.; Fritzsche, I.; Hauke, J.; Hahnen, E.; Coras, R.; et al. Deep sequencing reveals increased DNA methylation in chronic rat epilepsy. Acta Neuropathol. 2013, 126, 741–756. [Google Scholar] [CrossRef] [Green Version]
  17. Kobow, K.; Jabari, S.; Pieper, T.; Kudernatsch, M.; Polster, T.; Woermann, F.G.; Kalbhenn, T.; Hamer, H.; Rossler, K.; Muhlebner, A.; et al. Mosaic trisomy of chromosome 1q in human brain tissue associates with unilateral polymicrogyria, very early-onset focal epilepsy, and severe developmental delay. Acta Neuropathol. 2020, 140, 881–891. [Google Scholar] [CrossRef]
  18. Kobow, K.; Ziemann, M.; Kaipananickal, H.; Khurana, I.; Muhlebner, A.; Feucht, M.; Hainfellner, J.A.; Czech, T.; Aronica, E.; Pieper, T.; et al. Genomic DNA methylation distinguishes subtypes of human focal cortical dysplasia. Epilepsia 2019, 60, 1091–1103. [Google Scholar] [CrossRef] [Green Version]
  19. Fisher, R.S.; Acevedo, C.; Arzimanoglou, A.; Bogacz, A.; Cross, J.H.; Elger, C.E.; Engel, J., Jr.; Forsgren, L.; French, J.A.; Glynn, M.; et al. ILAE official report: A practical clinical definition of epilepsy. Epilepsia 2014, 55, 475–482. [Google Scholar] [CrossRef] [Green Version]
  20. Loscher, W. Animal Models of Seizures and Epilepsy: Past, Present, and Future Role for the Discovery of Antiseizure Drugs. Neurochem. Res. 2017, 42, 1873–1888. [Google Scholar] [CrossRef]
  21. Holtzman, D.; Obana, K.; Olson, J. Hyperthermia-Induced Seizures in the Rat Pup—a Model for Febrile Convulsions in Children. Science 1981, 213, 1034–1036. [Google Scholar] [CrossRef]
  22. Dube, C.; Chen, K.; Eghbal-Ahmadi, M.; Brunson, K.; Soltesz, I.; Baram, T.Z. Prolonged febrile seizures in the immature rat model enhance hippocampal excitability long term. Ann. Neurol. 2000, 47, 336–344. [Google Scholar] [CrossRef]
  23. Bergamasco, B.; Benna, P.; Ferrero, P.; Gavinelli, R. Neonatal hypoxia and epileptic risk: A clinical prospective study. Epilepsia 1984, 25, 131–136. [Google Scholar] [CrossRef]
  24. Chiba, S. Long-term effect of postnatal hypoxia on the seizure susceptibility in rats. Life Sci. 1985, 37, 1597–1604. [Google Scholar] [CrossRef]
  25. Jensen, F.E.; Wang, C. Hypoxia-induced hyperexcitability in vivo and in vitro in the immature hippocampus. Epilepsy Res. 1996, 26, 131–140. [Google Scholar] [CrossRef]
  26. Jensen, F.E.; Wang, C.; Stafstrom, C.E.; Liu, Z.; Geary, C.; Stevens, M.C. Acute and chronic increases in excitability in rat hippocampal slices after perinatal hypoxia In vivo. J. Neurophysiol. 1998, 79, 73–81. [Google Scholar] [CrossRef] [Green Version]
  27. Owens, J., Jr.; Robbins, C.A.; Wenzel, H.J.; Schwartzkroin, P.A. Acute and chronic effects of hypoxia on the developing hippocampus. Ann. Neurol. 1997, 41, 187–199. [Google Scholar] [CrossRef]
  28. Coulter, D.A.; Rafiq, A.; Shumate, M.; Gong, Q.Z.; DeLorenzo, R.J.; Lyeth, B.G. Brain injury-induced enhanced limbic epileptogenesis: Anatomical and physiological parallels to an animal model of temporal lobe epilepsy. Epilepsy Res. 1996, 26, 81–91. [Google Scholar] [CrossRef]
  29. Golarai, G.; Greenwood, A.C.; Feeney, D.M.; Connor, J.A. Physiological and structural evidence for hippocampal involvement in persistent seizure susceptibility after traumatic brain injury. J. Neurosci. 2001, 21, 8523–8537. [Google Scholar] [CrossRef]
  30. McIntosh, T.K.; Vink, R.; Noble, L.; Yamakami, I.; Fernyak, S.; Soares, H.; Faden, A.L. Traumatic brain injury in the rat: Characterization of a lateral fluid-percussion model. Neuroscience 1989, 28, 233–244. [Google Scholar] [CrossRef]
  31. Anderson, A.E.; Hrachovy, R.A.; Antalffy, B.A.; Armstrong, D.L.; Swann, J.W. A chronic focal epilepsy with mossy fiber sprouting follows recurrent seizures induced by intrahippocampal tetanus toxin injection in infant rats. Neuroscience 1999, 92, 73–82. [Google Scholar] [CrossRef]
  32. Jefferys, J.G.; Williams, S.F. Physiological and behavioural consequences of seizures induced in the rat by intrahippocampal tetanus toxin. Brain 1987, 110, 517–532. [Google Scholar] [CrossRef] [PubMed]
  33. Mitchell, J.; Gatherer, M.; Sundstrom, L.E. Loss of hilar somatostatin neurons following tetanus toxin-induced seizures. Acta Neuropathol. 1995, 89, 425–430. [Google Scholar] [CrossRef] [PubMed]
  34. Fujikawa, D.G. The temporal evolution of neuronal damage from pilocarpine-induced status epilepticus. Brain Res. 1996, 725, 11–22. [Google Scholar] [CrossRef]
  35. Mello, L.E.; Cavalheiro, E.A.; Tan, A.M.; Kupfer, W.R.; Pretorius, J.K.; Babb, T.L.; Finch, D.M. Circuit mechanisms of seizures in the pilocarpine model of chronic epilepsy: Cell loss and mossy fiber sprouting. Epilepsia 1993, 34, 985–995. [Google Scholar] [CrossRef]
  36. Sutula, T.; Cavazos, J.; Golarai, G. Alteration of long-lasting structural and functional effects of kainic acid in the hippocampus by brief treatment with phenobarbital. J. Neurosci. 1992, 12, 4173–4187. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. McNamara, J.O.; Constant Byrne, M.; Dasheiff, R.M.; Gregory Fitz, J. The kindling model of epilepsy: A review. Prog. Neurobiol. 1980, 15, 139–159. [Google Scholar] [CrossRef]
  38. Mckinney, R.A.; Debanne, D.; Gähwiler, B.H.; Thompson, S.M. Lesion-induced axonal sprouting and hyperexcitability in the hippocampus in vitro: Implications for the genesis of posttraumatic epilepsy. Nat. Med. 1997, 3, 990–996. [Google Scholar] [CrossRef]
  39. Costamagna, G.; Andreoli, L.; Corti, S.; Faravelli, I. iPSCs-Based Neural 3D Systems: A Multidimensional Approach for Disease Modeling and Drug Discovery. Cells 2019, 8, 1438. [Google Scholar] [CrossRef] [Green Version]
  40. Bradley, J.A.; Luithardt, H.H.; Metea, M.R.; Strock, C.J. In Vitro Screening for Seizure Liability Using Microelectrode Array Technology. Toxicol. Sci. 2018, 163, 240–253. [Google Scholar] [CrossRef]
  41. Giordano, G.; Costa, L.G. Primary neurons in culture and neuronal cell lines for in vitro neurotoxicological studies. Methods Mol. Biol. 2011, 758, 13–27. [Google Scholar]
  42. Otto, F.; Görtz, P.; Fleischer, W.; Siebler, M. Cryopreserved rat cortical cells develop functional neuronal networks on microelectrode arrays. J. Neurosci. Methods 2003, 128, 173–181. [Google Scholar] [CrossRef]
  43. Kiese, K.; Jablonski, J.; Hackenbracht, J.; Wrosch, J.K.; Groemer, T.W.; Kornhuber, J.; Blumcke, I.; Kobow, K. Epigenetic control of epilepsy target genes contributes to a cellular memory of epileptogenesis in cultured rat hippocampal neurons. Acta Neuropathol. Commun. 2017, 5, 79. [Google Scholar] [CrossRef] [Green Version]
  44. Basarsky, T.A.; Parpura, V.; Haydon, P.G. Hippocampal synaptogenesis in cell culture: Developmental time course of synapse formation, calcium influx, and synaptic protein distribution. J. Neurosci. 1994, 14, 6402–6411. [Google Scholar] [CrossRef] [Green Version]
  45. Srinivas, K.V.; Jain, R.; Saurav, S.; Sikdar, S.K. Small-world network topology of hippocampal neuronal network is lost, in an in vitro glutamate injury model of epilepsy. Eur. J. Neurosci. 2007, 25, 3276–3286. [Google Scholar] [CrossRef]
  46. Ferreira, T.A.; Blackman, A.V.; Oyrer, J.; Jayabal, S.; Chung, A.J.; Watt, A.J.; Sjöström, P.J.; van Meyel, D.J. Neuronal morphometry directly from bitmap images. Nat. Methods 2014, 11, 982–984. [Google Scholar] [CrossRef] [Green Version]
  47. Schindelin, J.; Arganda-Carreras, I.; Frise, E.; Kaynig, V.; Longair, M.; Pietzsch, T.; Preibisch, S.; Rueden, C.; Saalfeld, S.; Schmid, B.; et al. Fiji: An open-source platform for biological-image analysis. Nat. Methods 2012, 9, 676–682. [Google Scholar] [CrossRef] [Green Version]
  48. Cohen, E.; Ivenshitz, M.; Amor-Baroukh, V.; Greenberger, V.; Segal, M. Determinants of spontaneous activity in networks of cultured hippocampus. Brain Res. 2008, 1235, 21–30. [Google Scholar] [CrossRef]
  49. Gnatkovsky, V.; Francione, S.; Cardinale, F.; Mai, R.; Tassi, L.; Lo Russo, G.; de Curtis, M. Identification of reproducible ictal patterns based on quantified frequency analysis of intracranial EEG signals. Epilepsia 2011, 52, 477–488. [Google Scholar] [CrossRef]
  50. Krueger, F.; Andrews, S.R. Bismark: A flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics 2011, 27, 1571–1572. [Google Scholar] [CrossRef]
  51. Wu, H.; Wang, C.; Wu, Z. A new shrinkage estimator for dispersion improves differential expression detection in RNA-seq data. Biostatistics 2013, 14, 232–243. [Google Scholar] [CrossRef] [Green Version]
  52. Jaffe, A.E.; Murakami, P.; Lee, H.; Leek, J.T.; Fallin, M.D.; Feinberg, A.P.; Irizarry, R.A. Bump hunting to identify differentially methylated regions in epigenetic epidemiology studies. Int. J. Epidemiol. 2012, 41, 200–209. [Google Scholar] [CrossRef] [Green Version]
  53. Martinowich, K.; Hattori, D.; Wu, H.; Fouse, S.; He, F.; Hu, Y.; Fan, G.; Sun, Y.E. DNA methylation-related chromatin remodeling in activity-dependent BDNF gene regulation. Science 2003, 302, 890–893. [Google Scholar] [CrossRef] [Green Version]
  54. Guo, J.U.; Ma, D.K.; Mo, H.; Ball, M.P.; Jang, M.H.; Bonaguidi, M.A.; Balazer, J.A.; Eaves, H.L.; Xie, B.; Ford, E.; et al. Neuronal activity modifies the DNA methylation landscape in the adult brain. Nat. Neurosci. 2011, 14, 1345–1351. [Google Scholar] [CrossRef]
  55. Blumcke, I.; Zuschratter, W.; Schewe, J.C.; Suter, B.; Lie, A.A.; Riederer, B.M.; Meyer, B.; Schramm, J.; Elger, C.E.; Wiestler, O.D. Cellular pathology of hilar neurons in Ammon’s horn sclerosis. J. Comp. Neurol. 1999, 414, 437–453. [Google Scholar] [CrossRef]
  56. Miller-Delaney, S.F.; Bryan, K.; Das, S.; McKiernan, R.C.; Bray, I.M.; Reynolds, J.P.; Gwinn, R.; Stallings, R.L.; Henshall, D.C. Differential DNA methylation profiles of coding and non-coding genes define hippocampal sclerosis in human temporal lobe epilepsy. Brain 2015, 138, 616–631. [Google Scholar] [CrossRef] [Green Version]
  57. Miller-Delaney, S.F.; Das, S.; Sano, T.; Jimenez-Mateos, E.M.; Bryan, K.; Buckley, P.G.; Stallings, R.L.; Henshall, D.C. Differential DNA methylation patterns define status epilepticus and epileptic tolerance. J. Neurosci. 2012, 32, 1577–1588. [Google Scholar] [CrossRef] [Green Version]
  58. Williams-Karnesky, R.L.; Sandau, U.S.; Lusardi, T.A.; Lytle, N.K.; Farrell, J.M.; Pritchard, E.M.; Kaplan, D.L.; Boison, D. Epigenetic changes induced by adenosine augmentation therapy prevent epileptogenesis. J. Clin. Investig. 2013, 123, 3552–3563. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Calcium (Ca2+) imaging identified seizure-like activity in vitro. Glutamate-treated rat primary neuronal cultures showed temporal epileptogenic network activity starting after a short latency period with synchronized spiking and later development of seizure-like events (SLE; (a,b)), which progressed and became longer, showing a duration up to two minutes (c). SLEs were characterized by a sharp-onset/sharp-offset transient superimposed on low-voltage fast activity, high-amplitude rhythmic bursts, and postictal depression, thereby showing strong similarities with human stereo-EEG (sEEG) seizure. Traces from individual cells and sEEG contacts are displayed for comparison. Relative power spectrograms summarizing network activity in vitro and in vivo show typical down-chirping (i.e., linear decrease in frequency) (d). Ctrl—control (blue), Glu—glutamate (red), sEEG—intracranial stereo-electroencephalogram recording, SLE—seizure-like event. Asterisk (*) indicates significance.
Figure 1. Calcium (Ca2+) imaging identified seizure-like activity in vitro. Glutamate-treated rat primary neuronal cultures showed temporal epileptogenic network activity starting after a short latency period with synchronized spiking and later development of seizure-like events (SLE; (a,b)), which progressed and became longer, showing a duration up to two minutes (c). SLEs were characterized by a sharp-onset/sharp-offset transient superimposed on low-voltage fast activity, high-amplitude rhythmic bursts, and postictal depression, thereby showing strong similarities with human stereo-EEG (sEEG) seizure. Traces from individual cells and sEEG contacts are displayed for comparison. Relative power spectrograms summarizing network activity in vitro and in vivo show typical down-chirping (i.e., linear decrease in frequency) (d). Ctrl—control (blue), Glu—glutamate (red), sEEG—intracranial stereo-electroencephalogram recording, SLE—seizure-like event. Asterisk (*) indicates significance.
Cells 10 03004 g001
Figure 2. Neuronal morphology in epileptogenesis in vitro. (a) Representative images of individual neurons at different time points and sham or glutamate treatment. (b) Sholl analysis of 30 captured neurons per time point and treatment. The number of neurite branches or intersections was plotted against the distance from the soma. Schematic of Sholl analysis shown in the middle. Asterisks (*) indicate significance. d—day.
Figure 2. Neuronal morphology in epileptogenesis in vitro. (a) Representative images of individual neurons at different time points and sham or glutamate treatment. (b) Sholl analysis of 30 captured neurons per time point and treatment. The number of neurite branches or intersections was plotted against the distance from the soma. Schematic of Sholl analysis shown in the middle. Asterisks (*) indicate significance. d—day.
Cells 10 03004 g002
Figure 3. Relative loss of GABAergic synapses and Calbindin-immunoreactive cells. (a) Glutamatergic and GABAergic synapses were stained with Vglut1 (blue, white arrow) and Vgat (green, asterisk), respectively, for manual and automated quantification. Only synapses showing co-localization with Bsn (red) were counted. Confocal microscopy at 63×. (b) Processed immunofluorescence images for pre- and post-synapses of Fiji plug-in “Synapse Counter” and resulting co-localizations. Scale bars are 500 nm. (c) Relative quantification of GABAergic over glutamatergic synapses at 3 and 7 days after glutamate treatment as compared with time-matched sham-treated controls. (d,e) Relative amount of Calb1-positive neurons compared with Calb1-negative neurons in glutamate and sham-treated cultures after 7 days. Asterisks (*) indicate significance. ctrl—control, d—day, glu—glutamate.
Figure 3. Relative loss of GABAergic synapses and Calbindin-immunoreactive cells. (a) Glutamatergic and GABAergic synapses were stained with Vglut1 (blue, white arrow) and Vgat (green, asterisk), respectively, for manual and automated quantification. Only synapses showing co-localization with Bsn (red) were counted. Confocal microscopy at 63×. (b) Processed immunofluorescence images for pre- and post-synapses of Fiji plug-in “Synapse Counter” and resulting co-localizations. Scale bars are 500 nm. (c) Relative quantification of GABAergic over glutamatergic synapses at 3 and 7 days after glutamate treatment as compared with time-matched sham-treated controls. (d,e) Relative amount of Calb1-positive neurons compared with Calb1-negative neurons in glutamate and sham-treated cultures after 7 days. Asterisks (*) indicate significance. ctrl—control, d—day, glu—glutamate.
Cells 10 03004 g003
Figure 4. (a,b) PCA plot and heatmap showing distinct DNA methylation patterns in control (blue) and epileptogenic networks (red). Clustered heatmap of methylation percentages (ranging from 0% to 100 %) of all differentially methylated positions with an F-value of 0.05 or lower. Colors are relative, with darker colors representing higher methylation percentages. (c) Functional relationships between differentially methylated genes and their encoded proteins (5mC up = green, 5mC down red), whereby the spatial location of nodes was shifted into their corresponding subcellular locations. Proteins encoded by genes targeted by differential DNA methylation in our analysis mapped to the extracellular space, plasma membrane, cytosol, nucleus, and other. Direct and indirect interactions between DMGs and their respective encoding proteins are shown, illustrating a highly complex interaction network contributing to the epileptic phenotype in glutamate-treated cultures. A simplified overview for orientation is displayed in blue (lower left corner). 5mC—5-methyl-Cytosine, PI3K/AKT—Phosphatidylinositol—4,5—Bisphosphate 3—Kinase / AKT Serine/Threonine Kinase 1, ECM—extracellular matrix, MAPK—Mitogen-Activated Protein Kinase 1, PCA—principal component analysis.
Figure 4. (a,b) PCA plot and heatmap showing distinct DNA methylation patterns in control (blue) and epileptogenic networks (red). Clustered heatmap of methylation percentages (ranging from 0% to 100 %) of all differentially methylated positions with an F-value of 0.05 or lower. Colors are relative, with darker colors representing higher methylation percentages. (c) Functional relationships between differentially methylated genes and their encoded proteins (5mC up = green, 5mC down red), whereby the spatial location of nodes was shifted into their corresponding subcellular locations. Proteins encoded by genes targeted by differential DNA methylation in our analysis mapped to the extracellular space, plasma membrane, cytosol, nucleus, and other. Direct and indirect interactions between DMGs and their respective encoding proteins are shown, illustrating a highly complex interaction network contributing to the epileptic phenotype in glutamate-treated cultures. A simplified overview for orientation is displayed in blue (lower left corner). 5mC—5-methyl-Cytosine, PI3K/AKT—Phosphatidylinositol—4,5—Bisphosphate 3—Kinase / AKT Serine/Threonine Kinase 1, ECM—extracellular matrix, MAPK—Mitogen-Activated Protein Kinase 1, PCA—principal component analysis.
Cells 10 03004 g004
Table 1. DMGs with a role in brain development, function, and disease. 5mC—5-methyl-Cytosin/ DNA methylation, ADHD—attention deficit hyperactivity disorder, BDNF—brain derived neurotrophic factor, DMG—differentially methylated gene, ECM—extracellular matrix, EIEE—early infantile epileptic encephalopathy, ERK/MAPK—extracellular signal-regulated kinases/ mitogen—activated protein kinase, FAME—familial adult myoclonic epilepsy, GPCR—G—protein coupled receptor, IDD—intellectual developmental disorder, MCD—malformation of cortical development, mTOR—mammalian target of rapamycin.
Table 1. DMGs with a role in brain development, function, and disease. 5mC—5-methyl-Cytosin/ DNA methylation, ADHD—attention deficit hyperactivity disorder, BDNF—brain derived neurotrophic factor, DMG—differentially methylated gene, ECM—extracellular matrix, EIEE—early infantile epileptic encephalopathy, ERK/MAPK—extracellular signal-regulated kinases/ mitogen—activated protein kinase, FAME—familial adult myoclonic epilepsy, GPCR—G—protein coupled receptor, IDD—intellectual developmental disorder, MCD—malformation of cortical development, mTOR—mammalian target of rapamycin.
Gene NameSymbol5mCCompartmentFunctionDisease Association
Collagen XIV a 1Col14a1Extracellular spaceECM componentEpilepsy (FAME)
Metalloproteinase 13Mmp13ECM component, involved in neuroprotection and neurorepairAlzheimer’s disease, stroke
Neural EGFL Like 1Nell1Growth factor, control of cell growth and differentiation, nervous system developmentNeuroblastoma, craniosynostosis
Transforming Growth Factor AlphaTgfaGrowth factor, brain developmentHypothalamic hamartoma
Angiomotin Like 1Amotl1Plasma membraneControls paracellular permeability and maintains cell polarity, tight junctionsADHD
CUB and Sushi Multiple Domains 2Csmd2Synaptic transmembrane protein required for neuronal maturation, regulates the development, and maintenance of dendrites and synapsesPsychiatric disease
Catenin Alpha Like 1Ctnnal1Cadherin binding, scaffolding proteinHirschsprung disease
G Protein Subunit Beta 5Gnb5Neuronal signalingIDD, ADHD
Neuronal Membrane Glycoprotein M6BGpm6bNeuronal glycoprotein, involved in neuronal differentiation, myelination, maintenance of actin cytoskeleton, role in neuroplasticity Depression
5-Hydroxytryptamine Receptor 6Htr6GPCR, regulates cholinergic neuronal transmission in the brain, binds antidepressantsAlzheimer’s disease, schizophrenia
Potassium Voltage-Gated Channel Subfamily H Member 5Kcnh5Outward-rectifying, non-inactivating channel regulating neurotransmitter and hormone releaseEpilepsy (EIEE, Otahara Syndrome)
Semaphorin 6dSema6dAxon pathfinding
Sorbin and SH3 Domain Containing 2Sorbs2Part of the actin cytoskeleton, role in dendritic development, memory, stiffness sensing, and contractile force generationIntellectual disability
Sortilin Related VPS10 Domain Containing Receptor 2Sorcs2Neuropeptide receptor, binds precursor forms of NGF (proNGF) and BDNF (proBDNF), regulation of dendritic spine density, required for normal neurite branching and extension in response to BDNF, mediates BDNF-dependent synaptic plasticity, long-term depression, and long-term potentiationHuntington’s disease
Trafficking Protein Particle Complex 9Trappc9Activator of NF-kappa-B, role in neuronal differentiationIntellectual disability, MCD
WD Repeat Domain 1Wdr1Cofilin cofactor, involved in actin cytoskeletal dynamics, activator of NF-kappa-B, role in neuronal differentiationGlioblastoma, intellectual disability, MCD
BCL2 Apoptosis RegulatorBcl2Cytoplasm Suppresses apoptosis in neural cellsGlioma
DoublecortinDcxMicrotubule organization, neuronal migration, brain developmentMCD, epilepsy
Growth Arrest Specific 7Gas7Expressed primarily in terminally differentiated brain cells, role in neuronal development
G-protein Coupled Receptor 137cGpr137cRegulates Rag and mTORC1 localization and activity
Synaptotagmin 7Syt7Ca2+ sensor involved in exocytosis of secretory and synaptic vesicles; short-term synaptic potentiationBipolar disorder
Synaptotagmin 17Syt17Role in regulating fusion of intracellular vesicles with the plasma membrane, controls neurite outgrowth and synaptic plasticity
Thymosin Beta 10Tmsb10Important role in the organization of the cytoskeleton, binds to and sequesters actin monomers and therefore inhibits actin polymerization
WD Repeat Domain, Phosphoinositide Interacting 1Wipi1Scaffolding proteinNeurodegeneration
WW And C2 Domain Containing 1Wwc1Regulator of Hippo signaling, regulates collagen-stimulated activation of the ERK/MAPK cascade, plays a role in cognition and memory performanceIntellectual disability
AF4/FMR2 Family Member 1Aff1Nucleus Transcription factorFragile X intellectual disability, ataxia
AF4/FMR2 Family Member 2Aff2Transcription factorFragile X intellectual disability
ID Complex Subunit 4 HomologGid4Mediator complex component, is required for activation of RNA Pol II transcription by DNA-bound transcription factorsSpeech delay, intellectual disability, MCD
H2A Clustered Histone 21H2ac21Chromatin organization, core component of nucleosome
Nuclear Factor I BNfibTranscription factor, essential for proper brain developmentMacrocephaly and mental retardation
RNA Binding Motif Protein X-LinkedRbmxRegulates pre-mRNA alternative splice site selection, e.g., for Tau protein
SET Domain Containing 5Setd5Histone-Lysine Methyltransferase, mediates H3K36me3Autosomal dominant intellectual disability
Small Nuclear Ribonucleoprotein Polypeptide NSnrpnpre-mRNA processing; may contribute to tissue-specific alternative splicingPrader–Willi Syndrome
SNRPN Upstream Reading FrameSnurf Prader–Willi Syndrome, Angelman Syndrome
Sterile Alpha Motif Domain Containing 12Samd12 DNA-binding molecule of unclear functionFAME
TATA-Box Binding Protein Like 2Tbpl2 Transcription factorRolandic epilepsy
R3H Domain Containing 1CUB and Sushi Multiple Domains 2R3hdm1Csmd2 Required for neuronal maturation, regulates development and maintenance of dendrites and synapsesMild intellectual disability, psychiatric disease
Willebrand Factor A Domain-Containing Protein 3BR3H Domain Containing 1Vwa3bR3hdm1Other Role in brain development, involved in apoptotic signaling in neuronal cells. Role in neuron dendritic growth and branching Spinocerebellar ataxia with intellectual disability, mild intellectual disability
Willebrand Factor A Domain-Containing Protein 3BVwa3b Role in brain development, involved in apoptotic signaling in neuronal cellsSpinocerebellar ataxia with intellectual disability
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jablonski, J.; Hoffmann, L.; Blümcke, I.; Fejtová, A.; Uebe, S.; Ekici, A.B.; Gnatkovsky, V.; Kobow, K. Experimental Epileptogenesis in a Cell Culture Model of Primary Neurons from Rat Brain: A Temporal Multi-Scale Study. Cells 2021, 10, 3004. https://doi.org/10.3390/cells10113004

AMA Style

Jablonski J, Hoffmann L, Blümcke I, Fejtová A, Uebe S, Ekici AB, Gnatkovsky V, Kobow K. Experimental Epileptogenesis in a Cell Culture Model of Primary Neurons from Rat Brain: A Temporal Multi-Scale Study. Cells. 2021; 10(11):3004. https://doi.org/10.3390/cells10113004

Chicago/Turabian Style

Jablonski, Janos, Lucas Hoffmann, Ingmar Blümcke, Anna Fejtová, Steffen Uebe, Arif B. Ekici, Vadym Gnatkovsky, and Katja Kobow. 2021. "Experimental Epileptogenesis in a Cell Culture Model of Primary Neurons from Rat Brain: A Temporal Multi-Scale Study" Cells 10, no. 11: 3004. https://doi.org/10.3390/cells10113004

APA Style

Jablonski, J., Hoffmann, L., Blümcke, I., Fejtová, A., Uebe, S., Ekici, A. B., Gnatkovsky, V., & Kobow, K. (2021). Experimental Epileptogenesis in a Cell Culture Model of Primary Neurons from Rat Brain: A Temporal Multi-Scale Study. Cells, 10(11), 3004. https://doi.org/10.3390/cells10113004

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop