Epileptic seizures from abnormal networks: why some seizures defy predictability

Seizure prediction - ready for a new era

Epilepsy is a common disorder characterized by recurrent seizures. An overwhelming majority of people with epilepsy regard the unpredictability of seizures as a major issue. More than 30 years of international effort have been devoted to the prediction of seizures, aiming to remove the burden of unpredictability and to couple novel, time-specific treatment to seizure prediction technology. A highly influential review published in 2007 concluded that insufficient evidence indicated that seizures could be predicted. Since then, several advances have been made, including successful prospective seizure prediction using intracranial EEG in a small number of people in a trial of a real-time seizure prediction device. In this Review, we examine advances in the field, including EEG databases, seizure prediction competitions, the prospective trial mentioned and advances in our understanding of the mechanisms of seizures. We argue that these advances, together with statistical evaluations, set the stage for a resurgence in efforts towards the development of seizure prediction methodologies. We propose new avenues of investigation involving a synergy between mechanisms, models, data, devices and algorithms and refine the existing guidelines for the development of seizure prediction technology to instigate development of a solution that removes the burden of the unpredictability of seizures.

Coupled Oscillators Model of Hyperexcitable Neuroglial Networks

Glial populations within neuronal networks of the brain have recently gained much interest in the context of hyperexcitability and epilepsy. In this paper, we present an oscillator-based neuroglial model capable of generating Spontaneous Electrical Discharges (SEDs) in hyperexcitable conditions. The network is composed of 16 coupled Cognitive Rhythm Generators (CRGs), which are oscillator-based mathematical constructs previously described by our research team. CRGs are well-suited for modeling assemblies of excitable cells, and in this network, each represents one of the following populations: excitatory pyramidal cells, inhibitory interneurons, astrocytes, and microglia. We investigated various pathways leading to hyperexcitability, and our results suggest an important role for astrocytes and microglia in the generation of SEDs of various durations. Analysis of the resultant SEDs revealed two underlying duration distributions with differing properties. Particularly, short and long SEDs are associated with deterministic and random underlying processes, respectively. The mesoscale of this model makes it well-suited for (a) the elucidation of glia-related hypotheses in hyperexcitable conditions, (b) use as a testing platform for neuromodulation purposes, and (c) a hardware implementation for closed-loop neuromodulation.

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.

Impact of Neuronal Membrane Damage on the Local Field Potential in a Large-Scale Simulation of Cerebral Cortex

Within multiscale brain dynamics, the structure–function relationship between cellular changes at a lower scale and coordinated oscillations at a higher scale is not well understood. This relationship may be particularly relevant for understanding functional impairments after a mild traumatic brain injury (mTBI) when current neuroimaging methods do not reveal morphological changes to the brain common in moderate to severe TBI such as diffuse axonal injury or gray matter lesions. Here, we created a physiology-based model of cerebral cortex using a publicly released modeling framework (GEneral NEural SImulation System) to explore the possibility that performance deficits characteristic of blast-induced mTBI may reflect dysfunctional, local network activity influenced by microscale neuronal damage at the cellular level. We operationalized microscale damage to neurons as the formation of pores on the neuronal membrane based on research using blast paradigms, and in our model, pores were simulated by a change in membrane conductance. We then tracked changes in simulated electrical activity. Our model contained 585 simulated neurons, comprised of 14 types of cortical and thalamic neurons each with its own compartmental morphology and electrophysiological properties. Comparing the functional activity of neurons before and after simulated damage, we found that simulated pores in the membrane reduced both action potential generation and local field potential (LFP) power in the 1–40 Hz range of the power spectrum. Furthermore, the location of damage modulated the strength of these effects: pore formation on simulated axons reduced LFP power more strongly than did pore formation on the soma and the dendrites. These results indicate that even small amounts of cellular damage can negatively impact functional activity of larger scale oscillations, and our findings suggest that multiscale modeling provides a promising avenue to elucidate these relationships.

A Phase-Locked Loop Epilepsy Network Emulator

Most seizure forecasting employs statistical learning techniques that lack a representation of the network interactions that give rise to seizures. We present an epilepsy network emulator (ENE) that uses a network of interconnected phase-locked loops (PLLs) to model synchronous, circuit-level oscillations between electrocorticography (ECoG) electrodes. Using ECoG data from a canine-epilepsy model (Davis et al. 2011) and a physiological entropy measure (approximate entropy or ApEn, Pincus 1995), we demonstrate the entropy of the emulator phases increases dramatically during ictal periods across all ECoG recording sites and across all animals in the sample. Further, this increase precedes the observable voltage spikes that characterize seizure activity in the ECoG data. These results suggest that the ENE is sensitive to phase-domain information in the neural circuits measured by ECoG and that an increase in the entropy of this measure coincides with increasing likelihood of seizure activity. Understanding this unpredictable phase-domain electrical activity present in ECoG recordings may provide a target for seizure detection and feedback control.

Interaction between synaptic inhibition and glial-potassium dynamics leads to diverse seizure transition modes in biophysical models of human focal seizures

How focal seizures initiate and evolve in human neocortex remains a fundamental problem in neuroscience. Here, we use biophysical neuronal network models of neocortical patches to study how the interaction between inhibition and extracellular potassium ([K (+)] o ) dynamics may contribute to different types of focal seizures. Three main types of propagated focal seizures observed in recent intracortical microelectrode recordings in humans were modelled: seizures characterized by sustained (∼30-60 Hz) gamma local field potential (LFP) oscillations; seizures where the onset in the propagated site consisted of LFP spikes that later evolved into rhythmic (∼2-3 Hz) spike-wave complexes (SWCs); and seizures where a brief stage of low-amplitude fast-oscillation (∼10-20 Hz) LFPs preceded the SWC activity. Our findings are fourfold: (1) The interaction between elevated [K (+)] o (due to abnormal potassium buffering by glial cells) and the strength of synaptic inhibition plays a predominant role in shaping these three types of seizures. (2) Strengthening of inhibition leads to the onset of sustained narrowband gamma seizures. (3) Transition into SWC seizures is obtained either by the weakening of inhibitory synapses, or by a transient strengthening followed by an inhibitory breakdown (e.g. GABA depletion). This reduction or breakdown of inhibition among fast-spiking (FS) inhibitory interneurons increases their spiking activity and leads them eventually into depolarization block. Ictal spike-wave discharges in the model are then sustained solely by pyramidal neurons. (4) FS cell dynamics are also critical for seizures where the evolution into SWC activity is preceded by low-amplitude fast oscillations. Different levels of elevated [K (+)] o were important for transitions into and maintenance of sustained gamma oscillations and SWC discharges. Overall, our modelling study predicts that the interaction between inhibitory interneurons and [K (+)] o glial buffering under abnormal conditions may explain different types of ictal transitions and dynamics during propagated seizures in human focal epilepsy.

Plic-1, a new target in repressing epileptic seizure by regulation of GABAAR function in patients and a rat model of epilepsy

Plic-1 regulates GABAAR expression at synaptic sites during epileptic seizure. Plic-1 prolongs the seizure latency and reduces the seizure severity in epileptic rats. Plic-1 affects the inhibitory function by changing the mIPSCs and evoked IPSCs of the phasic GABA-ergic synaptic current.

Classifying normal and abnormal status based on video recordings of epileptic patients

Based on video recordings of the movement of the patients with epilepsy, this paper proposed a human action recognition scheme to detect distinct motion patterns and to distinguish the normal status from the abnormal status of epileptic patients. The scheme first extracts local features and holistic features, which are complementary to each other. Afterwards, a support vector machine is applied to classification. Based on the experimental results, this scheme obtains a satisfactory classification result and provides a fundamental analysis towards the human-robot interaction with socially assistive robots in caring the patients with epilepsy (or other patients with brain disorders) in order to protect them from injury.

Characterization of seizure-like events recorded in vivo in a mouse model of Rett syndrome

Rett syndrome is a neurodevelopmental disorder caused by mutations in the X-linked gene encoding methyl-CpG-binding protein 2 (MECP2). Spontaneous recurrent discharge episodes are displayed in Rett-related seizures as in other types of epilepsies. The aim of this paper is to investigate the seizure-like event (SLE) and inter-SLE states in a female MeCP2-deficient mouse model of Rett syndrome and compare them to those found in other spontaneous recurrent epilepsy models. The study was performed on a small population of female MeCP2-deficient mice using telemetric local field potential (LFP) recordings over a 24 h period. Durations of SLEs and inter-SLEs were extracted using a rule-based automated SLE detection system for both daytime and nighttime, as well as high and low power levels of the delta frequency range (0.5-4 Hz) of the recorded LFPs. The results suggest SLE occurrences are not influenced by circadian rhythms, but had a significantly greater association with delta power. Investigating inter-SLE and SLE states by fitting duration histograms to the gamma distribution showed that SLE initiation and termination were associated with random and deterministic mechanisms, respectively. These findings when compared to reported studies on epilepsy suggest that Rett-related seizures share many similarities with absence epilepsy.

Statistical analysis of sleep spindle occurrences

Spindles - a hallmark of stage II sleep - are a transient oscillatory phenomenon in the EEG believed to reflect thalamocortical activity contributing to unresponsiveness during sleep. Currently spindles are often classified into two classes: fast spindles, with a frequency of around 14 Hz, occurring in the centro-parietal region; and slow spindles, with a frequency of around 12 Hz, prevalent in the frontal region. Here we aim to establish whether the spindle generation process also exhibits spatial heterogeneity. Electroencephalographic recordings from 20 subjects were automatically scanned to detect spindles and the time occurrences of spindles were used for statistical analysis. Gamma distribution parameters were fit to each inter-spindle interval distribution, and a modified Wald-Wolfowitz lag-1 correlation test was applied. Results indicate that not all spindles are generated by the same statistical process, but this dissociation is not spindle-type specific. Although this dissociation is not topographically specific, a single generator for all spindle types appears unlikely.

Predicting single-neuron activity in locally connected networks

The characterization of coordinated activity in neuronal populations has received renewed interest in the light of advancing experimental techniques that allow recordings from multiple units simultaneously. Across both in vitro and in vivo preparations, nearby neurons show coordinated responses when spontaneously active and when subject to external stimuli. Recent work (Truccolo, Hochberg, & Donoghue, 2010 ) has connected these coordinated responses to behavior, showing that small ensembles of neurons in arm-related areas of sensorimotor cortex can reliably predict single-neuron spikes in behaving monkeys and humans. We investigate this phenomenon using an analogous point process model, showing that in the case of a computational model of cortex responding to random background inputs, one is similarly able to predict the future state of a single neuron by considering its own spiking history, together with the spiking histories of randomly sampled ensembles of nearby neurons. This model exhibits realistic cortical architecture and displays bursting episodes in the two distinct connectivity schemes studied. We conjecture that the baseline predictability we find in these instances is characteristic of locally connected networks more broadly considered.