Multivariate regression methods for estimating velocity of ictal discharges from human microelectrode recordings

Seizure event detection using intravital two-photon calcium imaging data

Significance

Intravital cellular calcium imaging has emerged as a powerful tool to investigate how different types of neurons interact at the microcircuit level to produce seizure activity, with newfound potential to understand epilepsy. Although many methods exist to measure seizure-related activity in traditional electrophysiology, few yet exist for calcium imaging.

Aim

To demonstrate an automated algorithmic framework to detect seizure-related events using calcium imaging-including the detection of pre-ictal spike events, propagation of the seizure wavefront, and terminal spreading waves for both population-level activity and that of individual cells.

Approach

We developed an algorithm for precise recruitment detection of population and individual cells during seizure-associated events, which broadly leverages averaged population activity and high-magnitude slope features to detect single-cell pre-ictal spike and seizure recruitment. We applied this method to data recorded using awake in vivo two-photon calcium imaging during pentylenetetrazol-induced seizures in mice.

Results

We demonstrate that our detected recruitment times are concordant with visually identified labels provided by an expert reviewer and are sufficiently accurate to model the spatiotemporal progression of seizure-associated traveling waves.

Conclusions

Our algorithm enables accurate cell recruitment detection and will serve as a useful tool for researchers investigating seizure dynamics using calcium imaging.

Identifying sources of human interictal discharges with travelling wave and white matter propagation

Interictal epileptiform discharges have been shown to propagate from focal epileptogenic sources as travelling waves or through more rapid white matter conduction. We hypothesize that both modes of propagation are necessary to explain interictal discharge timing delays. We propose a method that, for the first time, incorporates both propagation modes to identify unique potential sources of interictal activity. We retrospectively analysed 38 focal epilepsy patients who underwent intracranial EEG recordings and diffusion-weighted imaging for epilepsy surgery evaluation. Interictal discharges were detected and localized to the most likely source based on relative delays in time of arrival across electrodes, incorporating travelling waves and white matter propagation. We assessed the influence of white matter propagation on distance of spread, timing and clinical interpretation of interictal activity. To evaluate accuracy, we compared our source localization results to earliest spiking regions to predict seizure outcomes. White matter propagation helps to explain the timing delays observed in interictal discharge sequences, underlying rapid and distant propagation. Sources identified based on differences in time of receipt of interictal discharges are often distinct from the leading electrode location. Receipt of activity propagating rapidly via white matter can occur earlier than more local activity propagating via slower cortical travelling waves. In our cohort, our source localization approach was more accurate in predicting seizure outcomes than the leading electrode location. Inclusion of white matter in addition to travelling wave propagation in our model of discharge spread did not improve overall accuracy but allowed for identification of unique and at times distant potential sources of activity, particularly in patients with persistent postoperative seizures. Since distant white matter propagation can occur more rapidly than local travelling wave propagation, combined modes of propagation within an interictal discharge sequence can decouple the commonly assumed relationship between spike timing and distance from the source. Our findings thus highlight the clinical importance of recognizing the presence of dual modes of propagation during interictal discharges, as this may be a cause of clinical mislocalization.

Seizure Event Detection Using Intravital Two-Photon Calcium Imaging Data

Significance

Genetic cellular calcium imaging has emerged as a powerful tool to investigate how different types of neurons interact at the microcircuit level to produce seizure activity, with newfound potential to understand epilepsy. Although many methods exist to measure seizure-related activity in traditional electrophysiology, few yet exist for calcium imaging.

Aim

To demonstrate an automated algorithmic framework to detect seizure-related events using calcium imaging - including the detection of pre-ictal spike events, propagation of the seizure wavefront, and terminal spreading waves for both population-level activity and that of individual cells.

Approach

We developed an algorithm for precise recruitment detection of population and individual cells during seizure-associated events, which broadly leverages averaged population activity and high-magnitude slope features to detect single-cell pre-ictal spike and seizure recruitment. We applied this method to data recorded using awake in vivo two-photon calcium imaging during pentylenetetrazol induced seizures in mice.

Results

We demonstrate that our detected recruitment times are concordant with visually identified labels provided by an expert reviewer and are sufficiently accurate to model the spatiotemporal progression of seizure-associated traveling waves.

Conclusions

Our algorithm enables accurate cell recruitment detection and will serve as a useful tool for researchers investigating seizure dynamics using calcium imaging.

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.

Human acute microelectrode array recordings with broad cortical access, single-unit resolution, and parallel behavioral monitoring

There are vast gaps in our understanding of the organization and operation of the human nervous system at the level of individual neurons and their networks. Here, we report reliable and robust acute multichannel recordings using planar microelectrode arrays (MEAs) implanted intracortically in awake brain surgery with open craniotomies that grant access to large parts of the cortical hemisphere. We obtained high-quality extracellular neuronal activity at the microcircuit, local field potential level and at the cellular, single-unit level. Recording from the parietal association cortex, a region rarely explored in human single-unit studies, we demonstrate applications on these complementary spatial scales and describe traveling waves of oscillatory activity as well as single-neuron and neuronal population responses during numerical cognition, including operations with uniquely human number symbols. Intraoperative MEA recordings are practicable and can be scaled up to explore cellular and microcircuit mechanisms of a wide range of human brain functions.

Subsets of cortico-cortical evoked potentials propagate as traveling waves

Emerging evidence suggests that the temporal dynamics of cortico-cortical evoked potentials (CCEPs) may be used to characterize the patterns of information flow between and within brain networks. At present, however, the spatiotemporal dynamics of CCEP propagation cortically and subcortically are incompletely understood. We hypothesized that CCEPs propagate as an evoked traveling wave emanating from the site of stimulation. To elicit CCEPs, we applied single-pulse stimulation to stereoelectroencephalography (SEEG) electrodes implanted in 21 adult patients with intractable epilepsy. For each robust CCEP, we measured the timing of the maximal descent in evoked local field potentials and broadband high-gamma power (70-150 Hz) envelopes relative to the distance between the recording and stimulation contacts using three different metrics (i.e., Euclidean distance, path length, geodesic distance), representing direct, subcortical, and transcortical propagation, respectively. Many evoked responses to single-pulse electrical stimulation appear to propagate as traveling waves (~17-30%), even in the sparsely sampled, three-dimensional SEEG space. These results provide new insights into the spatiotemporal dynamics of CCEP propagation.

Tracking Multisite Seizure Propagation Using Ictal High-Gamma Activity

PURPOSE

Spatial patterns of long-range seizure propagation in epileptic networks have not been well characterized. Here, we use ictal high-gamma activity (HGA) as a proxy of intense neuronal population firing to map the spatial evolution of seizure recruitment.

METHODS

Ictal HGA (80-150 Hz) was analyzed in 13 patients with 72 seizures recorded by stereotactic depth electrodes, using previously validated methods. Distinct spatial clusters of channels with the ictal high-gamma signature were identified, and seizure hubs were defined as stereotypically recruited nonoverlapping clusters. Clusters correlated with asynchronous seizure terminations to provide supportive evidence for independent seizure activity at these sites. The spatial overlap between seizure hubs and interictal ripples was compared.

RESULTS

Ictal HGA was detected in 71% of seizures and 10% of implanted contacts, enabling tracking of contiguous and noncontiguous seizure recruitment. Multiple seizure hubs were identified in 54% of cases, including 43% of patients thought preoperatively to have unifocal epilepsy. Noncontiguous recruitment was associated with asynchronous seizure termination (odds ratio = 19.7; p = 0.029). Interictal ripples demonstrated greater spatial overlap with ictal HGA in cases with single seizure hubs compared with those with multiple hubs (100% vs. 66% per patient; p = 0.03).

CONCLUSIONS

Ictal HGA may serve as a useful adjunctive biomarker to distinguish contiguous seizure spread from propagation to remote seizure sites. High-gamma sites were found to cluster in stereotyped seizure hubs rather than being broadly distributed. Multiple hubs were common even in cases that were considered unifocal.

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.

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.

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.

Accuracy of omni-planar and surface casting of epileptiform activity for intracranial seizure localization

OBJECTIVE

Intracranial electroencephalography (ICEEG) recordings are performed for seizure localization in medically refractory epilepsy. Signal quantifications such as frequency power can be projected as heatmaps on personalized three-dimensional (3D) reconstructed cortical surfaces to distill these complex recordings into intuitive cinematic visualizations. However, simultaneously reconciling deep recording locations and reliably tracking evolving ictal patterns remain significant challenges.

METHODS

We fused oblique magnetic resonance imaging (MRI) slices along depth probe trajectories with cortical surface reconstructions and projected dynamic heatmaps using a simple mathematical metric of epileptiform activity (line-length). This omni-planar and surface casting of epileptiform activity approach (OPSCEA) thus illustrated seizure onset and spread among both deep and superficial locations simultaneously with minimal need for signal processing supervision. We utilized the approach on 41 patients at our center implanted with grid, strip, and/or depth electrodes for localizing medically refractory seizures. Peri-ictal data were converted into OPSCEA videos with multiple 3D brain views illustrating all electrode locations. Five people of varying expertise in epilepsy (medical student through epilepsy attending level) attempted to localize the seizure-onset zones.

RESULTS

We retrospectively compared this approach with the original ICEEG study reports for validation. Accuracy ranged from 73.2% to 97.6% for complete or overlapping onset lobe(s), respectively, and ~56.1% to 95.1% for the specific focus (or foci). Higher answer certainty for a given case predicted better accuracy, and scorers had similar accuracy across different training levels.

SIGNIFICANCE

In an era of increasing stereo-EEG use, cinematic visualizations fusing omni-planar and surface functional projections appear to provide a useful adjunct for interpreting complex intracranial recordings and subsequent surgery planning.

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.

Isoflurane-Induced Burst Suppression Is a Thalamus-Modulated, Focal-Onset Rhythm With Persistent Local Asynchrony and Variable Propagation Patterns in Rats

Background: Inhalational anesthetic-induced burst suppression (BS) is classically considered a bilaterally synchronous rhythm. However, local asynchrony has been predicted in theoretical studies and reported in patients with pre-existing focal pathology. Method: We used high-speed widefield calcium imaging to study the spatiotemporal dynamics of isoflurane-induced BS in rats. Results: We found that isoflurane-induced BS is not a globally synchronous rhythm. In the neocortex, neural activity first emerged in a spatially shifting, variably localized focus. Subsequent propagation across the whole cortex was rapid, typically within <100 milliseconds, giving the superficial resemblance to global synchrony. Neural activity remained locally asynchronous during the bursts, forming complex recurrent propagating waves. Despite propagation variability, spatial sequences of burst propagation were largely preserved between the hemispheres, and neural activity was highly correlated between the homotopic areas. The critical role of the thalamus in cortical burst initiation was demonstrated by using unilateral thalamic tetrodotoxin injection. Conclusion: The classical impression that anesthetics-induced BS is a state of global brain synchrony is inaccurate. Bursts are a series of shifting local cortical events facilitated by thalamic projection that unfold as rapid, bilaterally asynchronous propagating waves.

Burst suppression uncovers rapid widespread alterations in network excitability caused by an acute seizure focus

Burst suppression is an electroencephalogram pattern of globally symmetric alternating high amplitude activity and isoelectricity that can be induced by general anaesthetics. There is scattered evidence that burst suppression may become spatially non-uniform in the setting of underlying pathology. Here, we induced burst suppression with isoflurane in rodents and then created a neocortical acute seizure focus with injection of 4-aminopyridine (4-AP) in somatosensory cortex. Burst suppression events were recorded before and after creation of the focus using bihemispheric wide-field calcium imaging and multielectrode arrays. We find that the seizure focus elicits a rapid alteration in triggering, initiation, and propagation of burst suppression events. Compared with the non-seizing brain, bursts are triggered from the thalamus, initiate in regions uniquely outside the epileptic focus, elicit marked increases of multiunit activity and propagate towards the seizure focus. These findings support the rapid, widespread impact of focal epilepsy on the extended brain network.

How technology is driving the landscape of epilepsy surgery

This article emphasizes the role of the technological progress in changing the landscape of epilepsy surgery and provides a critical appraisal of robotic applications, laser interstitial thermal therapy, intraoperative imaging, wireless recording, new neuromodulation techniques, and high-intensity focused ultrasound. Specifically, (a) it relativizes the current hype in using robots for stereo-electroencephalography (SEEG) to increase the accuracy of depth electrode placement and save operating time; (b) discusses the drawback of laser interstitial thermal therapy (LITT) when it comes to the need for adequate histopathologic specimen and the fact that the concept of stereotactic disconnection is not new; (c) addresses the ratio between the benefits and expenditure of using intraoperative magnetic resonance imaging (MRI), that is, the high technical and personnel expertise needed that might restrict its use to centers with a high case load, including those unrelated to epilepsy; (d) soberly reviews the advantages, disadvantages, and future potentials of neuromodulation techniques with special emphasis on the differences between closed and open-loop systems; and (e) provides a critical outlook on the clinical implications of focused ultrasound, wireless recording, and multipurpose electrodes that are already on the horizon. This outlook shows that although current ultrasonic systems do have some limitations in delivering the acoustic energy, further advance of this technique may lead to novel treatment paradigms. Furthermore, it highlights that new data streams from multipurpose electrodes and wireless transmission of intracranial recordings will become available soon once some critical developments will be achieved such as electrode fidelity, data processing and storage, heat conduction as well as rechargeable technology. A better understanding of modern epilepsy surgery will help to demystify epilepsy surgery for the patients and the treating physicians and thereby reduce the surgical treatment gap.

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.

Role of paroxysmal depolarization in focal seizure activity

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

Multiscale recordings reveal the dynamic spatial structure of human seizures

The cellular activity underlying human focal seizures, and its relationship to key signatures in the EEG recordings used for therapeutic purposes, has not been well characterized despite many years of investigation both in laboratory and clinical settings. The increasing use of microelectrodes in epilepsy surgery patients has made it possible to apply principles derived from laboratory research to the problem of mapping the spatiotemporal structure of human focal seizures, and characterizing the corresponding EEG signatures. In this review, we describe results from human microelectrode studies, discuss some data interpretation pitfalls, and explain the current understanding of the key mechanisms of ictogenesis and seizure spread.

Theta and Alpha Oscillations Are Traveling Waves in the Human Neocortex

Human cognition requires the coordination of neural activity across widespread brain networks. Here, we describe a new mechanism for large-scale coordination in the human brain: traveling waves of theta and alpha oscillations. Examining direct brain recordings from neurosurgical patients performing a memory task, we found contiguous clusters of cortex in individual patients with oscillations at specific frequencies within 2 to 15 Hz. These oscillatory clusters displayed spatial phase gradients, indicating that they formed traveling waves that propagated at ∼0.25-0.75 m/s. Traveling waves were relevant behaviorally because their propagation correlated with task events and was more consistent when subjects performed the task well. Human traveling theta and alpha waves can be modeled by a network of coupled oscillators because the direction of wave propagation correlated with the spatial orientation of local frequency gradients. Our findings suggest that oscillations support brain connectivity by organizing neural processes across space and time.