Cortical activation mapping of epileptiform activity derived from interictal ECoG spikes

Spike propagation mapping reveals effective connectivity and predicts surgical outcome in epilepsy

Neurosurgical intervention is the best available treatment for selected patients with drug resistant epilepsy. For these patients, surgical planning requires biomarkers that delineate the epileptogenic zone, the brain area that is indispensable for the generation of seizures. Interictal spikes recorded with electrophysiological techniques are considered key biomarkers of epilepsy. Yet, they lack specificity, mostly because they propagate across brain areas forming networks. Understanding the relationship between interictal spike propagation and functional connections among the involved brain areas may help develop novel biomarkers that can delineate the epileptogenic zone with high precision. Here, we reveal the relationship between spike propagation and effective connectivity among onset and areas of spread and assess the prognostic value of resecting these areas. We analysed intracranial EEG data from 43 children with drug resistant epilepsy who underwent invasive monitoring for neurosurgical planning. Using electric source imaging, we mapped spike propagation in the source domain and identified three zones: onset, early-spread and late-spread. For each zone, we calculated the overlap and distance from surgical resection. We then estimated a virtual sensor for each zone and the direction of information flow among them via Granger causality. Finally, we compared the prognostic value of resecting these zones, the clinically-defined seizure onset zone and the spike onset on intracranial EEG channels by estimating their overlap with resection. We observed a spike propagation in source space for 37 patients with a median duration of 95 ms (interquartile range: 34-206), a spatial displacement of 14 cm (7.5-22 cm) and a velocity of 0.5 m/s (0.3-0.8 m/s). In patients with good surgical outcome (25 patients, Engel I), the onset had higher overlap with resection [96% (40-100%)] than early-spread [86% (34-100%), P = 0.01] and late-spread [59% (12-100%), P = 0.002], and it was also closer to resection than late-spread [5 mm versus 9 mm, P = 0.007]. We found an information flow from onset to early-spread in 66% of patients with good outcomes, and from early-spread to onset in 50% of patients with poor outcome. Finally, resection of spike onset, but not area of spike spread or the seizure onset zone, predicted outcome with positive predictive value of 79% and negative predictive value of 56% (P = 0.04). Spatiotemporal mapping of spike propagation reveals information flow from onset to areas of spread in epilepsy brain. Surgical resection of the spike onset disrupts the epileptogenic network and may render patients with drug resistant epilepsy seizure-free without having to wait for a seizure to occur during intracranial monitoring.

Reproducibility of interictal spike propagation in children with refractory epilepsy

OBJECTIVE

Interictal spikes are a characteristic feature of invasive electroencephalography (EEG) recordings in children with refractory epilepsy. Spikes frequently co-occur across multiple brain regions with discernable latencies, suggesting that spikes can propagate through distributed neural networks. The purpose of this study was to examine the long-term reproducibility of spike propagation patterns over hours to days of interictal recording.

METHODS

Twelve children (mean age 13.1 years) were retrospectively studied. A mean ± standard deviation (SD) of 47.2 ± 40.1 hours of interictal EEG recordings were examined per patient (range 17.5-166.5 hours). Interictal recordings were divided into 30-minute segments. Networks were extracted based on the frequency of spike coactivation between pairs of electrodes. For each 30-minute segment, electrodes were assigned a "Degree Preference (DP)" based on the tendency to appear upstream or downstream within propagation sequences. The consistency of DPs across segments ("DP-Stability") was quantified using the Spearman rank correlation.

RESULTS

Regions exhibited highly stable preferences to appear upstream, intermediate, or downstream in spike propagation sequences. Across networks, the mean ± SD DP-Stability was 0.88 ± 0.07, indicating that propagation patterns observed in 30-minute segments were representative of the patterns observed in the full interictal window. At the group level, regions involved in seizure generation appeared more upstream in spike propagation sequences.

SIGNIFICANCE

Interictal spike propagation is a highly reproducible output of epileptic networks. These findings shed new light on the spatiotemporal dynamics that may constrain the network mechanisms of refractory epilepsy.

Surgical resection of ripple onset predicts outcome in pediatric epilepsy

OBJECTIVE

In patients with medically refractory epilepsy (MRE), interictal ripples (80-250Hz) are observed in large brain areas whose resection may be unnecessary for seizure freedom. This limits their utility as epilepsy biomarkers for surgery. We assessed the spatiotemporal propagation of interictal ripples on intracranial electroencephalography (iEEG) in children with MRE, compared it with the propagation of spikes, identified ripples that initiated propagation (onset-ripples), and evaluated their clinical value as epilepsy biomarkers.

METHODS

Twenty-seven children who underwent epilepsy surgery were studied. We identified propagation sequences of ripples and spikes across multiple iEEG contacts and calculated each ripple or spike latency from the propagation onset. We classified ripples and spikes into categories (ie, onset, spread, and isolated) based on their spatiotemporal characteristics and correlated their mean rate inside and outside resection with outcome (good outcome, Engel 1 versus poor outcome, Engel≥2). We determined, as onset-zone, spread-zone, and isolated-zone, the areas generating the corresponding ripple or spike category and evaluated the predictive value of their resection.

RESULTS

We observed ripple propagation in all patients and spike propagation in 25 patients. Mean rate of onset-ripples inside resection predicted the outcome (odds ratio = 5.37; p = 0.02) and correlated with Engel class (rho = -0.55; p = 0.003). Resection of the onset-ripple-zone was associated with good outcome (p = 0.047). No association was found for the spread-ripple-zone, isolated-ripple-zone, or any spike-zone.

INTERPRETATION

Interictal ripples propagate across iEEG contacts in children with MRE. The association between the onset-ripple-zone resection and good outcome indicates that onset-ripples are promising epilepsy biomarkers, which estimate the epileptogenic tissue better than spread-ripples or onset-spikes. Ann Neurol 2018;84:331-346.

Recurring Functional Interactions Predict Network Architecture of Interictal and Ictal States in Neocortical Epilepsy

Human epilepsy patients suffer from spontaneous seizures, which originate in brain regions that also subserve normal function. Prior studies demonstrate focal, neocortical epilepsy is associated with dysfunction, several hours before seizures. How does the epileptic network perpetuate dysfunction during baseline periods? To address this question, we developed an unsupervised machine learning technique to disentangle patterns of functional interactions between brain regions, or subgraphs, from dynamic functional networks constructed from approximately 100 h of intracranial recordings in each of 22 neocortical epilepsy patients. Using this approach, we found: (1) subgraphs from ictal (seizure) and interictal (baseline) epochs are topologically similar, (2) interictal subgraph topology and dynamics can predict brain regions that generate seizures, and (3) subgraphs undergo slower and more coordinated fluctuations during ictal epochs compared to interictal epochs. Our observations suggest that seizures mark a critical shift away from interictal states that is driven by changes in the dynamical expression of strongly interacting components of the epileptic network.

Spatiotemporal Mapping of Interictal Spike Propagation: A Novel Methodology Applied to Pediatric Intracranial EEG Recordings

Synchronized cortical activity is implicated in both normative cognitive functioning and many neurologic disorders. For epilepsy patients with intractable seizures, irregular synchronization within the epileptogenic zone (EZ) is believed to provide the network substrate through which seizures initiate and propagate. Mapping the EZ prior to epilepsy surgery is critical for detecting seizure networks in order to achieve postsurgical seizure control. However, automated techniques for characterizing epileptic networks have yet to gain traction in the clinical setting. Recent advances in signal processing and spike detection have made it possible to examine the spatiotemporal propagation of interictal spike discharges across the epileptic cortex. In this study, we present a novel methodology for detecting, extracting, and visualizing spike propagation and demonstrate its potential utility as a biomarker for the EZ. Eighteen presurgical intracranial EEG recordings were obtained from pediatric patients ultimately experiencing favorable (i.e., seizure-free, n = 9) or unfavorable (i.e., seizure-persistent, n = 9) surgical outcomes. Novel algorithms were applied to extract multichannel spike discharges and visualize their spatiotemporal propagation. Quantitative analysis of spike propagation was performed using trajectory clustering and spatial autocorrelation techniques. Comparison of interictal propagation patterns revealed an increase in trajectory organization (i.e., spatial autocorrelation) among Sz-Free patients compared with Sz-Persist patients. The pathophysiological basis and clinical implications of these findings are considered.

Simultaneous Recording of Single-Neuron Activities and Broad-Area Intracranial Electroencephalography: Electrode Design and Implantation Procedure

BACKGROUND:

There has been growing interest in clinical single-neuron recording to better understand epileptogenicity and brain function. It is crucial to compare this new information, single-neuronal activity, with that obtained from conventional intracranial electroencephalography during simultaneous recording. However, it is difficult to implant microwires and subdural electrodes during a single surgical operation because the stereotactic frame hampers flexible craniotomy.

OBJECTIVE:

To describe newly designed electrodes and surgical techniques for implanting them with subdural electrodes that enable simultaneous recording from hippocampal neurons and broad areas of the cortical surface.

METHODS:

We designed a depth electrode that does not protrude into the dura and pulsates naturally with the brain. The length and tract of the depth electrode were determined preoperatively between the lateral subiculum and the lateral surface of the temporal lobe. A frameless navigation system was used to insert the depth electrode. Surface grids and ventral strips were placed before and after the insertion of the depth electrodes, respectively. Finally, a microwire bundle was inserted into the lumen of the depth electrode. We evaluated the precision of implantation, the recording stability, and the recording rate with microwire electrodes.

RESULTS:

Depth-microwire electrodes were placed with a precision of 3.6 mm. The mean successful recording rate of single- or multiple-unit activity was 14.8%, which was maintained throughout the entire recording period.

CONCLUSION:

We achieved simultaneous implantation of microwires, depth electrodes, and broad-area subdural electrodes. Our method enabled simultaneous and stable recording of hippocampal single-neuron activities and multichannel intracranial electroencephalography.

Automatic detection of prominent interictal spikes in intracranial EEG: validation of an algorithm and relationsip to the seizure onset zone

OBJECTIVE

To develop an algorithm for the automatic quantitative description and detection of spikes in the intracranial EEG and quantify the relationship between prominent spikes and the seizure onset zone.

METHODS

An algorithm was developed for the quantification of time-frequency properties of spikes (upslope, instantaneous energy, downslope) and their statistical representation in a univariate generalized extreme value distribution. Its performance was evaluated in comparison to expert detection of spikes in intracranial EEG recordings from 10 patients. It was subsequently used in 18 patients to detect prominent spikes and quantify their spatial relationship to the seizure onset area.

RESULTS

The algorithm displayed an average sensitivity of 63.4% with a false detection rate of 3.2 per minute for the detection of individual spikes and an average sensitivity of 88.6% with a false detection rate of 1.4% for the detection of intracranial EEG contacts containing the most prominent spikes. Prominent spikes occurred closer to the seizure onset area than less prominent spikes but they overlapped with it only in a minority of cases (3/18).

CONCLUSIONS

Automatic detection and quantification of the morphology of spikes increases their utility to localize the seizure onset area. Prominent spikes tend to originate mostly from contacts located in the close vicinity of the seizure onset area rather than from within it.

SIGNIFICANCE

Quantitative analysis of time-frequency characteristics and spatial distribution of intracranial spikes provides complementary information that may be useful for the localization of the seizure-onset zone.

Dynamic imaging of seizure activity in pediatric epilepsy patients

OBJECTIVE

To investigate the feasibility of using noninvasive EEG source imaging approach to image continuous seizure activity in pediatric epilepsy patients.

METHODS

Nine pediatric patients with medically intractable epilepsy were included in this study. Eight of the patients had extratemporal lobe epilepsy and one had temporal lobe epilepsy. All of the patients underwent resective surgery and seven of them underwent intracranial EEG (iEEG) monitoring. The ictal EEG was analyzed using a noninvasive dynamic seizure imaging (DSI) approach. The DSI approach separates scalp EEGs into independent components and extracts the spatio-temporal ictal features to achieve dynamic imaging of seizure sources. Surgical resection and intracranial recordings were used to validate the noninvasive imaging results.

RESULTS

The DSI determined seizure onset zones (SOZs) in these patients were localized within or in close vicinity to the surgically resected region. In the seven patients with intracranial monitoring, the estimated seizure onset sources were concordant with the seizure onset zones of iEEG. The DSI also localized the multiple foci involved in the later seizure propagation, which were confirmed by the iEEG recordings.

CONCLUSIONS

Dynamic seizure imaging can noninvasively image the seizure activations in pediatric patients with both temporal and extratemporal lobe epilepsy.

SIGNIFICANCE

EEG seizure imaging can potentially be used to noninvasively image the SOZs and aid the pre-surgical planning in pediatric epilepsy patients.

Source connectivity analysis from MEG and its application to epilepsy source localization

We report an approach to perform source connectivity analysis from MEG, and initially evaluate this approach to interictal MEG to localize epileptogenic foci and analyze interictal discharge propagations in patients with medically intractable epilepsy. Cortical activities were reconstructed from MEG using individual realistic geometry boundary element method head models. Directional connectivity among cortical regions of interest was then estimated using directed transfer function. The MEG source connectivity analysis method was implemented in the eConnectome software, which is open-source and freely available at http://econnectome.umn.edu . As an initial evaluation, the method was applied to study MEG interictal spikes from five epilepsy patients. Estimated primary epileptiform sources were consistent with surgically resected regions, suggesting the feasibility of using cortical source connectivity analysis from interictal MEG for potential localization of epileptiform activities.

Comparison of five directed graph measures for identification of leading interictal epileptic regions

Directed graphs (digraphs) derived from interictal periods of intracerebral EEG (iEEG) recordings can be used to estimate the leading interictal epileptic regions for presurgery evaluations. For this purpose, quantification of the emittance contribution of each node to the rest of digraph is important. However, the usual digraph measures are not very well suited for this quantification. Here, we compare the efficiency of recently introduced local information (LI) measure and a new measure called total global efficiency with classical measures like global efficiency, local efficiency and node degree. For evaluation, the estimated leading interictal epileptic regions based on five measures are compared with seizure onset zones obtained by visual inspection of epileptologists for five patients. The comparison revealed the superior performance of the LI measure. We showed efficiency of different digraph measures for the purpose of source and sink node identification.

Noninvasive cortical imaging of epileptiform activities from interictal spikes in pediatric patients

Improved non-invasive localization of the epileptogenic foci prior to epilepsy surgery would improve surgical outcome in patients with partial seizure disorders. A critical component for the identification of the epileptogenic brain is the analysis of electrophysiological data obtained during ictal activity from prolonged intracranial recordings. The development of a noninvasive means to identify the seizure onset zone (SOZ) would thus play an important role in treating patients with intractable epilepsy. In the present study, we have investigated non-invasive imaging of epileptiform activity in patients with medically intractable epilepsy by means of a cortical potential imaging (CPI) technique. Eight pediatric patients (1M/7F, ages 4-14 years) with intractable partial epilepsy were studied. Each patient had multiple (6 to 14) interictal spikes (IIS) subjected to the CPI analysis. Realistic geometry boundary element head models were built using each individual's MRI in order to maximize the imaging precision. CPI analysis was performed on the IISs, and extrema in the estimated CPI images were compared with SOZs as determined from the ictal electrocorticogram (ECoG) recordings, as well as the resected areas in the patients and surgical outcomes. The distances between the maximum cortical activities of the IISs reflected by the estimated cortical potential distributions and the SOZs were determined to quantitatively evaluate the performance of the CPI in localizing the epileptogenic zone. Ictal ECoG recordings revealed that six patients exhibited a single epileptogenic focus while two patients had multiple foci. In each patient, the CPI results revealed an area of activity overlapping with the SOZs as identified by ictal ECoG. The distance from the extreme of the CPI images at the peak of IIS to the nearest intracranial electrode associated with the onset of the ictal activity was evaluated for each patient and the averaged distance was 4.6mm. In the group of patients studied, the CPI imaged epileptogenic foci were within the resected areas. According to the follow-up of the eight patients included, two were seizure free and six had substantial reduction in seizure frequency. These promising results demonstrate the potential for noninvasive localization of the epileptogenic focus from interictal scalp EEG recordings. Confirmation of our results may have a significant impact on the process of presurgical planning in pediatric patients with intractable epilepsy by dramatically reducing or potentially eliminating the use of intracranial recording.

Identification of epileptogenic foci from causal analysis of ECoG interictal spike activity

OBJECTIVE

In patients with intractable epilepsy, the use of interictal spikes as surrogate markers of the epileptogenic cortex has generated significant interest. Previous studies have suggested that the cortical generators of the interictal spikes are correlated with the epileptogenic cortex as identified from the ictal recordings. We hypothesize that causal analysis of the functional brain networks during interictal spikes are correlated with the clinically-defined epileptogenic zone.

METHODS

We employed a time-varying causality measure, the adaptive directed transfer function (ADTF), to identify the cortical sources of the interictal spike activity in eight patients with medically intractable neocortical-onset epilepsy. The results were then compared to the foci identified by the epileptologists.

RESULTS

In all eight patients, the majority of the ADTF-calculated source activity was observed within the clinically-defined SOZs. Furthermore, in three of the five patients with two separate epileptogenic foci, the calculated source activity was correlated with both cortical sites.

CONCLUSIONS

The ADTF method identified the cortical sources of the interictal spike activity as originating from the same cortical locations as the recorded ictal activity.

SIGNIFICANCE

Evaluation of the sources of the cortical networks obtained during interictal spikes may provide information as to the generators underlying the ictal activity.