Dynamic Network Connectivity Analysis to Identify Epileptogenic Zones Based on Stereo-Electroencephalography

High frequency oscillation network dynamics predict outcome in non-palliative epilepsy surgery

High frequency oscillations are a promising biomarker of outcome in intractable epilepsy. Prior high frequency oscillation work focused on counting high frequency oscillations on individual channels, and it is still unclear how to translate those results into clinical care. We show that high frequency oscillations arise as network discharges that have valuable properties as predictive biomarkers. Here, we develop a tool to predict patient outcome before surgical resection is performed, based on only prospective information. In addition to determining high frequency oscillation rate on every channel, we performed a correlational analysis to evaluate the functional connectivity of high frequency oscillations in 28 patients with intracranial electrodes. We found that high frequency oscillations were often not solitary events on a single channel, but part of a local network discharge. Eigenvector and outcloseness centrality were used to rank channel importance within the connectivity network, then used to compare patient outcome by comparison with the seizure onset zone or a proportion within the proposed resected channels (critical resection percentage). Combining the knowledge of each patient's seizure onset zone resection plan along with our computed high frequency oscillation network centralities and high frequency oscillation rate, we develop a Naïve Bayes model that predicts outcome (positive predictive value: 100%) better than predicting based upon fully resecting the seizure onset zone (positive predictive value: 71%). Surgical margins had a large effect on outcomes: non-palliative patients in whom most of the seizure onset zone was resected ('definitive surgery', ≥ 80% resected) had predictable outcomes, whereas palliative surgeries (<80% resected) were not predictable. These results suggest that the addition of network properties of high frequency oscillations is more accurate in predicting patient outcome than seizure onset zone alone in patients with most of the seizure onset zone removed and offer great promise for informing clinical decisions in surgery for refractory epilepsy.

SEEG-based epileptic seizure network modeling and analysis for pre-surgery evaluation

Stereo-electroencephalography is a minimally invasive technique for patients with refractory epilepsy pursuing surgery to reduce or control seizures. Electrodes are implanted based on pre-surgery evaluations and can collect deep brain activities for surgery decisions. This paper presents a methodology to analyze stereo-electroencephalography and assist clinicians by recommending the optimal surgical option and target areas for focal epilepsy patients. A seizure network (graph) model is proposed to characterize the spatial distribution and temporal changes of ictal events. The network nodes and edges correspond to specific epileptogenic regions and propagation/impact pathways (weighted by directed transfer function), respectively. We then employ a K-means clustering strategy to group nodes into a few clusters, from which the target surgical areas can be identified. Ten patients with different types of focal seizures were thoroughly analyzed. Promising consistency between results of our method's recommendations, clinical decisions and surgery outcomes were observed.

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.

Neurobehavioural comorbidities of epilepsy: towards a network-based precision taxonomy

Cognitive and behavioural comorbidities are prevalent in childhood and adult epilepsies and impose a substantial human and economic burden. Over the past century, the classic approach to understanding the aetiology and course of these comorbidities has been through the prism of the medical taxonomy of epilepsy, including its causes, course, characteristics and syndromes. Although this 'lesion model' has long served as the organizing paradigm for the field, substantial challenges to this model have accumulated from diverse sources, including neuroimaging, neuropathology, neuropsychology and network science. Advances in patient stratification and phenotyping point towards a new taxonomy for the cognitive and behavioural comorbidities of epilepsy, which reflects the heterogeneity of their clinical presentation and raises the possibility of a precision medicine approach. As we discuss in this Review, these advances are informing the development of a revised aetiological paradigm that incorporates sophisticated neurobiological measures, genomics, comorbid disease, diversity and adversity, and resilience factors. We describe modifiable risk factors that could guide early identification, treatment and, ultimately, prevention of cognitive and broader neurobehavioural comorbidities in epilepsy and propose a road map to guide future research.

Efficacy of cathodal transcranial direct current stimulation on electroencephalographic functional networks in patients with focal epilepsy: Preliminary findings

BACKGROUND

Neuromodulation is a promising therapeutic alternative for epilepsy. We aimed to explore the efficacy and safety of cathodal transcranial current direct stimulation (ctDCS) on electroencephalographic functional networks in focal epilepsy.

METHODS

A sham-controlled, double-blinded, randomized study was conducted on 25 participants with focal epilepsy who underwent a 5-day, -1.0 mA, 20 min ctDCS, which targeted at the most active interictal epileptiform discharge (IED) region. We examined the electroencephalograms (EEGs) at baseline, immediately and at 4 weeks following ctDCS. The graph theory-based brain networks were established through time-variant partial directed coherence (TVPDC), and were calculated between each pair of EEG signals. The functional networks were characterized using average clustering coefficient, characteristic path length, and small-worldness index. The seizure frequencies, IEDs, graph-theory metrics and cognitive tests were compared.

RESULTS

Preliminary findings indicated an IED reduction of 30.2% at the end of 5-day active ctDCS compared to baseline (p < 0.10) and a significant IED reduction of 33.4% 4 weeks later (p < 0.05). In terms of the EEG functional network, the small-worldness index significantly reduced by 3.5% (p < 0.05) and the characteristic path length increased by 1.8% (p < 0.10) at the end of the session compared to the baseline. No obvious change was found in the seizure frequency during follow-up (p > 0.05). The Mini-Mental State Examination (MMSE) showed no difference between the active and sham groups (p > 0.05). No severe adverse reactions were observed.

CONCLUSIONS

In focal epilepsy, the 5-day consecutive ctDCS may potentially decrease the IEDs and ameliorate the EEG functional network, proposing a novel personalized therapeutic scenario for epilepsy.

Regional and global resting-state functional MR connectivity in temporal lobe epilepsy: Results from the Epilepsy Connectome Project

Temporal lobe epilepsy (TLE) has been conceptualized as focal disease with a discrete neurobiological focus and can respond well to targeted resection or ablation. In contrast, the neuro-cognitive deficits resulting from TLE can be widespread involving regions beyond the primary epileptic network. We hypothesize that this seemingly paradoxical findings can be explained by differences in connectivity between the primary epileptic region which is hyper-connected and its secondary influence on global connectome organization. This hypothesis is tested using regional and global graph theory metrics where we anticipate that regional mesial-temporal hyperconnectivity will be found and correlate with seizure frequency while global networks will be disorganized and be more closely associated with neuro-cognitive deficits. Resting-state fMRI was used to examine temporal lobe regional connectivity and global functional connectivity from 102 patients with TLE and 55 controls. Connectivity matrices were calculated for subcortical volumes and cortical parcellations. Graph theory metrics (global clustering coefficient (GCC), degree, closeness) were compared between groups and in relation to neuropsychological profiles and disease covariates using permutation testing and causal analysis. In TLE there was a decrease in GCC (p = 0.0345) associated with a worse neuropsychological profile (p = 0.0134). There was increased connectivity in the left hippocampus/amygdala (degree p = 0.0103, closeness p = 0.0104) and a decrease in connectivity in the right lateral temporal lobe (degree p = 0.0186, closeness p = 0.0122). A ratio between the hippocampus/amygdala and lateral temporal lobe-temporal lobe connectivity ratio (TLCR) revealed differences between TLE and controls for closeness (left p = 0.00149, right p = 0.0494) and for degree on left p = 0.00169; with trend on right p = 0.0567. Causal analysis suggested that "Epilepsy Activity" (seizure frequency, anti-seizure medications) was associated with increase in TLCR but not in GCC, while cognitive decline was associated with decreased GCC. These findings support the hypothesis that in TLE there is hyperconnectivity in the hippocampus/amygdala and hypoconnectivity in the lateral temporal lobe associated with "Epilepsy Activity." While, global connectome disorganization was associated with worse neuropsychological phenotype.

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.

Seizure-onset regions demonstrate high inward directed connectivity during resting-state: An SEEG study in focal epilepsy

OBJECTIVE

In patients with medically refractory focal epilepsy, stereotactic-electroencephalography (SEEG) can aid in localizing epileptogenic regions for surgical treatment. SEEG, however, requires long hospitalizations to record seizures, and ictal interpretation can be incomplete or inaccurate. Our recent work showed that non-directed resting-state analyses may identify brain regions as epileptogenic or uninvolved. Our present objective is to map epileptogenic networks in greater detail and more accurately identify seizure-onset regions using directed resting-state SEEG connectivity.

METHODS

In 25 patients with focal epilepsy who underwent SEEG, 2 minutes of resting-state, artifact-free, SEEG data were selected and functional connectivity was estimated. Using standard clinical interpretation, brain regions were classified into four categories: ictogenic, early propagation, irritative, or uninvolved. Three non-directed connectivity measures (mutual information [MI] strength, and imaginary coherence between and within regions) and four directed measures (partial directed coherence [PDC] and directed transfer function [DTF], inward and outward strength) were calculated. Logistic regression was used to generate a predictive model of ictogenicity.

RESULTS

Ictogenic regions had the highest and uninvolved regions had the lowest MI strength. Although both PDC and DTF inward strengths were highest in ictogenic regions, outward strengths did not differ among categories. A model incorporating directed and nondirected connectivity measures demonstrated an area under the receiver-operating characteristic (ROC) curve (AUC) of 0.88 in predicting ictogenicity of individual regions. The AUC of this model was 0.93 when restricted to patients with favorable postsurgical seizure outcomes.

SIGNIFICANCE

Directed connectivity measures may help identify epileptogenic networks without requiring ictal recordings. Greater inward but not outward connectivity in ictogenic regions at rest may represent broad inhibitory input to prevent seizure generation.

Model-based design for seizure control by stimulation

OBJECTIVE

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

APPROACH

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

MAIN RESULTS

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

SIGNIFICANCE

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

[Isolated effective coherence analysis of epileptogenic networks in temporal lobe epilepsy using stereo-electroencephalography]

Stereo-electroencephalography (SEEG) is widely used to record the electrical activity of patients' brain in clinical. The SEEG-based epileptogenic network can better describe the origin and the spreading of seizures, which makes it an important measure to localize epileptogenic zone (EZ). SEEG data from six patients with refractory epilepsy are used in this study. Five of them are with temporal lobe epilepsy, and the other is with extratemporal lobe epilepsy. The node outflow (out-degree) and inflow (in-degree) of information are calculated in each node of epileptic network, and the overlay between selected nodes and resected nodes is analyzed. In this study, SEEG data is transformed to bipolar montage, and then the epileptic network is established by using independent effective coherence (iCoh) method. The SEEG segments at onset, middle and termination of seizures in Delta, Theta, Alpha, Beta, and Gamma rhythms are used respectively. Finally, the K-means clustering algorithm is applied on the node values of out-degree and in-degree respectively. The nodes in the cluster with high value are compared with the resected regions. The final results show that the accuracy of selected nodes in resected region in the Delta, Alpha and Beta rhythm are 0.90, 0.88 and 0.89 based on out-degree values in temporal lobe epilepsy patients respectively, while the in-degree values cannot differentiate them. In contrast, the out-degree values are higher outside the temporal lobe in the patient with extratemporal lobe epilepsy. Based on the out-degree feature in low-frequency epileptic network, this study provides a potential quantitative measure for identifying patients with temporal lobe epilepsy in clinical.

Precursors of seizures due to specific spatial-temporal modifications of evolving large-scale epileptic brain networks

Knowing when, where, and how seizures are initiated in large-scale epileptic brain networks remains a widely unsolved problem. Seizure precursors – changes in brain dynamics predictive of an impending seizure – can now be identified well ahead of clinical manifestations, but either the seizure onset zone or remote brain areas are reported as network nodes from which seizure precursors emerge. We aimed to shed more light on the role of constituents of evolving epileptic networks that recurrently transit into and out of seizures. We constructed such networks from more than 3200 hours of continuous intracranial electroencephalograms recorded in 38 patients with medication refractory epilepsy. We succeeded in singling out predictive edges and predictive nodes. Their particular characteristics, namely edge weight respectively node centrality (a fundamental concept of network theory), from the pre-ictal periods of 78 out of 97 seizures differed significantly from the characteristics seen during inter-ictal periods. The vast majority of predictive nodes were connected by most of the predictive edges, but these nodes never played a central role in the evolving epileptic networks. Interestingly, predictive nodes were entirely associated with brain regions deemed unaffected by the focal epileptic process. We propose a network mechanism for a transition into the pre-seizure state, which puts into perspective the role of the seizure onset zone in this transition and highlights the necessity to reassess current concepts for seizure generation and seizure prevention.

Phase Synchronization Dynamics of Neural Network during Seizures

Epilepsy has been considered as a network-level disorder characterized by recurrent seizures, which result from network reorganization with evolution of synchronization. In this study, the brain networks were established by calculating phase synchronization based on electrocorticogram (ECoG) signals from eleven refractory epilepsy patients. Results showed that there was a significant increase of synchronization prior to seizure termination and no significant difference of the transitions of network states among the preseizure, seizure, and postseizure periods. Those results indicated that synchronization might participate in termination of seizures, and the network states transitions might not dominate the seizure evolution.

Localization of the Epileptogenic Zone Using Interictal MEG and Machine Learning in a Large Cohort of Drug-Resistant Epilepsy Patients

Objective: Epilepsy surgery results in seizure freedom in the majority of drug-resistant patients. To improve surgery outcome we studied whether MEG metrics combined with machine learning can improve localization of the epileptogenic zone, thereby enhancing the chance of seizure freedom.

Methods: Presurgical interictal MEG recordings of 94 patients (64 seizure-free >1y post-surgery) were analyzed to extract four metrics in source space: delta power, low-to-high-frequency power ratio, functional connectivity (phase lag index), and minimum spanning tree betweenness centrality. At the group level, we estimated the overlap of the resection area with the five highest values for each metric and determined whether this overlap differed between surgery outcomes. At the individual level, those metrics were used in machine learning classifiers (linear support vector machine (SVM) and random forest) to distinguish between resection and non-resection areas and between surgery outcome groups.

Results: The highest values, for all metrics, overlapped with the resection area in more than half of the patients, but the overlap did not differ between surgery outcome groups. The classifiers distinguished the resection areas from non-resection areas with 59.94% accuracy (95% confidence interval: 59.67–60.22%) for SVM and 60.34% (59.98–60.71%) for random forest, but could not differentiate seizure-free from not seizure-free patients [43.77% accuracy (42.08–45.45%) for SVM and 49.03% (47.25–50.82%) for random forest].

Significance: All four metrics localized the resection area but did not distinguish between surgery outcome groups, demonstrating that metrics derived from interictal MEG correspond to expert consensus based on several presurgical evaluation modalities, but do not yet localize the epileptogenic zone. Metrics should be improved such that they correspond to the resection area in seizure-free patients but not in patients with persistent seizures. It is important to test such localization strategies at an individual level, for example by using machine learning or individualized models, since surgery is individually tailored.

A Time-Varying Connectivity Analysis from Distributed EEG Sources: A Simulation Study

Time-varying connectivity analysis based on sources reconstructed using inverse modeling of electroencephalographic (EEG) data is important to understand the dynamic behaviour of the brain. We simulated cortical data from a visual spatial attention network with a time-varying connectivity structure, and then simulated the propagation to the scalp to obtain EEG data. Distributed EEG source modeling using sLORETA was applied. We compared different dipole (representing a source) selection strategies based on their time series in a region of interest. Next, we estimated multivariate autoregressive (MVAR) parameters using classical Kalman filter and general linear Kalman filter approaches followed by the calculation of partial directed coherence (PDC). MVAR parameters and PDC values for the selected sources were compared with the ground-truth. We found that the best strategy to extract the time series of a region of interest was to select a dipole with time series showing the highest correlation with the average time series in the region of interest. Dipole selection based on power or based on the largest singular value offer comparable alternatives. Among the different Kalman filter approaches, the use of a general linear Kalman filter was preferred to estimate PDC based connectivity except when only a small number of trials are available. In the latter case, the classical Kalman filter can be an alternative.