Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Blog
    • Collections
    • Podcast
  • TOPICS
    • Cognition and Behavior
    • Development
    • Disorders of the Nervous System
    • History, Teaching and Public Awareness
    • Integrative Systems
    • Neuronal Excitability
    • Novel Tools and Methods
    • Sensory and Motor Systems
  • ALERTS
  • FOR AUTHORS
  • ABOUT
    • Overview
    • Editorial Board
    • For the Media
    • Privacy Policy
    • Contact Us
    • Feedback
  • SUBMIT

User menu

Search

  • Advanced search
eNeuro

eNeuro

Advanced Search

 

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Blog
    • Collections
    • Podcast
  • TOPICS
    • Cognition and Behavior
    • Development
    • Disorders of the Nervous System
    • History, Teaching and Public Awareness
    • Integrative Systems
    • Neuronal Excitability
    • Novel Tools and Methods
    • Sensory and Motor Systems
  • ALERTS
  • FOR AUTHORS
  • ABOUT
    • Overview
    • Editorial Board
    • For the Media
    • Privacy Policy
    • Contact Us
    • Feedback
  • SUBMIT
PreviousNext
Research ArticleNew Research, Disorders of the Nervous System

Network Properties Revealed during Multi-Scale Calcium Imaging of Seizure Activity in Zebrafish

Jing Liu and Scott C. Baraban
eNeuro 25 February 2019, 6 (1) ENEURO.0041-19.2019; DOI: https://doi.org/10.1523/ENEURO.0041-19.2019
Jing Liu
Department of Neurological Surgery and Weill Institute for Neuroscience, University of California, San Francisco, CA 94143
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Jing Liu
Scott C. Baraban
Department of Neurological Surgery and Weill Institute for Neuroscience, University of California, San Francisco, CA 94143
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Scott C. Baraban
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Seizures are characterized by hypersynchronization of neuronal networks. Understanding these networks could provide a critical window for therapeutic control of recurrent seizure activity, i.e., epilepsy. However, imaging seizure networks has largely been limited to microcircuits in vitro or small “windows” in vivo. Here, we combine fast confocal imaging of genetically encoded calcium indicator (GCaMP)-expressing larval zebrafish with local field potential (LFP) recordings to study epileptiform events at whole-brain and single-neuron levels in vivo. Using an acute seizure model (pentylenetetrazole, PTZ), we reliably observed recurrent electrographic ictal-like events associated with generalized activation of all major brain regions and uncovered a well-preserved anterior-to-posterior seizure propagation pattern. We also examined brain-wide network synchronization and spatiotemporal patterns of neuronal activity in the optic tectum microcircuit. Brain-wide and single-neuronal level analysis of PTZ-exposed and 4-aminopyridine (4-AP)-exposed zebrafish revealed distinct network dynamics associated with seizure and non-seizure hyperexcitable states, respectively. Neuronal ensembles, comprised of coactive neurons, were also uncovered during interictal-like periods. Taken together, these results demonstrate that macro- and micro-network calcium motifs in zebrafish may provide a greater understanding of epilepsy.

  • epilepsy
  • fast confocal
  • neuronal networks
  • synchronization
  • whole-brain imaging
  • zebrafish

Significance Statement

Monitoring the dynamic activities in large-scale neuronal networks is critical to understanding seizure initiation and propagation. Here, we used well-established larval zebrafish seizure protocols and fast confocal imaging of genetically encoded calcium indicator (GCaMP)-expressing fish to investigate the epileptic network properties at brain-wide and single-cell levels. We revealed the rapid propagation of seizure activity from anterior-to-posterior brain regions in zebrafish CNS. We also showed that micro-ensembles of neuronal subpopulations are active during interictal-like periods in a manner similar to that seen in human electrophysiology data sets. Our findings demonstrate that these non-invasive optical imaging approaches will advance our understanding of the network basis underlying seizures and facilitate the development of methods to suppress these events.

Introduction

Synchronization of neuronal networks is a commonly observed event in the brain and underlies pathologic generation and propagation of epileptic seizures. Although electrophysiological features of an individual seizure event are well characterized using local field potential (LFP) and/or single-unit electrophysiological recordings (Grasse et al., 2013; Park et al., 2013; Schmidt et al., 2014; Truccolo et al., 2014; Weiss et al., 2016), these approaches do not capture brain-wide spatiotemporal network dynamics underlying these events. Efforts to image these networks have primarily been limited to small microcircuits contained within brain slice preparations or imaging “windows” in head-restrained mice (Trevelyan et al., 2006; Muldoon et al., 2013, 2015; Grosser et al., 2014; Kuzum et al., 2014; Hongo et al., 2015; Lillis et al., 2015; Rossi et al., 2017; Wenzel et al., 2017; Liou et al., 2018). Recent advances in imaging technology combined with stable expression of genetically encoded calcium indicators (GCaMPs; Tian et al., 2009; Akerboom et al., 2012; Chen et al., 2013) have led to a revolution in brain-wide activity mapping. GCaMPs are employed in a range of model organisms; among these, larval zebrafish (Danio rerio) offer an unparalleled combination of higher vertebrate relevance (including human), optical transparency, and a small brain size (comprising ∼100,000 neurons) allowing spatial coverage of multiple brain regions. Moreover, the availability of transgenic zebrafish, in which GCaMPs are under the control of neuronal elavl3 or NeuroD promoters has greatly facilitated the potential for assessing neural networks. Indeed, brain-wide imaging in GCaMP-expressing zebrafish has already been used to monitor spatiotemporal patterns of neural activity during fictive swimming (Ahrens et al., 2013a; Thiele et al., 2014), motor adaptation (Ahrens et al., 2012), visual processing (Akerboom et al., 2012; Muto et al., 2013), and prey capture (Muto and Kawakami, 2013; Semmelhack et al., 2014). During this time, zebrafish have also emerged as a valuable vertebrate model for neurologic disease, particularly for epilepsy (Hortopan et al., 2010; Ramirez et al., 2012; Stewart et al., 2014; Grone and Baraban, 2015; Martín-Jiménez et al., 2015; Griffin et al., 2017). Combining these advances, it is now possible to non-invasively monitor activity throughout the zebrafish nervous system during an epileptic seizure event.

Brain-wide imaging during seizures, in contrast to invasive electrophysiological recordings, increases the likelihood that network patterns will be discovered. To this end, we combined well-established acute seizure (pentylenetetrazole; PTZ; Baraban et al., 2005; Baxendale et al., 2012; Afrikanova et al., 2013; Buenafe et al., 2013; Orellana-Paucar et al., 2013; Rahn et al., 2014; Siebel et al., 2015; Barbalho et al., 2016; Copmans et al., 2018) and non-seizure hyperexcitability (4-aminopyridine; 4-AP) models (Ellis et al., 2012; Kumar et al., 2016; Cassar et al., 2017; Winter et al., 2017) with spinning disk confocal microscopy, enabling fast (20–30 fps) in vivo imaging of network events. Imaging studies were performed at brain-wide and single-cell resolution levels. As anesthetics can modify neuronal activity (Greenberg et al., 2008), experiments were done using “awake” agarose-embedded larval zebrafish paralyzed with pancuronium to reduce movement artifact. Simultaneous measurements of LFPs were used to confirm seizure activity. By analyzing brain-wide changes in calcium dynamics in an acute seizure model, we revealed neuronal populations in telencephalon initiating ictal-like seizure events that rapidly generalized to all CNS regions. The spatiotemporal patterns of calcium signal showed a rapid propagation of seizure activity from anterior-to-posterior brain regions. Subsequent analysis of spontaneous activity patterns established a clear pattern of synchronization during ictal-like seizure events consistent with a long-standing notion that seizures are a manifestation of hypersynchronous neural networks. Automated calcium signal processing and spike detection confirmed this high degree of synchronization at a single-neuron level. Finally, we showed that neuronal ensembles (Tao et al., 2011; Truccolo et al., 2011, 2014; Miller et al., 2014) are active during interictal-like periods. We suggest that these non-invasive optical imaging approaches could provide a more comprehensive understanding of complex neuronal networks involved in the generation and propagation of seizures.

Materials and Methods

Zebrafish lines

All procedures followed National Institute of Health and the University of California, San Francisco guidelines and were approved by the Institutional Animal Care and Use Committee. Calcium imaging experiments were performed on the Tg(neurod1:GCaMP6f) line (Rupprecht et al., 2016) in the nacre (mitfa-/-) background (kindly provided by Dr. C. Wyart from Institut du Cerveau et de la Moelle Épinière, Paris) on 5–6 d postfertilization (dpf). Adult zebrafish were maintained at 28°C on a 14/10 h light/dark cycle following standard methods. Larvae were raised in embryo media consisting of 0.03% Instant Ocean (Aquarium Systems, Inc.) and 0.0002% methylene blue in reverse osmosis-distilled water. At 5–6 dpf, zebrafish have not yet experienced sexual differentiation (Liew and Orbán, 2014).

Calcium imaging and electrophysiology

Zebrafish larvae were paralyzed in 300 µM pancuronium (Abcam), and then immobilized in a drop of 2% low-melting point agarose in a recording chamber (Fig. 1A). The recording chamber was placed on the stage of a Zeiss Axiocam upright microscope using a 5× or 20× objective (for whole-brain and neuron-level imaging, respectively), Yokogawa CSU-X1 Spinning Disk Confocal and a 470-nm laser light source (3i LaserStack). The recording chamber was filled with embryo media containing pancuronium (300 µM). Epileptiform activities were induced by bath application of PTZ (10 mM; Sigma-Aldrich) or 4-AP (4 mM; Sigma-Aldrich), and the fish were allowed to establish consistent epileptiform activities for at least 40 min before imaging studies were initiated (Fig. 1C). Images of 512 × 512 pixels were acquired at 20–30 Hz with an EMCCD camera (Photometrics Evolve). Multiple 5-min recordings (6000–9000 frames) were acquired for each experiment using SlideBook software (3i Intelligent Imaging Innovations). To record LFPs, a glass microelectrode was placed under visual guidance in the optic tectum or cerebellum. Electrodes were filled with 2 M NaCl, and electrical activity was recorded using an Axopatch 1D amplifier (Molecular Devices). Voltage records were low-pass filtered at 1 kHz (–3 dB, 8-pole Bessel), high-pass filtered at 0.1 Hz, digitized at 10 kHz using a Digidata 1520 A/D interface, and stored on a PC computer running Axoclamp software (Molecular Devices).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

High-speed confocal calcium imaging at different scales. A, Experimental setup. Simultaneous LFP recording and fast confocal imaging (20–30 fps) in agar-embedded larval zebrafish exposed to pancuronium (300 µM) with PTZ (10 mM) or 4-AP (4 mM). B, Representative high-resolution imaging of neurod1:GCaMP6f expressed larval zebrafish on 5–6 dpf. C, Experimental workflow. Recordings were obtained ∼40 min after drug application, and 5× and 20× objectives were used for whole-brain and neuron-level imaging, respectively. D, Schematic illustration depicting sub-regions of the larval zebrafish brain. E, ROIs and representative calcium traces (ΔF/F) of PTZ-induced ictal-like events; 1, 2: pallium; 3, 4: habenula; 5, 6: neuropil; 7, 8: SPV; 9, 10: cerebellum; 11, 12: hindbrain. F, Neuronal microcircuits within optic tectum and representative calcium traces of individual SPV neurons with PTZ-induced ictal-like event break-ins (underlined in red). Cerebellum was included as an indicator of ictal-like events that involve all brain regions. Scale bars as indicated in figure.

Image analysis

The raw images were first processed to correct motion drift using the Template Matching image registration plugin for ImageJ. Regions of interest (ROIs) for brain regions and single neurons were manually segmented with the ROI manager of ImageJ (Fig. 1E,F). Corrected image stacks and ROI segmentation file were then imported to MATLAB (MathWorks) for calcium fluorescence signal extraction and analysis using the FluroroSNNAP software (Patel et al., 2015). The calcium fluorescence signal for each brain region/neuron was obtained by averaging all pixels within the ROI. The fluorescence changes (ΔF/F) were calculated by subtracting each data point with the mean of the lower 50% of values within previous 10-s sliding window and normalized to the mean of the lower 50% of values within previous 10-s sliding window. Automated event detection was performed with the template-matching algorithm (Schultz et al., 2009) included in the FluroroSNNAP. A time-varying correlation coefficient between fluorescence trace and calcium transient templates (from the event wave form library) was calculated. Fluorescence transients with amplitude ΔF/F > 0.1 and correlation coefficient >0.85 were identified as events. Local maxima of the correlation coefficient trace identify the timestamps of calcium events. Since the algorithm may sometimes give errors, the event train was then manually corrected by deleting falsely detected events and adding events missed by the algorithm.

To determine the correlation between simultaneous calcium activity and LFP, the amplitude and duration of ictal-like seizure events in PTZ-treated zebrafish were evaluated. The amplitude is the absolute value of the maximum event peak, and the duration is the interval between the first and last points of the trace leaving and returning the baseline.

Network synchronization

Whole-brain network synchronization was evaluated with the network analysis module of FluroroSNNAP.

Based on the instantaneous phase of the fluorescence trace of each brain region, a pair-wise phase synchronization matrix was generated to identify patterns of network activity (Patel et al., 2015). According to the matrix, an eigenvalue-based algorithm was applied to calculate the global synchronization index which range from 0 (non-coordinated activity) to 1 (completely synchronized activity). For mathematical details see (Li et al., 2007). Briefly, the eigenvalues of the synchronization matrix were calculated and then normalized to generate an array of synchronization indices. The largest entry in this array is the global synchronization index. If the time series are fully uncorrelated, all elements in the synchronization matrix would be 0, and the eigenvalue would be 0, resulting a global synchronization index equal to 0. Alternatively, if the network was completely synchronized, all elements in the synchronization matrix would be 1, and the maximal eigenvalue would be equal to the number of time series with other eigenvalues falling to zero, resulting a global synchronization index equal to 1.

Seizure propagation

To determine the recruitment timing of each brain region to an ictal-like event, we took the time derivative (slope) of the fluorescence trace (smoothed by taking moving average of 50 data points), and used the local maxima (the fastest change in fluorescence) to identify the onset time point (Wenzel et al., 2017). Each brain region has two data points corresponding to the left and right counterparts. We then calculated the relative lag of each brain region with respect to the earliest onset within an ictal-like event and normalized it to the interval between the earliest onset and the latest onset. Normalized relative lags within an ictal-like event were divided into three quantile groups (first third, early; second third, middle; third third, late), and the probability of each brain region falling into each quantile group was calculated.

Detection of ensemble events

Ensemble events were defined as coactivation of a group of neurons in which a statistically significant number of neurons are active compared with surrogate datasets. We used a sliding window to generate a time series of coactivation of neurons by counting the number of events within a 0.5-s (10 frames) window. The binary event data were shuffled 2000 times within neurons, and the sliding window counting was performed. Frames with an observed number of coactive neurons >99.9% of all surrogate values (p < 0.001) were identified as highly active frames with ensemble event. To analyze the spatial extent of neuronal ensembles, the distances between the ensemble centroid and individual neurons within the ensemble were calculated using standard Euclidian distances.

Experimental design and statistical analysis

We used one-way ANOVA with Holm–Sidak post hoc analysis for multiple variable comparison, and χ2 test for onset timing distribution analysis, which were performed in Prism (GraphPad Software). Correlation coefficients (R 2) between calcium and electrical signals were calculated in MATLAB (MathWorks). Detection of neuronal ensembles was performed in MATLAB by comparing with surrogate datasets (see Materials and Methods). Individual analyses are described in Results.

Results

To study synchronized neural networks, we used two well-characterized pharmacological models employing bath application of PTZ (10 mM) or 4-AP (4 mM). These acute chemoconvulsant models work through different mechanisms: (1) PTZ acts as a GABAA receptor antagonist to disinhibit network activity (Baraban et al., 2005; Baxendale et al., 2012; Afrikanova et al., 2013; Buenafe et al., 2013; Orellana-Paucar et al., 2013; Rahn et al., 2014; Siebel et al., 2015; Barbalho et al., 2016; Copmans et al., 2018) and (2) 4-AP is a potassium channel blocker that enhances neuronal firing activity (Ellis et al., 2012; Kumar et al., 2016; Cassar et al., 2017; Winter et al., 2017). Together, these acute models provide a strategy to study recurrent spontaneous seizure events (PTZ) and non-seizure hyperexcitability (4-AP), respectively.

Multi-scale imaging of epileptiform activity

High-speed confocal microscopy was used to monitor ictal-like activity in agar-embedded larval zebrafish (Fig. 1A) with neuronal expression of GCaMP (neurod1:GCaMP6f; Fig. 1B). Larvae were immobilized in agar and imaged using a 5× (brain-wide) or 20× (single-cell) objective. For brain-wide studies, the confocal plane of view was adjusted to include forebrain (contains pallium and habenula), optic tectum (contains neuropil and stratum periventriculare; SPV), cerebellum and hindbrain (Fig. 1D,E). ROIs were drawn manually over anatomic areas corresponding to major brain regions, and fluorescence changes over time (ΔF/F) were extracted from time-lapse movies acquired at 20–30 Hz (Fig. 1E,F). In all PTZ-exposed larvae (n = 43), we observed large calcium transients at a frequency of one to two events per min, which occurred synchronously across all brain ROIs consistent with a generalized seizure classification (Fig. 1E; Movie 3). Large calcium transients were not seen in Control larvae bathed in embryo media (n = 9; Movie 1). As seizures are commonly associated with increased synchrony in neuronal populations (Penfield and Jasper, 1954; Jiruska et al., 2013), we also monitored calcium activity in SPV neurons within optic tectum microcircuits at single-cell resolution, which is facilitated by the mosaic neurod1-driven expression of GCaMP6f in this brain region (Fig. 1F). Exposure to PTZ elicited waves of calcium activity that showed a high degree of synchrony across all neurons (Movie 6).

Movie 1.

Representative whole-brain imaging with 5× objective from control fish. Movie played at 5× speed.

To determine whether synchronous calcium events correspond to electrographic discharges, we simultaneously monitored LFPs with an electrode placed in optic tectum SPV or cerebellum (as noted in Fig. 2B) in PTZ-treated larvae. Long-duration multi-spike ictal-like discharges coincided with large-amplitude calcium transients in every recording (Fig. 2A). The correlation between whole-brain calcium transient amplitude and LFP amplitude was low (R 2 = 0.208, n = 54 events from 10 PTZ-treated larvae; Fig. 2C). The relationship between whole-brain calcium transient duration and electrographic event duration was approximately linear (R 2 = 0.839, n = 54 events from 10 PTZ-treated larvae; Fig. 2D).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Correlation between LFP and calcium transients. A, Representative simultaneous calcium traces (green) from the cerebellum and LFP (black) recorded from optic tectum/cerebellum with recurrent PTZ-induced ictal-like seizures. Scale bars as indicated in figure. B, Fluorescence images of calcium activity during baseline and peak ictal-like event. The intensity of fluorescence is color coded as shown in the color bar. Events 1 and 2 as noted in A. Asterisk indicates the LFP recording site. C, D, Correlation between LFP and calcium transients in ictal-like event amplitude (C) and duration (D). Corresponding correlation coefficient R 2 is indicated in figure; n = 54 events from 10 PTZ-treated larvae.

Brain-wide network synchronization

To investigate characteristics of synchronization, we analyzed brain-wide network dynamics in hyperexcitable 4-AP-treated and epileptic PTZ-treated zebrafish. Age-matched untreated larvae were studied as Controls. 4-AP and PTZ-treated larvae showed distinct activity patterns in both whole-brain average fluorescence changes and LFP (Fig. 3A). 4-AP-treated larvae were characterized by high-frequency short-duration spikes in LFP and frequent low-amplitude calcium transients. In contrast, PTZ-treated larvae consistently exhibited recurrent high-amplitude long-duration ictal-like events in both LFP and calcium traces, and significantly increased area under curve (AUC) in fluorescence changes compared with control and 4-AP-treated larvae (p < 0.001, one-way ANOVA; n = 4, 7, and 7 fish for control, 4-AP, and PTZ condition, respectively; Fig. 3B).

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Brain-wide network synchronization. A, Representative whole-brain average fluorescence changes (ΔF/F; green) and LFP (black) recorded from control, 4-AP, and PTZ-treated larvae. B, Comparison of normalized AUC of fluorescence changes (norm. ΔF/F AUC). C, Representative fluorescence changes of each brain region under different conditions. ROIs refer to Figure 1E. D, Correlation matrices of whole-brain network activity for representative fish. E, Comparison of global synchronization (sync.) index. Scale bars as indicated in figure; n = 4, 7, and 7 fish for control, 4-AP, and PTZ condition, respectively. Error bars indicate SD. Statistical significance is indicated as ***p < 0.001.

Next, we analyzed the activity of subregions within the larval whole-brain network. Representative fluorescence traces of each brain region (ROI numbering refers to Fig. 1E) and corresponding correlation matrices are shown in Figure 3C,D. Larvae exposed to 4-AP showed frequent short-duration and low-amplitude calcium transients in cerebellum and hindbrain regions only. In contrast, PTZ-treated larvae showed widespread high-amplitude and long-duration calcium transients across the entire nervous system, and significantly increased brain-wide synchronization compared with Control and 4-AP-treated larvae (p < 0.001, one-way ANOVA; n = 4, 7, and 7 fish for control, 4-AP, and PTZ condition, respectively; Fig. 3E). Whole-brain activity of control, 4-AP, and PTZ-treated fish are exemplified in Movies 1, 2, 3, respectively.

Movie 2.

Representative whole-brain imaging with 5× objective from 4-AP fish. Movie played at 5× speed.

Movie 3.

Representative whole-brain imaging with 5× objective from PTZ fish. Movie played at 5× speed.

Seizure propagation

To determine brain regions recruited during seizure initiation and propagation, we analyzed calcium changes in multiple brain regions (Fig. 1E) during PTZ exposure (n = 30 events from 7 PTZ-treated larvae). Figure 4A shows representative fluorescence traces of each brain region during an ictal-like event break-in. The red trace is the time derivative of ΔF/F moving average of 50 nearby data points, and the local maxima position corresponds to the onset time point of the brain region. We calculated the relative lag of each brain region and normalized it to the earliest-to-latest onsets interval. Arranging the normalized relative lag of brain regions revealed a coarsely preserved temporal ordering across ictal-like seizure events (Fig. 4B). Figure 4C shows a color-coded onset time mapping from a representative PTZ-induced ictal-like seizure events in one zebrafish (average of four ictal-like events). We then divided the normalized relative lags within one ictal-like event into three quantile groups: early, middle and late (Fig. 4D), and calculated the probability of each brain region falling into each temporal group. This analysis revealed a conserved propagation pattern from the anterior-to-posterior brain (χ2 test X2(10) = 375.1, p < 0.001; n = 30 events from seven PTZ-treated larvae; Fig. 4E).

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Seizure propagation across the whole-brain. A, Representative calcium trace for each brain region within an PTZ-induced ictal-like event. Event onset time points were determined by the local maxima of the time derivative of moving average (of 50 data points) of ΔF/F. Scale bars as indicated in figure. B, Plot of normalized (norm.) relative lag from the earliest onset. C, Color-coded onset time mapping in a representative fish brain (average of four ictal events). D, Normalized relative lag for three quantile groups (first third early onset, second third middle onset, and third third late onset). E, Probability of brain regions falling into each quantile group. χ2 test X2(10) = 375.1, p < 0.001. Error bars indicate SD; n = 30 events from 7 PTZ-treated larvae.

Neuronal ensembles in optic tectum microcircuits

Next, we investigated network dynamics of the optic tectum microcircuit at single-neuron resolution. We analyzed individual SPV neuronal activity during interictal periods in representative control, 4-AP, and PTZ-treated fish, since ictal-like seizure events involve over 90% of the neuronal population in the field of view. Cerebellum activity, where baseline neurod1:GCaMP6f expression is prominent (Fig. 1B), was also extracted as a reference of ictal-like events that generalize to the entire nervous system. Based on binary event data generated from fluorescence changes with an automatic event detection algorithm (Patel et al., 2015), we constructed raster plots of neuronal activity, and then used a sliding window technique to generate a coactive neuron number time series (Fig. 5A). Ensemble events, namely a statistically significant number of coactive neurons compared with surrogate datasets, are marked by red arrows (Fig. 5A, bottom panel), and corresponding coactive neurons are colored red in the raster plot (Fig. 5A, middle panel). Chemoconvulsant (4-AP or PTZ) exposed larvae showed more frequent ensemble event occurrence compared with Controls. The neural firing rate (i.e., number of events in a 5-min recording epoch) distribution shows that both drug exposures increase excitability within the optic tectum network (Fig. 5B). For control fish, most neurons were inactive and firing at a low rate during the recording epoch. For 4-AP-treated fish, about half of the neurons exhibit firing rates between 1 and 10 events per 5 min, with ∼20% of neurons firing ≥10 events per 5 min. For PTZ-treated fish, no neurons were silent and firing rates were between 1 and 8 events per 5 min. Neurons in PTZ-treated fish showed significantly higher average firing rates than control fish, and neurons in 4-AP-treated fish had significantly higher average firing rate than both control and PTZ-treated fish (p < 0.001, one-way ANOVA; n = 75, 76, and 68 neurons for control, 4-AP, and PTZ-treated fish, respectively; Fig. 5C).

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

Neuronal ensembles in optic tectum microcircuits. A, Representative calcium traces (gray, individual SPV neuron; green, cerebellum), neuron activation raster plot (red, ensemble neurons) and number of coactive neurons (ensembles were indicated by red arrow) of control, 4-AP, and PTZ-treated fish. Cerebellum activity was shown as an indicator of ictal-like events that involve all brain regions. B, C, Histogram (B) and distributions (C) of the neural firing rate (event count per 5 min); n = 75, 76, and 68 neurons for control, 4-AP, and PTZ-treated fish, respectively. D, Spatial mapping of ensembles depicted in A. Dots in the same color represent coactive neurons, and lines are distances from ensemble neurons to the corresponding ensemble centroid. Each ensemble is indicated by a color. Scale bars as indicated in figure. E, Distribution of coactive neuron numbers as percentage of frames (6000 frames for each recording). F, G, Histogram (F) and distributions (G) of distances from ensemble neurons to the corresponding ensemble centroid; n = 17, 36, and 42 ensemble neurons for control, 4-AP, and PTZ-treated fish, respectively. Error bars indicate SD. Statistical significance is indicated as ***p < 0.001.

We then analyzed the spatial properties of neuronal ensembles. Representative spatial distribution plots of these ensembles (as noted in Fig. 5A) onto the optic tectum are shown in Figure 5D. Ensembles were detected in all control, 4-AP, and PTZ groups. Figure 5E shows the number of coactive neuron distributions within these ensembles. Control fish data showed few frames with coactive neurons (6000 frames for each recording). While 4-AP-treated fish showed more frames with small coactive neuron groups (number of coactive neurons <5), consistent with overall network hyperexcitability. PTZ-treated fish had more frames with larger coactive neuron groups (number of coactive neurons ≥5) than either Control or 4-AP-treated fish. In an epileptic network that is able to generate ictal-like events (i.e., PTZ-treated fish), spatially distributed ensembles are composed of larger groups of neurons.

The distance between ensemble centroid and each neuron within an ensemble also showed differences between groups (Fig. 5F,G). Ensembles observed in hyperexcitable 4-AP-treated fish network showed significantly larger spatial extent than in control or PTZ-treated fish (p < 0.001, one-way ANOVA; n = 17, 36, and 42 ensemble neurons for control, 4-AP, and PTZ-treated fish, respectively). Network dynamics of optic microcircuits are exemplified in Movie 4, 5, 6 for control, 4-AP, and PTZ-treated fish, respectively.

Movie 4.

Representative imaging in optic tectum microcircuits with 20× objective from control fish. Movie played at 5× speed.

Movie 5.

Representative imaging in optic tectum microcircuits with 20× objective from 4-AP fish. Movie played at 5× speed.

Movie 6.

Representative imaging in optic tectum microcircuits with 20× objective from PTZ fish. Movie played at 5× speed.

Discussion

Here, we combined sensitive GCaMPs and high-speed multi-scale monitoring of brain activity to characterize network dynamics in the CNS of larval zebrafish, following chemoconvulsant treatment. Our experiments revealed that ictal-like seizures are characterized by synchronous activation of the whole-brain network involving all major brain regions and propagate rapidly in a preserved temporal ordering: from anterior-to-posterior brain. PTZ-induced seizure and 4-AP-induced non-seizure hyperexcitability models showed distinct network dynamics in both whole-brain network and optic tectum microcircuits. Analysis of individual neuronal activities during interictal periods revealed neuronal ensembles in optic tectum, which demonstrated distinct properties in epileptic versus normal brain.

Optical monitoring in larval zebrafish provides a unique opportunity to map neuronal activity in a brain-wide manner. These measurements are not possible with electrophysiological approaches and led to spatiotemporal analysis of seizure initiation and propagation across the zebrafish CNS. Previous imaging studies using an acute PTZ seizure model in zebrafish focused on optical mapping of seizure propagation (Tao et al., 2011) or brain activities coincident with tail movements (Turrini et al., 2017), using mathematical modeling to reveal network connectivity changes during seizures (Rosch et al., 2018), and functional profiling of drug action in multiple brain regions (Winter et al., 2017). Generalized brain-wide activation during ictal-like seizure events was reported in these studies using acquisition speed ranges from 1 to 20 Hz. Here, using a fast spinning-disk confocal we observed similar patterns of brain activity with acquisition speeds up to 30 Hz and discovered a reproducible anterior-to-posterior propagation pattern of rapidly generalized optical seizures (Fig. 4).

To confirm the correlation between electrographic and optical activity, we acquired LFP and calcium imaging data simultaneously (Fig. 2). A high degree of correlation was observed in event duration but not in event amplitude (Fig. 2C,D). This is not surprising, as LFP event amplitude is dependent on electrode tip position in the brain relative to sources and sinks of current (Buzsáki et al., 2012), and calcium signal amplitude is linked to focal plane and GCaMP expression intensity variabilities between fish. Although both methods report generalized ictal-like seizure events consistently, LFP and calcium imaging do not measure identical activity. For example, single and multi-unit background activity shown in LFP recording in Control fish was not seen in corresponding calcium trace (Fig. 3A, left); 4-AP fish showed frequent high-amplitude spikes in LFP (Fig. 3A, middle), but little to no fluorescence change was observed in the recording region (optic tectum/cerebellum; ROIs 7–10 in Fig. 3C, middle). These observations reveal the relative advantages and disadvantages of each approach. Compared with LFP recording, GCaMP-based imaging provides superior spatial resolution but inferior detection sensitivity and speed (Chen et al., 2013), which may not effectively detect sub-threshold synaptic transmission or capture fast spikes seen in LFP recordings. Conversely, LFP recordings are limited to a single recording site in the brain.

4-AP has been reported to induce ictal-like discharges in vitro (Galvan et al., 1982; Avoli et al., 2002) and in vivo (Szente and Pongrácz, 1979; Peña and Tapia, 2000), but it is not a robust state in all conditions (Perreault and Avoli, 1991, 1992; Mattia et al., 1993; Zhu et al., 2008; Avoli and Jefferys, 2016). In this work, ictal-like seizure events were not consistently observed in 4-AP-exposed zebrafish larvae at a concentration of 4 mM, except a fatal long-lasting burst after prolonged exposures (data not shown), which is consistent with other zebrafish-based imaging studies (Winter et al., 2017). 4-AP-induced neuronal hyperexcitability revealed by LFP recordings (Fig. 3) and imaging of microcircuit dynamics (Fig. 5) was also in line with previous reported behavioral manifestations in 4-AP-treated zebrafish (Kumar et al., 2016) and provided us with an hyperexcitable state that was not complicated by the presence of ictal-like activity.

When examining neuronal dynamics in optic tectum microcircuits, we detected ensembles of coactive neurons during interictal periods and showed that, in PTZ-treated fish, ensembles were spatially localized, while in 4-AP-treated fish, ensembles were smaller (contained less neurons) and dispersed (Fig. 5). Neuronal ensembles could also be detected in Control fish, but the ensemble event occurrence was rare. Similar observations of functional clusters were reported in calcium imaging studies from hippocampal slices prepared from chronically epileptic mice (Muldoon et al., 2013). Interestingly, “microseizures,” which are discharges from highly confined neuron groups within ictal onset zones, have been observed in intracranial recordings from microelectrode arrays in human patients with partial epilepsy (Schevon et al., 2008, 2010; Stead et al., 2010). Although precise spatiotemporal information of micro-discharge clusters is absent in these studies, it raises an intriguing hypothesis that ictal networks are composed of sparse pathologic microdomains that can generate microseizures, and the spreading of microseizures results in ictal seizures (Stead et al., 2010). Here we employed a well-established acute zebrafish seizure model and high-speed calcium imaging with single-cell resolution to map neural activations. This strategy revealed neuronal ensembles in an ictal network (the optic tectum microcircuit), which possibly correspond to human microseizure domains. Future studies using three-dimensional light sheet imaging (Ahrens et al., 2013b; Portugues et al., 2014; Dunn et al., 2016) with enhanced scanning speed could improve the neuron sampling and provide an even better understanding of the role of neuronal ensembles in seizure generation and propagation.

In summary, we established a zebrafish-based calcium imaging platform to evaluate network dynamics underlying recurrent seizures, including network synchronization, seizure propagation and neuronal ensembles within local microcircuits. Future work with improved imaging techniques to explore the interplay between micro-ensembles and macro-seizures, combined with zebrafish models mimicking genetic forms of epilepsy, will provide insights into seizure prediction and suppression.

Acknowledgments

Acknowledgements: We thank Ehud Isacoff and members of the Isacoff laboratory (Claire Oldfield and Carlos Pantoja) at University of California Berkeley for early training and guidance using 3i imaging calcium imaging strategies in larval zebrafish.

Footnotes

  • S.C.B. is a co-Founder and Scientific Advisor for EpyGenix Therapeutics. J.L. declares no competing financial interests.

  • This work was supported by National Institutes of Health/National Institute of Neurological Disorders and Stroke Grants R01-NS096976 and R01-NS103139 and by a University of California San Francisco - University of California Berkeley (UCSF-UCB) Sackler Family Exchange sabbatical fellowship (to S.C.B.).

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

References

  1. ↵
    Afrikanova T, Serruys A-SK, Buenafe OEM, Clinckers R, Smolders I, de Witte PAM, Crawford AD, Esguerra CV (2013) Validation of the zebrafish pentylenetetrazol seizure model: locomotor versus electrographic responses to antiepileptic drugs. PLoS One 8:e54166. doi:10.1371/journal.pone.0054166 pmid:23342097
    OpenUrlCrossRefPubMed
  2. ↵
    Ahrens MB, Li JM, Orger MB, Robson DN, Schier AF, Engert F, Portugues R (2012) Brain-wide neuronal dynamics during motor adaptation in zebrafish. Nature 485:471–477. doi:10.1038/nature11057 pmid:22622571
    OpenUrlCrossRefPubMed
  3. ↵
    Ahrens MB, Huang K-H, Narayan S, Mensh BD, Engert F (2013a) Two-photon calcium imaging during fictive navigation in virtual environments. Front Neural Circuits 7:104. doi:10.3389/fncir.2013.00104 pmid:23761738
    OpenUrlCrossRefPubMed
  4. ↵
    Ahrens MB, Orger MB, Robson DN, Li JM, Keller PJ (2013b) Whole-brain functional imaging at cellular resolution using light-sheet microscopy. Nat Methods 10:413–420. doi:10.1038/nmeth.2434 pmid:23524393
    OpenUrlCrossRefPubMed
  5. ↵
    Akerboom J, Chen TW, Wardill TJ, Tian L, Marvin JS, Mutlu S, Calderón NC, Esposti F, Borghuis BG, Sun XR, Gordus A, Orger MB, Portugues R, Engert F, Macklin JJ, Filosa A, Aggarwal A, Kerr RA, Takagi R, Kracun S, et al. (2012) Optimization of a GCaMP calcium indicator for neural activity imaging. J Neurosci 32:13819–13840. doi:10.1523/JNEUROSCI.2601-12.2012 pmid:23035093
    OpenUrlAbstract/FREE Full Text
  6. ↵
    Avoli M, Jefferys JGR (2016) Models of drug-induced epileptiform synchronization in vitro. J Neurosci Methods 260:26–32. doi:10.1016/j.jneumeth.2015.10.006 pmid:26484784
    OpenUrlCrossRefPubMed
  7. ↵
    Avoli M, D’Antuono M, Louvel J, Köhling R, Biagini G, Pumain R, D’Arcangelo G, Tancredi V (2002) Network and pharmacological mechanisms leading to epileptiform synchronization in the limbic system in vitro. Prog Neurobiol 68:167–207. doi:10.1016/S0301-0082(02)00077-1 pmid:12450487
    OpenUrlCrossRefPubMed
  8. ↵
    Baraban SC, Taylor MR, Castro PA, Baier H (2005) Pentylenetetrazole induced changes in zebrafish behavior, neural activity and c-fos expression. Neuroscience 131:759–768. doi:10.1016/j.neuroscience.2004.11.031 pmid:15730879
    OpenUrlCrossRefPubMed
  9. ↵
    Barbalho PG, Lopes-Cendes I, Maurer-Morelli CV (2016) Indomethacin treatment prior to pentylenetetrazole-induced seizures downregulates the expression of il1b and cox2 and decreases seizure-like behavior in zebrafish larvae. BMC Neurosci 17:12. doi:10.1186/s12868-016-0246-y pmid:26961169
    OpenUrlCrossRefPubMed
  10. ↵
    Baxendale S, Holdsworth CJ, Santoscoy PLM, Harrison MRM, Fox J, Parkin CA, Ingham PW, Cunliffe VT (2012) Identification of compounds with anti-convulsant properties in a zebrafish model of epileptic seizures. Dis Model Mech 5:773–784. doi:10.1242/dmm.010090 pmid:22730455
    OpenUrlAbstract/FREE Full Text
  11. ↵
    Buenafe OE, Orellana-Paucar A, Maes J, Huang H, Ying X, De Borggraeve W, Crawford AD, Luyten W, Esguerra CV, de Witte P (2013) Tanshinone IIA exhibits anticonvulsant activity in zebrafish and mouse seizure models. ACS Chem Neurosci 4:1479–1487. doi:10.1021/cn400140e pmid:23937066
    OpenUrlCrossRefPubMed
  12. ↵
    Buzsáki G, Anastassiou CA, Koch C (2012) The origin of extracellular fields and currents—EEG, ECoG, LFP and spikes. Nat Rev Neurosci 13:407–420. doi:10.1038/nrn3241 pmid:22595786
    OpenUrlCrossRefPubMed
  13. ↵
    Cassar S, Breidenbach L, Olson A, Huang X, Britton H, Woody C, Sancheti P, Stolarik D, Wicke K, Hempel K, LeRoy B (2017) Measuring drug absorption improves interpretation of behavioral responses in a larval zebrafish locomotor assay for predicting seizure liability. J Pharmacol Toxicol Methods 88:56–63. doi:10.1016/j.vascn.2017.07.002 pmid:28712933
    OpenUrlCrossRefPubMed
  14. ↵
    Chen TW, Wardill TJ, Sun Y, Pulver SR, Renninger SL, Baohan A, Schreiter ER, Kerr RA, Orger MB, Jayaraman V, Looger LL, Svoboda K, Kim DS (2013) Ultrasensitive fluorescent proteins for imaging neuronal activity. Nature 499:295–300. doi:10.1038/nature12354 pmid:23868258
    OpenUrlCrossRefPubMed
  15. ↵
    Copmans D, Orellana-Paucar AM, Steurs G, Zhang Y, Ny A, Foubert K, Exarchou V, Siekierska A, Kim Y, De Borggraeve W, Dehaen W, Pieters L, de Witte PAM (2018) Methylated flavonoids as anti-seizure agents: naringenin 4′,7-dimethyl ether attenuates epileptic seizures in zebrafish and mouse models. Neurochem Int 112:124–133. doi:10.1016/j.neuint.2017.11.011 pmid:29174382
    OpenUrlCrossRefPubMed
  16. ↵
    Dunn TW, Mu Y, Narayan S, Randlett O, Naumann EA, Yang C-T, Schier AF, Freeman J, Engert F, Ahrens MB (2016) Brain-wide mapping of neural activity controlling zebrafish exploratory locomotion. Elife 5:e12741. doi:10.7554/elife.12741 pmid:27003593
    OpenUrlCrossRefPubMed
  17. ↵
    Ellis LD, Seibert J, Soanes KH (2012) Distinct models of induced hyperactivity in zebrafish larvae. Brain Res 1449:46–59. doi:10.1016/j.brainres.2012.02.022 pmid:22386495
    OpenUrlCrossRefPubMed
  18. ↵
    Galvan M, Grafe P, Bruggencate GT (1982) Convulsant actions of 4-aminopyridine on the guinea-pig olfactory cortex slice. Brain Res 241:75–86. doi:10.1016/0006-8993(82)91230-6 pmid:7104708
    OpenUrlCrossRefPubMed
  19. ↵
    Grasse DW, Karunakaran S, Moxon KA (2013) Neuronal synchrony and the transition to spontaneous seizures. Exp Neurol 248:72–84. doi:10.1016/j.expneurol.2013.05.004 pmid:23707218
    OpenUrlCrossRefPubMed
  20. ↵
    Greenberg DS, Houweling AR, Kerr JND (2008) Population imaging of ongoing neuronal activity in the visual cortex of awake rats. Nat Neurosci 11:749–751. doi:10.1038/nn.2140 pmid:18552841
    OpenUrlCrossRefPubMed
  21. ↵
    Griffin A, Hamling KR, Knupp K, Hong S, Lee LP, Baraban SC (2017) Clemizole and modulators of serotonin signalling suppress seizures in Dravet syndrome. Brain 140:669–683. doi:10.1093/brain/aww342 pmid:28073790
    OpenUrlCrossRefPubMed
  22. ↵
    Grone BP, Baraban SC (2015) Animal models in epilepsy research: legacies and new directions. Nat Neurosci 18:339–343. doi:10.1038/nn.3934 pmid:25710835
    OpenUrlCrossRefPubMed
  23. ↵
    Grosser S, Queenan BN, Lalchandani RR, Vicini S (2014) Hilar somatostatin interneurons contribute to synchronized GABA activity in an in vitro epilepsy model. PLoS One 9:e86250. doi:10.1371/journal.pone.0086250 pmid:24465989
    OpenUrlCrossRefPubMed
  24. ↵
    Hongo Y, Takasu K, Ikegaya Y, Hasegawa M, Sakaguchi G, Ogawa K (2015) Heterogeneous effects of antiepileptic drugs in an in vitro epilepsy model – a functional multineuron calcium imaging study. Eur J Neurosci 42:1818–1829. doi:10.1111/ejn.12945 pmid:25967117
    OpenUrlCrossRefPubMed
  25. ↵
    Hortopan GA, Dinday MT, Baraban SC (2010) Zebrafish as a model for studying genetic aspects of epilepsy. Dis Model Mech 3:144–148. doi:10.1242/dmm.002139 pmid:20212082
    OpenUrlAbstract/FREE Full Text
  26. ↵
    Jiruska P, de Curtis M, Jefferys JGR, Schevon CA, Schiff SJ, Schindler K (2013) Synchronization and desynchronization in epilepsy: controversies and hypotheses. J Physiol 591:787–797. doi:10.1113/jphysiol.2012.239590 pmid:23184516
    OpenUrlCrossRefPubMed
  27. ↵
    Kumar MG, Rowley S, Fulton R, Dinday MT, Baraban SC, Patel M (2016) Altered glycolysis and mitochondrial respiration in a zebrafish model of Dravet syndrome. eNeuro 3: ENEURO.0008-16.2016. doi:10.1523/eneuro.0008-16.2016 pmid:27066534
    OpenUrlCrossRefPubMed
  28. ↵
    Kuzum D, Takano H, Shim E, Reed JC, Juul H, Richardson AG, de Vries J, Bink H, Dichter MA, Lucas TH, Coulter DA, Cubukcu E, Litt B (2014) Transparent and flexible low noise graphene electrodes for simultaneous electrophysiology and neuroimaging. Nat Commun 5:5259. doi:10.1038/ncomms6259 pmid:25327632
    OpenUrlCrossRefPubMed
  29. ↵
    Li X, Cui D, Jiruska P, Fox JE, Yao X, Jefferys JGR (2007) Synchronization measurement of multiple neuronal populations. J Neurophysiol 98:3341–3348. doi:10.1152/jn.00977.2007 pmid:17913983
    OpenUrlCrossRefPubMed
  30. ↵
    Liew WC, Orbán L (2014) Zebrafish sex: a complicated affair. Brief Funct Genomics 13:172–187. doi:10.1093/bfgp/elt041 pmid:24148942
    OpenUrlCrossRefPubMed
  31. ↵
    Lillis KP, Wang Z, Mail M, Zhao GQ, Berdichevsky Y, Bacskai B, Staley KJ (2015) Evolution of network synchronization during early epileptogenesis parallels synaptic circuit alterations. J Neurosci 35:9920–9934. doi:10.1523/JNEUROSCI.4007-14.2015 pmid:26156993
    OpenUrlAbstract/FREE Full Text
  32. ↵
    Liou J, Ma H, Wenzel M, Zhao M, Baird-Daniel E, Smith EH, Daniel A, Emerson R, Yuste R, Schwartz TH, Schevon CA (2018) Role of inhibitory control in modulating focal seizure spread. Brain 141:2083–2097. doi:10.1093/brain/awy116 pmid:29757347
    OpenUrlCrossRefPubMed
  33. ↵
    Martín-Jiménez R, Campanella M, Russell C (2015) New zebrafish models of neurodegeneration. Curr Neurol Neurosci Rep 15:33. doi:10.1007/s11910-015-0555-z pmid:25903297
    OpenUrlCrossRefPubMed
  34. ↵
    Mattia D, Hwa GGC, Avoli M (1993) Epileptiform activity induced by 4-aminopyridine in guinea-pig and rat neocortices. Neurosci Lett 154:157–160. doi:10.1016/0304-3940(93)90195-q pmid:8103197
    OpenUrlCrossRefPubMed
  35. ↵
    Miller JK, Ayzenshtat I, Carrillo-Reid L, Yuste R (2014) Visual stimuli recruit intrinsically generated cortical ensembles. Proc Natl Acad Sci USA 111:E4053–E4061. doi:10.1073/pnas.1406077111 pmid:25201983
    OpenUrlAbstract/FREE Full Text
  36. ↵
    Muldoon SF, Soltesz I, Cossart R (2013) Spatially clustered neuronal assemblies comprise the microstructure of synchrony in chronically epileptic networks. Proc Natl Acad Sci USA 110:3567–3572. doi:10.1073/pnas.1216958110 pmid:23401510
    OpenUrlAbstract/FREE Full Text
  37. ↵
    Muldoon SF, Villette V, Tressard T, Malvache A, Reichinnek S, Bartolomei F, Cossart R (2015) GABAergic inhibition shapes interictal dynamics in awake epileptic mice. Brain 138:2875–2890. doi:10.1093/brain/awv227 pmid:26280596
    OpenUrlCrossRefPubMed
  38. ↵
    Muto A, Kawakami K (2013) Prey capture in zebrafish larvae serves as a model to study cognitive functions. Front Neural Circuits 7:110. doi:10.3389/fncir.2013.00110 pmid:23781176
    OpenUrlCrossRefPubMed
  39. ↵
    Muto A, Ohkura M, Abe G, Nakai J, Kawakami K (2013) Real-time visualization of neuronal activity during perception. Curr Biol 23:307–311. doi:10.1016/j.cub.2012.12.040 pmid:23375894
    OpenUrlCrossRefPubMed
  40. ↵
    Orellana-Paucar AM, Afrikanova T, Thomas J, Aibuldinov YK, Dehaen W, de Witte PAM, Esguerra CV (2013) Insights from zebrafish and mouse models on the activity and safety of ar-turmerone as a potential drug candidate for the treatment of epilepsy. PLoS One 8:e81634. doi:10.1371/journal.pone.0081634 pmid:24349101
    OpenUrlCrossRefPubMed
  41. ↵
    Park YS, Hochberg LR, Eskandar EN, Cash SS, Truccolo W (2013) Early detection of human epileptic seizures based on intracortical local field potentials. Int IEEE EMBS Conf Neural Eng 2013:323–326.
    OpenUrl
  42. ↵
    Patel TP, Man K, Firestein BL, Meaney DF (2015) Automated quantification of neuronal networks and single-cell calcium dynamics using calcium imaging. J Neurosci Methods 243:26–38. doi:10.1016/j.jneumeth.2015.01.020 pmid:25629800
    OpenUrlCrossRefPubMed
  43. ↵
    Peña F, Tapia R (2000) Seizures and neurodegeneration induced by 4-aminopyridine in rat hippocampus in vivo: role of glutamate- and GABA-mediated neurotransmission and of ion channels. Neuroscience 101:547–561. doi:10.1016/s0306-4522(00)00400-0 pmid:11113304
    OpenUrlCrossRefPubMed
  44. ↵
    Penfield W, Jasper H (1954) Epilepsy and the functional anatomy of the human brain. Oxford, UK: Little, Brown & Co.
  45. ↵
    Perreault P, Avoli M (1991) Physiology and pharmacology of epileptiform activity induced by 4-aminopyridine in rat hippocampal slices. J Neurophysiol 65:771–785. doi:10.1152/jn.1991.65.4.771 pmid:1675671
    OpenUrlCrossRefPubMed
  46. ↵
    Perreault P, Avoli M (1992) 4-aminopyridine-induced epileptiform activity and a GABA-mediated long-lasting depolarization in the rat hippocampus. J Neurosci 12:104–115. doi:10.1523/JNEUROSCI.12-01-00104.1992 pmid:1309571
    OpenUrlAbstract/FREE Full Text
  47. ↵
    Portugues R, Feierstein CE, Engert F, Orger MB (2014) Whole-brain activity maps reveal stereotyped, distributed networks for visuomotor behavior. Neuron 81:1328–1343. doi:10.1016/j.neuron.2014.01.019 pmid:24656252
    OpenUrlCrossRefPubMed
  48. ↵
    Rahn JJ, Bestman JE, Josey BJ, Inks ES, Stackley KD, Rogers CE, Chou CJ, Chan SSL (2014) Novel Vitamin K analogs suppress seizures in zebrafish and mouse models of epilepsy. Neuroscience 259:142–154. doi:10.1016/j.neuroscience.2013.11.040 pmid:24291671
    OpenUrlCrossRefPubMed
  49. ↵
    Ramirez IBR, Pietka G, Jones DR, Divecha N, Alia A, Baraban SC, Hurlstone AFL, Lowe M (2012) Impaired neural development in a zebrafish model for Lowe syndrome. Hum Mol Genet 21:1744–1759. doi:10.1093/hmg/ddr608 pmid:22210625
    OpenUrlCrossRefPubMed
  50. ↵
    Rosch RE, Hunter PR, Baldeweg T, Friston KJ, Meyer MP (2018) Calcium imaging and dynamic causal modelling reveal brain-wide changes in effective connectivity and synaptic dynamics during epileptic seizures. PLoS Comput Biol 14:e1006375. doi:10.1371/journal.pcbi.1006375 pmid:30138336
    OpenUrlCrossRefPubMed
  51. ↵
    Rossi LF, Wykes RC, Kullmann DM, Carandini M (2017) Focal cortical seizures start as standing waves and propagate respecting homotopic connectivity. Nat Commun 8:217. doi:10.1038/s41467-017-00159-6 pmid:28794407
    OpenUrlCrossRefPubMed
  52. ↵
    Rupprecht P, Prendergast A, Wyart C, Friedrich RW (2016) Remote z-scanning with a macroscopic voice coil motor for fast 3D multiphoton laser scanning microscopy. Biomed Opt Express 7:1656–1671. doi:10.1364/BOE.7.001656 pmid:27231612
    OpenUrlCrossRefPubMed
  53. ↵
    Schevon CA, Ng SK, Cappell J, Goodman RR, McKhann G, Waziri A, Branner A, Sosunov A, Schroeder CE, Emerson RG (2008) Microphysiology of epileptiform activity in human neocortex. J Clin Neurophysiol off Publ Am Electroencephalogr Soc 25:321–330. doi:10.1097/WNP.0b013e31818e8010 pmid:18997628
    OpenUrlCrossRefPubMed
  54. ↵
    Schevon CA, Goodman RR, McKhann G, Emerson RG (2010) Propagation of epileptiform activity on a submillimeter scale. J Clin Neurophysiol 27:406–411. doi:10.1097/WNP.0b013e3181fdf8a1 pmid:21076338
    OpenUrlCrossRefPubMed
  55. ↵
    Schmidt H, Petkov G, Richardson MP, Terry JR (2014) Dynamics on networks: the role of local dynamics and global networks on the emergence of hypersynchronous neural activity. PLoS Comput Biol 10:e1003947. doi:10.1371/journal.pcbi.1003947 pmid:25393751
    OpenUrlCrossRefPubMed
  56. ↵
    Schultz SR, Kitamura K, Post-Uiterweer A, Krupic J, Häusser M (2009) Spatial pattern coding of sensory information by climbing fiber-evoked calcium signals in networks of neighboring cerebellar Purkinje cells. J Neurosci 29:8005–8015. doi:10.1523/JNEUROSCI.4919-08.2009 pmid:19553440
    OpenUrlAbstract/FREE Full Text
  57. ↵
    Semmelhack JL, Donovan JC, Thiele TR, Kuehn E, Laurell E, Baier H (2014) A dedicated visual pathway for prey detection in larval zebrafish. Elife 3:e04878. doi:10.7554/elife.04878 pmid:25490154
    OpenUrlCrossRefPubMed
  58. ↵
    Siebel AM, Menezes FP, da Costa Schaefer I, Petersen BD, Bonan CD (2015) Rapamycin suppresses PTZ-induced seizures at different developmental stages of zebrafish. Pharmacol Biochem Behav 139:163–168. doi:10.1016/j.pbb.2015.05.022 pmid:26051026
    OpenUrlCrossRefPubMed
  59. ↵
    Stead M, Bower M, Brinkmann BH, Lee K, Marsh WR, Meyer FB, Litt B, Van Gompel J, Worrell GA (2010) Microseizures and the spatiotemporal scales of human partial epilepsy. Brain 133:2789–2797. doi:10.1093/brain/awq190 pmid:20685804
    OpenUrlCrossRefPubMed
  60. ↵
    Stewart AM, Braubach O, Spitsbergen J, Gerlai R, Kalueff AV (2014) Zebrafish models for translational neuroscience research: from tank to bedside. Trends Neurosci 37:264–278. doi:10.1016/j.tins.2014.02.011 pmid:24726051
    OpenUrlCrossRefPubMed
  61. ↵
    Szente M, Pongrácz F (1979) Aminopyridine-induced seizure activity. Electroencephalogr Clin Neurophysiol 46:605–608. doi:10.1016/0013-4694(79)90014-2 pmid:88349
    OpenUrlCrossRefPubMed
  62. ↵
    Tao L, Lauderdale JD, Sornborger AT (2011) Mapping functional connectivity between neuronal ensembles with larval zebrafish transgenic for a ratiometric calcium indicator. Front Neural Circuits 5:2. doi:10.3389/fncir.2011.00002 pmid:21373259
    OpenUrlCrossRefPubMed
  63. ↵
    Thiele TR, Donovan JC, Baier H (2014) Descending control of swim posture by a midbrain nucleus in zebrafish. Neuron 83:679–691. doi:10.1016/j.neuron.2014.04.018 pmid:25066082
    OpenUrlCrossRefPubMed
  64. ↵
    Tian L, Hires SA, Mao T, Huber D, Chiappe ME, Chalasani SH, Petreanu L, Akerboom J, McKinney SA, Schreiter ER, Bargmann CI, Jayaraman V, Svoboda K, Looger LL (2009) Imaging neural activity in worms, flies and mice with improved GCaMP calcium indicators. Nat Methods 6:875–881. doi:10.1038/nmeth.1398 pmid:19898485
    OpenUrlCrossRefPubMed
  65. ↵
    Trevelyan AJ, Sussillo D, Watson BO, Yuste R (2006) Modular propagation of epileptiform activity: evidence for an inhibitory veto in neocortex. J Neurosci 26:12447–12455. doi:10.1523/JNEUROSCI.2787-06.2006 pmid:17135406
    OpenUrlAbstract/FREE Full Text
  66. ↵
    Truccolo W, Donoghue JA, Hochberg LR, Eskandar EN, Madsen JR, Anderson WS, Brown EN, Halgren E, Cash SS (2011) Single-neuron dynamics in human focal epilepsy. Nat Neurosci 14:635–641. doi:10.1038/nn.2782 pmid:21441925
    OpenUrlCrossRefPubMed
  67. ↵
    Truccolo W, Ahmed OJ, Harrison MT, Eskandar EN, Cosgrove GR, Madsen JR, Blum AS, Potter NS, Hochberg LR, Cash SS (2014) Neuronal ensemble synchrony during human focal seizures. J Neurosci 34:9927–9944. doi:10.1523/JNEUROSCI.4567-13.2014 pmid:25057195
    OpenUrlAbstract/FREE Full Text
  68. ↵
    Turrini L, Fornetto C, Marchetto G, Müllenbroich MC, Tiso N, Vettori A, Resta F, Masi A, Mannaioni G, Pavone FS, Vanzi F (2017) Optical mapping of neuronal activity during seizures in zebrafish. Sci Rep 7:3025. doi:10.1038/s41598-017-03087-z pmid:28596596
    OpenUrlCrossRefPubMed
  69. ↵
    Weiss SA, Alvarado-Rojas C, Bragin A, Behnke E, Fields T, Fried I, Engel J, Staba R (2016) Ictal onset patterns of local field potentials, high frequency oscillations, and unit activity in human mesial temporal lobe epilepsy. Epilepsia 57:111–121. doi:10.1111/epi.13251 pmid:26611159
    OpenUrlCrossRefPubMed
  70. ↵
    Wenzel M, Hamm JP, Peterka DS, Yuste R (2017) Reliable and elastic propagation of cortical seizures in vivo. Cell Rep 19:2681–2693. doi:10.1016/j.celrep.2017.05.090 pmid:28658617
    OpenUrlCrossRefPubMed
  71. ↵
    Winter MJ, Windell D, Metz J, Matthews P, Pinion J, Brown JT, Hetheridge MJ, Ball JS, Owen SF, Redfern WS, Moger J, Randall AD, Tyler CR (2017) 4-dimensional functional profiling in the convulsant-treated larval zebrafish brain. Sci Rep 7:6581. doi:10.1038/s41598-017-06646-6 pmid:28747660
    OpenUrlCrossRefPubMed
  72. ↵
    Zhu L, Polley N, Mathews GC, Delpire E (2008) NKCC1 and KCC2 prevent hyperexcitability in the mouse hippocampus. Epilepsy Res 79:201–212. doi:10.1016/j.eplepsyres.2008.02.005 pmid:18394864
    OpenUrlCrossRefPubMed

Synthesis

Reviewing Editor: Liset Menendez de la Prida, Instituto Cajal CSIC

Decisions are customarily a result of the Reviewing Editor and the peer reviewers coming together and discussing their recommendations until a consensus is reached. When revisions are invited, a fact-based synthesis statement explaining their decision and outlining what is needed to prepare a revision will be listed below. Note: If this manuscript was transferred from JNeurosci and a decision was made to accept the manuscript without peer review, a brief statement to this effect will instead be what is listed below.

Manuscript transfered from JN after reviewing. Authors have made a good job in addressing issues raised by the previous reviewers at JN. I have personally checked their answers to their points and found the current version suitable for eNeuron. This is an interesting paper that combines genetically-encoded calcium indicators with whole brain imaging in zebrafish larvae to monitor pharmacologically induced activity patterns. The paper is relevant with potential application in anti-seizure drug screening. I therefore accept this paper for publication.

For the records, here is a synthesis of the issues raised in the reviewing version at JN.

Major comments:

From looking at the data, I cannot agree that the current method for detecting the regional onset of an ictal-like event relates to the actual regional onset of the event. Looking at Fig 4, there is a clear upward “notch” in the df/f signal at the onset of the event, and the since the 50 point averaging smooths the signal, the point of steepest ascent does not depict the “onset” that one would naturally ascribe. It seems that the region is slow to fully participate in the event however, given the soft shape of the df/f rise, which I assume is due to the large spatial coverage of the region and a non uniform start of the event within the region? Perhaps it would be more accurate to look at the peak point of participation (max amp of df/f within the event) and assess seizure spreading based on this? Given the shape of the signals, I cannot imagine a biological interpretation of the maximum of the time derivative as is currently calculated.

I find the use of the phrase “neuronal ensemble” to refer both to the time point, and neurons involved in time points containing a statistically significant increase in neuronal firing to be rather confusing. It would be more clear to define “events” in time and then discuss the neuronal population of said events.

Was any sort of cluster analysis applied to determine if “neuronal ensembles” reflected the co-activation of smaller ensembles as in Muldoon et al, 2013? At least in the PTZ data, it looks like certain “neuronal ensembles” contain the same neurons in different events (small ensembles in Fig 5A); do multiple small ensembles combine in other larger events?

Is there a relationship between the neuronal populations involved in events in 4-AP data and PTZ data? For example, if neurons 3,4,5 fired together in a neuronal ensemble in 4-AP events, did they also fire together in PTZ events? This is somewhat related to my comment above about the larger PTZ events potentially being comprised of multiple smaller ensembles. Given the relationship between spatial properties of the ensembles under the two conditions this seems not likely, but it would be worth investigating in more detail.

I noted with interest that the large electrographic discharges in 4-AP were generally not associated with very large Ca signal, consistent with others' observations in rodent tissue that this, at least initially, triggers primarily interneuronal activation (eg various works by Avoli & de Curtis, mainly of electrophysiological recordings, but relevant is Parrish et al 2018, J.Neurophys, which also examined 4AP activity using GCaMP expressing mice) - this relationship to the rodent work should probably be discussed. In fact, 4-AP did not appear to trigger full seizures, which was perhaps counter to expectation - full seizures were only really seen with the PTZ. They described a consistent anterior to posterior propagation pattern in these PTZ induced seizures. The 4-AP only showed “ensemble” activity, albeit at a higher rate than in untreated control preparations, and no description of any propagation was provided of these - I'm guessing that they did not propagate, but were “standing waves”, but could this be clarified, please?

The lack of seizures in the 4-AP is a shame - I would have expected these. It means that the only seizures are essentially from a disinhibition model, which is likely to show the least interesting neuronal dynamics. Disinhibition models in rodent tissue also propagate widely and with great speed, whereas far more interesting dynamics generally arise from the interplay between inhibitory elements in the network. Having a model that showed both the interesting ensembles and also transitions into seizures may have had a lot more substance, since it would allow one to examine pathways into the seizures. One benefit of being able to image the entire brain is that it does allow one to hone in on where the seizures start. This appears to be the pallium, yet the analysis of ensembles focussed largely on a different brain area, namely the tectum.

A final major concern is that this felt like it was primarily a presentation of new technique, and there was not a whole lot of new neurobiology. I also thought that the background pitch was a little disingenuous in its claims that previous imaging / recordings failed to capture widespread network dynamics. It should be noted that the MEA studies use arrays that are larger than the entire larval zebrafish brain. Furthermore, the Carandini lab (Rossi et al, 2017, Nat Comm) has used a similar approach to visualise epileptic propagation across wide fields of view in mice. There were no references generally about GECIs or GCaMP.

Minor comments:

How are the calcium event and electrographic event amplitude and duration defined?

The reference for calculating the global synchronization index in the “Network Synchronization” section appears to be incorrect. Please better define this measure and include the proper reference.

It seems that the “correlation matrices” of Figure 3E actually pairwise synchronization matrices computed using the phase synchronization metric.

In terms of presentation, the manuscript was generally well written, and I did not find many mistakes. A handful of points might be clarified. In figure 2C, were these mixed data from the two models, or all PTZ (I'm guessing, since 4-AP only produced low level activation). Likewise, in figure 4, is all the data from PTZ treated animals too (again, I'm guessing because these are ictal)?

Figure 5 - the legend indicates that these are analyses of optic tectum activity, yet panel A is labelled as “cerebellum”.

I was slightly surprised that the 4-AP appeared to show larger ensembles than the PTZ, but was this because larger PTZ events invariably escalated to be full seizures?

The terminology did not seem quite right to me, with the 4-AP data being described as a model of hyperexcitability. I'd have said that the PTZ also fell under that umbrella term. This is stated clearly early in the results section, but elsewhere, when the authors refer to hyperexcitable states, it seems odd not to include the PTZ data in these, and I know that is not what is intended here. I just felt that this introduced potential for confusion, and the authors might care to rethink their terminology.

Back to top

In this issue

eneuro: 6 (1)
eNeuro
Vol. 6, Issue 1
January/February 2019
  • Table of Contents
  • Index by author
  • Ed Board (PDF)
Email

Thank you for sharing this eNeuro article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Network Properties Revealed during Multi-Scale Calcium Imaging of Seizure Activity in Zebrafish
(Your Name) has forwarded a page to you from eNeuro
(Your Name) thought you would be interested in this article in eNeuro.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Network Properties Revealed during Multi-Scale Calcium Imaging of Seizure Activity in Zebrafish
Jing Liu, Scott C. Baraban
eNeuro 25 February 2019, 6 (1) ENEURO.0041-19.2019; DOI: 10.1523/ENEURO.0041-19.2019

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Share
Network Properties Revealed during Multi-Scale Calcium Imaging of Seizure Activity in Zebrafish
Jing Liu, Scott C. Baraban
eNeuro 25 February 2019, 6 (1) ENEURO.0041-19.2019; DOI: 10.1523/ENEURO.0041-19.2019
del.icio.us logo Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Significance Statement
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Acknowledgments
    • Footnotes
    • References
    • Synthesis
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • epilepsy
  • fast confocal
  • neuronal networks
  • synchronization
  • whole-brain imaging
  • zebrafish

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

New Research

  • Impaired AMPARs translocation into dendritic spines with motor skill learning in the Fragile X mouse model
  • Characterization of the Tau Interactome in Human Brain Reveals Isoform-Dependent Interaction with 14-3-3 Family Proteins
  • Taste-odor association learning alters the dynamics of intra-oral odor responses in the posterior piriform cortex of awake rats
Show more New Research

Disorders of the Nervous System

  • Characterization of the Tau Interactome in Human Brain Reveals Isoform-Dependent Interaction with 14-3-3 Family Proteins
  • Impaired AMPARs translocation into dendritic spines with motor skill learning in the Fragile X mouse model
  • Glycolytic System in Axons Supplement Decreased ATP Levels after Axotomy of the Peripheral Nerve
Show more Disorders of the Nervous System

Subjects

  • Disorders of the Nervous System

  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Latest Articles
  • Issue Archive
  • Blog
  • Browse by Topic

Information

  • For Authors
  • For the Media

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
  • Feedback
(eNeuro logo)
(SfN logo)

Copyright © 2023 by the Society for Neuroscience.
eNeuro eISSN: 2373-2822

The ideas and opinions expressed in eNeuro do not necessarily reflect those of SfN or the eNeuro Editorial Board. Publication of an advertisement or other product mention in eNeuro should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in eNeuro.