Neural fragility as an EEG marker of the seizure onset zone

Neuronal synchrony and critical bistability: Mechanistic biomarkers for localizing the epileptogenic network

OBJECTIVE

Postsurgical seizure freedom in drug-resistant epilepsy (DRE) patients varies from 30% to 80%, implying that in many cases the current approaches fail to fully map the epileptogenic zone (EZ). We aimed to advance a novel approach to better characterize epileptogenicity and investigate whether the EZ encompasses a broader epileptogenic network (EpiNet) beyond the seizure zone (SZ) that exhibits seizure activity.

METHODS

We first used computational modeling to test putative complex systems-driven and systems neuroscience-driven mechanistic biomarkers for epileptogenicity. We then used these biomarkers to extract features from resting-state stereoelectroencephalograms recorded from DRE patients and trained supervised classifiers to localize the SZ against gold standard clinical localization. To further explore the prevalence of pathological features in an extended brain network outside of the clinically identified SZ, we also used unsupervised classification.

RESULTS

Supervised SZ classification trained on individual features achieved accuracies of .6-.7 area under the receiver operating characteristic curve (AUC). Combining all criticality and synchrony features further improved the AUC to .85. Unsupervised classification discovered an EpiNet-like cluster of brain regions, in which 51% of brain regions were outside of the SZ. Brain regions in the EpiNet-like cluster engaged in interareal hypersynchrony and locally exhibited high-amplitude bistability and excessive inhibition, which was strikingly similar to the high seizure risk regime revealed by our computational modeling.

SIGNIFICANCE

The finding that combining biomarkers improves SZ localization accuracy indicates that the novel mechanistic biomarkers for epileptogenicity employed here yield synergistic information. On the other hand, the discovery of SZ-like brain dynamics outside of the clinically defined SZ provides empirical evidence of an extended pathophysiological EpiNet.

Simulated resections and RNS placement can optimize post-operative seizure outcomes when guided by fast ripple networks

In medication-resistant epilepsy, the goal of epilepsy surgery is to make a patient seizure free with a resection/ablation that is as small as possible to minimize morbidity. The standard of care in planning the margins of epilepsy surgery involves electroclinical delineation of the seizure onset zone (SOZ) and incorporation of neuroimaging findings from MRI, PET, SPECT, and MEG modalities. Resecting cortical tissue generating high-frequency oscillations (HFOs) has been investigated as a more efficacious alternative to targeting the SOZ. In this study, we used a support vector machine (SVM), with four distinct fast ripple (FR: 350-600 Hz on oscillations, 200-600 Hz on spikes) metrics as factors. These metrics included the FR resection ratio (RR), a spatial FR network measure, and two temporal FR network measures. The SVM was trained by the value of these four factors with respect to the actual resection boundaries and actual seizure free labels of 18 patients with medically refractory focal epilepsy. Leave one out cross-validation of the trained SVM in this training set had an accuracy of 0.78. We next used a simulated iterative virtual resection targeting the FR sites that were highest rate and showed most temporal autonomy. The trained SVM utilized the four virtual FR metrics to predict virtual seizure freedom. In all but one of the nine patients seizure free after surgery, we found that the virtual resections sufficient for virtual seizure freedom were larger in volume (p<0.05). In nine patients who were not seizure free, a larger virtual resection made five virtually seizure free. We also examined 10 medically refractory focal epilepsy patients implanted with the responsive neurostimulator system (RNS) and virtually targeted the RNS stimulation contacts proximal to sites generating FR at highest rates to determine if the simulated value of the stimulated SOZ and stimulated FR metrics would trend toward those patients with a better seizure outcome. Our results suggest: 1) FR measures can accurately predict whether a resection, defined by the standard of care, will result in seizure freedom; 2) utilizing FR alone for planning an efficacious surgery can be associated with larger resections; 3) when FR metrics predict the standard of care resection will fail, amending the boundaries of the planned resection with certain FR generating sites may improve outcome; and 4) more work is required to determine if targeting RNS stimulation contact proximal to FR generating sites will improve seizure outcome.

A spatial perturbation framework to validate implantation of the epileptogenic zone

Stereo-electroencephalography (SEEG) is the gold standard to delineate surgical targets in focal drug-resistant epilepsy. SEEG uses electrodes placed directly into the brain to identify the seizure-onset zone (SOZ). However, its major constraint is limited brain coverage, potentially leading to misidentification of the 'true' SOZ. Here, we propose a framework to assess adequate SEEG sampling by coupling epileptic biomarkers with their spatial distribution and measuring the system's response to a perturbation of this coupling. We demonstrate that the system's response is strongest in well-sampled patients when virtually removing the measured SOZ. We then introduce the spatial perturbation map, a tool that enables qualitative assessment of the implantation coverage. Probability modelling reveals a higher likelihood of well-implanted SOZs in seizure-free patients or non-seizure free patients with incomplete SOZ resections, compared to non-seizure-free patients with complete resections. This highlights the framework's value in sparing patients from unsuccessful surgeries resulting from poor SEEG coverage.

Utility of intracranial EEG networks depends on re-referencing and connectivity choice

Studies of intracranial EEG networks have been used to reveal seizure generators in patients with drug-resistant epilepsy. Intracranial EEG is implanted to capture the epileptic network, the collection of brain tissue that forms a substrate for seizures to start and spread. Interictal intracranial EEG measures brain activity at baseline, and networks computed during this state can reveal aberrant brain tissue without requiring seizure recordings. Intracranial EEG network analyses require choosing a reference and applying statistical measures of functional connectivity. Approaches to these technical choices vary widely across studies, and the impact of these technical choices on downstream analyses is poorly understood. Our objective was to examine the effects of different re-referencing and connectivity approaches on connectivity results and on the ability to lateralize the seizure onset zone in patients with drug-resistant epilepsy. We applied 48 pre-processing pipelines to a cohort of 125 patients with drug-resistant epilepsy recorded with interictal intracranial EEG across two epilepsy centres to generate intracranial EEG functional connectivity networks. Twenty-four functional connectivity measures across time and frequency domains were applied in combination with common average re-referencing or bipolar re-referencing. We applied an unsupervised clustering algorithm to identify groups of pre-processing pipelines. We subjected each pre-processing approach to three quality tests: (i) the introduction of spurious correlations; (ii) robustness to incomplete spatial sampling; and (iii) the ability to lateralize the clinician-defined seizure onset zone. Three groups of similar pre-processing pipelines emerged: common average re-referencing pipelines, bipolar re-referencing pipelines and relative entropy-based connectivity pipelines. Relative entropy and common average re-referencing networks were more robust to incomplete electrode sampling than bipolar re-referencing and other connectivity methods (Friedman test, Dunn-Šidák test P < 0.0001). Bipolar re-referencing reduced spurious correlations at non-adjacent channels better than common average re-referencing (Δ mean from machine ref = -0.36 versus -0.22) and worse in adjacent channels (Δ mean from machine ref = -0.14 versus -0.40). Relative entropy-based network measures lateralized the seizure onset hemisphere better than other measures in patients with temporal lobe epilepsy (Benjamini-Hochberg-corrected P < 0.05, Cohen's d: 0.60-0.76). Finally, we present an interface where users can rapidly evaluate intracranial EEG pre-processing choices to select the optimal pre-processing methods tailored to specific research questions. The choice of pre-processing methods affects downstream network analyses. Choosing a single method among highly correlated approaches can reduce redundancy in processing. Relative entropy outperforms other connectivity methods in multiple quality tests. We present a method and interface for researchers to optimize their pre-processing methods for deriving intracranial EEG brain networks.

Generalisability of epileptiform patterns across time and patients

The complexity of localising the epileptogenic zone (EZ) contributes to surgical resection failures in achieving seizure freedom. The distinct patterns of epileptiform activity during interictal and ictal phases, varying across patients, often lead to suboptimal localisation using electroencephalography (EEG) features. We posed two key questions: whether neural signals reflecting epileptogenicity generalise from interictal to ictal time windows within each patient, and whether epileptiform patterns generalise across patients. Utilising an intracranial EEG dataset from 55 patients, we extracted a large battery of simple to complex features from stereo-EEG (SEEG) and electrocorticographic (ECoG) neural signals during interictal and ictal windows. Our features (n = 34) quantified many aspects of the signals including statistical moments, complexities, frequency-domain and cross-channel network attributes. Decision tree classifiers were then trained and tested on distinct time windows and patients to evaluate the generalisability of epileptogenic patterns across time and patients, respectively. Evidence strongly supported generalisability from interictal to ictal time windows across patients, particularly in signal power and high-frequency network-based features. Consistent patterns of epileptogenicity were observed across time windows within most patients, and signal features of epileptogenic regions generalised across patients, with higher generalisability in the ictal window. Signal complexity features were particularly contributory in cross-patient generalisation across patients. These findings offer insights into generalisable features of epileptic neural activity across time and patients, with implications for future automated approaches to supplement other EZ localisation methods.

High-performance prediction of epilepsy surgical outcomes based on the genetic neural networks and hybrid iEEG marker

Accurately identification of the seizure onset zone (SOZ) is pivotal for successful surgery in patients with medically refractory epilepsy. The purpose of this study is to improve the performance of model predicting the epilepsy surgery outcomes using genetic neural network (GNN) model based on a hybrid intracranial electroencephalography (iEEG) marker. We extracted 21 SOZ related markers based on iEEG data from 79 epilepsy patients. The least absolute shrinkage and selection operator (LASSO) regression was employed to integrated seven markers, selected after testing in pairs with all 21 biomarkers and 7 machine learning models, into a hybrid marker. Based on the hybrid marker, we devised a GNN model and compared its predictive performance for surgical outcomes with six other mainstream machine-learning models. Compared to the mainstream models, underpinning the GNN with the hybrid iEEG marker resulted in a better prediction of surgical outcomes, showing a significant increase of the prediction accuracy from approximately 87% to 94.3% (P = 0.0412). This study suggests that the hybrid iEEG marker can improve the performance of model predicting the epilepsy surgical outcomes, and validates the effectiveness of the GNN in characterizing and analyzing complex relationships between clinical data variables.

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.

Advanced Electrode Technologies for Noninvasive Brain-Computer Interfaces

Brain-computer interfaces (BCIs) have garnered significant attention in recent years due to their potential applications in medical, assistive, and communication technologies. Building on this, noninvasive BCIs stand out as they provide a safe and user-friendly method for interacting with the human brain. In this work, we provide a comprehensive overview of the latest developments and advancements in material, design, and application of noninvasive BCIs electrode technology. We also explore the challenges and limitations currently faced by noninvasive BCI electrode technology and sketch out the technological roadmap from three dimensions: Materials and Design; Performances; Mode and Function. We aim to unite research efforts within the field of noninvasive BCI electrode technology, focusing on the consolidation of shared goals and fostering integrated development strategies among a diverse array of multidisciplinary researchers.

Soft Bioelectronics for Neuroengineering: New Horizons in the Treatment of Brain Tumor and Epilepsy

Soft bioelectronic technologies for neuroengineering have shown remarkable progress, which include novel soft material technologies and device design strategies. Such technological advances that are initiated from fundamental brain science are applied to clinical neuroscience and provided meaningful promises for significant improvement in the diagnosis efficiency and therapeutic efficacy of various brain diseases recently. System-level integration strategies in consideration of specific disease circumstances can enhance treatment effects further. Here, recent advances in soft implantable bioelectronics for neuroengineering, focusing on materials and device designs optimized for the treatment of intracranial disease environments, are reviewed. Various types of soft bioelectronics for neuroengineering are categorized and exemplified first, and then details for the sensing and stimulating device components are explained. Next, application examples of soft implantable bioelectronics to clinical neuroscience, particularly focusing on the treatment of brain tumor and epilepsy are reviewed. Finally, an ideal system of soft intracranial bioelectronics such as closed-loop-type fully-integrated systems is presented, and the remaining challenges for their clinical translation are discussed.

Dynamic evolution of the anterior cingulate-insula network during seizures

OBJECTIVES

In physiological situations, the anterior cingulate cortex (ACC) and anterior insular cortex (AIC) are prone to coactivation. The functional connectivity and interaction between ACC and AIC in the context of epilepsy remain unclear. This study aimed to investigate the dynamic coupling between these two brain regions during seizures.

METHODS

Patients who underwent stereoelectroencephalography (SEEG) recording were included in this study. The SEEG data were visually inspected and quantitatively analyzed. The narrowband oscillations and aperiodic components at seizure onset were parameterized. The frequency-specific non-linear correlation analysis was applied to the functional connectivity. The excitation/inhibition ratio (E:I ratio) reflected by the aperiodic slope was performed to evaluate the excitability.

RESULTS

Twenty patients were included in the study, with 10 diagnosed with anterior cingulate epilepsy and 10 with anterior insular epilepsy. In both types of epilepsy, the correlation coefficient (h2 ) between the ACC and AIC at seizure onset exhibited a significantly higher value than that during interictal and preictal periods (p < 0.05). The direction index (D) showed a significant increase at seizure onset, serving as an indicator for the direction of information flow between these two brain regions with up to 90% accuracy. The E:I ratio increased significantly at seizure onset, with the seizure-onset zone (SOZ) demonstrating a more pronounced increase compared to non-SOZ (p < 0.05). For seizures originating from AIC, the E:I ratio was significantly higher in the AIC than in the ACC (p = 0.0364).

CONCLUSIONS

In the context of epilepsy, the ACC and AIC are dynamically coupled during seizures. The functional connectivity and excitability exhibit a significant increase at seizure onset. By analyzing connectivity and excitability, the SOZ in ACC and AIC can be identified. The direction index (D) serves as an indicator for the direction of information flow from SOZ to non-SOZ. Notably, the excitability of SOZ changes more significantly than that of non-SOZ.

Passive and active markers of cortical excitability in epilepsy

Electroencephalography (EEG) has been the primary diagnostic tool in clinical epilepsy for nearly a century. Its review is performed using qualitative clinical methods that have changed little over time. However, the intersection of higher resolution digital EEG and analytical tools developed in the past decade invites a re-exploration of relevant methodology. In addition to the established spatial and temporal markers of spikes and high-frequency oscillations, novel markers involving advanced postprocessing and active probing of the interictal EEG are gaining ground. This review provides an overview of the EEG-based passive and active markers of cortical excitability in epilepsy and of the techniques developed to facilitate their identification. Several different emerging tools are discussed in the context of specific EEG applications and the barriers we must overcome to translate these tools into clinical practice.

Complementary structural and functional abnormalities to localise epileptogenic tissue

BACKGROUND

When investigating suitability for epilepsy surgery, people with drug-refractory focal epilepsy may have intracranial EEG (iEEG) electrodes implanted to localise seizure onset. Diffusion-weighted magnetic resonance imaging (dMRI) may be acquired to identify key white matter tracts for surgical avoidance. Here, we investigate whether structural connectivity abnormalities, inferred from dMRI, may be used in conjunction with functional iEEG abnormalities to aid localisation of the epileptogenic zone (EZ), improving surgical outcomes in epilepsy.

METHODS

We retrospectively investigated data from 43 patients (42% female) with epilepsy who had surgery following iEEG. Twenty-five patients (58%) were free from disabling seizures (ILAE 1 or 2) at one year. Interictal iEEG functional, and dMRI structural connectivity abnormalities were quantified by comparison to a normative map and healthy controls. We explored whether the resection of maximal abnormalities related to improved surgical outcomes, in both modalities individually and concurrently. Additionally, we suggest how connectivity abnormalities may inform the placement of iEEG electrodes pre-surgically using a patient case study.

FINDINGS

Seizure freedom was 15 times more likely in patients with resection of maximal connectivity and iEEG abnormalities (p = 0.008). Both modalities separately distinguished patient surgical outcome groups and when used simultaneously, a decision tree correctly separated 36 of 43 (84%) patients.

INTERPRETATION

Our results suggest that both connectivity and iEEG abnormalities may localise epileptogenic tissue, and that these two modalities may provide complementary information in pre-surgical evaluations.

FUNDING

This research was funded by UKRI, CDT in Cloud Computing for Big Data, NIH, MRC, Wellcome Trust and Epilepsy Research UK.

[Analysis of neural fragility in epileptic zone based on stereoelectroencephalography]

There are some limitations in the localization of epileptogenic zone commonly used by human eyes to identify abnormal discharges of intracranial electroencephalography in epilepsy. However, at present, the accuracy of the localization of epileptogenic zone by extracting intracranial electroencephalography features needs to be further improved. As a new method using dynamic network model, neural fragility has potential application value in the localization of epileptogenic zone. In this paper, the neural fragility analysis method was used to analyze the stereoelectroencephalography signals of 35 seizures in 20 patients, and then the epileptogenic zone electrodes were classified using the random forest model, and the classification results were compared with the time-frequency characteristics of six different frequency bands extracted by short-time Fourier transform. The results showed that the area under curve (AUC) of epileptic focus electrodes based on time-frequency analysis was 0.870 (delta) to 0.956 (high gamma), and its classification accuracy increased with the increase of frequency band, while the AUC by using neural fragility could reach 0.957. After fusing the neural fragility and the time-frequency characteristics of the γ and high γ band, the AUC could be further increased to 0.969, which was improved on the original basis. This paper verifies the effectiveness of neural fragility in identifying epileptogenic zone, and provides a theoretical reference for its further clinical application.

Complementary structural and functional abnormalities to localise epileptogenic tissue

Background

When investigating suitability for epilepsy surgery, people with drug-refractory focal epilepsy may have intracranial EEG (iEEG) electrodes implanted to localise seizure onset. Diffusion-weighted magnetic resonance imaging (dMRI) may be acquired to identify key white matter tracts for surgical avoidance. Here, we investigate whether structural connectivity abnormalities, inferred from dMRI, may be used in conjunction with functional iEEG abnormalities to aid localisation of the epileptogenic zone (EZ), improving surgical outcomes in epilepsy.

Methods

We retrospectively investigated data from 43 patients with epilepsy who had surgery following iEEG. Twenty-five patients (58%) were free from disabling seizures (ILAE 1 or 2) at one year. Interictal iEEG functional, and dMRI structural connectivity abnormalities were quantified by comparison to a normative map and healthy controls. We explored whether the resection of maximal abnormalities related to improved surgical outcomes, in both modalities individually and concurrently. Additionally, we suggest how connectivity abnormalities may inform the placement of iEEG electrodes pre-surgically using a patient case study.

Findings

Seizure freedom was 15 times more likely in patients with resection of maximal connectivity and iEEG abnormalities (p=0.008). Both modalities separately distinguished patient surgical outcome groups and when used simultaneously, a decision tree correctly separated 36 of 43 (84%) patients.

Interpretation

Our results suggest that both connectivity and iEEG abnormalities may localise epileptogenic tissue, and that these two modalities may provide complementary information in pre-surgical evaluations.

Funding

This research was funded by UKRI, CDT in Cloud Computing for Big Data, NIH, MRC, Wellcome Trust and Epilepsy Research UK.

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.

Intracranial EEG Structure-Function Coupling and Seizure Outcomes After Epilepsy Surgery

BACKGROUND AND OBJECTIVES

Surgery is an effective treatment for drug-resistant epilepsy, which modifies the brain's structure and networks to regulate seizure activity. Our objective was to examine the relationship between brain structure and function to determine the extent to which this relationship affects the success of the surgery in controlling seizures. We hypothesized that a stronger association between brain structure and function would lead to improved seizure control after surgery.

METHODS

We constructed functional and structural brain networks in patients with drug-resistant focal epilepsy by using presurgery functional data from intracranial EEG (iEEG) recordings, presurgery and postsurgery structural data from T1-weighted MRI, and presurgery diffusion-weighted MRI. We quantified the relationship (coupling) between structural and functional connectivity by using the Spearman rank correlation and analyzed this structure-function coupling at 2 spatial scales: (1) global iEEG network level and (2) individual iEEG electrode contacts using virtual surgeries. We retrospectively predicted postoperative seizure freedom by incorporating the structure-function connectivity coupling metrics and routine clinical variables into a cross-validated predictive model.

RESULTS

We conducted a retrospective analysis on data from 39 patients who met our inclusion criteria. Brain areas implanted with iEEG electrodes had stronger structure-function coupling in seizure-free patients compared with those with seizure recurrence (p = 0.002, d = 0.76, area under the receiver operating characteristic curve [AUC] = 0.78 [95% CI 0.62-0.93]). Virtual surgeries on brain areas that resulted in stronger structure-function coupling of the remaining network were associated with seizure-free outcomes (p = 0.007, d = 0.96, AUC = 0.73 [95% CI 0.58-0.89]). The combination of global and local structure-function coupling measures accurately predicted seizure outcomes with a cross-validated AUC of 0.81 (95% CI 0.67-0.94). These measures were complementary to other clinical variables and, when included for prediction, resulted in a cross-validated AUC of 0.91 (95% CI 0.82-1.0), accuracy of 92%, sensitivity of 93%, and specificity of 91%.

DISCUSSION

Our study showed that the strength of structure-function connectivity coupling may play a crucial role in determining the success of epilepsy surgery. By quantitatively incorporating structure-function coupling measures and standard-of-care clinical variables into presurgical evaluations, we may be able to better localize epileptogenic tissue and select patients for epilepsy surgery.

CLASSIFICATION OF EVIDENCE

This is a Class IV retrospective case series showing that structure-function mapping may help determine the outcome from surgical resection for treatment-resistant focal epilepsy.

Optogenetic stimulation reveals a latent tipping point in cortical networks during ictogenesis

Brain-state transitions are readily apparent from changes in brain rhythms,1 but are difficult to predict, suggestive that the underlying cause is latent to passive recording methods. Among the most important transitions, clinically, are the starts of seizures. We here show that an 'active probing' approach may have several important benefits for epileptic management, including by helping predict these transitions. We used mice expressing the optogenetic actuator, channelrhodopsin, in pyramidal cells, allowing this population to be stimulated in isolation. Intermittent stimulation at frequencies as low as 0.033 Hz (period = 30 s) delayed the onset of seizure-like events in an acute brain slice model of ictogenesis, but the effect was lost if stimulation was delivered at even lower frequencies (1/min). Notably, active probing additionally provides advance indication of when seizure-like activity is imminent, revealed by monitoring the postsynaptic response to stimulation. The postsynaptic response, recorded extracellularly, showed an all-or-nothing change in both amplitude and duration, a few hundred seconds before seizure-like activity began-a sufficient length of time to provide a helpful warning of an impending seizure. The change in the postsynaptic response then persisted for the remainder of the recording, indicative of a state change from a pre-epileptic to a pro-epileptic network. This occurred in parallel with a large increase in the stimulation-triggered Ca2+ entry into pyramidal dendrites, and a step increase in the number of evoked postsynaptic action potentials, both consistent with a reduction in the threshold for dendritic action potentials. In 0 Mg2+ bathing media, the reduced threshold was not associated with changes in glutamatergic synaptic function, nor of GABAergic release from either parvalbumin or somatostatin interneurons, but simulations indicate that the step change in the optogenetic response can instead arise from incremental increases in intracellular [Cl-]. The change in the response to stimulation was replicated by artificially raising intracellular [Cl-], using the optogenetic chloride pump, halorhodopsin. By contrast, increases in extracellular [K+] cannot account for the firing patterns in the response to stimulation, although this, and other cellular changes, may contribute to ictal initiation in other circumstances. We describe how these various cellular changes form a synergistic network of positive feedback mechanisms, which may explain the precipitous nature of seizure onset. This model of seizure initiation draws together several major lines of epilepsy research as well as providing an important proof-of-principle regarding the utility of open-loop brain stimulation for clinical management of the condition.

Pre-ictal fluctuation of EEG functional connectivity discriminates seizure phenotypes in mesial temporal lobe epilepsy

OBJECTIVE

We explored whether quantifiable differences between clinical seizures (CSs) and subclinical seizures (SCSs) occur in the pre-ictal state.

METHODS

We analyzed pre-ictal stereo-electroencephalography (SEEG) retrospectively across mesial temporal lobe epilepsy patients with recorded CSs and SCSs. Power spectral density and functional connectivity (FC) were quantified within and between the seizure onset zone (SOZ) and the early propagation zone (PZ), respectively. To evaluate the fluctuation of neural connectivity, FC variability was computed. Measures were further verified by a logistic regression model to evaluate their classification potentiality through the area under the receiver-operating-characteristics curve (AUC).

RESULTS

Fifty-four pre-ictal SEEG epochs (27 CSs and 27 SCSs) were selected among 14 patients. Within the SOZ, pre-ictal FC variability of CSs was larger than SCSs in 1-45 Hz during 30 seconds before seizure onset. Pre-ictal FC variability between the SOZ and PZ was larger in SCSs than CSs in 55-80 Hz within 1 minute before onset. Using these two variables, the logistic regression model achieved an AUC of 0.79 when classifying CSs and SCSs.

CONCLUSIONS

Pre-ictal FC variability within/between epileptic zones, not signal power or FC value, distinguished SCSs from CSs.

SIGNIFICANCE

Pre-ictal epileptic network stability possibly marks seizure phenotypes, contributing insights into ictogenesis and potentially helping seizure prediction.

Autoregressive models for biomedical signal processing

Autoregressive models are ubiquitous tools for the analysis of time series in many domains such as computational neuroscience and biomedical engineering. In these domains, data is, for example, collected from measurements of brain activity. Crucially, this data is subject to measurement errors as well as uncertainties in the underlying system model. As a result, standard signal processing using autoregressive model estimators may be biased. We present a framework for autoregressive modelling that incorporates these uncertainties explicitly via an overparameterised loss function. To optimise this loss, we derive an algorithm that alternates between state and parameter estimation. Our work shows that the procedure is able to successfully denoise time series and successfully reconstruct system parameters.Clinical relevance- This new paradigm can be used in a multitude of applications in neuroscience such as brain-computer interface data analysis and better understanding of brain dynamics in diseases such as epilepsy.

Quantitative approaches to guide epilepsy surgery from intracranial EEG

Over the past 10 years, the drive to improve outcomes from epilepsy surgery has stimulated widespread interest in methods to quantitatively guide epilepsy surgery from intracranial EEG (iEEG). Many patients fail to achieve seizure freedom, in part due to the challenges in subjective iEEG interpretation. To address this clinical need, quantitative iEEG analytics have been developed using a variety of approaches, spanning studies of seizures, interictal periods, and their transitions, and encompass a range of techniques including electrographic signal analysis, dynamical systems modeling, machine learning and graph theory. Unfortunately, many methods fail to generalize to new data and are sensitive to differences in pathology and electrode placement. Here, we critically review selected literature on computational methods of identifying the epileptogenic zone from iEEG. We highlight shared methodological challenges common to many studies in this field and propose ways that they can be addressed. One fundamental common pitfall is a lack of open-source, high-quality data, which we specifically address by sharing a centralized high-quality, well-annotated, multicentre dataset consisting of >100 patients to support larger and more rigorous studies. Ultimately, we provide a road map to help these tools reach clinical trials and hope to improve the lives of future patients.

Personalised virtual brain models in epilepsy

Individuals with drug-resistant focal epilepsy are candidates for surgical treatment as a curative option. Before surgery can take place, the patient must have a presurgical evaluation to establish whether and how surgical treatment might stop their seizures without causing neurological deficits. Virtual brains are a new digital modelling technology that map the brain network of a person with epilepsy, using data derived from MRI. This technique produces a computer simulation of seizures and brain imaging signals, such as those that would be recorded with intracranial EEG. When combined with machine learning, virtual brains can be used to estimate the extent and organisation of the epileptogenic zone (ie, the brain regions related to seizure generation and the spatiotemporal dynamics during seizure onset). Virtual brains could, in the future, be used for clinical decision making, to improve precision in localisation of seizure activity, and for surgical planning, but at the moment these models have some limitations, such as low spatial resolution. As evidence accumulates in support of the predictive power of personalised virtual brain models, and as methods are tested in clinical trials, virtual brains might inform clinical practice in the near future.

Six-Center Assessment of CNN-Transformer with Belief Matching Loss for Patient-Independent Seizure Detection in EEG

Neurologists typically identify epileptic seizures from electroencephalograms (EEGs) by visual inspection. This process is often time-consuming, especially for EEG recordings that last hours or days. To expedite the process, a reliable, automated, and patient-independent seizure detector is essential. However, developing a patient-independent seizure detector is challenging as seizures exhibit diverse characteristics across patients and recording devices. In this study, we propose a patient-independent seizure detector to automatically detect seizures in both scalp EEG and intracranial EEG (iEEG). First, we deploy a convolutional neural network with transformers and belief matching loss to detect seizures in single-channel EEG segments. Next, we extract regional features from the channel-level outputs to detect seizures in multi-channel EEG segments. At last, we apply post-processing filters to the segment-level outputs to determine seizures' start and end points in multi-channel EEGs. Finally, we introduce the minimum overlap evaluation scoring as an evaluation metric that accounts for minimum overlap between the detection and seizure, improving upon existing assessment metrics. We trained the seizure detector on the Temple University Hospital Seizure (TUH-SZ) dataset and evaluated it on five independent EEG datasets. We evaluate the systems with the following metrics: sensitivity (SEN), precision (PRE), and average and median false positive rate per hour (aFPR/h and mFPR/h). Across four adult scalp EEG and iEEG datasets, we obtained SEN of 0.617-1.00, PRE of 0.534-1.00, aFPR/h of 0.425-2.002, and mFPR/h of 0-1.003. The proposed seizure detector can detect seizures in adult EEGs and takes less than 15[Formula: see text]s for a 30[Formula: see text]min EEG. Hence, this system could aid clinicians in reliably identifying seizures expeditiously, allocating more time for devising proper treatment.

Prediction of epilepsy surgery outcome using foramen ovale EEG - A machine learning approach

INTRODUCTION

Patients with drug-resistant focal epilepsy may benefit from ablative or resective surgery. In presurgical work-up, intracranial EEG markers have been shown to be useful in identification of the seizure onset zone and prediction of post-surgical seizure freedom. However, in most cases, implantation of depth or subdural electrodes is performed, exposing patients to increased risks of complications.

METHODS

We analysed EEG data recorded from a minimally invasive approach utilizing foramen ovale (FO) and epidural peg electrodes using a supervised machine learning approach to predict post-surgical seizure freedom. Power-spectral EEG features were incorporated in a logistic regression model predicting one-year post-surgical seizure freedom. The prediction model was validated using repeated 5-fold cross-validation and compared to outcome prediction based on clinical and scalp EEG variables.

RESULTS

Forty-seven patients (26 patients with post-surgical 1-year seizure freedom) were included in the study, with 31 having FO and 27 patients having peg onset seizures. The area under the receiver-operating curve for post-surgical seizure freedom (Engel 1A) prediction in patients with FO onset seizures was 0.74 ± 0.23 using electrophysiology features, compared to 0.66 ± 0.22 for predictions based on clinical and scalp EEG variables (p < 0.001). The most important features for prediction were spectral power in the gamma and high gamma ranges. EEG data from peg electrodes was not informative in predicting post-surgical outcomes.

CONCLUSION

In this hypothesis-generating study, a data-driven approach based on EEG features derived from FO electrodes recordings outperformed the predictive ability based solely on clinical and scalp EEG variables. Pending validation in future studies, this method may provide valuable post-surgical prognostic information while minimizing risks of more invasive diagnostic approaches.

Delineating epileptogenic networks using brain imaging data and personalized modeling in drug-resistant epilepsy

Precise estimates of epileptogenic zone networks (EZNs) are crucial for planning intervention strategies to treat drug-resistant focal epilepsy. Here, we present the virtual epileptic patient (VEP), a workflow that uses personalized brain models and machine learning methods to estimate EZNs and to aid surgical strategies. The structural scaffold of the patient-specific whole-brain network model is constructed from anatomical T1 and diffusion-weighted magnetic resonance imaging. Each network node is equipped with a mathematical dynamical model to simulate seizure activity. Bayesian inference methods sample and optimize key parameters of the personalized model using functional stereoelectroencephalography recordings of patients' seizures. These key parameters together with their personalized model determine a given patient's EZN. Personalized models were further used to predict the outcome of surgical intervention using virtual surgeries. We evaluated the VEP workflow retrospectively using 53 patients with drug-resistant focal epilepsy. VEPs reproduced the clinically defined EZNs with a precision of 0.6, where the physical distance between epileptogenic regions identified by VEP and the clinically defined EZNs was small. Compared with the resected brain regions of 25 patients who underwent surgery, VEP showed lower false discovery rates in seizure-free patients (mean, 0.028) than in non-seizure-free patients (mean, 0.407). VEP is now being evaluated in an ongoing clinical trial (EPINOV) with an expected 356 prospective patients with epilepsy.

Graph theoretical measures of fast ripple networks improve the accuracy of post-operative seizure outcome prediction

Fast ripples (FR) are a biomarker of epileptogenic brain, but when larger portions of FR generating regions are resected seizure freedom is not always achieved. To evaluate and improve the diagnostic accuracy of FR resection for predicting seizure freedom we compared the FR resection ratio (RR) with FR network graph theoretical measures. In 23 patients FR were semi-automatically detected and quantified in stereo EEG recordings during sleep. MRI normalization and co-registration localized contacts and relation to resection margins. The number of FR, and graph theoretical measures, which were spatial (i.e., FR rate-distance radius) or temporal correlational (i.e., FR mutual information), were compared with the resection margins and with seizure outcome We found that the FR RR did not correlate with seizure-outcome (p > 0.05). In contrast, the FR rate-distance radius resected difference and the FR MI mean characteristic path length RR did correlate with seizure-outcome (p < 0.05). Retesting of positive FR RR patients using either FR rate-distance radius resected difference or the FR MI mean characteristic path length RR reduced seizure-free misclassifications from 44 to 22% and 17%, respectively. These results indicate that graph theoretical measures of FR networks can improve the diagnostic accuracy of the resection of FR events for predicting seizure freedom.

Synergistic Positive Feedback Mechanisms Underlying Seizure Initiation

Investigations into seizure initiation, in recent years, have focused almost entirely upon alterations of interneuronal function, chloride homeostasis, and extracellular potassium levels. In contrast, little attention has been directed toward a possible role of dendritic plateau potentials in the actual ictogenic transition, despite a substantial literature dating back 40 years regarding its importance generally in epilepsy. Here, we argue that an increase in dendritic excitability, coordinated across the population of pyramidal cells, is a key stage in ictogenesis.

Cortico-cortical evoked potentials in response to varying stimulation intensity improves seizure localization

OBJECTIVE

As single pulse electrical stimulation (SPES) is increasingly utilized to help localize the seizure onset zone (SOZ), it is important to understand how stimulation intensity can affect the ability to use cortico-cortical evoked potentials (CCEPs) to delineate epileptogenic regions.

METHODS

We studied 15 drug-resistant epilepsy patients undergoing intracranial EEG monitoring and SPES with titrations of stimulation intensity. The N1 amplitude and distribution of CCEPs elicited in the SOZ and non-seizure onset zone (nSOZ) were quantified at each intensity. The separability of the SOZ and nSOZ using N1 amplitudes was compared between models using responses to titrations, responses to one maximal intensity, or both.

RESULTS

At 2 mA and above, the increase in N1 amplitude with current intensity was greater for responses within the SOZ, and SOZ response distribution was maximized by 4-6 mA. Models incorporating titrations achieved better separability of SOZ and nSOZ compared to those using one maximal intensity.

CONCLUSIONS

We demonstrated that differences in CCEP amplitude over a range of current intensities can improve discriminability of SOZ regions.

SIGNIFICANCE

This study provides insight into the underlying excitability of the SOZ and how differences in current-dependent amplitudes of CCEPs may be used to help localize epileptogenic sites.

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.

Graph-generative neural network for EEG-based epileptic seizure detection via discovery of dynamic brain functional connectivity

Dynamic complexity in brain functional connectivity has hindered the effective use of signal processing or machine learning methods to diagnose neurological disorders such as epilepsy. This paper proposed a new graph-generative neural network (GGN) model for the dynamic discovery of brain functional connectivity via deep analysis of scalp electroencephalogram (EEG) signals recorded from various regions of a patient's scalp. Brain functional connectivity graphs are generated for the extraction of spatial-temporal resolution of various onset epilepsy seizure patterns. Our supervised GGN model was substantiated by seizure detection and classification experiments. We train the GGN model using a clinically proven dataset of over 3047 epileptic seizure cases. The GGN model achieved a 91% accuracy in classifying seven types of epileptic seizure attacks, which outperformed the 65%, 74%, and 82% accuracy in using the convolutional neural network (CNN), graph neural networks (GNN), and transformer models, respectively. We present the GGN model architecture and operational steps to assist neuroscientists or brain specialists in using dynamic functional connectivity information to detect neurological disorders. Furthermore, we suggest to merge our spatial-temporal graph generator design in upgrading the conventional CNN and GNN models with dynamic convolutional kernels for accuracy enhancement.

Towards network-guided neuromodulation for epilepsy

Epilepsy is well-recognized as a disorder of brain networks. There is a growing body of research to identify critical nodes within dynamic epileptic networks with the aim to target therapies that halt the onset and propagation of seizures. In parallel, intracranial neuromodulation, including deep brain stimulation and responsive neurostimulation, are well-established and expanding as therapies to reduce seizures in adults with focal-onset epilepsy; and there is emerging evidence for their efficacy in children and generalized-onset seizure disorders. The convergence of these advancing fields is driving an era of 'network-guided neuromodulation' for epilepsy. In this review, we distil the current literature on network mechanisms underlying neurostimulation for epilepsy. We discuss the modulation of key 'propagation points' in the epileptogenic network, focusing primarily on thalamic nuclei targeted in current clinical practice. These include (i) the anterior nucleus of thalamus, now a clinically approved and targeted site for open loop stimulation, and increasingly targeted for responsive neurostimulation; and (ii) the centromedian nucleus of the thalamus, a target for both deep brain stimulation and responsive neurostimulation in generalized-onset epilepsies. We discuss briefly the networks associated with other emerging neuromodulation targets, such as the pulvinar of the thalamus, piriform cortex, septal area, subthalamic nucleus, cerebellum and others. We report synergistic findings garnered from multiple modalities of investigation that have revealed structural and functional networks associated with these propagation points - including scalp and invasive EEG, and diffusion and functional MRI. We also report on intracranial recordings from implanted devices which provide us data on the dynamic networks we are aiming to modulate. Finally, we review the continuing evolution of network-guided neuromodulation for epilepsy to accelerate progress towards two translational goals: (i) to use pre-surgical network analyses to determine patient candidacy for neurostimulation for epilepsy by providing network biomarkers that predict efficacy; and (ii) to deliver precise, personalized and effective antiepileptic stimulation to prevent and arrest seizure propagation through mapping and modulation of each patients' individual epileptogenic networks.

Computerized application for epilepsy in China: Does the era of artificial intelligence comes?

Epilepsy, one of the most common neurological diseases in China, is notorious for its spontaneous, unprovoked and recurrent seizures. The etiology of epilepsy varies among individual patients, including congenital gene mutation, traumatic injury, infections, etc. This heterogeneity partly hampered the accurate diagnosis and choice of appropriate treatments. Encouragingly, great achievements have been achieved in computational science, making it become a key player in medical fields gradually and bringing new hope for rapid and accurate diagnosis as well as targeted therapies in epilepsy. Here, we historically review the advances of computerized applications in epilepsy-especially those tremendous findings achieved in China-for different purposes including seizure prediction, localization of epileptogenic zone, post-surgical prognosis, etc. Special attentions are paid to the great progress based on artificial intelligence (AI), which is more "sensitive", "smart" and "in-depth" than human capacities. At last, we give a comprehensive discussion about the disadvantages and limitations of current computerized applications for epilepsy and propose some future directions as further stepping stones to embrace "the era of AI" in epilepsy.

Multiple Sources of Fast Traveling Waves during Human Seizures: Resolving a Controversy

During human seizures, organized waves of voltage activity rapidly sweep across the cortex. Two contradictory theories describe the source of these fast traveling waves: either a slowly advancing narrow region of multiunit activity (an ictal wavefront) or a fixed cortical location. Limited observations and different analyses prevent resolution of these incompatible theories. Here we address this disagreement by combining the methods and microelectrode array recordings (N = 11 patients, 2 females, N = 31 seizures) from previous human studies to analyze the traveling wave source. We find, inconsistent with both existing theories, a transient relationship between the ictal wavefront and traveling waves, and multiple stable directions of traveling waves in many seizures. Using a computational model that combines elements of both existing theories, we show that interactions between an ictal wavefront and fixed source reproduce the traveling wave dynamics observed in vivo We conclude that combining both existing theories can generate the diversity of ictal traveling waves.SIGNIFICANCE STATEMENT The source of voltage discharges that propagate across cortex during human seizures remains unknown. Two candidate theories exist, each proposing a different discharge source. Support for each theory consists of observations from a small number of human subject recordings, analyzed with separately developed methods. How the different, limited data and different analysis methods impact the evidence for each theory is unclear. To resolve these differences, we combine the unique, human microelectrode array recordings collected separately for each theory and analyze these combined data with a unified approach. We show that neither existing theory adequately describes the data. We then propose a new theory that unifies existing proposals and successfully reproduces the voltage discharge dynamics observed in vivo.

Critical dynamics in the spread of focal epileptic seizures: Network connectivity, neural excitability and phase transitions

Focal epileptic seizures can remain localized or, alternatively, spread across brain areas, often resulting in impairment of cognitive function and loss of consciousness. Understanding the factors that promote spread is important for developing better therapeutic approaches. Here, we show that: (1) seizure spread undergoes “critical” phase transitions in models (epileptor-networks) that capture the neural dynamics of spontaneous seizures while incorporating patient-specific brain network connectivity, axonal delays and identified epileptogenic zones (EZs). We define a collective variable for the spreading dynamics as the spread size, i.e. the number of areas or nodes in the network to which a seizure has spread. Global connectivity strength and excitability in the surrounding non-epileptic areas work as phase-transition control parameters for this collective variable. (2) Phase diagrams are predicted by stability analysis of the network dynamics. (3) In addition, the components of the Jacobian’s leading eigenvector, which tend to reflect the connectivity strength and path lengths from the EZ to surrounding areas, predict the temporal order of network-node recruitment into seizure. (4) However, stochastic fluctuations in spread size in a near-criticality region make predictability more challenging. Overall, our findings support the view that within-patient seizure-spread variability can be characterized by phase-transition dynamics under transient variations in network connectivity strength and excitability across brain areas. Furthermore, they point to the potential use and limitations of model-based prediction of seizure spread in closed-loop interventions for seizure control.

Reconstructing distant interactions of multiple paths between perceptible nodes in dark networks

Quantitative research of interdisciplinary fields, including biological and social systems, has attracted great attention in recent years. Complex networks are popular and important tools for the investigations. Explosively increasing data are created by practical networks, from which useful information about dynamic networks can be extracted. From data to network structure, i.e., network reconstruction, is a crucial task. There are many difficulties in fulfilling network reconstruction, including data shortage (existence of hidden nodes) and time delay for signal propagation between adjacent nodes. In this paper a deep network reconstruction method is proposed, which can work in the conditions that even only two nodes (say A and B) are perceptible and all other network nodes are hidden. With a well-designed stochastic driving on node A, this method can reconstruct multiple interaction paths from A to B based on measured data. The distance, effective intensity, and transmission time delay of each path can be inferred accurately.

GATSMOTE: Improving Imbalanced Node Classification on Graphs via Attention and Homophily

In recent decades, non-invasive neuroimaging techniques and graph theories have enabled a better understanding of the structural patterns of the human brain at a macroscopic level. As one of the most widely used non-invasive techniques, an electroencephalogram (EEG) may collect non-neuronal signals from “bad channels”. Automatically detecting these bad channels represents an imbalanced classification task; research on the topic is rather limited. Because the human brain can be naturally modeled as a complex graph network based on its structural and functional characteristics, we seek to extend previous imbalanced node classification techniques to the bad-channel detection task. We specifically propose a novel edge generator considering the prominent small-world organization of the human brain network. We leverage the attention mechanism to adaptively calculate the weighted edge connections between each node and its neighboring nodes. Moreover, we follow the homophily assumption in graph theory to add edges between similar nodes. Adding new edges between nodes sharing identical labels shortens the path length, thus facilitating low-cost information messaging.

Graph theoretical measures of fast ripples support the epileptic network hypothesis

The epileptic network hypothesis and epileptogenic zone hypothesis are two theories of ictogenesis. The network hypothesis posits that coordinated activity among interconnected nodes produces seizures. The epileptogenic zone hypothesis posits that distinct regions are necessary and sufficient for seizure generation. High-frequency oscillations, and particularly fast ripples, are thought to be biomarkers of the epileptogenic zone. We sought to test these theories by comparing high-frequency oscillation rates and networks in surgical responders and non-responders, with no appreciable change in seizure frequency or severity, within a retrospective cohort of 48 patients implanted with stereo-EEG electrodes. We recorded inter-ictal activity during non-rapid eye movement sleep and semi-automatically detected and quantified high-frequency oscillations. Each electrode contact was localized in normalized coordinates. We found that the accuracy of seizure onset zone electrode contact classification using high-frequency oscillation rates was not significantly different in surgical responders and non-responders, suggesting that in non-responders the epileptogenic zone partially encompassed the seizure onset zone(s) (P > 0.05). We also found that in the responders, fast ripple on oscillations exhibited a higher spectral content in the seizure onset zone compared with the non-seizure onset zone (P < 1 × 10⁻⁵). By contrast, in the non-responders, fast ripple had a lower spectral content in the seizure onset zone (P < 1 × 10⁻⁵). We constructed two different networks of fast ripple with a spectral content >350 Hz. The first was a rate–distance network that multiplied the Euclidian distance between fast ripple-generating contacts by the average rate of fast ripple in the two contacts. The radius of the rate–distance network, which excluded seizure onset zone nodes, discriminated non-responders, including patients not offered resection or responsive neurostimulation due to diffuse multifocal onsets, with an accuracy of 0.77 [95% confidence interval (CI) 0.56–0.98]. The second fast ripple network was constructed using the mutual information between the timing of the events to measure functional connectivity. For most non-responders, this network had a longer characteristic path length, lower mean local efficiency in the non-seizure onset zone, and a higher nodal strength among non-seizure onset zone nodes relative to seizure onset zone nodes. The graphical theoretical measures from the rate–distance and mutual information networks of 22 non- responsive neurostimulation treated patients was used to train a support vector machine, which when tested on 13 distinct patients classified non-responders with an accuracy of 0.92 (95% CI 0.75–1). These results indicate patients who do not respond to surgery or those not selected for resection or responsive neurostimulation can be explained by the epileptic network hypothesis that is a decentralized network consisting of widely distributed, hyperexcitable fast ripple-generating nodes.

Dynamical Network Models From EEG and MEG for Epilepsy Surgery—A Quantitative Approach

There is an urgent need for more informative quantitative techniques that non-invasively and objectively assess strategies for epilepsy surgery. Invasive intracranial electroencephalography (iEEG) remains the clinical gold standard to investigate the nature of the epileptogenic zone (EZ) before surgical resection. However, there are major limitations of iEEG, such as the limited spatial sampling and the degree of subjectivity inherent in the analysis and clinical interpretation of iEEG data. Recent advances in network analysis and dynamical network modeling provide a novel aspect toward a more objective assessment of the EZ. The advantage of such approaches is that they are data-driven and require less or no human input. Multiple studies have demonstrated success using these approaches when applied to iEEG data in characterizing the EZ and predicting surgical outcomes. However, the limitations of iEEG recordings equally apply to these studies—limited spatial sampling and the implicit assumption that iEEG electrodes, whether strip, grid, depth or stereo EEG (sEEG) arrays, are placed in the correct location. Therefore, it is of interest to clinicians and scientists to see whether the same analysis and modeling techniques can be applied to whole-brain, non-invasive neuroimaging data (from MRI-based techniques) and neurophysiological data (from MEG and scalp EEG recordings), thus removing the limitation of spatial sampling, while safely and objectively characterizing the EZ. This review aims to summarize current state of the art non-invasive methods that inform epilepsy surgery using network analysis and dynamical network models. We also present perspectives on future directions and clinical applications of these promising approaches.