Abstract
In secondary epilepsy, a seizure-prone neural network evolves during the latent period between brain injury and the onset of spontaneous seizures. The nature of the evolution is largely unknown, and even its completeness at the onset of seizures has recently been challenged by measures of gradually decreasing intervals between subsequent seizures. Sequential calcium imaging of neuronal activity, in the pyramidal cell layer of mouse hippocampal in vitro preparations, during early post-traumatic epileptogenesis demonstrated rapid increases in the fraction of neurons that participate in interictal activity. This was followed by more gradual increases in the rate at which individual neurons join each developing seizure, the pairwise correlation of neuronal activities as a function of the distance separating the pair, and network-wide measures of functional connectivity. These data support the continued evolution of synaptic connectivity in epileptic networks beyond the latent period: early seizures occur when recurrent excitatory pathways are largely polysynaptic, while ongoing synaptic remodeling after the onset of epilepsy enhances intranetwork connectivity as well as the onset and spread of seizure activity.
Introduction
Post-traumatic epilepsy is increasingly being recognized as an important public health problem (Ferguson et al., 2010; Raymont et al., 2010). We know very little about the processes that underlie the development of epilepsy after brain injury (Lowenstein, 2009; Noebels et al., 2012; Pitkänen and Immonen, 2014), and thus cannot develop rational treatments for this particularly refractory form of secondary epilepsy (Temkin, 2009; Prince et al., 2012; Goldberg and Coulter, 2013).
Computer modeling studies (Houweling et al., 2005) and hippocampal slice culture experiments (Bausch et al., 2006) have proposed that epileptogenesis could be driven by neuronal homeostasis, the process by which neurons normalize their input in the face of experience-dependent changes in synaptic weighting (Turrigiano et al., 1998; Pozo and Goda, 2010). The resultant networks may have different architecture than the networks built during normal development. For example, small-world networks, in which a minority of neurons are connected to an outsized number of other neurons (Watts and Strogatz, 1998), are a robust means by which to engender seizures (Netoff et al., 2004; Dyhrfjeld-Johnsen et al., 2007).
Anatomical studies demonstrate loss of principal cells, interneurons, and astrogliosis in experimental and human epilepsy (Buckmaster and Dudek, 1997; Maglóczky and Freund, 2005; Clasadonte and Haydon, 2012; Huusko et al., 2015). Axon sprouting is strongly correlated with epileptogenesis, although whether sprouting is necessary for epilepsy has been questioned (Heng et al., 2013). Similarly, loss of interneurons correlates poorly with subsequent epilepsy (Huusko et al., 2015), and in at least some models inhibition recovers by the time seizures commence (Thind et al., 2010).
The time course of epileptogenesis has been quantified by studying the frequency of spontaneous seizures after brain injury induced by status epilepticus (Williams et al., 2009). Seizure probability continues to increase after the first spontaneous seizure, implying that the processes underlying epileptogenesis continue after the first spontaneous seizures. Human data support the ongoing evolution of epilepsy (Elwes et al., 1988), including progressive regional atrophy (Bernhardt et al., 2013); the gradual emergence of drug resistance (Berg et al., 2006); and recrudescence of seizures years after seizure surgery (McIntosh et al., 2004, 2012).
The development of pathological synchronization during epileptogenesis is unexplored. Spikes precede seizures, but spikes do not increase in duration to become seizures (Williams et al., 2009, White et al., 2010). The topology of network synchronization remains hypothetical. For example, seizures could coalesce from smaller scattered nodes of physiological (Ikegaya et al., 2004; Bonifazi et al., 2009) or pathological (Huberfeld et al., 2011; Sabolek et al., 2012) synchrony. Alternatively, the fraction of neurons in the network participating could increase regardless of spatial constraints.
The organotypic hippocampal slice preparation provides an experimentally accessible preparation in which to study post-traumatic epileptogenesis (McBain et al., 1989; Berdichevsky et al., 2012; Albus et al., 2013). By visualizing neural activity during epileptogenesis, we can address questions regarding network structure and the evolution of pathological synchronization. We sequentially imaged, in the same slice cultures, the calcium activity of neurons throughout the early course of epilepsy.
Materials and Methods
Organotypic slice cultures.
Organotypic slice culture protocols were approved by the Massachusetts General Hospital Subcommittee on Research Animal Care. For imaging experiments, membrane interface type organotypic slice cultures were prepared as described previously (Stoppini et al., 1991) from mice expressing Yellow Cameleon 3.6 in all Nestin-positive cells. Briefly, isolated hippocampi from P5–P8 pups of either sex were cut into 400 μm slices on a McIlwain tissue chopper (Mickle Laboratory Engineering). Slices were plated on membrane culture inserts (PICMORG50; Millipore) and placed in 50 mm coverslip glass-bottom Petri dishes (Fig. 1A; HBST-5040; WillCo Wells). For electrophysiological experiments, to facilitate transfer of sample to an electrode-accessible recording chamber, roller tube cultures were prepared, from C57BL/6 mice (patch-clamp experiments) or Sprague Dawley rats (kynurenate washout experiment) of either sex, on glass coverslips as previously described (Gähwiler, 1981; Dyhrfjeld-Johnsen et al., 2010). Cultures were incubated at 5% CO2, 36°C in 1.5 ml Neurobasal-A growth medium with 2% B27, 500 μm GlutaMAX, and 0.03 mg/ml gentamycin added (all from Invitrogen). Growth medium was changed every 3–4 d. All imaging and loose-patch recordings were made in growth medium.
Imaging epileptogenesis. A, B, To image physiological changes associated with the emergence of an epileptic network, slices were grown in an optically accessible environment and imaged in a sterile, humidified, CO2-controlled stage-mounted incubator. C, Slices were made from animals expressing calcium indicator YC3.6 in all brain cells and the CA1 pyramidal cell layer was imaged using two-photon targeted path scanning. D, Simultaneous calcium imaging and loose-patch electrophysiological recordings determined that we could fluorometrically detect single action potentials with significant spike-triggered averaging (D2), but could not significantly differentiate them from baseline on a single trial basis (D3, D4). Small bursts of action potentials, corresponding to a change in YFP/CFP ratio of 0.2, were distinguishable from baseline. The light pink trace in D1 was low-pass filtered at 5 Hz, as were all analyzed calcium signals. The dark red trace was low-pass filtered at 1 Hz and is shown for reference only. E, The orientation of the pyramidal cell layer (yellow curve) and notable features such as blood vessels (red arrowheads) were used as fiducial markers in returning to the same subregion of CA1 on consecutive days. F, G, Astrocytes were readily distinguishable from neurons based on both calcium dynamics and spatial distribution.
Electrophysiology.
Organotypic slice cultures were transferred to a submerged recording chamber and perfused (2.5 ml/min) with prewarmed (36°C) oxygenated (95% O2 and 5% CO2) ACSF containing the following (in mm): 119 NaCl, 1.25 NaH2PO4, 3 KCl, 26 NaHCO3, 11 glucose, 2.5 CaCl2, and 2 MgSO4.7H2O, pH 7.4, for 30 min before electrophysiological recording. Recording electrodes were pulled from borosilicate glass capillaries (Sutter Instruments) using a micropipette puller (Model P-97; Sutter Instruments) with resistance 4–7 MΩ when filled with internal solution containing the following (in mm): 120 Cs-Gluconate, 5 MgCl2, 0.6 EGTA, 30 HEPES, 4 Mg-ATP, 0.4 Na2-GTP, 10 phosphocreatine-Tris, and 5 QX-314; 290 mOsm and pH was adjusted at 7.2 with CSOH. Dual whole-cell patch-clamp recordings were made from a pair of CA1 pyramidal neurons with visualization using infrared differential interference contrast camera on an upright Nikon Eclipse FN1 microscope. Signal acquisition was performed using a Multiclamp amplifier (Multiclamp 700B; Molecular Devices) with Clampex 10 software (Molecular Devices). Signals were sampled at 10 kHz and filtered at 2 kHz. Data were stored on a PC for off-line analysis after digitization by A/D converter (Digidata 1440A; Molecular Devices). The resting membrane potential of neurons was measured in current-clamp mode within 5 s of establishing whole-cell recordings to minimize the effect of internal solution on intrinsic properties of neurons. Membrane potential was corrected for a liquid junction potential of 13.7 mV. Neurons with resting membrane potential less than −60 mV were discarded. Neuronal input and access resistances were repetitively monitored during the experiment. Cells with changes of input and access resistance changes of >20% were excluded from analysis.
Spontaneous EPSCs (sEPSCs) were collected at a holding potential of −70 mV, which is close to the calculated (in whole-cell configuration, EGABA is heavily influenced by electrode Cl− concentration) reverse potential of GABA, which is −68 mV. sIPSCs were also measured in the same neuron by raising the holding potential to 0–10 mV, a potential close to the reverse potential of excitatory input, to minimize the amplitude of concurrent sEPSCs. mEPSCs and mIPSCs were recorded from the same neurons after TTX (1 μm) application for 30 min. Postsynaptic currents were detected using custom software (DClamp; www.ieeg.org/?q=node/34).
Loose-patch recordings were performed by sealing a patch electrode to the cell membrane and current clamping with zero current (Fig. 1D). To identify the recorded cell for two-photon imaging, and to verify that the membrane was not ruptured, electrodes were filled with ACSF containing Alexa Fluor 594.
Imaging in culture.
All imaging was performed in a microscope stage-mounted miniature incubator (Fig. 1B; TC-MIS-30; Bioscience Tools) set to 37°C, and equipped with a humidified, sterilely filtered, CO2 controller set to 5% (CO2-MI, CO2-500ML; Bioscience Tools). Slices were imaged through the coverslip glass-bottomed Petri dish using an objective inverter (LSM Tech) to effectively use an upright two-photon microscope as an inverted microscope. To maintain sterility, the Petri dishes remained covered throughout the imaging sessions and were handled with sterile gloves. A mark was placed on the dish to ensure that the slice was oriented in the same direction on consecutive days. The objective was positioned at a depth of 50 μm up from the bottom (membrane side) of the slice. The same region of the CA1 pyramidal cell layer was imaged over several days by manually positioning the slice such that the angle of stratum pyramidale and (when possible) any landmarks (e.g., blood vessels) matched those in the image acquired on the previous day (Fig. 1E).
Two-photon images were acquired using custom-designed software and the scan head from a Radiance 2000MP (Bio-Rad), equipped with a 20× 0.95 NA water-immersion objective (Olympus), and photomultiplier tubes with appropriate filters for YFP (545/30) and CFP (450/80). A Spectra-Physics Mai Tai laser (Newport), set to 860 nm, was used for two-photon excitation. All images of calcium dynamics were acquired using targeted path scanning (TPS; Lillis et al., 2008). Astrocytes were excluded from analysis based on three characteristics. (1) Patterns of activity show that cells with suprathreshold events that had a >2 s and 10–90% rise time were excluded (Fig. 1E; Dani et al., 1992). (2) Morphology describes astrocytes with a “spongiform” appearance with highly ramified processes (Bushong et al., 2004). (3) Spatial distribution demonstrates that astrocytic soma tended to be primarily located outside the stratum pyramidale (Fig. 1F). To compute “calcium” traces, YFP and CFP signals were first low-pass filtered at 5 Hz. Next, a first-order polynomial was fit to the YFP/CFP ratio, using least-squares minimization and was subtracted from the YFP/CFP ratio to correct for the different rates of bleaching in YFP and CFP.
Statistics.
All p values given were calculated using Pearson's correlation coefficient to test for significant changes in measured parameters over time. All values are presented as mean ± SE in text and figures.
Epileptiform event detection.
Epileptiform events were identified using the mean, baseline-subtracted YFP/CFP ratio: mean ΔR. Peaks in mean ΔR represent large amplitude and/or highly synchronous activity in the recorded cells and approximately match a field potential recording (Lillis et al., 2012). Excursions in mean ΔR exceeding 0.15 were deemed epileptiform events. The event peak was determined by finding the maximum of mean ΔR between the initial threshold crossing and ending after 5 s or when mean ΔR drops below 0.15, whichever came first. Criteria for distinguishing interictal spikes from seizures in animal models of epilepsy are only recently beginning to be developed (Fisher et al., 2014) and often cannot be described by a single parameter. In the data presented here, there were two distinct types of events: brief calcium (mean ΔR) transients separated by >1 s and tonic clonic-like seizures in which calcium remained elevated for several seconds, then resolved with synchronous bursting (Figs. 3⇓–5). These events were well separated using duration alone: those lasting longer than 10 s were categorized as seizure-like events, while all other events were interictal epileptiform discharges. Although these are not clinical seizures as the terms imply, for the sake of being succinct we henceforth refer to these events as “seizures” and “interictal bursts.” For the purposes of characterizing interictal burst and seizure onset dynamics, individual cellular threshold crossings (threshold: Δr = 0.2; Fig. 1D) were used to calculate an onset time for each cell relative to the interictal burst/seizure peak.
Network analysis.
Zero-lag correlation was used as a measure of functional connectivity in the network analyses presented here. For all cross-correlation computations, data were first standardized by subtracting the mean and dividing by the SD. Correlation coefficients were computed between all pairs standardized calcium traces. To determine significance, correlation analyses were performed on 200 sets of randomized surrogate calcium traces with amplitude spectra matching the recorded data. Surrogate time series were generated using the iterative amplitude adjusted Fourier transform algorithm (Schreiber and Schmitz, 1996). For each pair of neurons, a distribution of correlation strengths was computed using the 200 pairs of surrogate traces and the 99th percentile was chosen as a threshold for significance.
The MatlabBGL package (Gleich, 2009) was used to compute the graph analytical statistics mean path length and betweenness centrality (graphically defined in Fig. 7). For these analyses, each neuron was treated as a node on the graph. Mean path length was computed using the inverse of correlation strength as the “distance” between two nodes. Betweenness centrality was computed as the average fraction of shortest paths on which a given node falls (i.e., the degree to which a neuron acts as a “hub” between other neurons; Freeman, 1977). Windowed statistics were calculated using a 500 sample sliding window, shifted in 50 sample steps, corresponding to ∼13 s windows with a 1.3 s step size. Sample period is inherently recording specific with TPS recordings, but did not significantly change with time in culture (p = 0.24). Thus, seizure-triggered correlation averages were computed across slices using correlation windows as a time base and time axes for Figure 7 are approximate.
Results
The natural history of epileptogenesis in hippocampal organotypic slice cultures has been characterized using both acutely sampled and chronic field potential recordings (Dyhrfjeld-Johnsen et al., 2010; Berdichevsky et al., 2012). Here we describe changes in the spatial organization of epileptiform activity in the CA1 pyramidal cell layer as interictal bursts and seizures emerge.
We first conducted electrophysiological recordings to test whether epilepsy evolves in this preparation as a consequence of alterations in intrinsic neuronal properties or imbalances in the rate of excitatory versus inhibitory synaptogenesis. Dual whole-cell patch-clamp recordings of spontaneous and miniature postsynaptic currents were made from CA1 pyramidal neurons at 1, 5, and 8 DIV. mPSCs were recorded during the blockade of spiking activity and represent the underlying network structure, while sPSCs are driven by action potential-dependent synaptic activity and are thus a measure of spontaneous network dynamics (which presumably, at least partially result from network structure). Any epileptiform discharges present during recording of sPSCs were excluded from analysis.
Electrophysiological properties measured are summarized in Table 1. Briefly, the frequency of recorded excitatory and inhibitory spontaneous (Fig. 2C; sEPSC: p = 2.4 × 10−35, n = 193 cells; sIPSC: p = 1.95 × 10−14, n = 153 cells) and miniature (Fig. 2D; mEPSC: p = 0.0015, n = 78 cells; mIPSC: p = 1.05 × 10−5, n = 81) postsynaptic currents increased over the first 8 DIV. The sEPSC/sIPSC amplitude ratio did not significantly change (DIV 1, 0.96 ± 0.06; DIV 5, 1.00 ± 0.11; DIV 8, 1.14 ± 0.06; p = 0.154; n = 137), while the sEPSC/sIPSC frequency ratio increased (DIV 1, 0.87 ± 0.14; DIV 5, 1.19 ± 0.09; DIV 8, 1.84 ± 0.16; p = 1.35 × 10−6; n = 137), reflecting changes in network activity. However, neither the mEPSC/mIPSC amplitude ratio (DIV 1, 0.64 ± 0.06; DIV 5, 1.47 ± 0.17; DIV 8, 1.03 ± 0.09; p = 0.3617; n = 59 cells) nor the mEPSC/mIPSC frequency ratio (DIV 1, 0.31 ± 0.08; DIV 5, 0.89 ± 0.34; DIV 8, 0.62 ± 0.33; p = 0.2848; n = 59 cells) significantly changed over time. The resting membrane potential also did not change with time in culture (DIV 1, 69.7 ± 4.4 mV; DIV 5, 71.1 ± 5.3 mV; DIV 8, 69.8 ± 4.3 mV; p = 0.75, n = 147 cells). However, the input resistance significantly dropped (from DIV 1, 205 ± 71 MΩ to DIV 5, 159 ± 52 MΩ and DIV 8, 148 ± 43 MΩ; p = 1.3 × 10−5, n = 148 cells). Together, these electrophysiological results are consistent with extensive synaptogenesis and axonal outgrowth seen during normal development (Gomez-Di Cesare et al., 1997; Tyzio et al., 2003), but steady-state membrane properties were not commensurate with changes in intrinsic excitability and miniature postsynaptic currents did not reflect an imbalance of inhibitory versus excitatory synaptogenesis.
Summary of electrophysiological changes during epileptogenesis
Dual whole-cell patch recordings. A, B, Pairs of CA1 pyramidal neurons were recorded in voltage clamp to measure both spontaneous (excluding epileptiform discharges) and miniature synaptic currents. C, D, Spontaneous and miniature excitatory and IPSC frequency increased throughout the first 8 DIV. E, The correlation strength of (nonictal) spontaneous excitatory and inhibitory activity significantly increased as epileptogenesis progressed. F, The trend toward increased mEPSC correlation did not reach significance, but mIPSC correlation slightly increased by DIV 8. (**p < 0.01, ***p < 0.001).
Epilepsy could also develop as a consequence of excessive recurrent synaptic connectivity. To begin to assess this, we measured the degree of correlation of synaptic activity in paired whole-cell recordings. The cross-correlation of spontaneous and miniature EPSCs between the two recorded neurons was computed for epochs of time containing no epileptiform activity. There was a steady increase in sEPSC (p = 1.54 × 10−7, n = 78 pairs) and sIPSC (p = 9.46 × 10−6, n = 56 pairs) correlation (Fig. 2A,E), while the mEPSC correlation did not significantly change (p = 0.15, n = 23 pairs) and mIPSC correlation increased slightly (Fig. 2B,F; p = 0.0017, n = 29 pairs).
Imaging network changes in epileptogenesis
To track network correlation changes on a larger scale, we developed a method for sequential sterile imaging of cellular calcium dynamics in hippocampal organotypic slice cultures prepared from mice expressing the genetically encoded calcium indicator Yellow Cameleon 3.6 (YC3.6) in Nestin-expressing cells (Fig. 1A–C). Two-photon TPS imaging (Lillis et al., 2008) was used to image YC3.6-expressing cells in CA1 stratum pyramidale. Cell-attached patch recordings revealed that single action potentials could be fluorometrically detected with spike-triggered averaging (Fig. 1D2), but could not be reliably distinguished from noise on a single trial basis (Fig. 1D3,D4). Brief bursts of action potentials resulted in a change in the YFP/CFP ratio of >0.2 that was discernible in single trials (Fig. 1D1, light pink trace; t = 173). This threshold, Δr = 0.2, was used for displaying raster plot data in Figures 3⇓–5, computing the relative amount of nonictal activity (Fig. 4D), and computing “cell onset times” (Figs. 4, 5), but correlation analyses were performed on unthresholded data. As described below, we first imaged with hourly resolution the advent and progression of interictal bursting (Fig. 3). Next we characterized, with daily recordings and correlation analysis, the development of seizures and the evolution of seizure onset dynamics (Figs. 4⇓⇓–7).
Imaging the emergence of interictal bursts. To examine the earliest emergence of synchronous activity, three slices were imaged hourly for the first 24 h in vitro (n = 3 slices, 72 recordings, 4912 detected bursts). A, Raster plots for interictal bursts highlighted in pink box are shown in the second column. The first signs of correlated firing appeared 2–3 h after slicing, when subsets of the imaged cells began to exhibit weakly correlated activity (top two rows, left two columns). By the end of the DIV 1, cells were recruited more efficiently into the population burst (bottom two rows, left two columns). This apparent buildup of activity was quantified by computing cell onset times relative to interictal burst peak. Burst peak times were calculated using the mean calcium trace, while individual cellular calcium threshold crossings (threshold: YFP/CFP ΔR = 0.2) were used to calculate an onset time for each cell relative to the burst peak. The increasing slopes and rightward shifts in the cumulative distribution of onset times (third and fourth columns) demonstrate that population bursts became more rapid and synchronous. B, Burst onset time, defined as 80th to 20th percentile of CDFs, steadily decreased.
Interictal activity
To examine the earliest emergence of synchronous activity, three slices were imaged hourly for the first 24 h in vitro (n = 3 slices, 72 recordings, 4912 detected bursts). The first signs of correlated firing appeared 2–3 h after slicing, when subsets of the imaged cells began to exhibit weakly correlated bouts of activity (Fig. 3A), which were termed interictal bursts. Interictal burst peak times were calculated using the mean, baseline-subtracted YFP/CFP fluorescence ratios, computed across all recorded neurons (Fig. 3, left column, mean ΔR). Individual cellular calcium threshold crossings (threshold: YFP/CFP Δr = 0.2; Fig. 1D) were used to calculate an onset time for each cell relative to the burst peak. The increasing slope and rightward shift in the cumulative distribution of onset times suggests that bursts become more synchronous over the first 24 h in culture (more rapid onsets in Fig. 3A, right column, B; p = 2.4 × 10−29; n = 3 slices, 72 recordings, 4912 detected bursts).
Ictal activity
To study the emergence and evolution of seizures, daily recordings were made during the first week in vitro from 30 cells in CA1 stratum pyramidale (n = 5 slices, 40 recordings, 113 detected seizures). Tonic clonic-like seizures were apparent in the recorded calcium traces as early as 1 DIV (Fig. 4A,B). We quantified changes in “functional connectivity” of the network by computing the cross-correlation between calcium traces for all pairs of neurons. Plots of correlation coefficient versus physical distance between cells reveal predominantly local high correlations early in epileptogenesis (negative slope of correlation vs distance lines; Fig. 4C). On later DIV, correlation coefficients were more uniformly high, nearly independent of distance, as evidenced by the correlation versus distance slope approaching zero (Fig. 4C,D).
Evolution of ictal activity. Representative recordings from different phases of epileptogenesis are shown as raster plots of the 30 imaged cells (A) and mean change in YFP/CFP ratio (B; i.e., mean of 30 cells shown in A). Tonic clonic-like seizure activity was apparent as early as 1 DIV. There was substantial uncorrelated activity between seizures on DIV 0–4. After DIV 4, dynamics more closely resembled status epilepticus, with essentially no “baseline activity” between postictal depression and the subsequent seizures. C, To examine the spatial dependence of functional connectivity, correlation (corr) strength for each pair of calcium traces (pink dots) was plotted as a function of the physical distance between the cells. Pairwise correlations (pink dots) were averaged in 50 μm bins (black stars) and fit with a line (solid green). Early in epileptogenesis, the strongest correlations were predominantly local. At later DIV, correlation strength was uniformly high and nearly independent of distance between cells. D, The slope of the correlation versus distance plots shown in C approached zero as epileptogenesis progresses. E, The fraction of time during which there was detectable neural activity decreased over the first 8 DIV, corresponding to reduced neuronal activity between seizures.
Raster plots of all 30 calcium traces (Figs. 4A, 5A) reveal both synchronous activation and uncorrelated activity between seizures on 0–4 DIV. After 4 DIV, dynamics more closely resembled status epilepticus, with essentially no “baseline activity” between postictal depression and the subsequent seizures. This decrease in detectable (ΔR > 0.2) interictal activity (Fig. 4E; p = 5.8531 × 10−6; n = 5 slices, 40 recordings, 113 detected seizures) contributed to the increase in correlation shown in Figure 4, C and D, because at later dates the correlation reflected a larger fraction of highly synchronous ictal activity versus less synchronous interictal activity.
Seizure onset becomes more instantaneous. A, Mean ΔR (left column) and cellular calcium traces (second column) were used to quantify the rate at which neurons joined seizures. Seizure peak times were calculated as described for burst peaks in Figure 3. CDFs of the onset times of individual neurons relative to the seizure peak (third and fourth column) became steeper, indicating a more rapid and step-like seizure onset. B, A plot of seizure onset time, defined as 80th to 20th percentile from the CDF data, reveals the trend that seizure onset became more instantaneous at later DIV.
In addition to the increased long-range correlation and decreased interictal activity, seizure onset became more instantaneous at later DIV. On early DIV, there was a gradual increase in pre-ictal activity (Figs. 4A,B, 5A). This apparent buildup of activity was quantified, as with burst activation, by computing an onset time, relative to seizure peak, for each cell. Seizure peak times and cell onset times were calculated from the mean ΔR and individual neuronal calcium traces, respectively (Fig. 5A, left two columns). Cumulative distribution functions (CDFs) of neuronal onset times relative to the seizure peak reveal that seizure onset became more rapid, synchronous, and step-like over the first 8 d in culture (Fig. 5A, right two columns). Seizure onset time, calculated as the time difference between the 80th and 20th percentiles of the CDF, revealed a trend toward more rapid seizure onset as epileptogenesis progressed (Fig. 5B; p = 1.03 × 10−16; n = 5 slices, 40 recordings, 113 detected seizures). Notably, these changes did not drastically alter the tonic clonic-like structure of the seizures (Figs. 4, 5) or the mean seizure duration (p = 0.3216).
Dynamic changes in functional connectivity
In this study, we used the correlation between cellular calcium activities as a measure of functional network connectivity. Since functional network structure depends on the measured activity, it inherently changes with the dynamic state of the network (Feldt et al., 2010). We quantified these transitions in functional network state using a 500 sample windowed, zero-lag cross-correlation between all pairs of neuronal calcium traces.
During the early phase of epileptogenesis, the gradual buildup of synchronous activity preceding seizure onset is reflected in the steady increase in mean correlation (Fig. 6A,B). Phase-randomized surrogates were generated for each pair of calcium segments and used (as described in Materials and Methods) to determine a threshold for significant correlations (Fig. 6C). This information was used to generate graphs of network connectivity, wherein each node represents a neuron and vertices are drawn between significantly correlated pairs (Fig. 6D).
Dynamic changes in functional connectivity. Windowed correlation analyses were used to quantify changes in functional network connectivity that correspond to different modes of activity. A, In early epileptogenesis (DIV 1), there is a slow pre-ictal buildup of neuronal calcium activity apparent in plots of mean ΔR. B, There is a corresponding increase in mean correlation. C, Significant correlation (corr) thresholds were determined using phase-randomized surrogates. Histograms represent the distribution of cell pairs' correlation minus the 99th percentile of phase-randomized surrogate correlation for each pair (values >0 represent significant correlations). D, When correlation exceeded threshold, vertices were drawn between the pair of neurons (nodes) on the graph of “functional connectivity.” E, The mean path length (PL; where “distance” of a vertex = 1/correlation) between all pairs of neurons steadily decreased during the pre-ictal buildup of activity. F, Maximum betweenness centrality (BC), a measure of small-world network structure, peaked during the times corresponding to transitions to (pink) and from (green) seizure activity, at which time levels of functional connectivity were moderate, but remained low when functional connectivity was low (yellow) or uniform (blue). G, H, Log-log histograms of degree, the number significant correlations per cell, for the sample shown in A–F (G) and for all samples (H) reveal an approximately power law distribution during times of peak betweenness centrality, suggesting transient scale-free network topologies.
Mean path length and betweenness centrality are higher order measures of connectivity that were computed to quantify changes in functional network structure. Mean path length is calculated by finding the shortest path (where distance = 1/correlation coefficient) between all pairs of neurons (Fig. 7B). During the pre-ictal increase in mean correlation, mean path length steadily dropped (Fig. 6E). Given its intimate connection to correlation (Newman, 2001; Opsahl et al., 2010), this finding is not unexpected. Betweenness centrality is calculated as the number of shortest paths on which each neuron falls (i.e., the extent to which each neuron acts as a “hub” in the network; Fig. 7C). Since the presence of strongly hub-like nodes is one of the defining characteristics of a small-world network, maximum betweenness centrality (MBC) represents an index of the extent to which the network can be described as having small-world structure (Fig. 6F). When there was little synchronous activity in the network (Fig. 6, yellow box), connectivity was low, thus there were no hub-like neurons, and MBC was low. In the middle of a seizure, connectivity was essentially all to all; thus, the shortest path between two nodes was a direct connection (i.e., there were no opportunities for hubs between directly connected neurons) and MBC was low (Fig. 6, blue box). In the transition between quiescence and seizure, when there was sparse synchronous activity, MBC peaked and small-world network structure transiently emerged (Fig. 6, pink and green boxes).
Correlation, mean path length (PL), and betweenness centrality (BC) at seizure onset. For each DIV, seizure-triggered averages were computed of network correlation (A), mean path length (B), and betweenness centrality (C). The slow buildup in mean correlation at DIV 1 transforms into a step-like increase by DIV 8. The correlation increase results in a congruent decrease in mean path length. On DIV 1–3, when there was significant pre-ictal activity, maximum betweenness centrality peaked. In late epileptogenesis, functional connectivity abruptly transitioned from low to uniformly high, so that maximum betweenness centrality remained low. Time 0 corresponds to seizure peak. Seizure-triggered averages were computed using the 10 correlation windows preceding seizure peak. Time axis is slightly variable for reasons described in Materials and Methods, but shown here in seconds for reference. Diagrams on the right describe how each network statistic was calculated.
Another way to test for small-world network structure is to plot the degree distribution (i.e., histogram of number of significant correlations per cell). A scale-free network, which is an “ultrasmall” world network (Cohen and Havlin, 2003), has a power law distribution. While the sampling size is relatively small in this study, the peaks in betweenness centrality correspond to a linear trend in log-log plots of degree distribution (Fig. 6G, pink and green boxes), supporting the idea that the network is transiently small world and perhaps also scale free. To characterize this transition across slices, seizures were analyzed from each slice on at least 2 d in which there was a pre-ictal buildup of activity (i.e., not an abrupt seizure onset). The connectivity patterns observed in the population data (Fig. 6H) mirror those observed in the illustrative example shown in Fig. 6G, suggesting that the state-dependent fluctuations in functional network architecture is not unique to the representative slice shown.
To quantify the evolution of functional network state transitions shown in Figure 6, seizure-triggered averages of correlation, mean path length, and betweenness centrality were computed with data from five slices followed over 8 DIV. During the first 3 d of seizures, there is a steady increase in correlation and decrease in mean path length (Fig. 7A,B). By 4 DIV there is a nearly step-like increase in correlation and decrease in mean path length at seizure onset. In early epileptogenesis, during the slow pre-ictal buildup of activity, maximum betweenness centrality peaked (Fig. 7C). Once seizures evolve into a more instantaneous onset, the network abruptly transitions from low connectivity to all-to-all connectivity; thus, maximum betweenness centrality remained low (Fig. 7C).
Discussion
Serial calcium imaging of network activation during early epileptogenesis revealed a progressive increase in the rate of synchronization that continuously evolves from the earliest interictal clustering of activity through the development of the first seizures continuing through the development of electrographic status epilepticus. Dual whole-cell recordings demonstrated that this evolution is associated with increases in the rates of spontaneous and action-potential independent transmitter release, but not with proconvulsant alterations in membrane properties or with the balance of excitatory versus inhibitory synaptogenesis, as estimated from the rates of miniature synaptic currents. Electrophysiological and calcium imaging measures of interneuronal correlation increased throughout epileptogenesis, as did estimates of network connectivity based on these correlations. Thus we conclude that progressive enhancements of recurrent synaptic circuitry are the most likely substrate for epileptogenesis, although the pattern of this increased connectivity remains to be determined.
What are the critical elements of epileptogenesis?
Our prior studies indicate that in this model the glutamate receptor antagonist kynurenate blocked seizure activity (Berdichevsky et al., 2012), but on washout the slice immediately began to seize repeatedly at a rate fivefold higher than untreated controls, suggesting that synaptic reorganization proceeded in the absence of synaptic activity (consistent with studies of homeostatic synaptic plasticity; Turrigiano et al., 1998). Inhibition of the mTOR pathway, which is critical for synaptogenesis and long-term potentiation of existing synapses (Tang et al., 2002; Dwyer and Duman, 2013), is the only treatment demonstrated (thus far) to inhibit epileptogenesis in this model (Berdichevsky et al., 2013).
Both inhibitory and excitatory synaptic activity levels and the strength of network coupling progressively increased during epileptogenesis. Prior studies have demonstrated increases in recurrent glutamatergic synaptic activity from 1 to 3% in acute slices (Miles and Wong, 1986) to 33% in the hippocampal organotypic slice preparation (Debanne et al., 1995; Pavlidis and Madison, 1999). Prolific and ectopic anatomical sprouting (Gutiérrez and Heinemann, 1999; De Simoni et al., 2003) and dendritic spine formation and stabilization (Dailey and Smith, 1996) have been demonstrated in organotypic slices over the time frame studied here. The most parsimonious explanation of these data is that increases in the number of recurrent excitatory synapses between pyramidal cells subserve the increases in network connectivity, seizure probability, and the rate of ictal network activation in this model. This mechanism is also supported by the reduction in activity between epileptiform events (Fig. 4E): as recurrent excitation is increased, any activity in the network is more likely to lead to pathologically high levels of network activation. This mechanism supports the neuronal homeostasis hypothesis of replacement of lost afferent and efferent connections with recurrent synapses contributes to epileptogenesis after brain injury (Houweling et al., 2005).
Evolution of the epileptic network
The incidence of epilepsy and rate of epileptogenesis following brain trauma in human or experimental epilepsy are correlated with a number of factors including developmental stage, nature of injury, and size of animal. In pediatric patients recovering from brain trauma, synaptic remodeling is driven by both normal brain development and repair of injured brain tissue (Johnston, 2009). Similarly, in organotypic slice cultures prepared from a wide age range of juvenile rodent pups, previous studies have reported increases in synaptic connectivity (Debanne et al., 1995; Dailey and Smith, 1996; Pavlidis and Madison, 1999; De Simoni et al., 2003); stabilization of initially dynamic, transient filopodial protrusions (Dailey and Smith, 1996) into dendritic spines; and epileptogenesis (Albus et al., 2013). Whether the same homeostatic mechanisms guide developmental synapse formation and recovery (or epileptogenesis) following injury remains an important, unanswered question (Avanzini et al., 2014).
If, as our data suggest, recurrent synaptogenesis underlies epileptogenesis, then a key determinant of the rate of epileptogenesis following injury is the fraction of recurrent versus nonrecurrent targets to which the injured neurons can reconnect. If neurons homeostatically scale their synapses to restore a target firing rate (Turrigiano et al., 1998), having a smaller pool of reachable synaptic targets would naturally lead to a more recurrently connected network. As a first estimate, Figure 8 shows that the rate of epileptogenesis is approximately logarithmically related to brain volume, which is used here as an estimate of nonrecurrent synaptic targets.
Latent period of post-traumatic epileptogenesis increases with brain volume. This descriptive plot demonstrates that in both experimental animal models and human epilepsy the rate of epileptogenesis is inversely proportional to the volume of brain tissue involved. Latent period data are described previously (Annegers et al., 1998; Williams et al., 2009; Glass et al., 2011; Heinrich et al., 2011; Berdichevsky et al., 2012; Steinmetz et al., 2013). Volume data from (Tajima et al., 1993; Rengachary, 2005; Badea et al., 2007; Knickmeyer et al., 2008; Kavoi and Jameela, 2011).
The rate at which neuronal activity synchronized at the onset of a seizure increased throughout early epileptogenesis (Figs. 3, 4, 5, 7). This corresponds to an increase in long-distance correlation (Fig. 4) and to the correlation of spontaneous EPSCs in dual whole-cell recordings (Fig. 2). We interpret these data to indicate that ongoing synaptogenesis, driven either developmentally or by recovery after trauma, led to progressive shortening in excitatory polysynaptic pathways between neurons, making them shorter and more direct. A prediction of such an increase in connectivity might be that seizure probability would increase with increased excitatory connectivity in the network, but this was not observed (p = 0.0856). One possible explanation for this is that the interseizure interval was already as short as the typical postictal period of depression even early in epileptogenesis, so that we were not able to resolve increases in seizure probability in this model.
Evolution of ictal network dynamics
Functional connectivity in neuronal networks is quite dynamic at both the macroscopic (Hutchison et al., 2013) and microscopic (Feldt et al., 2010) levels. Although prior research into functional network connectivity has focused predominantly on excitatory connections between neurons (Watts and Strogatz, 1998; Netoff et al., 2004; Perin et al., 2011; Markov et al., 2013), we hypothesize that changes in functional connectivity occur, at least in part, because interneurons can mask (Snyder et al., 2015) or enhance (Hamilton et al., 2013) functional connectivity among populations of principal neurons. Recent evidence suggests that there is a pre-ictal failure of inhibition (Trevelyan et al., 2006; Huberfeld et al., 2011; Lillis et al., 2012; Zhang et al., 2012), resulting in a predominantly excitatory functional network. Using windowed network analysis, we found small-world network structure during this time window in the presumed exclusively principal neuron population (Fig. 6). Furthermore, functional connectivity dynamics (i.e., the trajectory of functional network architecture as the network cycled through ictal, postictal, and interictal states) changed throughout epileptogenesis (Fig. 7). It will be interesting to see if similar measures of network connectivity dynamics can be extracted from clinical EEG data, and whether such analyses will provide useful information regarding the stage or severity of clinical epilepsy.
Clinical implications
Activation of epileptic networks
In this in vitro preparation, the rate of epileptic network activation correlated strongly with measures of network connectivity and inversely with the amount of relatively normal interictal, nonepileptic activity (Fig. 4E). One implication of these findings is that the rate of ictal activation might be useful as a gauge of seizure probability. Severity of brain injury correlates with the incidence and severity of the subsequent epilepsy (Ferguson et al., 2010; Raymont et al., 2010). The current data support the idea that these correlations may arise from loss of afferent and efferent connections and the consequent degree of recurrent connectivity in the epileptic focus. A further extrapolation might be that MRI measures of axonal sheer injury (Skandsen et al., 2010), as a biomarker of deafferentation, would correlate with later epilepsy; this has not been studied to date.
A variety of patterns of ictal onset have been demonstrated in human seizures (Geiger and Harner, 1978; Fisher et al., 1992; Yoo et al., 2014). Interhemispheric seizure propagation time has been shown to predict surgical outcome, with slow propagation corresponding to better surgical outcome (Lieb et al., 1986; Weinand et al., 1992). Similarly, the dynamics of ictal correlations determine whether seizures spontaneously terminate or continue, as in status epilepticus (Kramer et al., 2012). In line with these findings in human electrophysiology, our data demonstrate, at a cellular level, a continual decrease in seizure propagation time (Fig. 5) and increase in long-distance correlation (Fig. 4C,D) as epileptogenesis progresses. These data support a recently proposed computational model, which suggests that more gradual ictal onsets would be associated with lower seizure probability (Jirsa et al., 2014). This might be best addressed by an in vivo experimental study testing whether the rate of ictal activation corresponds to seizure probability.
If replacement of lost afferent and efferent connections with recurrent synapses is sufficient to induce epileptogenesis after brain injury, then inhibition of recurrent synaptogenesis should prevent this mechanism of epileptogenesis. Nonspecific inhibition of synaptogenesis with mTOR antagonists was an effective antiepileptogenic strategy in this model (Berdichevsky et al., 2013), but this strategy is likely to also inhibit functional recovery after brain injury (Fawcett et al., 2012). A more selective antiepileptic strategy would exploit the intense activation of the epileptic network, for example, by activity-dependent weakening of recurrent synapses that are activated by interictal activity (Bains et al., 1999). Although initial efforts in this regard have been unsuccessful (Hellier et al., 2009), continued efforts to identify more effective strategies of activity-dependent synaptic (or even neuronal) pruning thus might realize an antiepileptic strategy that spares functional recovery.
Limitations
The organotypic hippocampal slice preparation is a model of epileptogenesis driven by very severe brain injury, resulting in accelerated epileptogenesis. In the experiments described here, epileptogenesis was likely further accelerated, because seizure activity began at 1–2 DIV, compared with 5–9 DIV in prior studies (Berdichevsky et al., 2012, 2013). Phototoxicity associated with high-intensity multiphoton excitation of fluorescence, perturbations due to transferring slices between the CO2 incubator and the microscope stage, and transient hypocapnia may have contributed to the more rapid development of seizure activity (Albus et al., 2013). The cortical dimensions are also quite different from in vivo, with one dimension very foreshortened; this is likely to alter the pattern of recurrent synaptogenesis. Additionally, although for simplicity we have used the terms “interictal burst” and “seizure” to describe the epileptiform activity in slice cultures, slices do not have seizures in the clinical sense. Thus, while they are an experimentally accessible tool for studying neural circuit reorganization in post-traumatic epileptogenesis, organotypic slice cultures do not recapitulate all factors that may contribute to clinical epilepsy (e.g., hormonal variations, involvement of other organ systems, comorbidities, etc.).
The fluorophore used here reports primarily bursts of action potentials rather than single action potentials, and we focused on somatic calcium transients rather than dendritic activity. Thus the analyzed activity was high amplitude. The organotypic slice is prepared from young animals, so that epileptogenesis and normal developmental alterations in neuronal circuitry overlap; however, this is also the case in pediatric secondary epilepsy. Finally, the number of neurons that could be assayed was limited by the mirror servo speeds, and we were not able to identify the same neurons in sequential imaging sessions. This limited the resolution of network structure and the evolution of the role of individual neurons in spike initiation and ictogenesis.
Future studies
Our calcium imaging reflected relatively intense neuronal activity in the CA1 pyramidal cell layer, and thus was likely to miss the earliest elements of network activation. More sensitive fluorophores, larger population samples, cell type-specific reporter expression, and analysis of dendritic and axonal activity should provide additional information regarding the patterns of epileptiform network activation. Recently developed voltage-sensitive fluorescent proteins (VSFPs) hold promise for recording subthreshold activity (Cao et al., 2013; Hochbaum et al., 2014; Zou et al., 2014). This would enable experiments wherein a single neuron is stimulated and responding postsynaptic cells are mapped, rather than relying on calcium correlation-based measures of functional connectivity. The ability of high-quality VSFPs to image hyperpolarization would also facilitate mapping of functional connectivity in interneuron networks, which anatomical evidence suggests may have small-world or scale-free architecture (Buzsáki et al., 2004; Soltesz, 2006). Experiments could be further enhanced by performing chronic imaging in situ, inside the primary water-jacketed CO2 incubator, to reduce disturbance of slice cultures and enable tracking of individual neurons over many days.
Footnotes
This work was supported by the National Institutes of Health grants NS034700, NS074772, and NS077908 (to K.S.), and EB000768 (to B.B.).
The authors declare no competing financial interests.
- Correspondence should be addressed to Kevin J Staley, 114 16th Street, #2600, Charlestown, MA 02129. staley.kevin{at}mgh.harvard.edu
