Human seizures couple across spatial scales through travelling wave dynamics

Explaining slow seizure propagation with white matter tractography

Epileptic seizures recorded with stereo-EEG can take a fraction of a second or several seconds to propagate from one region to another. What explains such propagation patterns? We combine tractography and stereo-EEG to determine the relationship between seizure propagation and the white matter architecture and to describe seizure propagation mechanisms. Patient-specific spatiotemporal seizure propagation maps were combined with tractography from diffusion imaging of matched subjects from the Human Connectome Project. The onset of seizure activity was marked on a channel-by-channel basis by two board-certified neurologists for all channels involved in the seizure. We measured the tract connectivity (number of tracts) between regions-of-interest pairs among the seizure onset zone, regions of seizure spread and non-involved regions. We also investigated how tract-connected the seizure onset zone is to regions of early seizure spread compared with regions of late spread. Comparisons were made after correcting for differences in distance. Sixty-nine seizures were marked across 26 patients with drug-resistant epilepsy; 11 were seizure free after surgery (Engel IA) and 15 were not (Engel IB-Engel IV). The seizure onset zone was more tract-connected to regions of seizure spread than to non-involved regions (P < 0.0001); however, regions of seizure spread were not differentially tract-connected to other regions of seizure spread compared with non-involved regions. In seizure-free patients only, regions of seizure spread were more tract-connected to the seizure onset zone than to other regions of spread (P < 0.0001). Over the temporal evolution of a seizure, the seizure onset zone was significantly more tract-connected to regions of early spread compared with regions of late spread in seizure-free patients only (P < 0.0001). By integrating information on structure, we demonstrate that seizure propagation is likely to be mediated by white matter tracts. The pattern of connectivity between seizure onset zone, regions of spread and non-involved regions demonstrates that the onset zone might be largely responsible for seizures propagating throughout the brain, rather than seizures propagating to intermediate points, from which further propagation takes place. Our findings also suggest that seizure propagation over seconds might be the result of a continuous bombardment of action potentials from the seizure onset zone to regions of spread. In non-seizure-free patients, the paucity of tracts from the presumed seizure onset zone to regions of spread suggests that the onset zone was missed. Fully understanding the structure-propagation relationship might eventually provide insight into selecting the correct targets for epilepsy surgery.

Focal seizures induce spatiotemporally organized spiking activity in the human cortex

Epileptic seizures are debilitating because of the clinical symptoms they produce. These symptoms, in turn, may stem directly from disruptions in neural coding. Recent evidence has suggested that the specific temporal order, or sequence, of spiking across a population of cortical neurons may encode information. Here, we investigate how seizures disrupt neuronal spiking sequences in the human brain by recording multi-unit activity from the cerebral cortex in five male participants undergoing monitoring for seizures. We find that pathological discharges during seizures are associated with bursts of spiking activity across a population of cortical neurons. These bursts are organized into highly consistent and stereotyped temporal sequences. As the seizure evolves, spiking sequences diverge from the sequences observed at baseline and become more spatially organized. The direction of this spatial organization matches the direction of the ictal discharges, which spread over the cortex as traveling waves. Our data therefore suggest that seizures can entrain cortical spiking sequences by changing the spatial organization of neuronal firing, providing a possible mechanism by which seizures create symptoms.

Temperature effects on neuronal synchronization in seizures

We present a computational model of networked neurons developed to study the effect of temperature on neuronal synchronization in the brain in association with seizures. The network consists of a set of chaotic bursting neurons surrounding a core tonic neuron in a square lattice with periodic boundary conditions. Each neuron is reciprocally coupled to its four nearest neighbors via temperature dependent gap junctions. Incorporating temperature in the gap junctions makes the coupling stronger when temperature rises, resulting in higher likelihood for synchrony in the network. Raising the temperature eventually makes the network elicit waves of synchronization in circular ripples that propagate from the center outwardly. We suggest this process as a possible underlying mechanism for seizures induced by elevated brain temperatures.

Neuronal tuning and population representations of shape and category in human visual cortex

Object recognition and categorization are essential cognitive processes which engage considerable neural resources in the human ventral visual stream. However, the tuning properties of human ventral stream neurons for object shape and category are virtually unknown. We performed large-scale recordings of spiking activity in human Lateral Occipital Complex in response to stimuli in which the shape dimension was dissociated from the category dimension. Consistent with studies in nonhuman primates, the neuronal representations were primarily shape-based, although we also observed category-like encoding for images of animals. Surprisingly, linear decoders could reliably classify stimulus category even in data sets that were entirely shape-based. In addition, many recording sites showed an interaction between shape and category tuning. These results represent a detailed study on shape and category coding at the neuronal level in the human ventral visual stream, furnishing essential evidence that reconciles human imaging and macaque single-cell studies.

Virtual brain twins: from basic neuroscience to clinical use

Virtual brain twins are personalized, generative and adaptive brain models based on data from an individual's brain for scientific and clinical use. After a description of the key elements of virtual brain twins, we present the standard model for personalized whole-brain network models. The personalization is accomplished using a subject's brain imaging data by three means: (1) assemble cortical and subcortical areas in the subject-specific brain space; (2) directly map connectivity into the brain models, which can be generalized to other parameters; and (3) estimate relevant parameters through model inversion, typically using probabilistic machine learning. We present the use of personalized whole-brain network models in healthy ageing and five clinical diseases: epilepsy, Alzheimer's disease, multiple sclerosis, Parkinson's disease and psychiatric disorders. Specifically, we introduce spatial masks for relevant parameters and demonstrate their use based on the physiological and pathophysiological hypotheses. Finally, we pinpoint the key challenges and future directions.

Multi-scale modelling of the epileptic brain: advantages of computational therapy exploration

Objective: Epilepsy is a complex disease spanning across multiple scales, from ion channels in neurons to neuronal circuits across the entire brain. Over the past decades, computational models have been used to describe the pathophysiological activity of the epileptic brain from different aspects. Traditionally, each computational model can aid in optimizing therapeutic interventions, therefore, providing a particular view to design strategies for treating epilepsy. As a result, most studies are concerned with generating specific models of the epileptic brain that can help us understand the certain machinery of the pathological state. Those specific models vary in complexity and biological accuracy, with system-level models often lacking biological details.Approach: Here, we review various types of computational model of epilepsy and discuss their potential for different therapeutic approaches and scenarios, including drug discovery, surgical strategies, brain stimulation, and seizure prediction. We propose that we need to consider an integrated approach with a unified modelling framework across multiple scales to understand the epileptic brain. Our proposal is based on the recent increase in computational power, which has opened up the possibility of unifying those specific epileptic models into simulations with an unprecedented level of detail.Main results: A multi-scale epilepsy model can bridge the gap between biologically detailed models, used to address molecular and cellular questions, and brain-wide models based on abstract models which can account for complex neurological and behavioural observations.Significance: With these efforts, we move toward the next generation of epileptic brain models capable of connecting cellular features, such as ion channel properties, with standard clinical measures such as seizure severity.

Sequentially activated discrete modules appear as traveling waves in neuronal measurements with limited spatiotemporal sampling

Numerous studies have identified traveling waves in the cortex and suggested they play important roles in brain processing. These waves are most often measured using macroscopic methods that are unable to assess the local spiking activity underlying wave dynamics. Here, we investigated the possibility that waves may not be traveling at the single neuron scale. We first show that sequentially activating two discrete brain areas can appear as traveling waves in EEG simulations. We next reproduce these results using an analytical model of two sequentially activated regions. Using this model, we were able to generate wave-like activity with variable directions, velocities, and spatial patterns, and to map the discriminability limits between traveling waves and modular sequential activations. Finally, we investigated the link between field potentials and single neuron excitability using large-scale measurements from turtle cortex ex vivo. We found that while field potentials exhibit wave-like dynamics, the underlying spiking activity was better described by consecutively activated spatially adjacent groups of neurons. Taken together, this study suggests caution when interpreting phase delay measurements as continuously propagating wavefronts in two different spatial scales. A careful distinction between modular and wave excitability profiles across scales will be critical for understanding the nature of cortical computations.

Pan-cortical electrophysiologic changes underlying attention

We previously reported that pan-cortical effects occur when cognitive tasks end afterdischarges. For this report, we analyzed wavelet cross-coherence changes during cognitive tasks used to terminate afterdischarges studying multiple time segments and multiple groups of inter-electrode-con distances. We studied 12 patients with intractable epilepsy, with 970 implanted electrode contacts, and 39,871 electrode contact combinations. When cognitive tasks ended afterdischarges, coherence varied similarly across the cortex throughout the tasks, but there were gradations with time, distance, and frequency: (1) They tended to progressively decrease relative to baseline with time and then to increase toward baseline when afterdischarges ended. (2) During most time segments, decreases from baseline were largest for the closest inter-contact distances, moderate for intermediate inter-contact distances, and smallest for the greatest inter-contact distances. With respect to our patients' intractable epilepsy, the changes found suggest that future therapies might treat regions beyond those closest to regions of seizure onset and treat later in a seizure's evolution. Similar considerations might apply to other disorders. Our findings also suggest that cognitive tasks can result in pan-cortical coherence changes that participate in underlying attention, perhaps complementing the better-known regional mechanisms.

Shape-changing electrode array for minimally invasive large-scale intracranial brain activity mapping

Large-scale brain activity mapping is important for understanding the neural basis of behaviour. Electrocorticograms (ECoGs) have high spatiotemporal resolution, bandwidth, and signal quality. However, the invasiveness and surgical risks of electrode array implantation limit its application scope. We developed an ultrathin, flexible shape-changing electrode array (SCEA) for large-scale ECoG mapping with minimal invasiveness. SCEAs were inserted into cortical surfaces in compressed states through small openings in the skull or dura and fully expanded to cover large cortical areas. MRI and histological studies on rats proved the minimal invasiveness of the implantation process and the high chronic biocompatibility of the SCEAs. High-quality micro-ECoG activities mapped with SCEAs from male rodent brains during seizures and canine brains during the emergence period revealed the spatiotemporal organization of different brain states with resolution and bandwidth that cannot be achieved using existing noninvasive techniques. The biocompatibility and ability to map large-scale physiological and pathological cortical activities with high spatiotemporal resolution, bandwidth, and signal quality in a minimally invasive manner offer SCEAs as a superior tool for applications ranging from fundamental brain research to brain-machine interfaces.

Rogue Waves: A Unifying Model of Human Seizure Propagation

Multiple Sources of Fast Traveling Waves During Human Seizures: Resolving a Controversy Schlafly ED, Marshall FA, Merricks EM, Eden UT, Cash SS, Schevon CA, Kramer MA. J Neurosci . 2022;42(36):6966-6982. doi:10.1523/JNEUROSCI.0338-22.2022 During human seizures organized waves of voltage activity rapidly sweep across the cortex. Two contradictory theories describe the source of these fast traveling waves: either a slowly advancing narrow region of multiunit activity (an ictal wavefront) or a fixed cortical location. Limited observations and different analyses prevent resolution of these incompatible theories. Here we address this disagreement by combining the methods and microelectrode array recordings (N = 11 patients, 2 females, N = 31 seizures) from previous human studies to analyze the traveling wave source. We find—inconsistent with both existing theories—a transient relationship between the ictal wavefront and traveling waves, and multiple stable directions of traveling waves in many seizures. Using a computational model that combines elements of both existing theories, we show that interactions between an ictal wavefront and fixed source reproduce the traveling wave dynamics observed in vivo. We conclude that combining both existing theories can generate the diversity of ictal traveling waves. Significance Statement: The source of voltage discharges that propagate across cortex during human seizures remains unknown. Two candidate theories exist, each proposing a different discharge source. Support for each theory consists of observations from a small number of human subject recordings, analyzed with separately developed methods. How the different, limited data and different analysis methods impact the evidence for each theory is unclear. To resolve these differences, we combine the unique, human microelectrode array recordings collected separately for each theory and analyze these combined data with a unified approach. We show that neither existing theory adequately describes the data. We then propose a new theory that unifies existing proposals and successfully reproduces the voltage discharge dynamics observed in vivo.

Microscale Neuronal Activity Collectively Drives Chaotic and Inflexible Dynamics at the Macroscale in Seizures

Neuronal activity propagates through the network during seizures, engaging brain dynamics at multiple scales. Such propagating events can be described through the avalanches framework, which can relate spatiotemporal activity at the microscale with global network properties. Interestingly, propagating avalanches in healthy networks are indicative of critical dynamics, where the network is organized to a phase transition, which optimizes certain computational properties. Some have hypothesized that the pathologic brain dynamics of epileptic seizures are an emergent property of microscale neuronal networks collectively driving the brain away from criticality. Demonstrating this would provide a unifying mechanism linking microscale spatiotemporal activity with emergent brain dysfunction during seizures. Here, we investigated the effect of drug-induced seizures on critical avalanche dynamics, using in vivo whole-brain two-photon imaging of GCaMP6s larval zebrafish (males and females) at single neuron resolution. We demonstrate that single neuron activity across the whole brain exhibits a loss of critical statistics during seizures, suggesting that microscale activity collectively drives macroscale dynamics away from criticality. We also construct spiking network models at the scale of the larval zebrafish brain, to demonstrate that only densely connected networks can drive brain-wide seizure dynamics away from criticality. Importantly, such dense networks also disrupt the optimal computational capacities of critical networks, leading to chaotic dynamics, impaired network response properties and sticky states, thus helping to explain functional impairments during seizures. This study bridges the gap between microscale neuronal activity and emergent macroscale dynamics and cognitive dysfunction during seizures.SIGNIFICANCE STATEMENT Epileptic seizures are debilitating and impair normal brain function. It is unclear how the coordinated behavior of neurons collectively impairs brain function during seizures. To investigate this we perform fluorescence microscopy in larval zebrafish, which allows for the recording of whole-brain activity at single-neuron resolution. Using techniques from physics, we show that neuronal activity during seizures drives the brain away from criticality, a regime that enables both high and low activity states, into an inflexible regime that drives high activity states. Importantly, this change is caused by more connections in the network, which we show disrupts the ability of the brain to respond appropriately to its environment. Therefore, we identify key neuronal network mechanisms driving seizures and concurrent cognitive dysfunction.

Interictal discharges in the human brain are travelling waves arising from an epileptogenic source

While seizure activity may be electrographically widespread, increasing evidence has suggested that ictal discharges may in fact represent travelling waves propagated from a focal seizure source. Interictal epileptiform discharges (IEDs) are an electrographic manifestation of excessive hypersynchronization of cortical activity that occur between seizures and are considered a marker of potentially epileptogenic tissue. The precise relationship between brain regions demonstrating IEDs and those involved in seizure onset, however, remains poorly understood. Here, we hypothesize that IEDs likewise reflect the receipt of travelling waves propagated from the same regions which give rise to seizures. Forty patients from our institution who underwent invasive monitoring for epilepsy, proceeded to surgery and had at least one year of follow-up were included in our study. Interictal epileptiform discharges were detected using custom software, validated by a clinical epileptologist. We show that IEDs reach electrodes in sequences with a consistent temporal ordering, and this ordering matches the timing of receipt of ictal discharges, suggesting that both types of discharges spread as travelling waves. We use a novel approach for localization of ictal discharges, in which time differences of discharge receipt at nearby electrodes are used to compute source location; similar algorithms have been used in acoustics and geophysics. We find that interictal discharges co-localize with ictal discharges. Moreover, interictal discharges tend to localize to the resection territory in patients with good surgical outcome and outside of the resection territory in patients with poor outcome. The seizure source may originate at, and also travel to, spatially distinct IED foci. Our data provide evidence that interictal discharges may represent travelling waves of pathological activity that are similar to their ictal counterparts, and that both ictal and interictal discharges emerge from common epileptogenic brain regions. Our findings have important clinical implications, as they suggest that seizure source localizations may be derived from interictal discharges, which are much more frequent than seizures.

Multilayer brain networks can identify the epileptogenic zone and seizure dynamics

Seizure generation, propagation, and termination occur through spatiotemporal brain networks. In this paper, we demonstrate the significance of large-scale brain interactions in high-frequency (80-200Hz) for the identification of the epileptogenic zone (EZ) and seizure evolution. To incorporate the continuity of neural dynamics, here we have modeled brain connectivity constructed from stereoelectroencephalography (SEEG) data during seizures using multilayer networks. After introducing a new measure of brain connectivity for temporal networks, named multilayer eigenvector centrality (mlEVC), we applied a consensus hierarchical clustering on the developed model to identify the EZ as a cluster of nodes with distinctive brain connectivity in the ictal period. Our algorithm could successfully predict electrodes inside the resected volume as EZ for 88% of participants, who all were seizure-free for at least 12 months after surgery. Our findings illustrated significant and unique desynchronization between EZ and the rest of the brain in the early to mid-seizure. We showed that aging and the duration of epilepsy intensify this desynchronization, which can be the outcome of abnormal neuroplasticity. Additionally, we illustrated that seizures evolve with various network topologies, confirming the existence of different epileptogenic networks in each patient. Our findings suggest not only the importance of early intervention in epilepsy but possible factors that correlate with disease severity. Moreover, by analyzing the propagation patterns of different seizures, we demonstrate the necessity of collecting sufficient data for identifying epileptogenic networks.

Localized chemogenetic silencing of inhibitory neurons: a novel mouse model of focal cortical epileptic activity

Focal cortical epilepsies are frequently refractory to available anticonvulsant drug therapies. One key factor contributing to this state is the limited availability of animal models that allow to reliably study focal cortical seizures and how they recruit surrounding brain areas in vivo. In this study, we selectively expressed the inhibitory chemogenetic receptor, hM4D, in GABAergic neurons in focal cortical areas using viral gene transfer. GABAergic silencing using Clozapine-N-Oxide (CNO) demonstrated reliable induction of local epileptiform events in the electroencephalogram signal of awake freely moving mice. Anesthetized mice experiments showed consistent induction of focal epileptiform-events in both the barrel cortex (BC) and the medial prefrontal cortex (mPFC), accompanied by high-frequency oscillations, a known characteristic of human seizures. Epileptiform-events showed propagation indication with favored propagation pathways: from the BC on 1 hemisphere to its counterpart and from the BC to the mPFC, but not vice-versa. Lastly, sensory whisker-pad stimulation evoked BC epileptiform events post-CNO, highlighting the potential use of this model in studying sensory-evoked seizures. Combined, our results show that targeted chemogenetic inhibition of GABAergic neurons using hM4D can serve as a novel, versatile, and reliable model of focal cortical epileptic activity suitable for systematically studying cortical ictogenesis in different cortical areas.

Coupled Dynamics of Stimulus-Evoked Gustatory Cortical and Basolateral Amygdalar Activity

Gustatory cortical (GC) single-neuron taste responses reflect taste quality and palatability in successive epochs. Ensemble analyses reveal epoch-to-epoch firing-rate changes in these responses to be sudden, coherent transitions. Such nonlinear dynamics suggest that GC is part of a recurrent network, producing these dynamics in concert with other structures. Basolateral amygdala (BLA), which is reciprocally connected to GC and central to hedonic processing, is a strong candidate partner for GC, in that BLA taste responses evolve on the same general clock as GC and because inhibition of activity in the BLA→GC pathway degrades the sharpness of GC transitions. These facts motivate, but do not test, our overarching hypothesis that BLA and GC act as a single, comodulated network during taste processing. Here, we provide just this test of simultaneous (BLA and GC) extracellular taste responses in female rats, probing the multiregional dynamics of activity to directly test whether BLA and GC responses contain coupled dynamics. We show that BLA and GC response magnitudes covary across trials and within single responses, and that changes in BLA-GC local field potential phase coherence are epoch specific. Such classic coherence analyses, however, obscure the most salient facet of BLA-GC coupling: sudden transitions in and out of the epoch known to be involved in driving gaping behavior happen near simultaneously in the two regions, despite huge trial-to-trial variability in transition latencies. This novel form of inter-regional coupling, which we show is easily replicated in model networks, suggests collective processing in a distributed neural network.SIGNIFICANCE STATEMENT There has been little investigation into real-time communication between brain regions during taste processing, a fact reflecting the dominant belief that taste circuitry is largely feedforward. Here, we perform an in-depth analysis of real-time interactions between GC and BLA in response to passive taste deliveries, using both conventional coherence metrics and a novel methodology that explicitly considers trial-to-trial variability and fast single-trial dynamics in evoked responses. Our results demonstrate that BLA-GC coherence changes as the taste response unfolds, and that BLA and GC specifically couple for the sudden transition into (and out of) the behaviorally relevant neural response epoch, suggesting (although not proving) that: (1) recurrent interactions subserve the function of the dyad as (2) a putative attractor network.

Temporal-lobe Epilepsy: Harmonic and Anharmonic Periodicity in Microeletrode Voltage

Temporal-lobe epilepsy in humans is often associated with widespread, synchronized neuron firing that co-occurs with traveling waves in local field potential. These traveling waves generate stochastic oscillations in a time series of microelectrode voltage, and previous work has deemed it informative for traveling-wave analysis to study the mean periodicity. This manuscript reveals that: a) mean voltage (i.e., traveling-wave periodicity) adequately explains the observed voltage periodicity only for a select few time intervals during seizure; and b) mean voltage has a 7 Hz cosine-series representation indicative of a nonlinear system response given alpha-rhythm input. The a) result implies that residual noise should be modelled explicitly, while b) motivates a departure from the conventional plane-wave modeling regime in source-localization efforts. The 7 Hz fundamental frequency is unsurprising given the relative transparency of the brain to 14 Hz alpha rhythms in neurophysiological diseases (14 Hz being a subharmonic frequency of the 7 Hz signal).

Multiple Sources of Fast Traveling Waves during Human Seizures: Resolving a Controversy

During human seizures, organized waves of voltage activity rapidly sweep across the cortex. Two contradictory theories describe the source of these fast traveling waves: either a slowly advancing narrow region of multiunit activity (an ictal wavefront) or a fixed cortical location. Limited observations and different analyses prevent resolution of these incompatible theories. Here we address this disagreement by combining the methods and microelectrode array recordings (N = 11 patients, 2 females, N = 31 seizures) from previous human studies to analyze the traveling wave source. We find, inconsistent with both existing theories, a transient relationship between the ictal wavefront and traveling waves, and multiple stable directions of traveling waves in many seizures. Using a computational model that combines elements of both existing theories, we show that interactions between an ictal wavefront and fixed source reproduce the traveling wave dynamics observed in vivo We conclude that combining both existing theories can generate the diversity of ictal traveling waves.SIGNIFICANCE STATEMENT The source of voltage discharges that propagate across cortex during human seizures remains unknown. Two candidate theories exist, each proposing a different discharge source. Support for each theory consists of observations from a small number of human subject recordings, analyzed with separately developed methods. How the different, limited data and different analysis methods impact the evidence for each theory is unclear. To resolve these differences, we combine the unique, human microelectrode array recordings collected separately for each theory and analyze these combined data with a unified approach. We show that neither existing theory adequately describes the data. We then propose a new theory that unifies existing proposals and successfully reproduces the voltage discharge dynamics observed in vivo.

Splitting of the magnetic encephalogram into «brain» and «non-brain» physiological signals based on the joint analysis of frequency-pattern functional tomograms and magnetic resonance images

The article considers the problem of dividing the encephalography data into two time series, that generated by the brain and that generated by other electrical sources located in the human head. The magnetic encephalograms and magnetic resonance images of the head were recorded in the Center for Neuromagnetism at NYU Grossman School of Medicine. Data obtained at McGill University and Montreal University were also used. Recordings were made in a magnetically shielded room and the gradiometers were designed to suppress external noise, making it possible to eliminate them from the data analysis. Magnetic encephalograms were analyzed by the method of functional tomography, based on the Fourier transform and on the solution of inverse problem for all frequencies. In this method, one spatial position is assigned to each frequency component. Magnetic resonance images of the head were evaluated to annotate the space to be included in the analysis. The included space was divided into two parts: «brain» and «non-brain». The frequency components were classified by the feature of their inclusion in one or the other part. The set of frequencies, designated as «brain», represented the partial spectrum of the brain signal, while the set of frequencies designated as «non-brain», represented the partial spectrum of the physiological noise produced by the head. Both partial spectra shared the same frequency band. From the partial spectra, a time series of the «brain» area signal and «non-brain» area head noise were reconstructed. Summary spectral power of the signal was found to be ten times greater than the noise. The proposed method makes it possible to analyze in detail both the signal and the noise components of the encephalogram and to filter the magnetic encephalogram.

Focal seizures are organized by feedback between neural activity and ion concentration changes

Human and animal EEG data demonstrate that focal seizures start with low-voltage fast activity, evolve into rhythmic burst discharges and are followed by a period of suppressed background activity. This suggests that processes with dynamics in the range of tens of seconds govern focal seizure evolution. We investigate the processes associated with seizure dynamics by complementing the Hodgkin-Huxley mathematical model with the physical laws that dictate ion movement and maintain ionic gradients. Our biophysically realistic computational model closely replicates the electrographic pattern of a typical human focal seizure characterized by low voltage fast activity onset, tonic phase, clonic phase and postictal suppression. Our study demonstrates, for the first time in silico, the potential mechanism of seizure initiation by inhibitory interneurons via the initial build-up of extracellular K+ due to intense interneuronal spiking. The model also identifies ionic mechanisms that may underlie a key feature in seizure dynamics, i.e., progressive slowing down of ictal discharges towards the end of seizure. Our model prediction of specific scaling of inter-burst intervals is confirmed by seizure data recorded in the whole guinea pig brain in vitro and in humans, suggesting that the observed termination pattern may hold across different species. Our results emphasize ionic dynamics as elementary processes behind seizure generation and indicate targets for new therapeutic strategies.

From pattern to process? Dual travelling waves, with contrasting propagation speeds, best describe a self-organised spatio-temporal pattern in population growth of a cyclic rodent

The dynamics of cyclic populations distributed in space result from the relative strength of synchronising influences and the limited dispersal of destabilising factors (activators and inhibitors), known to cause multi‐annual population cycles. However, while each of these have been well studied in isolation, there is limited empirical evidence of how the processes of synchronisation and activation–inhibition act together, largely owing to the scarcity of datasets with sufficient spatial and temporal scale and resolution. We assessed a variety of models that could be underlying the spatio‐temporal pattern, designed to capture both theoretical and empirical understandings of travelling waves using large‐scale (>35,000 km²), multi‐year (2011–2017) field monitoring data on abundances of common vole (Microtus arvalis), a cyclic agricultural rodent pest. We found most support for a pattern formed from the summation of two radial travelling waves with contrasting speeds that together describe population growth rates across the region.

Seizure classification with selected frequency bands and EEG montages: a Natural Language Processing approach

Individualized treatment is crucial for epileptic patients with different types of seizures. The differences among patients impact the drug choice as well as the surgery procedure. With the advance in machine learning, automatic seizure detection can ease the manual time-consuming and labor-intensive procedure for diagnose seizure in the clinical setting. In this paper, we present an electroencephalography (EEG) frequency bands (sub-bands) and montages selection (sub-zones) method for classifier training that exploits Natural Language Processing from individual patients' clinical report. The proposed approach is targeting for individualized treatment. We integrated the prior knowledge from patient's reports into the classifier-building process, mimicking the authentic thinking process of experienced neurologist's when diagnosing seizure using EEG. The keywords from clinical documents are mapped to the EEG data in terms of frequency bands and scalp EEG electrodes. The data of experiments are from the Temple University Hospital EEG seizure corpus, and the dataset is divided based on each group of patients with same seizure type and same recording electrode references. The classifier includes Random Forest, Support Vector Machine and Multi-Layer Perceptron. The classification performance indicates that competitive results can be achieve with a small portion of EEG the data. Using the sub-zones selection for Generalized Seizures (GNSZ) on all three electrodes, data are reduced by nearly 50% while the performance metrics remain at the same level with the whole frequency and zones. Moreover, when selecting by sub-zones and sub-bands together for GNSZ with Linked Ears reference, the data range reduced to 0.3% of whole range, and the performance deviates less than 3% from the results with whole range of data. Results show that using proposed approach may lead to more efficient implementations of the seizure classifier to be executed on power-efficient devices for long lasting real-time seizures detection.

Multimodal, Multiscale Insights into Hippocampal Seizures Enabled by Transparent, Graphene-Based Microelectrode Arrays

Hippocampal seizures are a defining feature of mesial temporal lobe epilepsy (MTLE). Area CA1 of the hippocampus is commonly implicated in the generation of seizures, which may occur because of the activity of endogenous cell populations or of inputs from other regions within the hippocampal formation. Simultaneously observing activity at the cellular and network scales in vivo remains challenging. Here, we present a novel technology for simultaneous electrophysiology and multicellular calcium imaging of CA1 pyramidal cells (PCs) in mice enabled by a transparent graphene-based microelectrode array (Gr MEA). We examine PC firing at seizure onset, oscillatory coupling, and the dynamics of the seizure traveling wave as seizures evolve. Finally, we couple features derived from both modalities to predict the speed of the traveling wave using bootstrap aggregated regression trees. Analysis of the most important features in the regression trees suggests a transition among states in the evolution of hippocampal seizures.

Consciousness is supported by near-critical slow cortical electrodynamics

Mounting evidence suggests that during conscious states, the electrodynamics of the cortex are poised near a critical point or phase transition and that this near-critical behavior supports the vast flow of information through cortical networks during conscious states. Here, we empirically identify a mathematically specific critical point near which waking cortical oscillatory dynamics operate, which is known as the edge-of-chaos critical point, or the boundary between stability and chaos. We do so by applying the recently developed modified 0-1 chaos test to electrocorticography (ECoG) and magnetoencephalography (MEG) recordings from the cortices of humans and macaques across normal waking, generalized seizure, anesthesia, and psychedelic states. Our evidence suggests that cortical information processing is disrupted during unconscious states because of a transition of low-frequency cortical electric oscillations away from this critical point; conversely, we show that psychedelics may increase the information richness of cortical activity by tuning low-frequency cortical oscillations closer to this critical point. Finally, we analyze clinical electroencephalography (EEG) recordings from patients with disorders of consciousness (DOC) and show that assessing the proximity of slow cortical oscillatory electrodynamics to the edge-of-chaos critical point may be useful as an index of consciousness in the clinical setting.

Human interictal epileptiform discharges are bidirectional traveling waves echoing ictal discharges

Interictal epileptiform discharges (IEDs), also known as interictal spikes, are large intermittent electrophysiological events observed between seizures in patients with epilepsy. Although they occur far more often than seizures, IEDs are less studied, and their relationship to seizures remains unclear. To better understand this relationship, we examined multi-day recordings of microelectrode arrays implanted in human epilepsy patients, allowing us to precisely observe the spatiotemporal propagation of IEDs, spontaneous seizures, and how they relate. These recordings showed that the majority of IEDs are traveling waves, traversing the same path as ictal discharges during seizures, and with a fixed direction relative to seizure propagation. Moreover, the majority of IEDs, like ictal discharges, were bidirectional, with one predominant and a second, less frequent antipodal direction. These results reveal a fundamental spatiotemporal similarity between IEDs and ictal discharges. These results also imply that most IEDs arise in brain tissue outside the site of seizure onset and propagate toward it, indicating that the propagation of IEDs provides useful information for localizing the seizure focus.

Ictal wavefront propagation in slices and simulations with conductance-based refractory density model

The mechanisms determining ictal discharge (ID) propagation are still not clear. In the present study, we aimed to examine these mechanisms in animal and mathematical models of epileptiform activity. Using double-patch and extracellular potassium ion concentration recordings in rat hippocampal-cortical slices, we observed that IDs moved at a speed of about 1 mm/s or less. The mechanisms of such slow propagation have been studied with a mathematical, conductance-based refractory density (CBRD) model that describes the GABA- and glutamatergic neuronal populations’ interactions and ion dynamics in brain tissue. The modeling study reveals two main factors triggerring IDs: (i) increased interneuronal activity leading to chloride ion accumulation and a consequent depolarizing GABAergic effect and (ii) the elevation of extracellular potassium ion concentration. The local synaptic transmission followed by local potassium ion extrusion and GABA receptor-mediated chloride ion accumulation underlies the ID wavefront’s propagation. In contrast, potassium ion diffusion in the extracellular space is slower and does not affect ID’s speed. The short discharges, constituting the ID, propagate much faster than the ID front. The accumulation of sodium ions inside neurons due to their hyperactivity and glutamatergic currents boosts the Na⁺/K⁺ pump, which terminates the ID. Knowledge of the mechanism of ID generation and propagation contributes to the development of new treatments against epilepsy.

During an epileptic seizure, neuronal excitation spreads across the brain tissue and is accompanied by significant changes in ionic concentrations. Ictal discharge front spreads at low speeds, less than 1 mm/s. Mechanisms underlying this phenomenon are not yet well understood. We study these mechanisms using electrophysiological recordings in brain slices and computer simulations. Our detailed biophysical model describing neuronal populations’ interaction, spatial propagation, and ionic dynamics reproduces the generation and propagation of spontaneously repeating ictal discharges. The simulations are consistent with our recordings of the electrical activity and the extracellular potassium ion concentration. We distinguished between the two alternative mechanisms of the ictal wavefront propagation: (i) the diffusion of potassium ions released from excited neurons, which depolarizes distant neurons and thus supports excitation, and (ii) the axonal spread of excitation followed by the local extracellular potassium ion accumulation that supports the excitation. Our simulations provide evidence in favor of the latter mechanism. Our experiment-based modeling contributes to a mathematical description of brain tissue functioning and potentially contributes to developing new treatments against epilepsy.

LeGUI: A Fast and Accurate Graphical User Interface for Automated Detection and Anatomical Localization of Intracranial Electrodes

Accurate anatomical localization of intracranial electrodes is important for identifying the seizure foci in patients with epilepsy and for interpreting effects from cognitive studies employing intracranial electroencephalography. Localization is typically performed by coregistering postimplant computed tomography (CT) with preoperative magnetic resonance imaging (MRI). Electrodes are then detected in the CT, and the corresponding brain region is identified using the MRI. Many existing software packages for electrode localization chain together separate preexisting programs or rely on command line instructions to perform the various localization steps, making them difficult to install and operate for a typical user. Further, many packages provide solutions for some, but not all, of the steps needed for confident localization. We have developed software, Locate electrodes Graphical User Interface (LeGUI), that consists of a single interface to perform all steps needed to localize both surface and depth/penetrating intracranial electrodes, including coregistration of the CT to MRI, normalization of the MRI to the Montreal Neurological Institute template, automated electrode detection for multiple types of electrodes, electrode spacing correction and projection to the brain surface, electrode labeling, and anatomical targeting. The software is written in MATLAB, core image processing is performed using the Statistical Parametric Mapping toolbox, and standalone executable binaries are available for Windows, Mac, and Linux platforms. LeGUI was tested and validated on 51 datasets from two universities. The total user and computational time required to process a single dataset was approximately 1 h. Automatic electrode detection correctly identified 4362 of 4695 surface and depth electrodes with only 71 false positives. Anatomical targeting was verified by comparing electrode locations from LeGUI to locations that were assigned by an experienced neuroanatomist. LeGUI showed a 94% match with the 482 neuroanatomist-assigned locations. LeGUI combines all the features needed for fast and accurate anatomical localization of intracranial electrodes into a single interface, making it a valuable tool for intracranial electrophysiology research.

Computational modeling of seizure spread on a cortical surface

In the field of computational epilepsy, neural field models helped to understand some large-scale features of seizure dynamics. These insights however remain on general levels, without translation to the clinical settings via personalization of the model with the patient-specific structure. In particular, a link was suggested between epileptic seizures spreading across the cortical surface and the so-called theta-alpha activity (TAA) pattern seen on intracranial electrographic signals, yet this link was not demonstrated on a patient-specific level. Here we present a single patient computational study linking the seizure spreading across the patient-specific cortical surface with a specific instance of the TAA pattern recorded in the patient. Using the realistic geometry of the cortical surface we perform the simulations of seizure dynamics in The Virtual Brain platform, and we show that the simulated electrographic signals qualitatively agree with the recorded signals. Furthermore, the comparison with the simulations performed on surrogate surfaces reveals that the best quantitative fit is obtained for the real surface. The work illustrates how the patient-specific cortical geometry can be utilized in The Virtual Brain for personalized model building, and the importance of such approach.

The impact of a closed-loop thalamocortical model on the spatiotemporal dynamics of cortical and thalamic traveling waves

Propagation of activity in spatially structured neuronal networks has been observed in awake, anesthetized, and sleeping brains. How these wave patterns emerge and organize across brain structures, and how network connectivity affects spatiotemporal neural activity remains unclear. Here, we develop a computational model of a two-dimensional thalamocortical network, which gives rise to emergent traveling waves similar to those observed experimentally. We illustrate how spontaneous and evoked oscillatory activity in space and time emerge using a closed-loop thalamocortical architecture, sustaining smooth waves in the cortex and staggered waves in the thalamus. We further show that intracortical and thalamocortical network connectivity, cortical excitation/inhibition balance, and thalamocortical or corticothalamic delay can independently or jointly change the spatiotemporal patterns (radial, planar and rotating waves) and characteristics (speed, direction, and frequency) of cortical and thalamic traveling waves. Computer simulations predict that increased thalamic inhibition induces slower cortical frequencies and that enhanced cortical excitation increases traveling wave speed and frequency. Overall, our results provide insight into the genesis and sustainability of thalamocortical spatiotemporal patterns, showing how simple synaptic alterations cause varied spontaneous and evoked wave patterns. Our model and simulations highlight the need for spatially spread neural recordings to uncover critical circuit mechanisms for brain functions.

The Effect of Anti-seizure Medications on the Propagation of Epileptic Activity: A Review

The propagation of epileptiform events is a highly interesting phenomenon from the pathophysiological point of view, as it involves several mechanisms of recruitment of neural networks. Extensive in vivo and in vitro research has been performed, suggesting that multiple networks as well as cellular candidate mechanisms govern this process, including the co-existence of wave propagation, coupled oscillator dynamics, and more. The clinical importance of seizure propagation stems mainly from the fact that the epileptic manifestations cannot be attributed solely to the activity in the seizure focus itself, but rather to the propagation of epileptic activity to other brain structures. Propagation, especially when causing secondary generalizations, poses a risk to patients due to recurrent falls, traumatic injuries, and poor neurological outcome. Anti-seizure medications (ASMs) affect propagation in diverse ways and with different potencies. Importantly, for drug-resistant patients, targeting seizure propagation may improve the quality of life even without a major reduction in simple focal events. Motivated by the extensive impact of this phenomenon, we sought to review the literature regarding the propagation of epileptic activity and specifically the effect of commonly used ASMs on it. Based on this body of knowledge, we propose a novel classification of ASMs into three main categories: major, minor, and intermediate efficacy in reducing the propagation of epileptiform activity.

Mean-Field Models for EEG/MEG: From Oscillations to Waves

Neural mass models have been used since the 1970s to model the coarse-grained activity of large populations of neurons. They have proven especially fruitful for understanding brain rhythms. However, although motivated by neurobiological considerations they are phenomenological in nature, and cannot hope to recreate some of the rich repertoire of responses seen in real neuronal tissue. Here we consider a simple spiking neuron network model that has recently been shown to admit an exact mean-field description for both synaptic and gap-junction interactions. The mean-field model takes a similar form to a standard neural mass model, with an additional dynamical equation to describe the evolution of within-population synchrony. As well as reviewing the origins of this next generation mass model we discuss its extension to describe an idealised spatially extended planar cortex. To emphasise the usefulness of this model for EEG/MEG modelling we show how it can be used to uncover the role of local gap-junction coupling in shaping large scale synaptic waves.

Bidirectional propagation of low frequency oscillations over the human hippocampal surface

The hippocampus is diversely interconnected with other brain systems along its axis. Cycles of theta-frequency activity are believed to propagate from the septal to temporal pole, yet it is unclear how this one-way route supports the flexible cognitive capacities of this structure. We leveraged novel thin-film microgrid arrays conformed to the human hippocampal surface to track neural activity two-dimensionally in vivo. All oscillation frequencies identified between 1-15 Hz propagated across the tissue. Moreover, they dynamically shifted between two roughly opposite directions oblique to the long axis. This predominant propagation axis was mirrored across participants, hemispheres, and consciousness states. Directionality was modulated in a participant who performed a behavioral task, and it could be predicted by wave amplitude topography over the hippocampal surface. Our results show that propagation directions may thus represent distinct meso-scale network computations, operating along versatile spatiotemporal processing routes across the hippocampal body.

Travelling waves reveal a dynamic seizure source in human focal epilepsy

Treatment of patients with drug-resistant focal epilepsy relies upon accurate seizure localization. Ictal activity captured by intracranial EEG has traditionally been interpreted to suggest that the underlying cortex is actively involved in seizures. Here, we hypothesize that such activity instead reflects propagated activity from a relatively focal seizure source, even during later time points when ictal activity is more widespread. We used the time differences observed between ictal discharges in adjacent electrodes to estimate the location of the hypothesized focal source and demonstrated that the seizure source, localized in this manner, closely matches the clinically and neurophysiologically determined brain region giving rise to seizures. Moreover, we determined this focal source to be a dynamic entity that moves and evolves over the time course of a seizure. Our results offer an interpretation of ictal activity observed by intracranial EEG that challenges the traditional conceptualization of the seizure source.

Evidence for spreading seizure as a cause of theta-alpha activity electrographic pattern in stereo-EEG seizure recordings

Intracranial electroencephalography is a standard tool in clinical evaluation of patients with focal epilepsy. Various early electrographic seizure patterns differing in frequency, amplitude, and waveform of the oscillations are observed. The pattern most common in the areas of seizure propagation is the so-called theta-alpha activity (TAA), whose defining features are oscillations in the θ - α range and gradually increasing amplitude. A deeper understanding of the mechanism underlying the generation of the TAA pattern is however lacking. In this work we evaluate the hypothesis that the TAA patterns are caused by seizures spreading across the cortex. To do so, we perform simulations of seizure dynamics on detailed patient-derived cortical surfaces using the spreading seizure model as well as reference models with one or two homogeneous sources. We then detect the occurrences of the TAA patterns both in the simulated stereo-electroencephalographic signals and in the signals of recorded epileptic seizures from a cohort of fifty patients, and we compare the features of the groups of detected TAA patterns to assess the plausibility of the different models. Our results show that spreading seizure hypothesis is qualitatively consistent with the evidence available in the seizure recordings, and it can explain the features of the detected TAA groups best among the examined models.

Data-driven method to infer the seizure propagation patterns in an epileptic brain from intracranial electroencephalography

Surgical interventions in epileptic patients aimed at the removal of the epileptogenic zone have success rates at only 60-70%. This failure can be partly attributed to the insufficient spatial sampling by the implanted intracranial electrodes during the clinical evaluation, leading to an incomplete picture of spatio-temporal seizure organization in the regions that are not directly observed. Utilizing the partial observations of the seizure spreading through the brain network, complemented by the assumption that the epileptic seizures spread along the structural connections, we infer if and when are the unobserved regions recruited in the seizure. To this end we introduce a data-driven model of seizure recruitment and propagation across a weighted network, which we invert using the Bayesian inference framework. Using a leave-one-out cross-validation scheme on a cohort of 45 patients we demonstrate that the method can improve the predictions of the states of the unobserved regions compared to an empirical estimate that does not use the structural information, yet it is on the same level as the estimate that takes the structure into account. Furthermore, a comparison with the performed surgical resection and the surgery outcome indicates a link between the inferred excitable regions and the actual epileptogenic zone. The results emphasize the importance of the structural connectome in the large-scale spatio-temporal organization of epileptic seizures and introduce a novel way to integrate the patient-specific connectome and intracranial seizure recordings in a whole-brain computational model of seizure spread.

Measuring the Coupling Direction between Neural Oscillations with Weighted Symbolic Transfer Entropy

Neural oscillations reflect rhythmic fluctuations in the synchronization of neuronal populations and play a significant role in neural processing. To further understand the dynamic interactions between different regions in the brain, it is necessary to estimate the coupling direction between neural oscillations. Here, we developed a novel method, termed weighted symbolic transfer entropy (WSTE), that combines symbolic transfer entropy (STE) and weighted probability distribution to measure the directionality between two neuronal populations. The traditional STE ignores the degree of difference between the amplitude values of a time series. In our proposed WSTE method, this information is picked up by utilizing a weighted probability distribution. The simulation analysis shows that the WSTE method can effectively estimate the coupling direction between two neural oscillations. In comparison with STE, the new method is more sensitive to the coupling strength and is more robust against noise. When applied to epileptic electrocorticography data, a significant coupling direction from the anterior nucleus of thalamus (ANT) to the seizure onset zone (SOZ) was detected during seizures. Considering the superiorities of the WSTE method, it is greatly advantageous to measure the coupling direction between neural oscillations and consequently characterize the information flow between different brain regions.

Acoustoelectric imaging of deep dipoles in a human head phantom for guiding treatment of epilepsy

OBJECTIVE

This study employs a human head model with real skull to demonstrate the feasibility of transcranial acoustoelectric brain imaging (tABI) as a new modality for electrical mapping of deep dipole sources during treatment of epilepsy with much better resolution and accuracy than conventional mapping methods.

APPROACH

This technique exploits an interaction between a focused ultrasound (US) beam and tissue resistivity to localize current source densities as deep as 63 mm at high spatial resolution (1 to 4 mm) and resolve fast time-varying currents with sub-ms precision.

MAIN RESULTS

Detection thresholds through a thick segment of the human skull at biologically safe US intensities was below 0.5 mA and within range of strong currents generated by the human brain.

SIGNIFICANCE

This work suggests that 4D tABI may emerge as a revolutionary modality for real-time high-resolution mapping of neuronal currents for the purpose of monitoring, staging, and guiding treatment of epilepsy and other brain disorders characterized by abnormal rhythms.

Seizure Onset Zone Lateralization Using a Non-linear Analysis of Micro vs. Macro Electroencephalographic Recordings During Seizure-Free Stages of the Sleep-Wake Cycle From Epilepsy Patients

The application of non-linear signal analysis techniques to biomedical data is key to improve our knowledge about complex physiological and pathological processes. In particular, the use of non-linear techniques to study electroencephalographic (EEG) recordings can provide an advanced characterization of brain dynamics. In epilepsy these dynamics are altered at different spatial scales of neuronal organization. We therefore apply non-linear signal analysis to EEG recordings from epilepsy patients derived with intracranial hybrid electrodes, which are composed of classical macro contacts and micro wires. Thereby, these electrodes record EEG at two different spatial scales. Our aim is to test the degree to which the analysis of the EEG recorded at these different scales allows us to characterize the neuronal dynamics affected by epilepsy. For this purpose, we retrospectively analyzed long-term recordings performed during five nights in three patients during which no seizures took place. As a benchmark we used the accuracy with which this analysis allows determining the hemisphere that contains the seizure onset zone, which is the brain area where clinical seizures originate. We applied the surrogate-corrected non-linear predictability score (ψ), a non-linear signal analysis technique which was shown previously to be useful for the lateralization of the seizure onset zone from classical intracranial EEG macro contact recordings. Higher values of ψ were found predominantly for signals recorded from the hemisphere containing the seizure onset zone as compared to signals recorded from the opposite hemisphere. These differences were found not only for the EEG signals recorded with macro contacts, but also for those recorded with micro wires. In conclusion, the information obtained from the analysis of classical macro EEG contacts can be complemented by the one of micro wire EEG recordings. This combined approach may therefore help to further improve the degree to which quantitative EEG analysis can contribute to the diagnostics in epilepsy patients.

Intracranial EEG in the 21st Century

Intracranial electroencephalography (iEEG) has been the mainstay of identifying the seizure onset zone (SOZ), a key diagnostic procedure in addition to neuroimaging when considering epilepsy surgery. In many patients, iEEG has been the basis for resective epilepsy surgery, to date still the most successful treatment for drug-resistant epilepsy. Intracranial EEG determines the location and resectability of the SOZ. Advances in recording and implantation of iEEG provide multiple options in the 21st century. This not only includes the choice between subdural electrodes (SDE) and stereoelectroencephalography (SEEG) but also includes the implantation and recordings from microelectrodes. Before iEEG implantation, especially in magnetic resonance imaging -negative epilepsy, a clear hypothesis for seizure generation and propagation should be based on noninvasive methods. Intracranial EEG implantation should be planned by a multidisciplinary team considering epileptic networks. Recordings from SDE and SEEG have both their advantages and disadvantages. Stereo-EEG seems to have a lower rate of complications that are clinically significant, but has limitations in spatial sampling of the cortical surface. Stereo-EEG can sample deeper areas of the brain including deep sulci and hard to reach areas such as the insula. To determine the epileptogenic zone, interictal and ictal information should be taken into consideration. Interictal spiking, low frequency slowing, as well as high frequency oscillations may inform about the epileptogenic zone. Ictally, high frequency onsets in the beta/gamma range are usually associated with the SOZ, but specialized recordings with combined macro and microelectrodes may in the future educate us about onset in higher frequency bands. Stimulation of intracranial electrodes triggering habitual seizures can assist in identifying the SOZ. Advanced computational methods such as determining the epileptogenicity index and similar measures may enhance standard clinical interpretation. Improved techniques to record and interpret iEEG may in the future lead to a greater proportion of patients being seizure free after epilepsy surgery.

Seizure pathways change on circadian and slower timescales in individual patients with focal epilepsy

Personalized medicine requires that treatments adapt to not only the patient but also changing factors within each individual. Although epilepsy is a dynamic disorder characterized by pathological fluctuations in brain state, surprisingly little is known about whether and how seizures vary in the same patient. We quantitatively compared within-patient seizure network evolutions using intracranial electroencephalographic (iEEG) recordings of over 500 seizures from 31 patients with focal epilepsy (mean 16.5 seizures per patient). In all patients, we found variability in seizure paths through the space of possible network dynamics. Seizures with similar pathways tended to occur closer together in time, and a simple model suggested that seizure pathways change on circadian and/or slower timescales in the majority of patients. These temporal relationships occurred independent of whether the patient underwent antiepileptic medication reduction. Our results suggest that various modulatory processes, operating at different timescales, shape within-patient seizure evolutions, leading to variable seizure pathways that may require tailored treatment approaches.

A simple model of epileptic seizure propagation: Potassium diffusion versus axo-dendritic spread

The mechanisms of epileptic discharge generation and spread are not yet fully known. A recently proposed simple biophysical model of interictal and ictal discharges, Epileptor-2, reproduces well the main features of neuronal excitation and ionic dynamics during discharge generation. In order to distinguish between two hypothesized mechanisms of discharge propagation, we extend the model to the case of two-dimensional propagation along the cortical neural tissue. The first mechanism is based on extracellular potassium diffusion, and the second is the propagation of spikes and postsynaptic signals along axons and dendrites. Our simulations show that potassium diffusion is too slow to reproduce an experimentally observed speed of ictal wavefront propagation (tenths of mm/s). By contrast, the synaptic mechanism predicts well the speed and synchronization of the pre-ictal bursts before the ictal front and the afterdischarges in the ictal core. Though this fact diminishes the role of diffusion and electrodiffusion, the model nevertheless highlights the role of potassium extrusion during neuronal excitation, which provides a positive feedback that changes at the ictal wavefront the balance of excitation versus inhibition in favor of excitation. This finding may help to find a target for a treatment to prevent seizure propagation.

Model-based design for seizure control by stimulation

OBJECTIVE

Current brain stimulation paradigms are largely empirical rather than theoretical. An opportunity exists to improve upon their modest effectiveness in closed-loop control strategies with the development of theoretically grounded, model-based designs.

APPROACH

Inspired by this need, here we couple experimental data and mathematical modeling with a control-theoretic strategy for seizure termination. We begin by exercising a dynamical systems approach to model seizures (n = 94) recorded using intracranial EEG (iEEG) from 21 patients with medication-resistant, localization-related epilepsy.

MAIN RESULTS

Although each patient's seizures displayed unique spatial and temporal patterns, their evolution can be parsimoniously characterized by the same model form. Idiosyncracies of the model can inform individualized intervention strategies, specifically in iEEG samples with well-localized seizure onset zones. Temporal fluctuations in the spatial profiles of the oscillatory modes show that seizure onset marks a transition into a regime in which the underlying system supports prolonged rhythmic and focal activity. Based on these observations, we propose a control-theoretic strategy that aims to stabilize ictal activity using static output feedback for linear time-invariant switching systems. Finally, we demonstrate in silico that our proposed strategy allows us to dampen the emerging focal oscillatory sources using only a small set of electrodes.

SIGNIFICANCE

Our integrative study informs the development of modulation and control algorithms for neurostimulation that could improve the effectiveness of implantable, closed-loop anti-epileptic devices.

A model for focal seizure onset, propagation, evolution, and progression

We developed a neural network model that can account for major elements common to human focal seizures. These include the tonic-clonic transition, slow advance of clinical semiology and corresponding seizure territory expansion, widespread EEG synchronization, and slowing of the ictal rhythm as the seizure approaches termination. These were reproduced by incorporating usage-dependent exhaustion of inhibition in an adaptive neural network that receives global feedback inhibition in addition to local recurrent projections. Our model proposes mechanisms that may underline common EEG seizure onset patterns and status epilepticus, and postulates a role for synaptic plasticity in the emergence of epileptic foci. Complex patterns of seizure activity and bi-stable seizure end-points arise when stochastic noise is included. With the rapid advancement of clinical and experimental tools, we believe that this model can provide a roadmap and potentially an in silico testbed for future explorations of seizure mechanisms and clinical therapies.

Controversies on the network theory of epilepsy: Debates held during the ICTALS 2019 conference

Debates on six controversial topics on the network theory of epilepsy were held during two debate sessions, as part of the International Conference for Technology and Analysis of Seizures, 2019 (ICTALS 2019) convened at the University of Exeter, UK, September 2-5 2019. The debate topics were (1) From pathologic to physiologic: is the epileptic network part of an existing large-scale brain network? (2) Are micro scale recordings pertinent for defining the epileptic network? (3) From seconds to years: do we need all temporal scales to define an epileptic network? (4) Is it necessary to fully define the epileptic network to control it? (5) Is controlling seizures sufficient to control the epileptic network? (6) Does the epileptic network want to be controlled? This article, written by the organizing committee for the debate sessions and the debaters, summarizes the arguments presented during the debates on these six topics.

Transient neocortical gamma oscillations induced by neuronal response modulation

In this paper a mean field model of spatio-temporal electroencephalographic activity in the neocortex is used to computationally study the emergence of neocortical gamma oscillations as a result of neuronal response modulation. It is shown using a numerical bifurcation analysis that gamma oscillations emerge robustly in the solutions of the model and transition to beta oscillations through coordinated modulation of the responsiveness of inhibitory and excitatory neuronal populations. The spatio-temporal pattern of the propagation of these oscillations across the neocortex is illustrated by solving the equations of the model using a finite element software package. Thereby, it is shown that the gamma oscillations remain localized to the regions of neuronal modulation. Moreover, it is discussed that the inherent spatial averaging effect of commonly used electrocortical measurement techniques can significantly alter the amplitude and pattern of fast oscillations in neocortical recordings, and hence can potentially affect physiological interpretations of these recordings.

Investigating the role of gap junctions in seizure wave propagation

The effect of gap junctions as well as the biological mechanisms behind seizure wave propagation is not completely understood. In this work, we use a simple neural field model to study the possible influence of gap junctions specifically on cortical wave propagation that has been observed in vivo preceding seizure termination. We consider a voltage-based neural field model consisting of an excitatory and an inhibitory population as well as both chemical and gap junction-like synapses. We are able to approximate important properties of cortical wave propagation previously observed in vivo before seizure termination. This model adds support to existing evidence from models and clinical data suggesting a key role of gap junctions in seizure wave propagation. In particular, we found that in this model gap junction-like connectivity determines the propagation of one-bump or two-bump traveling wave solutions with features consistent with the clinical data. For sufficiently increased gap junction connectivity, wave solutions cease to exist. Moreover, gap junction connectivity needs to be sufficiently low or moderate to permit the existence of linearly stable solutions of interest.

Next-generation neural mass and field modeling

The Wilson-Cowan population model of neural activity has greatly influenced our understanding of the mechanisms for the generation of brain rhythms and the emergence of structured brain activity. As well as the many insights that have been obtained from its mathematical analysis, it is now widely used in the computational neuroscience community for building large-scale in silico brain networks that can incorporate the increasing amount of knowledge from the Human Connectome Project. Here, we consider a neural population model in the spirit of that originally developed by Wilson and Cowan, albeit with the added advantage that it can account for the phenomena of event-related synchronization and desynchronization. This derived mean-field model provides a dynamic description for the evolution of synchrony, as measured by the Kuramoto order parameter, in a large population of quadratic integrate-and-fire model neurons. As in the original Wilson-Cowan framework, the population firing rate is at the heart of our new model; however, in a significant departure from the sigmoidal firing rate function approach, the population firing rate is now obtained as a real-valued function of the complex-valued population synchrony measure. To highlight the usefulness of this next-generation Wilson-Cowan style model, we deploy it in a number of neurobiological contexts, providing understanding of the changes in power spectra observed in electro- and magnetoencephalography neuroimaging studies of motor cortex during movement, insights into patterns of functional connectivity observed during rest and their disruption by transcranial magnetic stimulation, and to describe wave propagation across cortex.

Neural field models for latent state inference: Application to large-scale neuronal recordings

Large-scale neural recording methods now allow us to observe large populations of identified single neurons simultaneously, opening a window into neural population dynamics in living organisms. However, distilling such large-scale recordings to build theories of emergent collective dynamics remains a fundamental statistical challenge. The neural field models of Wilson, Cowan, and colleagues remain the mainstay of mathematical population modeling owing to their interpretable, mechanistic parameters and amenability to mathematical analysis. Inspired by recent advances in biochemical modeling, we develop a method based on moment closure to interpret neural field models as latent state-space point-process models, making them amenable to statistical inference. With this approach we can infer the intrinsic states of neurons, such as active and refractory, solely from spiking activity in large populations. After validating this approach with synthetic data, we apply it to high-density recordings of spiking activity in the developing mouse retina. This confirms the essential role of a long lasting refractory state in shaping spatiotemporal properties of neonatal retinal waves. This conceptual and methodological advance opens up new theoretical connections between mathematical theory and point-process state-space models in neural data analysis.

Developing statistical tools to connect single-neuron activity to emergent collective dynamics is vital for building interpretable models of neural activity. Neural field models relate single-neuron activity to emergent collective dynamics in neural populations, but integrating them with data remains challenging. Recently, latent state-space models have emerged as a powerful tool for constructing phenomenological models of neural population activity. The advent of high-density multi-electrode array recordings now enables us to examine large-scale collective neural activity. We show that classical neural field approaches can yield latent state-space equations and demonstrate that this enables inference of the intrinsic states of neurons from recorded spike trains in large populations.

The generation and propagation of the human alpha rhythm

The alpha rhythm is the longest-studied brain oscillation and has been theorized to play a key role in cognition. Still, its physiology is poorly understood. In this study, we used microelectrodes and macroelectrodes in surgical epilepsy patients to measure the intracortical and thalamic generators of the alpha rhythm during quiet wakefulness. We first found that alpha in both visual and somatosensory cortex propagates from higher-order to lower-order areas. In posterior cortex, alpha propagates from higher-order anterosuperior areas toward the occipital pole, whereas alpha in somatosensory cortex propagates from associative regions toward primary cortex. Several analyses suggest that this cortical alpha leads pulvinar alpha, complicating prevailing theories of a thalamic pacemaker. Finally, alpha is dominated by currents and firing in supragranular cortical layers. Together, these results suggest that the alpha rhythm likely reflects short-range supragranular feedback, which propagates from higher- to lower-order cortex and cortex to thalamus. These physiological insights suggest how alpha could mediate feedback throughout the thalamocortical system.

A two-layered brain network model and its chimera state

Based on the data of cerebral cortex, we present a two-layered brain network model of coupled neurons where the two layers represent the left and right hemispheres of cerebral cortex, respectively, and the links between the two layers represent the inter-couplings through the corpus callosum. By this model we show that abundant patterns of synchronization can be observed, especially the chimera state, depending on the parameters of system such as the coupling strengths and coupling phase. Further, we extend the model to a more general two-layered network to better understand the mechanism of the observed patterns, where each hemisphere of cerebral cortex is replaced by a highly clustered subnetwork. We find that the number of inter-couplings is another key parameter for the emergence of chimera states. Thus, the chimera states come from a matching between the structure parameters such as the number of inter-couplings and clustering coefficient etc and the dynamics parameters such as the intra-, inter-coupling strengths and coupling phase etc. A brief theoretical analysis is provided to explain the borderline of synchronization. These findings may provide helpful clues to understand the mechanism of brain functions.