Measuring resetting of brain dynamics at epileptic seizures: application of global optimization and spatial synchronization techniques

Estimated EEG functional connectivity and aperiodic component induced by vagal nerve stimulation in patients with drug-resistant epilepsy

Background

Vagal nerve stimulation (VNS) improves seizure frequency and quality of life in patients with drug-resistant epilepsy (DRE), although the exact mechanism is not fully understood. Previous studies have evaluated the effect of VNS on functional connectivity using the phase lag index (PLI), but none has analyzed its effect on EEG aperiodic parameters (offset and exponent), which are highly conserved and related to physiological functions.

Objective

This study aimed to evaluate the effect of VNS on PLI and aperiodic parameters and infer whether these changes correlate with clinical responses in subjects with DRE.

Materials and methods

PLI, exponent, and offset were derived for each epoch (and each frequency band for PLI), on scalp-derived 64-channel EEG traces of 10 subjects with DRE, recorded before and 1 year after VNS. PLI, exponent, and offset were compared before and after VNS for each patient on a global basis, individual scalp regions, and channels and separately in responders and non-responders. A correlation analysis was performed between global changes in PLI and aperiodic parameters and clinical response.

Results

PLI (global and regional) decreased after VNS for gamma and delta bands and increased for an alpha band in responders, but it was not modified in non-responders. Aperiodic parameters after VNS showed an opposite trend in responders vs. non-responders: both were reduced in responders after VNS, but they were increased in non-responders. Changes in aperiodic parameters correlated with the clinical response.

Conclusion

This study explored the action of VNS therapy from a new perspective and identified EEG aperiodic parameters as a new and promising method to analyze the efficacy of neuromodulation.

Identification of Hidden Sources by Estimating Instantaneous Causality in High-Dimensional Biomedical Time Series

The study of connectivity patterns of a system’s variables, such as multi-channel electroencephalograms (EEG), is of utmost importance towards a better understanding of its internal evolutionary mechanisms. Here, the problem of estimating the connectivity network from multivariate time series in the presence of prominent unobserved variables is addressed. The causality measure of partial mutual information from mixed embedding (PMIME), designed to estimate direct lag-causal effects in the presence of many observed variables, is adapted to estimate also zero-lag effects, the so-called instantaneous causality. We term the proposed advanced method, PMIME0. The estimation of instantaneous causality by PMIME0 is a signature of the presence of hidden source in the observed system, as demonstrated analytically in a toy model. It is further demonstrated that the PMIME0 identifies the true instantaneous with great accuracy in a variety of high-dimensional dynamical systems. The method is applied to EEG data with epileptiform discharges (EDs), and the results imply a strong impact of unobserved confounders during the EDs. This finding comes as a possible explanation for the increased levels of causality during epileptic seizures estimated by some measures affected by the presence of a common source.

Multivariate Matching Pursuit Decomposition and Normalized Gabor Entropy for Quantification of Preictal Trends in Epilepsy

Quantification of the complexity of signals recorded concurrently from multivariate systems, such as the brain, plays an important role in the study and characterization of their state and state transitions. Multivariate analysis of the electroencephalographic signals (EEG) over time is conceptually most promising in unveiling the global dynamics of dynamical brain disorders such as epilepsy. We employed a novel methodology to study the global complexity of the epileptic brain en route to seizures. The developed measures of complexity were based on Multivariate Matching Pursuit (MMP) decomposition of signals in terms of time–frequency Gabor functions (atoms) and Shannon entropy. The measures were first validated on simulation data (Lorenz system) and then applied to EEGs from preictal (before seizure onsets) periods, recorded by intracranial electrodes from eight patients with temporal lobe epilepsy and a total of 42 seizures, in search of global trends of complexity before seizures onset. Out of five Gabor measures of complexity we tested, we found that our newly defined measure, the normalized Gabor entropy (NGE), was able to detect statistically significant (p < 0.05) nonlinear trends of the mean global complexity across all patients over 1 h periods prior to seizures’ onset. These trends pointed to a slow decrease of the epileptic brain’s global complexity over time accompanied by an increase of the variance of complexity closer to seizure onsets. These results show that the global complexity of the epileptic brain decreases at least 1 h prior to seizures and imply that the employed methodology and measures could be useful in identifying different brain states, monitoring of seizure susceptibility over time, and potentially in seizure prediction.

Predictability and Resetting in a Case of Convulsive Status Epilepticus

In this case study, we present evidence of resetting of brain dynamics following convulsive status epilepticus (SE) that was treated successfully with antiepileptic medications (AEDs). The measure of effective inflow (EI), a novel measure of network connectivity, was applied to the continuously recorded multichannel intracranial stereoelectroencephalographic (SEEG) signals before, during and after SE. Results from this analysis indicate trends of progressive reduction of EI over hours up to the onset of SE, mainly at sites of the epileptogenic focus with reversal of those trends upon successful treatment of SE by AEDs. The proposed analytical framework is promising for elucidation of the pathology of neuronal network dynamics that could lead to SE, evaluation of the efficacy of SE treatment strategies, as well as the development of biomarkers for susceptibility to SE.

Direct Causal Networks for the Study of Transcranial Magnetic Stimulation Effects on Focal Epileptiform Discharges

Background: Transcranial magnetic stimulation (TMS) can have inhibitory effects on epileptiform discharges (EDs) of patients with focal seizures. However, the brain connectivity before, during and after EDs, with or without the administration of TMS, has not been extensively explored. Objective: To investigate the brain network of effective connectivity during ED with and without TMS in patients with focal seizures. Methods: For the effective connectivity a direct causality measure is applied termed partial mutual information from mixed embedding (PMIME). TMS-EEG data from two patients with focal seizures were analyzed. Each EEG record contained a number of EDs in the majority of which TMS was administered over the epileptic focus. As a control condition, sham stimulation over the epileptogenic zone or real TMS at a distance from the epileptic focus was also performed. The change in brain connectivity structure was investigated from the causal networks formed at each sliding window. Conclusion: The PMIME could detect distinct changes in the network structure before, within, and after ED. The administration of real TMS over the epileptic focus, in contrast to sham stimulation, terminated the ED prematurely in a node-specific manner and regained the network structure as if it would have terminated spontaneously.

Epileptic focus localization based on resting state interictal MEG recordings is feasible irrespective of the presence or absence of spikes

OBJECTIVE

To investigate whether epileptogenic focus localization is possible based on resting state connectivity analysis of magnetoencephalographic (MEG) data.

METHODS

A multivariate autoregressive (MVAR) model was constructed using the sensor space data and was projected to the source space using lead field and inverse matrix. The generalized partial directed coherence was estimated from the MVAR model in the source space. The dipole with the maximum information inflow was hypothesized to be within the epileptogenic focus.

RESULTS

Applying the focus localization algorithm (FLA) to the interictal MEG recordings from five patients with neocortical epilepsy, who underwent presurgical evaluation for the identification of epileptogenic focus, we were able to correctly localize the focus, on the basis of maximum interictal information inflow in the presence or absence of interictal epileptic spikes in the data, with three out of five patients undergoing resective surgery and being seizure free since.

CONCLUSION

Our preliminary results suggest that accurate localization of the epileptogenic focus may be accomplished using noninvasive spontaneous "resting-state" recordings of relatively brief duration and without the need to capture definite interictal and/or ictal abnormalities.

SIGNIFICANCE

Epileptogenic focus localization is possible through connectivity analysis of resting state MEG data irrespective of the presence/absence of spikes.

A novel dynamic update framework for epileptic seizure prediction

Epileptic seizure prediction is a difficult problem in clinical applications, and it has the potential to significantly improve the patients' daily lives whose seizures cannot be controlled by either drugs or surgery. However, most current studies of epileptic seizure prediction focus on high sensitivity and low false-positive rate only and lack the flexibility for a variety of epileptic seizures and patients' physical conditions. Therefore, a novel dynamic update framework for epileptic seizure prediction is proposed in this paper. In this framework, two basic sample pools are constructed and updated dynamically. Furthermore, the prediction model can be updated to be the most appropriate one for the prediction of seizures' arrival. Mahalanobis distance is introduced in this part to solve the problem of side information, measuring the distance between two data sets. In addition, a multichannel feature extraction method based on Hilbert-Huang transform and extreme learning machine is utilized to extract the features of a patient's preseizure state against the normal state. At last, a dynamic update epileptic seizure prediction system is built up. Simulations on Freiburg database show that the proposed system has a better performance than the one without update. The research of this paper is significantly helpful for clinical applications, especially for the exploitation of online portable devices.

A novel spatiotemporal analysis of peri-ictal spiking to probe the relation of spikes and seizures in epilepsy

The relation between epileptic spikes and seizures is an important but still unresolved question in epilepsy research. Preclinical and clinical studies have produced inconclusive results on the causality or even on the existence of such a relation. We set to investigate this relation taking in consideration seizure severity and spatial extent of spike rate. We developed a novel automated spike detection algorithm based on morphological filtering techniques and then tested the hypothesis that there is a pre-ictal increase and post-ictal decrease of the spatial extent of spike rate. Peri-ictal (around seizures) spikes were detected from intracranial EEG recordings in 5 patients with temporal lobe epilepsy. The 94 recorded seizures were classified into two classes, based on the percentage of brain sites having higher or lower rate of spikes in the pre-ictal compared to post-ictal periods, with a classification accuracy of 87.4%. This seizure classification showed that seizures with increased pre-ictal spike rate and spatial extent compared to the post-ictal period were mostly (83%) clinical seizures, whereas no such statistically significant (α = 0.05) increase was observed peri-ictally in 93% of sub-clinical seizures. These consistent across patients results show the existence of a causal relation between spikes and clinical seizures, and imply resetting of the preceding spiking process by clinical seizures.

Resetting of brain dynamics: epileptic versus psychogenic nonepileptic seizures

We investigated the possibility of differential diagnosis of patients with epileptic seizures (ES) and patients with psychogenic nonepileptic seizures (PNES) through an advanced analysis of the dynamics of the patients' scalp EEGs. The underlying principle was the presence of resetting of brain's preictal spatiotemporal entrainment following onset of ES and the absence of resetting following PNES. Long-term (days) scalp EEGs recorded from five patients with ES and six patients with PNES were analyzed. It was found that: (1) Preictal entrainment of brain sites was reset at ES (P<0.05) in four of the five patients with ES, and not reset (P=0.28) in the fifth patient. (2) Resetting did not occur (p>0.1) in any of the six patients with PNES. These preliminary results in patients with ES are in agreement with our previous findings from intracranial EEG recordings on resetting of brain dynamics by ES and are expected to constitute the basis for the development of a reliable and supporting tool in the differential diagnosis between ES and PNES. Finally, we believe that these results shed light on the electrophysiology of PNES by showing that occurrence of PNES does not assist patients in overcoming a pathological entrainment of brain dynamics. This article is part of a Supplemental Special Issue entitled The Future of Automated Seizure Detection and Prediction.

ABNORMAL INTERICTAL GAMMA ACTIVITY MAY MANIFEST A SEIZURE ONSET ZONE IN TEMPORAL LOBE EPILEPSY

Even though recent studies have suggested that seizures do not occur suddenly and that before a seizure there is a period with an increased probability of seizure occurrence, neurophysiological mechanisms of interictal and pre-seizure states are unknown. The ability of mathematical methods to provide much more sensitive tools for the detection of subtle changes in the electrical activity of the brain gives promise that electrophysiological markers of enhanced seizure susceptibility can be found even during interictal periods when EEG of epilepsy patients often looks 'normal'. Previously, we demonstrated in animals that hippocampal and neocortical gamma-band rhythms (30–100 Hz) intensify long before seizures caused by systemic infusion of kainic acid. Other studies in recent years have also drawn attention to the fast activity (>30 Hz) as a possible marker of epileptogenic tissue. The current study quantified gamma-band activity during interictal periods and seizures in intracranial EEG (iEEG) in 5 patients implanted with subdural grids/intracranial electrodes during their pre-surgical evaluation. In all our patients, we found distinctive (abnormal) bursts of gamma activity with a 3 to 100 fold increase in power at gamma frequencies with respect to selected by clinicians, quiescent, artifact-free, 7–20 min "normal" background (interictal) iEEG epochs 1 to 14 hours prior to seizures. Increases in gamma activity were largest in those channels which later displayed the most intensive electrographic seizure discharges. Moreover, location of gamma-band bursts correlated (with high specificity, 96.4% and sensitivity, 83.8%) with seizure onset zone (SOZ) determined by clinicians. Spatial localization of interictal gamma rhythms within SOZ suggests that the persistent presence of abnormally intensified gamma rhythms in the EEG may be an important tool for focus localization and possibly a determinant of epileptogenesis.

TRANSCUTANEOUS FOCAL ELECTRICAL STIMULATION VIA CONCENTRIC RING ELECTRODES REDUCES SYNCHRONY INDUCED BY PENTYLENETETRAZOLE IN BETA AND GAMMA BANDS IN RATS

Epilepsy is a neurological disorder that affects approximately one percent of the world population. Anti-epileptic drugs are ineffective in 25~30% of cases. Electrical stimulation to control seizures may be an additive therapy. We applied noninvasive transcutaneous focal electrical stimulation (TFES) via concentric ring electrodes on the scalp of rats after inducing seizures with pentylenetetrazole. We found a significant increase in synchrony within the beta-gamma bands during seizures and that TFES significantly reduced the synchrony of the beta-gamma activity and increased synchrony in the delta band.

Seizure prediction and its applications

Epilepsy is characterized by intermittent, paroxysmal, hypersynchronous electrical activity that may remain localized and/or spread and severely disrupt the brain's normal multitask and multiprocessing function. Epileptic seizures are the hallmarks of such activity. The ability to issue warnings in real time of impending seizures may lead to novel diagnostic tools and treatments for epilepsy. Applications may range from a warning to the patient to avert seizure-associated injuries, to automatic timely administration of an appropriate stimulus. Seizure prediction could become an integral part of the treatment of epilepsy through neuromodulation, especially in the new generation of closed-loop seizure control systems.

Predicting temporal lobe epileptic seizures based on zero-crossing interval analysis in scalp EEG

A novel real-time patient-specific algorithm to predict epileptic seizures is proposed. The method is based on the analysis of the positive zero-crossing intervals in the scalp electroencephalogram (EEG), describing the brain dynamics. In a moving-window analysis, the histogram of these intervals in each EEG epoch is computed, and the distribution of the histogram value in specific bins, selected using interictal and preictal references, is estimated based on the values obtained from the current epoch and the epochs of the last 5 min. The resulting distribution for each selected bin is then compared to two reference distributions (interictal and preictal), and a seizure prediction index is developed. Comparing this index with a patient-specific threshold for all EEG channels, a seizure prediction alarm is finally generated. The algorithm was tested on approximately 15.5 hours of multichannel scalp EEG recordings from three patients with temporal lobe epilepsy, including 14 seizures. 86% of seizures were predicted with an average prediction time of 20.8 min and a false prediction rate of 0.12/hr.

Automated real-time epileptic seizure detection in scalp EEG recordings using an algorithm based on wavelet packet transform

A novel wavelet-based algorithm for real-time detection of epileptic seizures using scalp EEG is proposed. In a moving-window analysis, the EEG from each channel is decomposed by wavelet packet transform. Using wavelet coefficients from seizure and nonseizure references, a patient-specific measure is developed to quantify the separation between seizure and nonseizure states for the frequency range of 1-30 Hz. Utilizing this measure, a frequency band representing the maximum separation between the two states is determined and employed to develop a normalized index, called combined seizure index (CSI). CSI is derived for each epoch of every EEG channel based on both rhythmicity and relative energy of that epoch as well as consistency among different channels. Increasing significantly during ictal states, CSI is inspected using one-sided cumulative sum test to generate proper channel alarms. Analyzing alarms from all channels, a seizure alarm is finally generated. The algorithm was tested on scalp EEG recordings from 14 patients, totaling approximately 75.8 h with 63 seizures. Results revealed a high sensitivity of 90.5%, a false detection rate of 0.51 h(-1) and a median detection delay of 7 s. The algorithm could also lateralize the focus side for patients with temporal lobe epilepsy.

An entropy-based approach to predict seizures in temporal lobe epilepsy using scalp EEG

We describe a novel algorithm for the prediction of epileptic seizures using scalp EEG. The method is based on the analysis of the positive zero-crossing interval series of the EEG signal and its first and second derivatives as a measure of brain dynamics. In a moving-window analysis, we estimated the probability density of these intervals and computed the differential entropy. The resultant entropy time series were then inspected using the cumulative sum (CUSUM) procedure to detect decreases as precursors of upcoming seizures. In the next step, the alarm sequences resulting from analysis of the EEG waveform and its derivatives were combined. Finally, a seizure prediction index was generated based on the spatio-temporal processing of the combined CUSUM alarms. We evaluated our algorithm using a dataset of approximately 21.5 hours of multichannel scalp EEG recordings from four patients with temporal lobe epilepsy, resulting in 87.5% sensitivity, a false prediction rate of 0.28/hr, and an average prediction time of 25 min.

Advances in the application of technology to epilepsy: the CIMIT/NIO Epilepsy Innovation Summit

In 2008, a group of clinicians, scientists, engineers, and industry representatives met to discuss advances in the application of engineering technologies to the diagnosis and treatment of patients with epilepsy. The presentations also provided a guide for further technological development, specifically in the evaluation of patients for epilepsy surgery, seizure onset detection and seizure prediction, intracranial treatment systems, and extracranial treatment systems. This article summarizes the discussions and demonstrates that cross-disciplinary interactions can catalyze collaborations between physicians and engineers to address and solve many of the pressing unmet needs in epilepsy.