Localization of Epileptogenic Zone With the Correction of Pathological Networks

A Multivariate Approach to Quantifying Risk Factors Impacting Stereotactic Robotic-Guided Stereoelectroencephalography

BACKGROUND AND OBJECTIVES

Stereoelectroencephalography (SEEG) is an important method for invasive monitoring to establish surgical candidacy in approximately half of refractory epilepsy patients. Identifying factors affecting lead placement can mitigate potential surgical risks. This study applies multivariate analyses to identify perioperative factors affecting stereotactic electrode placement.

METHODS

We collected registration and accuracy data for consecutive patients undergoing SEEG implantation between May 2022 and November 2023. Stereotactic robotic guidance, using intraoperative imaging and a novel frame-based fiducial, was used for planning and SEEG implantation. Entry-point (EE), target-point (TE), and angular errors were measured, and statistical univariate and multivariate linear regression analyses were performed.

RESULTS

Twenty-seven refractory epilepsy patients (aged 15-57 years) undergoing SEEG were reviewed. Sixteen patients had unilateral implantation (10 left-sided, 6 right-sided); 11 patients underwent bilateral implantation. The mean number of electrodes per patient was 18 (SD = 3) with an average registration mean error of 0.768 mm (SD = 0.108). Overall, 486 electrodes were reviewed. Univariate analysis showed significant correlations of lead error with skull thickness (EE: P = .003; TE: P = .012); entry angle (EE: P < .001; TE: P < .001; angular error: P = .030); lead length (TE: P = .020); and order of electrode implantation (EE: P = .003; TE: P = .001). Three multiple linear regression models were used. All models featured predictors of implantation region (157 temporal, 241 frontal, 79 parietal, 9 occipital); skull thickness (mean = 5.80 mm, SD = 2.97 mm); order (range: 1-23); and entry angle in degrees (mean = 75.47, SD = 11.66). EE and TE error models additionally incorporated lead length (mean = 44.08 mm, SD = 13.90 mm) as a predictor. Implantation region and entry angle were significant predictors of error (P ≤ .05).

CONCLUSION

Our study identified 2 primary predictors of SEEG lead error, region of implantation and entry angle, with nonsignificant contributions from lead length or order of electrode placement. Future considerations for SEEG may consider varying regional approaches and angles for more optimal accuracy in lead placement.

Predicting radiofrequency thermocoagulation surgical outcomes in refractory focal epilepsy patients using functional coupled neural mass model

Objective

The success rate of achieving seizure freedom after radiofrequency thermocoagulation surgery for patients with refractory focal epilepsy is about 20-40%. This study aims to enhance the prediction of surgical outcomes based on preoperative decisions through network model simulation, providing a reference for clinicians to validate and optimize surgical plans.

Methods

Twelve patients with epilepsy who underwent radiofrequency thermocoagulation were retrospectively reviewed in this study. A coupled model based on model subsets of the neural mass model was constructed by calculating partial directed coherence as the coupling matrix from stereoelectroencephalography (SEEG) signals. Multi-channel time-varying model parameters of excitation and inhibitions were identified by fitting the real SEEG signals with the coupled model. Further incorporating these model parameters, the coupled model virtually removed contacts destroyed in radiofrequency thermocoagulation or selected randomly. Subsequently, the coupled model after virtual surgery was simulated.

Results

The identified excitatory and inhibitory parameters showed significant difference before and after seizure onset (p < 0.05), and the trends of parameter changes aligned with the seizure process. Additionally, excitatory parameters of epileptogenic contacts were higher than that of non-epileptogenic contacts, and opposite findings were noticed for inhibitory parameters. The simulated signals of postoperative models to predict surgical outcomes yielded an area under the curve (AUC) of 83.33% and an accuracy of 91.67%.

Conclusion

The multi-channel coupled model proposed in this study with physiological characteristics showed a desirable performance for preoperatively predicting patients' prognoses.

Bifurcations and bursting in the Epileptor

The Epileptor is a phenomenological model for seizure activity that is used in a personalized large-scale brain modeling framework, the Virtual Epileptic Patient, with the aim of improving surgery outcomes for drug-resistant epileptic patients. Transitions between interictal and ictal states are modeled as bifurcations, enabling the definition of seizure classes in terms of onset/offset bifurcations. This establishes a taxonomy of seizures grounded in their essential underlying dynamics and the Epileptor replicates the activity of the most common class, as observed in patients with focal epilepsy, which is characterized by square-wave bursting properties. The Epileptor also encodes an additional mechanism to account for interictal spikes and spike and wave discharges. Here we use insights from a more generic model for square-wave bursting, based on the Unfolding Theory approach, to guide the bifurcation analysis of the Epileptor and gain a deeper understanding of the model and the role of its parameters. We show how the Epileptor's parameters can be modified to produce activities for other seizures classes of the taxonomy, as observed in patients, so that the large-scale brain models could be further personalized. Some of these classes have already been described in the literature in the Epileptor, others, predicted by the generic model, are new. Finally, we unveil how the interaction with the additional mechanism for spike and wave discharges alters the bifurcation structure of the main burster.

Kindling in humans: Does secondary epileptogenesis occur?

The relevance of secondary epileptogenesis for human epilepsy remains a controversial subject decades after it was first described in animal models. Whether or not a previously normal brain region can become independently epileptogenic through a kindling-like process has not, and cannot, be definitely proven in humans. Rather than reliance on direct experimental evidence, attempts to answering this question must depend on observational data. In this review, observations based largely upon contemporary surgical series will advance the case for secondary epileptogenesis in humans. As will be argued, hypothalamic hamartoma-related epilepsy provides the strongest case for this process; all the stages of secondary epileptogenesis can be observed. Hippocampal sclerosis (HS) is another pathology where the question of secondary epileptogenesis frequently arises, and observations from bitemporal and dual pathology series are explored. The verdict here is far more difficult to reach, in large part because of the scarcity of longitudinal cohorts; moreover, recent experimental data have challenged the claim that HS is acquired consequent to recurrent seizures. Synaptic plasticity more than seizure-induced neuronal injury is the likely mechanism of secondary epileptogenesis. Postoperative running-down phenomenon provides the best evidence that a kindling-like process occurs in some patients, evidenced by its reversal. Finally, a network perspective of secondary epileptogenesis is considered, as well as the possible role for subcortical surgical interventions.

Seizure onset zone identification using phase-amplitude coupling and multiple machine learning approaches for interictal electrocorticogram

Automatic seizure onset zone (SOZ) localization using interictal electrocorticogram (ECoG) improves the diagnosis and treatment of patients with medically refractory epilepsy. This study aimed to investigate the characteristics of phase-amplitude coupling (PAC) extracted from interictal ECoG and the feasibility of PAC serving as a promising biomarker for SOZ identification. We employed the mean vector length modulation index approach on the 20-s ECoG window to calculate PAC features between low-frequency rhythms (0.5-24 Hz) and high frequency oscillations (HFOs) (80-560 Hz). We used statistical measures to test the significant difference in PAC between the SOZ and non-seizure onset zone (NSOZ). To overcome the drawback of handcraft feature engineering, we established novel machine learning models to learn automatically the characteristics of the obtained PAC features and classify them to identify the SOZ. Besides, to handle imbalanced dataset classification, we introduced novel feature-wise/class-wise re-weighting strategies in conjunction with classifiers. In addition, we proposed a time-series nest cross-validation to provide more accurate and unbiased evaluations for this model. Seven patients with focal cortical dysplasia were included in this study. The experiment results not only showed that a significant coupling at band pairs of slow waves and HFOs exists in the SOZ when compared with the NSOZ, but also indicated the effectiveness of the PAC features and the proposed models in achieving better classification performance .

Personalized dynamic network models of the human brain as a future tool for planning and optimizing epilepsy therapy

Epilepsy is a common neurological disorder, with one third of patients not responding to currently available antiepileptic drugs. The proportion of pharmacoresistant epilepsies has remained unchanged for many decades. To cure epilepsy and control seizures requires a paradigm shift in the development of new approaches to epilepsy diagnosis and treatment. Contemporary medicine has benefited from the exponential growth of computational modeling, and the application of network dynamics theory to understanding and treating human brain disorders. In epilepsy, the introduction of these approaches has led to personalized epileptic network modeling that can explore the patient's seizure genesis and predict the functional impact of resection on its individual network's propensity to seize. The application of the dynamic systems approach to neurostimulation therapy of epilepsy allows designing stimulation strategies that consider the patient's seizure dynamics and long-term fluctuations in the stability of their epileptic networks. In this article, we review, in a nontechnical fashion suitable for a broad neuroscientific audience, recent progress in personalized dynamic brain network modeling that is shaping the future approach to the diagnosis and treatment of epilepsy.

The association of magnetoencephalography high-frequency oscillations with epilepsy types and a ripple-based method with source-level connectivity for mapping epilepsy sources

OBJECTIVE

To explore the association between high-frequency oscillations (HFOs) and epilepsy types and to improve the accuracy of source localization.

METHODS

Magnetoencephalography (MEG) ripples of 63 drug-resistant epilepsy patients were detected. Ripple rates, distribution, spatial complexity, and the clustering coefficient of ripple channels were used for the preliminary classification of lateral temporal lobe epilepsy (LTLE), mesial temporal lobe epilepsy (MTLE), and nontemporal lobe epilepsy (NTLE), mainly frontal lobe epilepsy (FLE). Furthermore, the seizure site identification was improved using the Tucker LCMV method and source-level betweenness centrality.

RESULTS

Ripple rates were significantly higher in MTLE than in LTLE and NTLE (p < 0.05). The LTLE and MTLE were mainly distributed in the temporal lobe, followed by the parietal lobe, occipital lobe, and frontal lobe, whereas MTLE ripples were mainly distributed in the frontal lobe, then parietal lobe and occipital lobe. Nevertheless, the NTLE ripples were primarily in the frontal lobe and partially in the occipital lobe (p < 0.05). Meanwhile, the spatial complexity of NTLE was significantly higher than that of LTLE and MTLE and was lowest in MTLE (p < 0.01). However, an opposite trend was observed for the standardized clustering coefficient compared with spatial complexity (p < 0.01). Finally, the tucker algorithm showed a higher percentage of ripples at the surgical site when the betweenness centrality was added (p < 0.01).

CONCLUSION

This study demonstrated that HFO rates, distribution, spatial complexity, and clustering coefficient of ripple channels varied considerably among the three epilepsy types. Additionally, tucker MEG estimation combined with ripple rates based on the source-level functional connectivity is a promising approach for presurgical epilepsy evaluation.

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.

High-Frequency Oscillations and Epileptogenic Network

Epilepsy is a network disease caused by aberrant neocortical large-scale connectivity spanning regions on the scale of several centimeters. High-frequency oscillations, characterized by the 80-600 Hz signals in electroencephalography, have been proven to be a promising biomarker of epilepsy that can be used in assessing the severity and susceptibility of epilepsy as well as the location of the epileptogenic zone. However, the presence of a high-frequency oscillation network remains a topic of debate as high-frequency oscillations have been previously thought to be incapable of propagation, and the relationship between high-frequency oscillations and the epileptogenic network has rarely been discussed. Some recent studies reported that high-frequency oscillations may behave like networks that are closely relevant to the epileptogenic network. Pathological highfrequency oscillations are network-driven phenomena and elucidate epileptogenic network development; high-frequency oscillations show different characteristics coincident with the epileptogenic network dynamics, and cross-frequency coupling between high-frequency oscillations and other signals may mediate the generation and propagation of abnormal discharges across the network.

Epidemic models characterize seizure propagation and the effects of epilepsy surgery in individualized brain networks based on MEG and invasive EEG recordings

Epilepsy surgery is the treatment of choice for drug-resistant epilepsy patients. However, seizure-freedom is currently achieved in only 2/3 of the patients after surgery. In this study we have developed an individualized computational model based on MEG brain networks to explore seizure propagation and the efficacy of different virtual resections. Eventually, the goal is to obtain individualized models to optimize resection strategy and outcome. We have modelled seizure propagation as an epidemic process using the susceptible-infected (SI) model on individual brain networks derived from presurgical MEG. We included 10 patients who had received epilepsy surgery and for whom the surgery outcome at least one year after surgery was known. The model parameters were tuned in in order to reproduce the patient-specific seizure propagation patterns as recorded with invasive EEG. We defined a personalized search algorithm that combined structural and dynamical information to find resections that maximally decreased seizure propagation for a given resection size. The optimal resection for each patient was defined as the smallest resection leading to at least a 90% reduction in seizure propagation. The individualized model reproduced the basic aspects of seizure propagation for 9 out of 10 patients when using the resection area as the origin of epidemic spreading, and for 10 out of 10 patients with an alternative definition of the seed region. We found that, for 7 patients, the optimal resection was smaller than the resection area, and for 4 patients we also found that a resection smaller than the resection area could lead to a 100% decrease in propagation. Moreover, for two cases these alternative resections included nodes outside the resection area. Epidemic spreading models fitted with patient specific data can capture the fundamental aspects of clinically observed seizure propagation, and can be used to test virtual resections in silico. Combined with optimization algorithms, smaller or alternative resection strategies, that are individually targeted for each patient, can be determined with the ultimate goal to improve surgery outcome. MEG-based networks can provide a good approximation of structural connectivity for computational models of seizure propagation, and facilitate their clinical use.

Deep-layer motif method for estimating information flow between EEG signals

Accurate identification for the information flow between epileptic seizure signals is the key to construct the directional epileptic brain network which can be used to localize epileptic focus. In this paper, our concern is on how to improve the direction identification of information flow and also investigate how it can be cut off or weakened. In view of this, we propose the deep-layer motif method. Based on the directional index (DI) estimation using permutation conditional mutual information, the effectiveness of the proposed deep-layer motif method is numerically assessed with the coupled mass neural model. Furthermore, we investigate the robustness of this method in considering the interference of autaptic coupling, time delay and short-term plasticity. Results show that compared to the simple 1-layer motif method, the 2nd- and 3rd-layer motif methods have the dominant enhancement effects for the direction identification. In particular, deep-layer motif method possesses good anti-jamming performance and good robustness in calculating DI. In addition, we investigate the effect of deep brain stimulation (DBS) on the information flow. It is found that this deep-layer motif method is still superior to the single-layer motif method in direction identification and is robust to weak DBS. However, the high-frequency strong DBS can effectively decrease the DI suggesting the weakened information flow. These results may give new insights into the seizure detection and control.

Optimization of epilepsy surgery through virtual resections on individual structural brain networks

The success of epilepsy surgery in patients with refractory epilepsy depends upon correct identification of the epileptogenic zone (EZ) and an optimal choice of the resection area. In this study we developed individualized computational models based upon structural brain networks to explore the impact of different virtual resections on the propagation of seizures. The propagation of seizures was modelled as an epidemic process [susceptible-infected-recovered (SIR) model] on individual structural networks derived from presurgical diffusion tensor imaging in 19 patients. The candidate connections for the virtual resection were all connections from the clinically hypothesized EZ, from which the seizures were modelled to start, to other brain areas. As a computationally feasible surrogate for the SIR model, we also removed the connections that maximally reduced the eigenvector centrality (EC) (large values indicate network hubs) of the hypothesized EZ, with a large reduction meaning a large effect. The optimal combination of connections to be removed for a maximal effect were found using simulated annealing. For comparison, the same number of connections were removed randomly, or based on measures that quantify the importance of a node or connection within the network. We found that 90% of the effect (defined as reduction of EC of the hypothesized EZ) could already be obtained by removing substantially less than 90% of the connections. Thus, a smaller, optimized, virtual resection achieved almost the same effect as the actual surgery yet at a considerably smaller cost, sparing on average 27.49% (standard deviation: 4.65%) of the connections. Furthermore, the maximally effective connections linked the hypothesized EZ to hubs. Finally, the optimized resection was equally or more effective than removal based on structural network characteristics both regarding reducing the EC of the hypothesized EZ and seizure spreading. The approach of using reduced EC as a surrogate for simulating seizure propagation can suggest more restrictive resection strategies, whilst obtaining an almost optimal effect on reducing seizure propagation, by taking into account the unique topology of individual structural brain networks of patients.

Sharp decrease in the Laplacian matrix rank of phase-space graphs: a potential biomarker in epilepsy

In this paper, phase space reconstruction from stereo-electroencephalography data of ten patients with focal epilepsy forms a series of graphs. Those obtained graphs reflect the transition characteristics of brain dynamical system from pre-seizure to seizure of epilepsy. Interestingly, it is found that the rank of Laplacian matrix of these graphs has a sharp decrease when a seizure is close to happen, which thus might be viewed as a new potential biomarker in epilepsy. In addition, the reliability of this method is numerically verified with a coupled mass neural model. In particular, our simulation suggests that this potential biomarker can play the roles of predictive effect or delayed awareness, depending on the bias current of the Gaussian noise. These results may give new insights into the seizure detection.

Dynamic Temporospatial Patterns of Functional Connectivity and Alterations in Idiopathic Generalized Epilepsy

The dynamic profile of brain function has received much attention in recent years and is also a focus in the study of epilepsy. The present study aims to integrate the dynamics of temporal and spatial characteristics to provide comprehensive and novel understanding of epileptic dynamics. Resting state fMRI data were collected from eighty-three patients with idiopathic generalized epilepsy (IGE) and 87 healthy controls (HC). Specifically, we explored the temporal and spatial variation of functional connectivity density (tvFCD and svFCD) in the whole brain. Using a sliding-window approach, for a given region, the standard variation of the FCD series was calculated as the tvFCD and the variation of voxel-wise spatial distribution was calculated as the svFCD. We found primary, high-level, and sub-cortical networks demonstrated distinct tvFCD and svFCD patterns in HC. In general, the high-level networks showed the highest variation, the subcortical and primary networks showed moderate variation, and the limbic system showed the lowest variation. Relative to HC, the patients with IGE showed weaken temporal and enhanced spatial variation in the default mode network and weaken temporospatial variation in the subcortical network. Besides, enhanced temporospatial variation in sensorimotor and high-level networks was also observed in patients. The hyper-synchronization of specific brain networks was inferred to be associated with the phenomenon responsible for the intrinsic propensity of generation and propagation of epileptic activities. The disrupted dynamic characteristics of sensorimotor and high-level networks might potentially contribute to the driven motion and cognition phenotypes in patients. In all, presently provided evidence from the temporospatial variation of functional interaction shed light on the dynamics underlying neuropathological profiles of epilepsy.

Validation of virtual resection on intraoperative interictal data acquired during epilepsy surgery

OBJECTIVE

A 'Virtual resection' consists of computationally simulating the effect of an actual resection on the brain. We validated two functional connectivity based virtual resection methods with the actual connectivity measured using post-resection intraoperative recordings.

METHODS

A non-linear association index was applied to pre-resection recordings from 11 extra-temporal focal epilepsy patients. We computed two virtual resection strategies: first, a 'naive' one obtained by simply removing from the connectivity matrix the electrodes that were resected; second, a virtual resection with partialization accounting for the influence of resected electrodes on not-resected electrodes. We validated the virtual resections with two analysis: 1) We tested with a Kolmogorov-Smirnov test if the distributions of connectivity values after the virtual resections differed from the actual post-resection connectivity distribution; 2) we tested if the overall effect of the resection measured by contrasting pre-resection and post-resection connectivity values is detectable with the virtual resection approach using a Kolmogorv-Smirnov test.

RESULTS

The estimation of post-resection connectivity values did not succeed for both methods. In the second analysis, the naive method failed completely to detect the effect found between pre-resection and post-resection connectivity distributions, while the partialization method agreed with post-resection measurements in detecting a drop connectivity compared to pre-resection recordings.

CONCLUSION

Our findings suggest that the partialization technique is superior to the naive method in detecting the overall effect after the resection.

SIGNIFICANCE

We pointed out how a realistic validation based on actual post-resection recordings reveals that virtual resection methods are not yet mature to inform the clinical decision-making.

Removal of Interictal MEG-Derived Network Hubs Is Associated With Postoperative Seizure Freedom

Objective: To investigate whether MEG network connectivity was associated with epilepsy duration, to identify functional brain network hubs in patients with refractory focal epilepsy, and assess if their surgical removal was associated with post-operative seizure freedom. Methods: We studied 31 patients with drug refractory focal epilepsy who underwent resting state magnetoencephalography (MEG), and structural magnetic resonance imaging (MRI) as part of pre-surgical evaluation. Using the structural MRI, we generated 114 cortical regions of interest, performed surface reconstruction and MEG source localization. Representative source localized signals for each region were correlated with each other to generate a functional brain network. We repeated this procedure across three randomly chosen one-minute epochs. Network hubs were defined as those with the highest intra-hemispheric mean correlations. Post-operative MRI identified regions that were surgically removed. Results: Greater mean MEG network connectivity was associated with a longer duration of epilepsy. Patients who were seizure free after surgery had more hubs surgically removed than patients who were not seizure free (AUC = 0.76, p = 0.01) consistently across three randomly chosen time segments. Conclusion: Our results support a growing literature implicating network hub involvement in focal epilepsy, the removal of which by surgery is associated with greater chance of post-operative seizure freedom.

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.

Distinguishing Dependent-Stage Secondary Epileptogenesis in a Complex Case of Giant Hypothalamic Hamartoma With Assistance of a Computational Method

Besides gelastic seizures, hypothalamic hamartoma (HH) is also noted for its susceptibility to remote secondary epileptogenesis. Although clinical observations have demonstrated its existence, and a three-stage theory has been proposed, how to determine whether a remote symptom is spontaneous or dependent on epileptic activities of HH is difficult in some cases. Herein, we report a case of new non-gelastic seizures in a 9-year-old female associated with a postoperatively remaining HH. Electrophysiological examinations and stereo-electroencephalography (SEEG) demonstrated seizure onsets with slow-wave and fast activities on the outside of the HH. By using computational methodologies to calculate the network dynamic effective connectivities, the importance of HH in the epileptic network was revealed. After SEEG-guided thermal coagulation of the remaining HH, the patient finally was seizure-free at the 2-year follow-up. This case showed the ability of computational methods to reveal information underlying complex SEEG signals, and further demonstrated the dependent-stage secondary epileptogenesis, which has been rarely reported.

Interictal intracranial electroencephalography for predicting surgical success: The importance of space and time

Objective

Predicting postoperative seizure freedom using functional correlation networks derived from interictal intracranial electroencephalography (EEG) has shown some success. However, there are important challenges to consider: (1) electrodes physically closer to each other naturally tend to be more correlated, causing a spatial bias; (2) implantation location and number of electrodes differ between patients, making cross-subject comparisons difficult; and (3) functional correlation networks can vary over time but are currently assumed to be static.

Methods

In this study, we address these three challenges using intracranial EEG data from 55 patients with intractable focal epilepsy. Patients additionally underwent preoperative magnetic resonance imaging (MRI), intraoperative computed tomography, and postoperative MRI, allowing accurate localization of electrodes and delineation of the removed tissue.

Results

We show that normalizing for spatial proximity between nearby electrodes improves prediction of postsurgery seizure outcomes. Moreover, patients with more extensive electrode coverage were more likely to have their outcome predicted correctly (area under the receiver operating characteristic curve > 0.9, P « 0.05) but not necessarily more likely to have a better outcome. Finally, our predictions are robust regardless of the time segment analyzed.

Significance

Future studies should account for the spatial proximity of electrodes in functional network construction to improve prediction of postsurgical seizure outcomes. Greater coverage of both removed and spared tissue allows for predictions with higher accuracy.

Identifying the Epileptogenic Zone With the Relative Strength of High-Frequency Oscillation: A Stereoelectroencephalography Study

Background

High-frequency oscillation (HFO) represents a promising biomarker of epileptogenicity. However, the significant interindividual differences among patients limit its application in clinical practice. Here, we applied and evaluated an individualized, frequency-based approach of HFO analysis in stereoelectroencephalography (SEEG) data for localizing the epileptogenic zones (EZs).

Methods

Clinical and SEEG data of 19 patients with drug-resistant focal epilepsy were retrospectively analyzed. The individualized spectral power of all signals recorded by electrode array, i.e., the relative strength of HFO, was computed with a wavelet method for each patient. Subsequently, the clinical value of the relative strength of HFO for identifying the EZ was evaluated.

Results

Focal increase in the relative strength of HFO in SEEG recordings were identified in all 19 patients. HFOs identified inside the clinically identified seizure onset zone had more spectral power than those identified outside (p < 0.001), and HFOs in 250–500 Hz band (fast ripples) seemed to be more specific identifying the EZ than in those in 80–250 Hz band (ripples) (p < 0.01). The resection of brain regions generating HFOs resulted in a favorable seizure outcome in 17 patients (17/19; 89.5%), while in the cases of other patients with poor outcomes, the brain regions generating HFOs were not removed completely.

Conclusion

The relative strength of HFO, especially fast ripples, is a promising effective biomarker for identifying the EZ and can lead to a favorable seizure outcome if used to guide epilepsy surgery.

High interictal connectivity within the resection zone is associated with favorable post-surgical outcomes in focal epilepsy patients

Patients with drug-resistant focal epilepsy are often candidates for invasive surgical therapies. In these patients, it is necessary to accurately localize seizure generators to ensure seizure freedom following intervention. While intracranial electroencephalography (iEEG) is the gold standard for mapping networks for surgery, this approach requires inducing and recording seizures, which may cause patient morbidity. The goal of this study is to evaluate the utility of mapping interictal (non-seizure) iEEG networks to identify targets for surgical treatment. We analyze interictal iEEG recordings and neuroimaging from 27 focal epilepsy patients treated via surgical resection. We generate interictal functional networks by calculating pairwise correlation of iEEG signals across different frequency bands. Using image coregistration and segmentation, we identify electrodes falling within surgically resected tissue (i.e. the resection zone), and compute node-level and edge-level synchrony in relation to the resection zone. We further associate these metrics with post-surgical outcomes. Greater overlap between resected electrodes and highly synchronous electrodes is associated with favorable post-surgical outcomes. Additionally, good-outcome patients have significantly higher connectivity localized within the resection zone compared to those with poorer postoperative seizure control. This finding persists following normalization by a spatially-constrained null model. This study suggests that spatially-informed interictal network synchrony measures can distinguish between good and poor post-surgical outcomes. By capturing clinically-relevant information during interictal periods, our method may ultimately reduce the need for prolonged invasive implants and provide insights into the pathophysiology of an epileptic brain. We discuss next steps for translating these findings into a prospectively useful clinical tool.

•

We analyze interictal iEEG recordings and neuroimaging from epilepsy patients

•

We find that high interictal strength selectivity is associated with better outcomes

•

This effect appears to be driven largely by connectivity within the resection zone

•

Interictal recordings can guide identification of seizure-generating networks

Identification of epileptic brain states by dynamic functional connectivity analysis of simultaneous EEG-fMRI: a dictionary learning approach

Most fMRI studies of the brain's intrinsic functional connectivity (FC) have assumed that this is static; however, it is now clear that it changes over time. This is particularly relevant in epilepsy, which is characterized by a continuous interchange between epileptic and normal brain states associated with the occurrence of epileptic activity. Interestingly, recurrent states of dynamic FC (dFC) have been found in fMRI data using unsupervised learning techniques, assuming either their sparse or non-sparse combination. Here, we propose an l₁-norm regularized dictionary learning (l₁-DL) approach for dFC state estimation, which allows an intermediate and flexible degree of sparsity in time, and demonstrate its application in the identification of epilepsy-related dFC states using simultaneous EEG-fMRI data. With this l₁-DL approach, we aim to accommodate a potentially varying degree of sparsity upon the interchange between epileptic and non-epileptic dFC states. The simultaneous recording of the EEG is used to extract time courses representative of epileptic activity, which are incorporated into the fMRI dFC state analysis to inform the selection of epilepsy-related dFC states. We found that the proposed l₁-DL method performed best at identifying epilepsy-related dFC states, when compared with two alternative methods of extreme sparsity (k-means clustering, maximum; and principal component analysis, minimum), as well as an l₀-norm regularization framework (l₀-DL), with a fixed amount of temporal sparsity. We further showed that epilepsy-related dFC states provide novel insights into the dynamics of epileptic networks, which go beyond the information provided by more conventional EEG-correlated fMRI analysis, and which were concordant with the clinical profile of each patient. In addition to its application in epilepsy, our study provides a new dFC state identification method of potential relevance for studying brain functional connectivity dynamics in general.