Linear time-varying model characterizes invasive EEG signals generated from complex epileptic networks

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.

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.

Sources and Sinks in Interictal iEEG Networks: An iEEG Marker of the Epileptogenic Zone

Around 30% of epilepsy patients have seizures that cannot be controlled with medication. The most effective treatments for medically resistant epilepsy are interventions that surgically remove the epileptogenic zone (EZ), the regions of the brain that initiate seizure activity. A precise identification of the EZ is essential for surgical success but unfortunately, current success rates range from 20-80%. Localization of the EZ requires visual inspection of intracranial EEG (iEEG) recordings during seizure events. The need for seizure occurrence makes the process both costly and time-consuming and in the end, less than 1% of the data captured is used to assist in EZ localization. In this study, we aim to leverage interictal (between seizures) data to localize the EZ. We develop and test the source-sink index as an interictal iEEG marker by identifying two groups of network nodes from a patient's interictal iEEG network: those that inhibit a set of their neighboring nodes ("sources") and the inhibited nodes themselves ("sinks"). Specifically, we i) estimate patient-specific dynamical network models from interictal iEEG data and ii) compute a source-sink index for every network node (iEEG channel) to identify pathological nodes that correspond to the EZ. Our results suggest that in patients with successful surgical outcomes, the source-sink index clearly separates the clinically identified EZ (CA-EZ) channels from other channels whereas in patients with failed outcomes CA-EZ channels cannot be distinguished from the rest of the network.

Neural fragility as an EEG marker of the seizure onset zone

Over 15 million patients with epilepsy worldwide do not respond to drugs. Successful surgical treatment requires complete removal or disconnection of the seizure onset zone (SOZ), brain region(s) where seizures originate. Unfortunately, surgical success rates vary between 30 and 70% because no clinically validated biological marker of the SOZ exists. We develop and retrospectively validate a new electroencephalogram (EEG) marker-neural fragility-in a retrospective analysis of 91 patients by using neural fragility of the annotated SOZ as a metric to predict surgical outcomes. Fragility predicts 43 out of 47 surgical failures, with an overall prediction accuracy of 76% compared with the accuracy of clinicians at 48% (successful outcomes). In failed outcomes, we identify fragile regions that were untreated. When compared to 20 EEG features proposed as SOZ markers, fragility outperformed in predictive power and interpretability, which suggests neural fragility as an EEG biomarker of the SOZ.

Network Fragility for Seizure Genesis in an Acute in vivo Model of Epilepsy

Epilepsy affects over 50 million people worldwide and 30% of patients' seizures are medically refractory. The process of localizing and removing the epileptogenic zone is error-prone and ill-posed in part because we do not understand how epilepsy manifests. It has recently been proposed that the epileptic cortex is fragile in the sense that seizures manifest through small perturbations in the synaptic connections that render the entire cortical network unstable. If the fragility of the cortical network could be computed over a period in which seizure genesis occurs, then it might elucidate network mechanisms correlated to the epileptogenic zone. In this study, we used local field potentials (LFP) from neocortex by implementing an acute model of epilepsy in mice. These recordings were used to develop a dynamical network model that quantifies the fragility of the nodes from LFP epochs of baseline activity, preictal and ictal states. Fragility was quantified by the generation of a linear time-varying model to which we then applied a perturbation to determine the sensitivity of nodes in the network. Spatiotemporal fragility maps showed clear quantifiable changes in the epileptogenic network's properties throughout different states of seizure genesis. We quantified this difference over a baseline, preictal and ictal periods to show that network fragility is modulated in the manifestation of epilepsy.

State-space models of evoked potentials to localize the seizure onset zone

Surgical removal of the seizure onset zone (SOZ) in epilepsy patients is a potentially curative treatment, but the process heavily relies on accurate localization of the SOZ via visual inspection. SPES (Single-pulse electrical stimulation) is a method recently used to explore inter-areal connectivity in vivo to probe functional brain networks such as language and motor networks, and to a much lesser degree, seizure networks. We hypothesized that a dynamical quantification of the connectivity networks derived from the evoked responses induced by SPES could also be used to localize the SOZ. To test our hypothesis, we used an intracranial EEG (iEEG) data set in which five epilepsy patients underwent extensive SPES evaluation. For each patient, and for each dataset that stimulated a different pair of electrodes, we constructed a state-space model from the patient's data. Specifically, we simultaneously estimated model parameters under an exogenous pulse input to a dynamical system whose state vector consisted of the response iEEG signals. Then, the size of the reachable state space, as quantified by the maximum singular value of the reachability matrix, σ_max(R), was computed and denoted as the "largest" network response possible when stimulating the given pair. Our results suggest high agreement between σ_max(R) and clinically annotated SOZ for patients with localizable SOZs.Clinical Relevance- Our study applies dynamical systems theory to identify epileptogenic brain regions, creating a novel tool that clinicians may use in surgical planning for medically-refractory epilepsy patients.

Estimating Intracranial EEG Signals at Missing Electrodes in Epileptic Networks

Epilepsy can be controlled by targeted treatment of the epileptogenic zone (EZ), the region in the brain where seizures originate. Identification of the EZ often requires visual inspection of invasive EEG recordings and thus relies heavily on placement of electrodes, such that they cover the EZ. A dense brain coverage would be ideal to obtain accurate boundaries of the EZ but is not possible due to surgical limitations. This gives rise to the "missing electrode problem", where clinicians desire to know what neural activity looks like between implanted electrodes. In this paper, we compare two methods for time series estimation of missing stereotactic EEG (SEEG) recordings. Specifically, we represent SEEG data as a sequence of Linear Time-Invariant (LTI) models. We then remove one signal from the data set and apply two different algorithms to simultaneously estimate the LTI models and the "missing" signal: (i) a Reduced-Order Observer in combination with Least Squares Estimation and (ii) an Expectation Maximization based Kalman Filter. The performance of each approach is evaluated in terms of (i) estimation error, (ii) sensitivity to initial conditions, and (iii) algorithm run-time. We found that the EM approach has smaller estimation errors and is less sensitive to initial conditions. However, the reduced-order observer has a run-time that is orders of magnitude faster.

Virtual Cortical Stimulation Mapping of Epilepsy Networks to Localize the Epileptogenic Zone

Cortical stimulation mapping (CSM) is a common clinical procedure for mapping eloquent cortex in epilepsy patients. Electrical responses to the stimulation, or after-discharges (ADs), that occur in response to stimulation can point to unstable regions of cortex that are more prone to spontaneous seizures. Clinicians are interested in identifying regions that start seizures, i.e., the epileptogenic zone (EZ), so that they can target treatment. However, during CSM, not all regions are stimulated, as it would be time-consuming and potentially harmful to the patient. This limits the clinician's ability to fully explore ADs to reliably localize the EZ. In this paper, we develop a virtual CSM procedure that processes pre-seizure intracranial EEG recordings obtained from epilepsy patients being treated at three different epilepsy centers. First, we identify a linear time varying network (LTVN) model from electrocorticography (ECoG) and stereo-EEG (SEEG) data using sparse least squares estimation for each patient. We then construct an virtual CSM by applying impulse perturbations to each electrode contact in the LTVN model and then measuring the ADs of the network. We summarize the l2-norm of the responses in the form of a heatmap that shows the spatio-temporal evolution of the ADs before, during, and after seizures. Finally we compute an impulse response ratio (IRR) metric from each heatmap, that measures the ratio between the mean norm of ADs of clinically annotated EZ contacts and the mean norm of ADs of the remaining contacts. We find that the IRR is higher in maps derived from patients with successful surgical outcomes and lower in failed surgical outcomes. This suggests that virtual CSM may provide valuable information to clinicians regarding EZ location.

Emerging roles of network analysis for epilepsy

In recent years there has been increasing interest in applying network science tools to EEG data. At the 2018 American Epilepsy Society conference in New Orleans, LA, the yearly session of the Engineering and Neurostimulation Special Interest Group focused on emerging, translational technologies to analyze seizure networks. Each speaker demonstrated practical examples of how network tools can be utilized in clinical care and provide additional data to help care for patients with intractable epilepsy. The groups presented advances using tools from functional connectivity, control theory, and graph theory to analyze human EEG data. These tools have great potential to augment clinical interpretation of EEG signals.