Effective connectivity analysis of iEEG and accurate localization of the epileptogenic focus at the onset of operculo-insular seizures

Localizing epileptogenic zones with high-frequency oscillations and directed connectivity

PURPOSE

Precise localization of the epileptogenic zone (EZ) is essential for epilepsy surgery. Existing methods often fail to detect slow onset patterns or similar neural activities presented in the recorded signals. To address this issue, we propose a new measure to quantify epileptogenicity, i.e., the connectivity high-frequency epileptogenicity index (cHFEI).

METHODS

The cHFEI method combines directed connectivity and high-frequency oscillations (HFOs) to measure the epileptogenicity of regions involved in a brain network. By applying this method to stereoelectroencephalography (SEEG) recordings of 49 seizures in 20 patients, we calculated the accuracy, sensitivity, and precision with a visually identified epileptogenic zone as a reference. The performance was evaluated by the confusion matrix and the area under the receiver operating characteristic (ROC) curve.

RESULTS

Epileptic network estimation based on cHFEI successfully distinguished brain regions involved in seizure onset from the propagation network. Moreover, cHFEI outperformed other existing detection methods in the estimation of EZs in all patients, with an average area under the ROC curve of 0.88 and an accuracy of 0.85.

CONCLUSIONS

cHFEI can characterize EZ in a robust manner despite various seizure onset patterns and has potential application in epilepsy therapy.

Cross-bispectrum connectivity of intracranial EEG: A novel approach to seizure onset zone localization

Connectivity analyses of intracranial electroencephalography (iEEG) could guide surgical planning for epilepsy surgery by improving the delineation of the seizure onset zone. Traditional approaches fail to quantify important interactions between frequency components. To assess if effective connectivity based on cross-bispectrum -a measure of nonlinear multivariate cross-frequency coupling- can quantitatively identify generators of seizure activity, cross-bispectrum connectivity between channels was computed from iEEG recordings of 5 patients (34 seizures) with good postsurgical outcome. Personalized thresholds of 50% and 80% of the maximum coupling values were used to identify generating electrode channels. In all patients, outflow coupling between α (8-15 Hz) and β (16-31 Hz) frequencies identified at least one electrode inside the resected seizure onset zone. With the 50% and 80% thresholds respectively, an average of 5 (44.7%; specificity = 82.6%) and 2 (22.5%; specificity = 99.0%) resected electrodes were correctly identified. Results show promise for the automatic identification of the seizure onset zone based on cross-bispectrum connectivity analysis.

An application of dynamical directed connectivity of ictal intracranial EEG recordings in seizure onset zone localization

BACKGROUND

Identification of the seizure onset zone (SOZ) is a challenging task in epilepsy surgery. Patients with epilepsy have an altered brain network, allowing connectivity-based analyses to have a great potential in SOZ identification. We investigated a dynamical directed connectivity analysis utilizing ictal intracranial electroencephalographic (iEEG) recordings and proposed an algorithm for SOZ identification based on grouping iEEG contacts.

NEW METHODS

Granger Causality was used for directed connectivity analysis in this study. The intracranial contacts were grouped into visually detected contacts (VDCs), which were identified as SOZ by epileptologists, and non-resected contacts (NRCs). The intragroup and intergroup directed connectivity for VDCs and NRCs were calculated around seizure onset. We then proposed an algorithm for SOZ identification based on the cross-correlation of intragroup outflow and inflow of SOZ candidate contacts.

RESULTS

Our results revealed that the intragroup connectivity of VDCs (VDC→VDC) was significantly larger than the intragroup connectivity of NRCs (NRC→NRC) and the intergroup connectivity between NRCs and VDCs (NRC→VDC) around seizure onset. We found that the proposed algorithm had 90.1 % accuracy for SOZ identification in the seizure-free patients.

COMPARISON WITH EXISTING METHODS

The existing connectivity-based methods for SOZ identification often use either outflow or inflow. In this study, SOZ contacts were identified by integrating outflow and inflow based on the cross correlation between these two measures.

CONCLUSIONS

The proposed group-based dynamical connectivity analysis in this study can aid our understanding of underlying seizure network and may be used to assist in identifying the SOZ contacts before epilepsy surgery.

Source-sink connectivity: a novel interictal EEG marker for seizure localization

Over 15 million epilepsy patients worldwide have drug-resistant epilepsy. Successful surgery is a standard of care treatment but can only be achieved through complete resection or disconnection of the epileptogenic zone, the brain region(s) where seizures originate. Surgical success rates vary between 20% and 80%, because no clinically validated biological markers of the epileptogenic zone exist. Localizing the epileptogenic zone is a costly and time-consuming process, which often requires days to weeks of intracranial EEG (iEEG) monitoring. Clinicians visually inspect iEEG data to identify abnormal activity on individual channels occurring immediately before seizures or spikes that occur interictally (i.e. between seizures). In the end, the clinical standard mainly relies on a small proportion of the iEEG data captured to assist in epileptogenic zone localization (minutes of seizure data versus days of recordings), missing opportunities to leverage these largely ignored interictal data to better diagnose and treat patients. IEEG offers a unique opportunity to observe epileptic cortical network dynamics but waiting for seizures increases patient risks associated with invasive monitoring. In this study, we aimed to leverage interictal iEEG data by developing a new network-based interictal iEEG marker of the epileptogenic zone. We hypothesized that when a patient is not clinically seizing, it is because the epileptogenic zone is inhibited by other regions. We developed an algorithm that identifies two groups of nodes from the interictal iEEG network: those that are continuously inhibiting a set of neighbouring nodes ('sources') and the inhibited nodes themselves ('sinks'). Specifically, patient-specific dynamical network models were estimated from minutes of iEEG and their connectivity properties revealed top sources and sinks in the network, with each node being quantified by source-sink metrics. We validated the algorithm in a retrospective analysis of 65 patients. The source-sink metrics identified epileptogenic regions with 73% accuracy and clinicians agreed with the algorithm in 93% of seizure-free patients. The algorithm was further validated by using the metrics of the annotated epileptogenic zone to predict surgical outcomes. The source-sink metrics predicted outcomes with an accuracy of 79% compared to an accuracy of 43% for clinicians' predictions (surgical success rate of this dataset). In failed outcomes, we identified brain regions with high metrics that were untreated. When compared with high frequency oscillations, the most commonly proposed interictal iEEG feature for epileptogenic zone localization, source-sink metrics outperformed in predictive power (by a factor of 1.2), suggesting they may be an interictal iEEG fingerprint of the epileptogenic zone.

Decoding Intracranial EEG With Machine Learning: A Systematic Review

Advances in intracranial electroencephalography (iEEG) and neurophysiology have enabled the study of previously inaccessible brain regions with high fidelity temporal and spatial resolution. Studies of iEEG have revealed a rich neural code subserving healthy brain function and which fails in disease states. Machine learning (ML), a form of artificial intelligence, is a modern tool that may be able to better decode complex neural signals and enhance interpretation of these data. To date, a number of publications have applied ML to iEEG, but clinician awareness of these techniques and their relevance to neurosurgery, has been limited. The present work presents a review of existing applications of ML techniques in iEEG data, discusses the relative merits and limitations of the various approaches, and examines potential avenues for clinical translation in neurosurgery. One-hundred-seven articles examining artificial intelligence applications to iEEG were identified from 3 databases. Clinical applications of ML from these articles were categorized into 4 domains: i) seizure analysis, ii) motor tasks, iii) cognitive assessment, and iv) sleep staging. The review revealed that supervised algorithms were most commonly used across studies and often leveraged publicly available timeseries datasets. We conclude with recommendations for future work and potential clinical applications.

L1-norm based time-varying brain neural network and its application to dynamic analysis for motor imagery

OBJECTIVE

EEG-based motor imagery (MI) brain-computer interface offers a promising way to improve the efficiency of motor rehabilitation and motor skill learning. In recent years, the power of dynamic network analysis for MI classification has been proved. In fact, its usability mainly depends on the accurate estimation of brain connection. However, traditional dynamic network estimation strategies such as adaptive directed transfer function (ADTF) are designed in the L2-norm. Usually, they estimate a series of pseudo connections caused by outliers, which results in biased features and further limits its online application. Thus, how to accurately infer dynamic causal relationship under outlier influence is urgent.

APPROACH

In this work, we proposed a novel ADTF, which solves the dynamic system in the L1-norm space (L1-ADTF), so as to restrict the outlier influence. To enhance its convergence, we designed an iteration strategy with the alternating direction method of multipliers (ADMM), which could be used for the solution of the dynamic state-space model restricted in the L1-norm space. Furthermore, we compared L1-ADTF to traditional ADTF and its dual extension across both simulation and real EEG experiments.

MAIN RESULTS

A quantitative comparison between L1-ADTF and other ADTFs in simulation studies demonstrates that fewer bias errors and more desirable dynamic state transformation patterns can be captured by the L1-ADTF. Application to real MI EEG datasets seriously noised by ocular artifacts also reveals the efficiency of the proposed L1-ADTF approach to extract the time-varying brain neural network patterns, even when more complex noises are involved.

SIGNIFICANCE

The L1-ADTF may not only be capable of tracking time-varying brain network state drifts robustly but may also be useful in solving a wide range of dynamic systems such as trajectory tracking problems and dynamic neural networks.

Structural Connectivity Alterations in Operculo-Insular Epilepsy

Operculo-insular epilepsy (OIE) is an under-recognized condition that can mimic temporal and extratemporal epilepsies. Previous studies have revealed structural connectivity changes in the epileptic network of focal epilepsy. However, most reports use the debated streamline-count to quantify ‘connectivity strength’ and rely on standard tracking algorithms. We propose a sophisticated cutting-edge method that is robust to crossing fibers, optimizes cortical coverage, and assigns an accurate microstructure-reflecting quantitative conectivity marker, namely the COMMIT (Convex Optimization Modeling for Microstructure Informed Tractography)-weight. Using our pipeline, we report the connectivity alterations in OIE. COMMIT-weighted matrices were created in all participants (nine patients with OIE, eight patients with temporal lobe epilepsy (TLE), and 22 healthy controls (HC)). In the OIE group, widespread increases in ‘connectivity strength’ were observed bilaterally. In OIE patients, ‘hyperconnections’ were observed between the insula and the pregenual cingulate gyrus (OIE group vs. HC group) and between insular subregions (OIE vs. TLE). Graph theoretic analyses revealed higher connectivity within insular subregions of OIE patients (OIE vs. TLE). We reveal, for the first time, the structural connectivity distribution in OIE. The observed pattern of connectivity in OIE likely reflects a diffuse epileptic network incorporating insular-connected regions and may represent a structural signature and diagnostic biomarker.

Dynamic analysis on simultaneous iEEG-MEG data via hidden Markov model

BACKGROUND

Intracranial electroencephalography (iEEG) recordings are used for clinical evaluation prior to surgical resection of the focus of epileptic seizures and also provide a window into normal brain function. A major difficulty with interpreting iEEG results at the group level is inconsistent placement of electrodes between subjects making it difficult to select contacts that correspond to the same functional areas. Recent work using time delay embedded hidden Markov model (HMM) applied to magnetoencephalography (MEG) resting data revealed a distinct set of brain states with each state engaging a specific set of cortical regions. Here we use a rare group dataset with simultaneously acquired resting iEEG and MEG to test whether there is correspondence between HMM states and iEEG power changes that would allow classifying iEEG contacts into functional clusters.

METHODS

Simultaneous MEG-iEEG recordings were performed at rest on 11 patients with epilepsy whose intracranial electrodes were implanted for pre-surgical evaluation. Pre-processed MEG sensor data was projected to source space. Time delay embedded HMM was then applied to MEG time series. At the same time, iEEG time series were analyzed with time-frequency decomposition to obtain spectral power changes with time. To relate MEG and iEEG results, correlations were computed between HMM probability time courses of state activation and iEEG power time course from the mid contact pair for each electrode in equally spaced frequency bins and presented as correlation spectra for the respective states and iEEG channels. Association of iEEG electrodes with HMM states based on significant correlations was compared to that based on the distance to peaks in subject-specific state topographies.

RESULTS

Five HMM states were inferred from MEG. Two of them corresponded to the left and the right temporal activations and had a spectral signature primarily in the theta/alpha frequency band. All the electrodes had significant correlations with at least one of the states (p < 0.05 uncorrected) and for 27/50 electrodes these survived within-subject FDR correction (q < 0.05). These correlations peaked in the theta/alpha band. There was a highly significant dependence between the association of states and electrodes based on functional correlations and that based on spatial proximity (p = 5.6e-6,χ2 test for independence). Despite the potentially atypical functional anatomy and physiological abnormalities related to epilepsy, HMM model estimated from the patient group was very similar to that estimated from healthy subjects.

CONCLUSION

Epilepsy does not preclude HMM analysis of interictal data. The resulting group functional states are highly similar to those reported for healthy controls. Power changes recorded with iEEG correlate with HMM state time courses in the alpha-theta band and the presence of this correlation can be related to the spatial location of electrode contacts close to the individual peaks of the corresponding state topographies. Thus, the hypothesized relation between iEEG contacts and HMM states exists and HMM could be further explored as a method for identifying comparable iEEG channels across subjects for the purposes of group analysis.

Seizure onset detection based on detection of changes in brain activity quantified by evolutionary game theory model

BACKGROUND AND OBJECTIVE

Epilepsy is one of the most common diseases of the nervous system, affecting about 1% of the world's population. The unpredictable nature of the epilepsy seizures deprives the patients and those around them of living a normal life. Therefore, the development of new methods that can help these patients will increase the life quality of these people and can bring a lot of economic savings in the health sector.

METHODS

In this study, we introduced a new framework for seizure onset detection. Our framework provides a new modelling for brain activity using evolutionary game theory and Kalman filter. If the patterns in the electroencephalogram (EEG) signal violate the predicted patterns by the proposed model, using a novel detection algorithm that has been also introduced in this paper, it can be determined whether the observed violation is the result of the onset of an epileptic seizure or not.

RESULTS

The proposed approach was able to detect the onset of all the seizures in CHB-MIT dataset with an average delay of -0.8 s and a false alarm of 0.39 per hour. Also, our proposed approach is about 20 times faster compared to recent studies.

CONCLUSIONS

The experimental results of applying the proposed framework on the CHB-MIT dataset show that our framework not only performed well with respect to the sensitivity, delay, and false alarm metrics but also performed much better in terms of run time compared to recent studies. This appropriate run time, along with other suitable metrics, makes it possible to use this framework in many cases where processing power or energy is limited and to think about creating new and inexpensive solutions for the treatment and care of people diagnosed with epilepsy.

Large-Scale Desynchronization During Interictal Epileptic Discharges Recorded With Intracranial EEG

It is increasingly recognized that deep understanding of epileptic seizures requires both localizing and characterizing the functional network of the region where they are initiated, i. e., the epileptic focus. Previous investigations of the epileptogenic focus' functional connectivity have yielded contrasting results, reporting both pathological increases and decreases during resting periods and seizures. In this study, we shifted paradigm to investigate the time course of connectivity in relation to interictal epileptiform discharges. We recruited 35 epileptic patients undergoing intracranial EEG (iEEG) investigation as part of their presurgical evaluation. For each patient, 50 interictal epileptic discharges (IEDs) were marked and iEEG signals were epoched around those markers. Signals were narrow-band filtered and time resolved phase-locking values were computed to track the dynamics of functional connectivity during IEDs. Results show that IEDs are associated with a transient decrease in global functional connectivity, time-locked to the peak of the discharge and specific to the high range of the gamma frequency band. Disruption of the long-range connectivity between the epileptic focus and other brain areas might be an important process for the generation of epileptic activity. Transient desynchronization could be a potential biomarker of the epileptogenic focus since 1) the functional connectivity involving the focus decreases significantly more than the connectivity outside the focus and 2) patients with good surgical outcome appear to have a significantly more disconnected focus than patients with bad outcomes.

Localization of the epileptogenic zone based on ictal stereo-electroencephalogram: Brain network and single-channel signal feature analysis

Accurate localization of the epileptogenic zone (EZ) is crucial for refractory focal epilepsy patients to achieve freedom from seizures following epilepsy surgery. In this study, ictal stereo-electroencephalography data from 35 patients with refractory focal epilepsy were analyzed. Effective networks based on partial directed coherence were analyzed, and a gray level co-occurrence matrix was applied to extract the time-varying features of the in-degree. These features, combined with the single-channel signal time-frequency features, including approximate entropy and line length, were used to localize the EZ based on a cluster algorithm. For all seizure-free patients (n = 23), the proposed method was effective in identifying the clinical-EZ-contacts and clinical-EZ-blocks, with an F1-score of 62.47 % and 72.18 %, respectively. The sensitivity was 96.00 % for the clinical-EZ-block identification, which provided the information for the decision-making of clinicians, prompting clinicians to focus on the identified EZ-blocks and their nearby contacts. The agreement between the EZ identified by the proposed method and the clinical-EZ was worse for non-seizure-free patients (n = 12) than for seizure-free patients. Furthermore, our method provided better results than using only brain network or single-channel signal features. This suggests that combining these complementary features can facilitate more accurate localization of the EZ.

History of the Network Approach in Epilepsy Surgery

We provide a history and overview of the network approach to epilepsy surgery. Models of the epileptogenic zone (EZ) have evolved considerably over the years with more recent models accounting for the connectivity and network properties of epileptic foci. Next, we describe several examples of network phenotypes of focal epilepsy and how these have the potential to influence surgical decision-making and patient outcome. Future research will provide new insight into how network models of the EZ can determine optimal surgical interventions that improve seizure outcomes and optimize cognitive outcomes.