Aberrant fast spiking interneuronal activity precedes seizure transitions in humans.

There is active debate regarding how GABAergic function changes during seizure initiation and propagation, and whether interneuronal activity drives or impedes the pathophysiology. Here, we track cell-type specific firing during spontaneous human seizures to identify neocortical mechanisms of inhibitory failure. Fast-spiking interneuron activity was maximal over 1 second before equivalent excitatory increases, and showed transitions to out-of-phase firing prior to local tissue becoming incorporated into the seizure-driving territory. Using computational modeling, we linked this observation to transient saturation block as a precursor to seizure invasion, as supported by multiple lines of evidence in the patient data. We propose that transient blocking of inhibitory firing due to selective fast-spiking interneuron saturation-resulting from intense excitatory synaptic drive-is a novel mechanism that contributes to inhibitory failure, allowing seizure propagation.

A Novel Quantitative Metric Based on a Complete and Unique Characterization of Neural Network Activity: 4D Shannon's Entropy.

The human brain comprises an intricate web of connections that generate complex neural networks capable of storing and processing information. This information depends on multiple factors, including underlying network structure, connectivity, and interactions; and thus, methods to characterize neural networks typically aim to unravel and interpret a combination of these factors. Here, we present four-dimensional (4D) Shannon's entropy, a novel quantitative metric of network activity based on the Triple Correlation Uniqueness (TCU) theorem. Triple correlation, which provides a complete and unique characterization of the network, relates three nodes separated by up to four spatiotemporal lags. Here, we evaluate the 4D entropy from the spatiotemporal lag probability distribution function (PDF) of the network activity's triple correlation. Given a spike raster, we compute triple correlation by iterating over time and space. Summing the contributions to the triple correlation over each of the spatial and temporal lag combinations generates a unique 4D spatiotemporal lag distribution, from which we estimate a PDF and compute Shannon's entropy. To outline our approach, we first compute 4D Shannon's entropy from feedforward motif-class patterns in a simulated spike raster. We then apply this methodology to spiking activity recorded from rat cortical cultures to compare our results to previously published results of pairwise (2D) correlated spectral entropy over time. We find that while first- and second-order metrics of activity (spike rate and cross-correlation) show agreement with previously published results, our 4D entropy computation (which also includes third-order interactions) reveals a greater depth of underlying network organization compared to published pairwise entropy. Ultimately, because our approach is based on the TCU, we propose that 4D Shannon's entropy is a more complete tool for neural network characterization.

Emerging Activity Patterns and Synaptogenesis in Dissociated Hippocampal Cultures.

Cultures of dissociated hippocampal neurons display a stereotypical development of network activity patterns within the first three weeks of maturation. During this process, network connections develop and the associated spiking patterns range from increasing levels of activity in the first two weeks to regular bursting activity in the third week of maturation. Characterization of network structure is important to examine the mechanisms underlying the emergent functional organization of neural circuits. To accomplish this, confocal microscopy techniques have been used and several automated synapse quantification algorithms based on (co)localization of synaptic structures have been proposed recently. However, these approaches suffer from the arbitrary nature of intensity thresholding and the lack of correction for random-chance colocalization. To address this problem, we developed and validated an automated synapse quantification algorithm that requires minimal operator intervention. Next, we applied our approach to quantify excitatory and inhibitory synaptogenesis using confocal images of dissociated hippocampal neuronal cultures captured at 5, 8, 14 and 20 days in vitro, the time period associated with the development of distinct neuronal activity patterns. As expected, we found that synaptic density increased with maturation, coinciding with increasing spiking activity in the network. Interestingly, the third week of the maturation exhibited a reduction in excitatory synaptic density suggestive of synaptic pruning that coincided with the emergence of regular bursting activity in the network.

Spatiotemporal spike-centered averaging reveals symmetry of temporal and spatial components of the spike-LFP relationship during human focal seizures.

The electrographic manifestation of neural activity can reflect the relationship between the faster action potentials of individual neurons and the slower fluctuations of the local field potential (LFP). This relationship is typically examined in the temporal domain using the spike-triggered average. In this study, we add a spatial component to this relationship. Here we first derive a theoretical model of the spike-LFP relationship across a macroelectrode. This mathematical derivation showed a special symmetry in the spike-LFP relationship wherein a sinc function in the temporal domain predicts a sinc function in the spatial domain. We show that this theoretical result is observed in a real-world system by characterizing the spike-LFP relationship using microelectrode array (MEA) recordings of human focal seizures. To do this, we present a approach, termed the spatiotemporal spike-centered average (st-SCA), that allows for visualization of the spike-LFP relationship in both the temporal and spatial domains. We applied this method to 25 MEA recordings obtained from seven patients with pharmacoresistant focal epilepsy. Of the five patients with MEAs implanted in recruited territory, three exhibited spatiotemporal patterns consistent with a sinc function, and two exhibited spatiotemporal patterns resembling deep wells of excitation. These results suggest that in some cases characterization of the spike-LFP relationship in the temporal domain is sufficient to predict the underlying spatial pattern. Finally, we discuss the biological interpretation of these findings and propose that the sinc function may reflect the role of mid-range excitatory connections during seizure activity.

Third-order motifs are sufficient to fully and uniquely characterize spatiotemporal neural network activity.

Neuroscientific analyses balance between capturing the brain's complexity and expressing that complexity in meaningful and understandable ways. Here we present a novel approach that fully characterizes neural network activity and does so by uniquely transforming raw signals into easily interpretable and biologically relevant metrics of network behavior. We first prove that third-order (triple) correlation describes network activity in its entirety using the triple correlation uniqueness theorem. Triple correlation quantifies the relationships among three events separated by spatial and temporal lags, which are triplet motifs. Classifying these motifs by their event sequencing leads to fourteen qualitatively distinct motif classes that embody well-studied network behaviors including synchrony, feedback, feedforward, convergence, and divergence. Within these motif classes, the summed triple correlations provide novel metrics of network behavior, as well as being inclusive of commonly used analyses. We demonstrate the power of this approach on a range of networks with increasingly obscured signals, from ideal noiseless simulations to noisy experimental data. This approach can be easily applied to any recording modality, so existing neural datasets are ripe for reanalysis. Triple correlation is an accessible signal processing tool with a solid theoretical foundation capable of revealing previously elusive information within recordings of neural networks.

Manifestation of Hippocampal Interictal Discharges on Clinical Scalp EEG Recordings.

PURPOSE: Epileptiform activity limited to deep sources such as the hippocampus currently lacks reliable scalp correlates. Recent studies, however, have found that a subset of hippocampal interictal discharges may be associated with visible scalp signals, suggesting that some types of hippocampal activity may be monitored noninvasively. The purpose of this study is to characterize the relationship between these scalp waveforms and the underlying intracranial activity. METHODS: Paired intracranial and scalp EEG recordings obtained from 16 patients were used to identify hippocampal interictal discharges. Discharges were grouped by waveform shape, and spike-triggered averages of the intracranial and scalp signals were calculated for each group. Cross-correlation of intracranial and scalp spike-triggered averages was used to determine their temporal relationship, and topographic maps of the scalp were generated for each group. RESULTS: Cross-correlation of intracranial and scalp correlates resulted in two classes of scalp waveforms-those with and without time delays from the associated hippocampal discharges. Scalp signals with no delay showed topographies with a broad field with higher amplitudes on the side ipsilateral to the discharges and a left-right flip in polarity-observations consistent with the volume conduction of a single unilateral deep source. In contrast, scalp correlates with time lags showed rotational dynamics, suggesting synaptic propagation mechanisms. CONCLUSIONS: The temporal relationship between the intracranial and scalp signals suggests that both volume conduction and synaptic propagation contribute to these scalp manifestations. Furthermore, the topographic evolution of these scalp waveforms may be used to distinguish spikes that are limited to the hippocampus from those that travel to or engage other brain areas.

Digital reconstruction of infraslow activity in human intracranial ictal recordings using a deconvolution-based inverse filter.

Infraslow activity (ISA) is a biomarker that has recently become of interest in the characterization of seizure recordings. Recent data from a small number of studies have suggested that the epileptogenic zone may be identified by the presence of ISA. Investigation of low frequency activity in clinical seizure recordings, however, has been hampered by technical limitations. EEG systems necessarily include a high-pass filter early in the measurement chain to remove large artifactual drifts that can saturate recording elements such as the amplifier. This filter, unfortunately, attenuates legitimately seizure-related low frequencies, making ISA difficult to study in clinical EEG recordings. In this study, we present a deconvolution-based digital inverse filter that allows recovery of attenuated low frequency activity in intracranial recordings of temporal lobe epilepsy patients. First, we show that the unit impulse response (UIR) of an EEG system can be characterized by differentiation of the system's step response. As proof of method, we present several examples that show that the low frequency component of a high-pass filtered signal can be restored by deconvolution with the UIR. We then demonstrate that this method can be applied to biologically relevant signals including clinical EEG recordings obtained from seizure patients. Finally, we discuss how this method can be applied to study ISA to identify and assess the seizure onset zone.

DC shifts, high frequency oscillations, ripples and fast ripples in relation to the seizure onset zone.

Efforts to improve epilepsy surgery outcomes have led to increased interest in the study of electroencephalographic oscillations outside the conventional EEG bands. These include fast activity above the gamma band, known as high frequency oscillations (HFOs), and infraslow activity (ISA) below the delta band, sometimes referred to as direct current (DC) or ictal baseline shifts (IBS). HFOs in particular have been extensively studied as potential biomarkers for epileptogenic tissue in light of evidence showing that resection of brain tissue containing HFOs is associated with good surgical outcomes. Not all HFOs are conclusively pathological, however, as they can be recorded in nonepileptic tissue and induced by cognitive, visual, or motor tasks. Consequently, efforts to distinguish between pathological and physiological HFOs have identified several traits specific to pathological HFOs, such as coupling with interictal spikes, association with delta waves, and stereotypical morphologies. On the opposite end of the EEG spectrum, sub-delta oscillations have been shown to co-localize with the seizure onset zones (SOZ) and appear in a narrower spatial distribution than activity in the conventional EEG frequency bands. In this report, we review studies that implicate HFOs and ISA in ictogenesis and discuss current limitations such as inter-observer variability and poor standardization of recording techniques. Furthermore, we propose that HFOs and ISA should be analyzed in addition to activity in the conventional EEG band during intracranial presurgical EEG monitoring to identify the best possible surgical margin.

Role of paroxysmal depolarization in focal seizure activity.

We analyze the role of inhibition in sustaining focal epileptic seizure activity. We review ongoing seizure activity at the mesoscopic scale that can be observed with microelectrode arrays as well as at the macroscale of standard clinical EEG. We provide clinical, experimental, and modeling data to support the hypothesis that paroxysmal depolarization (PD) is a critical component of the ictal machinery. We present dual-patch recordings in cortical cultures showing reduced synaptic transmission associated with presynaptic occurrence of PD, and we find that the PD threshold is cell size related. We further find evidence that optically evoked PD activity in parvalbumin neurons can promote propagation of neuronal excitation in neocortical networks in vitro. Spike sorting results from microelectrode array measurements around ictal wave propagation in human focal seizures demonstrate a strong increase in putative inhibitory firing with an approaching excitatory wave, followed by a sudden reduction of firing at passage. At the macroscopic level, we summarize evidence that this excitatory ictal wave activity is strongly correlated with oscillatory activity across a centimeter-sized cortical network. We summarize Wilson-Cowan-type modeling showing how inhibitory function is crucial for this behavior. Our findings motivated us to develop a network motif of neurons in silico, governed by a reduced version of the Hodgkin-Huxley formalism, to show how feedforward, feedback, PD, and local failure of inhibition contribute to observed dynamics across network scales. The presented multidisciplinary evidence suggests that the PD not only is a cellular marker or epiphenomenon but actively contributes to seizure activity.NEW & NOTEWORTHY We present mechanisms of ongoing focal seizures across meso- and macroscales of microelectrode array and standard clinical recordings, respectively. We find modeling, experimental, and clinical evidence for a dual role of inhibition across these scales: local failure of inhibition allows propagation of a mesoscopic ictal wave, whereas inhibition elsewhere remains intact and sustains macroscopic oscillatory activity. We present evidence for paroxysmal depolarization as a mechanism behind this dual role of inhibition in shaping ictal activity.

The Relationship Between Ictal Multi-Unit Activity and the Electrocorticogram.

During neocortical seizures in patients with epilepsy, microelectrode array recordings from the ictal core show a strong correlation between the fast, cellular spiking activities and the low-frequency component of the potential field, reflected in the electrocorticogram (ECoG). Here, we model the relationship between the cellular spike activity and this low-frequency component as the input and output signals of a linear time invariant system. Our approach is based on the observation that this relationship can be characterized by a so-called sinc function, the unit impulse response of an ideal (brick-wall) filter. Accordingly, using a brick-wall filter, we are able to convert ictal cellular spike inputs into an output that significantly correlates with the observed seizure activity in the ECoG (r = 0.40 - 0.56,p < 0.01) , while ECoG recordings of subsequent seizures within patients also show significant, but lower, correlations (r = 0.10 - 0.30,p < 0.01) . Furthermore, we can produce seizure-like output signals using synthetic spike trains with ictal properties. We propose a possible physiological mechanism to explain the observed properties associated with an ideal filter, and discuss the potential use of our approach for the evaluation of anticonvulsant strategies.

Cross-scale effects of neural interactions during human neocortical seizure activity.

Small-scale neuronal networks may impose widespread effects on large network dynamics. To unravel this relationship, we analyzed eight multiscale recordings of spontaneous seizures from four patients with epilepsy. During seizures, multiunit spike activity organizes into a submillimeter-sized wavefront, and this activity correlates significantly with low-frequency rhythms from electrocorticographic recordings across a 10-cm-sized neocortical network. Notably, this correlation effect is specific to the ictal wavefront and is absent interictally or from action potential activity outside the wavefront territory. To examine the multiscale interactions, we created a model using a multiscale, nonlinear system and found evidence for a dual role for feedforward inhibition in seizures: while inhibition at the wavefront fails, allowing seizure propagation, feedforward inhibition of the surrounding centimeter-scale networks is activated via long-range excitatory connections. Bifurcation analysis revealed that distinct dynamical pathways for seizure termination depend on the surrounding inhibition strength. Using our model, we found that the mesoscopic, local wavefront acts as the forcing term of the ictal process, while the macroscopic, centimeter-sized network modulates the oscillatory seizure activity.

Multiscale Aspects of Generation of High-Gamma Activity during Seizures in Human Neocortex.

High-gamma (HG; 80-150 Hz) activity in macroscopic clinical records is considered a marker for critical brain regions involved in seizure initiation; it is correlated with pathological multiunit firing during neocortical seizures in the seizure core, an area identified by correlated multiunit spiking and low frequency seizure activity. However, the effects of the spatiotemporal dynamics of seizure on HG power generation are not well understood. Here, we studied HG generation and propagation, using a three-step, multiscale signal analysis and modeling approach. First, we analyzed concurrent neuronal and microscopic network HG activity in neocortical slices from seven intractable epilepsy patients. We found HG activity in these networks, especially when neurons displayed paroxysmal depolarization shifts and network activity was highly synchronized. Second, we examined HG activity acquired with microelectrode arrays recorded during human seizures (n = 8). We confirmed the presence of synchronized HG power across microelectrode records and the macroscale, both specifically associated with the core region of the seizure. Third, we used volume conduction-based modeling to relate HG activity and network synchrony at different network scales. We showed that local HG oscillations require high levels of synchrony to cross scales, and that this requirement is met at the microscopic scale, but not within macroscopic networks. Instead, we present evidence that HG power at the macroscale may result from harmonics of ongoing seizure activity. Ictal HG power marks the seizure core, but the generating mechanism can differ across spatial scales.

Network burst activity in hippocampal neuronal cultures: the role of synaptic and intrinsic currents.

The goal of this work was to define the contributions of intrinsic and synaptic mechanisms toward spontaneous network-wide bursting activity, observed in dissociated rat hippocampal cell cultures. This network behavior is typically characterized by short-duration bursts, separated by order of magnitude longer interburst intervals. We hypothesize that while short-timescale synaptic processes modulate spectro-temporal intraburst properties and network-wide burst propagation, much longer timescales of intrinsic membrane properties such as persistent sodium (Nap) currents govern burst onset during interburst intervals. To test this, we used synaptic receptor antagonists picrotoxin, 6-cyano-7-nitroquinoxaline-2,3-dione (CNQX), and 3-(2-carboxypiperazine-4-yl)propyl-1-phosphonate (CPP) to selectively block GABAA, AMPA, and NMDA receptors and riluzole to selectively block Nap channels. We systematically compared intracellular activity (recorded with patch clamp) and network activity (recorded with multielectrode arrays) in eight different synaptic connectivity conditions: GABAA + NMDA + AMPA, NMDA + AMPA, GABAA + AMPA, GABAA + NMDA, AMPA, NMDA, GABAA, and all receptors blocked. Furthermore, we used mixed-effects modeling to quantify the aforementioned independent and interactive synaptic receptor contributions toward spectro-temporal burst properties including intraburst spike rate, burst activity index, burst duration, power in the local field potential, network connectivity, and transmission delays. We found that blocking intrinsic Nap currents completely abolished bursting activity, demonstrating their critical role in burst onset within the network. On the other hand, blocking different combinations of synaptic receptors revealed that spectro-temporal burst properties are uniquely associated with synaptic functionality and that excitatory connectivity is necessary for the presence of network-wide bursting. In addition to confirming the critical contribution of direct excitatory effects, mixed-effects modeling also revealed distinct combined (nonlinear) contributions of excitatory and inhibitory synaptic activity to network bursting properties.

Wilson-Cowan Equations for Neocortical Dynamics.

In 1972-1973 Wilson and Cowan introduced a mathematical model of the population dynamics of synaptically coupled excitatory and inhibitory neurons in the neocortex. The model dealt only with the mean numbers of activated and quiescent excitatory and inhibitory neurons, and said nothing about fluctuations and correlations of such activity. However, in 1997 Ohira and Cowan, and then in 2007-2009 Buice and Cowan introduced Markov models of such activity that included fluctuation and correlation effects. Here we show how both models can be used to provide a quantitative account of the population dynamics of neocortical activity.We first describe how the Markov models account for many recent measurements of the resting or spontaneous activity of the neocortex. In particular we show that the power spectrum of large-scale neocortical activity has a Brownian motion baseline, and that the statistical structure of the random bursts of spiking activity found near the resting state indicates that such a state can be represented as a percolation process on a random graph, called directed percolation.Other data indicate that resting cortex exhibits pair correlations between neighboring populations of cells, the amplitudes of which decay slowly with distance, whereas stimulated cortex exhibits pair correlations which decay rapidly with distance. Here we show how the Markov model can account for the behavior of the pair correlations.Finally we show how the 1972-1973 Wilson-Cowan equations can account for recent data which indicates that there are at least two distinct modes of cortical responses to stimuli. In mode 1 a low intensity stimulus triggers a wave that propagates at a velocity of about 0.3 m/s, with an amplitude that decays exponentially. In mode 2 a high intensity stimulus triggers a larger response that remains local and does not propagate to neighboring regions.

Modeling focal epileptic activity in the Wilson-cowan model with depolarization block.

UNLABELLED: Measurements of neuronal signals during human seizure activity and evoked epileptic activity in experimental models suggest that, in these pathological states, the individual nerve cells experience an activity driven depolarization block, i.e. they saturate. We examined the effect of such a saturation in the Wilson-Cowan formalism by adapting the nonlinear activation function; we substituted the commonly applied sigmoid for a Gaussian function. We discuss experimental recordings during a seizure that support this substitution. Next we perform a bifurcation analysis on the Wilson-Cowan model with a Gaussian activation function. The main effect is an additional stable equilibrium with high excitatory and low inhibitory activity. Analysis of coupled local networks then shows that such high activity can stay localized or spread. Specifically, in a spatial continuum we show a wavefront with inhibition leading followed by excitatory activity. We relate our model simulations to observations of spreading activity during seizures. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s13408-015-0019-4) contains supplementary material 1.

Electrical control of epilepsy.

Epilepsy afflicts approximately 1-2% of the world's population. The mainstay therapy for treating the chronic recurrent seizures that are emblematic of epilepsy are drugs that manipulate levels of neuronal excitability in the brain. However, approximately one-third of all epilepsy patients get little to no clinical relief from this therapeutic regimen. The use of electrical stimulation in many forms to treat drug-refractory epilepsy has grown markedly over the past few decades, with some devices and protocols being increasingly used as standard clinical treatment. This article seeks to review the fundamental modes of applying electrical stimulation-from the noninvasive to the nominally invasive to deep brain stimulation-for the control of seizures in epileptic patients. Therapeutic practices from the commonly deployed clinically to the experimental are discussed to provide an overview of the innovative neural engineering approaches being explored to treat this difficult disease.

Large-scale modeling of epileptic seizures: scaling properties of two parallel neuronal network simulation algorithms.

Our limited understanding of the relationship between the behavior of individual neurons and large neuronal networks is an important limitation in current epilepsy research and may be one of the main causes of our inadequate ability to treat it. Addressing this problem directly via experiments is impossibly complex; thus, we have been developing and studying medium-large-scale simulations of detailed neuronal networks to guide us. Flexibility in the connection schemas and a complete description of the cortical tissue seem necessary for this purpose. In this paper we examine some of the basic issues encountered in these multiscale simulations. We have determined the detailed behavior of two such simulators on parallel computer systems. The observed memory and computation-time scaling behavior for a distributed memory implementation were very good over the range studied, both in terms of network sizes (2,000 to 400,000 neurons) and processor pool sizes (1 to 256 processors). Our simulations required between a few megabytes and about 150 gigabytes of RAM and lasted between a few minutes and about a week, well within the capability of most multinode clusters. Therefore, simulations of epileptic seizures on networks with millions of cells should be feasible on current supercomputers.

Resonance in neocortical neurons and networks.

Neocortical networks produce oscillations that often correspond to characteristic physiological or pathological patterns. However, the mechanisms underlying the generation of and the transitions between such oscillatory states remain poorly understood. In this study, we examined resonance in mouse layer V neocortical pyramidal neurons. To accomplish this, we employed standard electrophysiology to describe cellular resonance parameters. Bode plot analysis revealed a range of resonance magnitude values in layer V neurons and demonstrated that both magnitude and phase response characteristics of layer V neocortical pyramidal neurons are modulated by changes in the extracellular environment. Specifically, increased resonant frequencies and total inductive areas were observed at higher extracellular potassium concentrations and more hyperpolarised membrane potentials. Experiments using pharmacological agents suggested that current through hyperpolarization-activated cyclic nucleotide-gated channels (I(h) ) acts as the primary driver of resonance in these neurons, with other potassium currents, such as A-type potassium current and delayed-rectifier potassium current (Kv1.4 and Kv1.1, respectively), contributing auxiliary roles. The persistent sodium current was also shown to play a role in amplifying the magnitude of resonance without contributing significantly to the phase response. Although resonance effects in individual neurons are small, their properties embedded in large networks may significantly affect network behavior and may have potential implications for pathological processes.