Seizure pathways change on circadian and slower timescales in individual patients with focal epilepsy

The part and the whole: how single nodes contribute to large-scale phase-locking in functional EEG networks

OBJECTIVE

The application of signal analysis techniques to electroencephalographic (EEG) recordings from epilepsy patients shows that epilepsy involves not only altered neuronal synchronization but also the reorganization of functional EEG networks. This study aims to assess the large-scale phase-locking of such functional networks and how individual network nodes contribute to this collective dynamics.

METHODS

We analyze the EEG recorded before, during and after seizures from sixteen patients with pharmacoresistant focal-onset epilepsy. The data is filtered to low (4-30 Hz) and high (80-150 Hz) frequencies. We define the multivariate phase-locking measure and the univariate phase-locking contribution measure. Surrogate signals are used to estimate baseline results expected under the null hypothesis that the EEG is a correlated linear stochastic process.

RESULTS

On average, nodes from inside and outside the seizure onset zone (SOZ) increase and decrease, respectively, the large-scale phase-locking. This difference becomes most evident in a joint analysis of low and high frequencies.

CONCLUSIONS

Nodes inside and outside the SOZ play opposite roles for the large-scale phase-locking in functional EEG network in epilepsy patients.

SIGNIFICANCE

The application of the phase-locking contribution measure to EEG recordings from epilepsy patients can potentially help in localizing the SOZ.

ECG-based epileptic seizure prediction: Challenges of current data-driven models

OBJECTIVE

Up to a third of patients with epilepsy fail to achieve satisfactory seizure control. A reliable method of predicting seizures would alleviate psychological and physical impact. Dysregulation in heart rate variability (HRV) has been found to precede epileptic seizures and may serve as an extracerebral predictive biomarker. This study aims to identify the preictal HRV dynamics and unveil the factors impeding the clinical application of ECG-based seizure prediction.

METHODS

Thirty-nine adult patients (eight women; median age: 38, [IQR = 31, 56.5]) with 252 seizures were included. Each patient had more than three recorded epileptic seizures, each at least 2 hours apart. For each seizure, one hour of ECG prior to seizure onset was analyzed and 97 HRV features were extracted from overlapping three-minute windows with 10s stride. Two separate patient-specific experiments were performed using a support vector machine (SVM). Firstly, the separability of training data was examined in a non-causal trial. Secondly, the prediction was attempted in pseudo-prospective conditions. Finally, visualized HRV data, clinical metadata, and results were correlated.

RESULTS

The mean receiver operating characteristic (ROC) area under the curve (AUC) for the non-causal experiment was 0.823 (±0.12), with 208 (82.5%) seizures achieving an improvement over chance (IoC) classification score (p < 0.05, Hanley & McNeil test). In pseudo-prospective classification, the ROC-AUC was 0.569 (±0.17), and 86 (49.4%) seizures were classified with IoC. Off-sample optimized SVMs failed to improve performance. Major limiting factors identified include non-stationarity, variable preictal duration and dynamics. The latter is expressed as both inter-seizure onset zone (SOZ) and intra-SOZ variability.

SIGNIFICANCE

The pseudo-prospective preictal classification achieving IoC in approximately half of tested seizures suggests the presence of genuine preictal HRV dynamics, but the overall performance does not warrant clinical application at present. The limiting factors identified are often overlooked in non-causal study designs. While current deterministic prediction methods prove inadequate, probabilistic approaches may offer a promising alternative.

PLAIN LANGUAGE SUMMARY

Many patients with epilepsy suffer from uncontrollable seizures and would greatly benefit from a reliable seizure prediction method. Currently, no such system is available to meet this need. Previous studies suggest that changes in the electrocardiogram (ECG) precede seizures by several minutes. In our work, we evaluated whether variations in heart rate could be used to predict epileptic seizures. Our findings indicate that we are still far from achieving results suitable for clinical application and highlight several limiting factors of present seizure prediction approaches.

Diminished circadian and ultradian rhythms of human brain activity in pathological tissue in vivo

Chronobiological rhythms, such as the circadian rhythm, have long been linked to neurological disorders, but it is currently unknown how pathological processes affect the expression of biological rhythms in the brain. Here, we use the unique opportunity of long-term, continuous intracranially recorded EEG from 38 patients (totalling 6338 hours) to delineate circadian (daily) and ultradian (minute to hourly) rhythms in different brain regions. We show that functional circadian and ultradian rhythms are diminished in pathological tissue, independent of regional variations. We further demonstrate that these diminished rhythms are persistent in time, regardless of load or occurrence of pathological events. These findings provide evidence that brain pathology is functionally associated with persistently diminished chronobiological rhythms in vivo in humans, independent of regional variations or pathological events. Future work interacting with, and restoring, these modulatory chronobiological rhythms may allow for novel therapies.

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.

Clinical Correlation of Altered Molecular Signatures in Epileptic Human Hippocampus and Amygdala

Widespread alterations in the expression of various genes could contribute to the pathogenesis of epilepsy. The expression levels of various genes, including major inhibitory and excitatory receptors, ion channels, cell type-specific markers, and excitatory amino acid transporters, were assessed and compared between the human epileptic hippocampus and amygdala, and findings from autopsy controls. Moreover, the potential correlation between molecular alterations in epileptic brain tissues and the clinical characteristics of patients undergoing epilepsy surgery was evaluated. Our findings revealed significant and complex changes in the expression of several key regulatory genes in both the hippocampus and amygdala of patients with intractable epilepsy. The expression changes in various genes differed considerably between the epileptic hippocampus and amygdala. Different correlation patterns were observed between changes in gene expression and clinical characteristics, depending on whether the patients were considered as a whole or were subdivided. Altered molecular signatures in different groups of epileptic patients, defined within a given category, could be viewed as diagnostic biomarkers. Distinct patterns of molecular changes that distinguish these groups from each other appear to be associated with epilepsy-specific functional consequences.

The time-evolving epileptic brain network: concepts, definitions, accomplishments, perspectives

Epilepsy is now considered a network disease that affects the brain across multiple levels of spatial and temporal scales. The paradigm shift from an epileptic focus-a discrete cortical area from which seizures originate-to a widespread epileptic network-spanning lobes and hemispheres-considerably advanced our understanding of epilepsy and continues to influence both research and clinical treatment of this multi-faceted high-impact neurological disorder. The epileptic network, however, is not static but evolves in time which requires novel approaches for an in-depth characterization. In this review, we discuss conceptual basics of network theory and critically examine state-of-the-art recording techniques and analysis tools used to assess and characterize a time-evolving human epileptic brain network. We give an account on current shortcomings and highlight potential developments towards an improved clinical management of epilepsy.

Individualized Seizure Cluster Prediction Using Machine Learning and Chronic Ambulatory Intracranial EEG

Epilepsy patients often experience acute repetitive seizures, known as seizure clusters, which can progress to prolonged seizures or status epilepticus if left untreated. Predicting the onset of seizure clusters is crucial to enable patients to receive preventative treatments. Additionally, studying the patterns of seizure clusters can help predict the seizure type (isolated or cluster) after observing a just occurred seizure. This paper presents machine learning models that use bivariate intracranial EEG (iEEG) features to predict seizure clustering. Specifically, we utilized relative entropy (REN) as a bivariate feature to capture potential differences in brain region interactions underlying isolated and cluster seizures. We analyzed a large ambulatory iEEG dataset collected from 15 patients and spanned up to 2 years of recordings for each patient, consisting of 3341 cluster seizures (from 427 clusters) and 369 isolated seizures. The dataset's substantial number of seizures per patient enabled individualized analyses and predictions. We observed that REN was significantly different between isolated and cluster seizures in majority of the patients. Machine learning models based on REN: 1) predicted whether a seizure will occur soon after a given seizure with up to 69.5% Area under the ROC Curve (AUC), 2) predicted if a seizure is the first one in a cluster with up to 55.3% AUC, outperforming baseline techniques. Overall, our findings could be beneficial in addressing the clinical burden associated with seizure clusters, enabling patients to receive timely treatments and improving their quality of life.

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.

Computational modeling allows unsupervised classification of epileptic brain states across species

Current advances in epilepsy treatment aim to personalize and responsively adjust treatment parameters to overcome patient heterogeneity in treatment efficiency. For tailoring treatment to the individual and the current brain state, tools are required that help to identify the patient- and time-point-specific parameters of epilepsy. Computational modeling has long proven its utility in gaining mechanistic insight. Recently, the technique has been introduced as a diagnostic tool to predict individual treatment outcomes. In this article, the Wendling model, an established computational model of epilepsy dynamics, is used to automatically classify epileptic brain states in intracranial EEG from patients (n = 4) and local field potential recordings from in vitro rat data (high-potassium model of epilepsy, n = 3). Five-second signal segments are classified to four types of brain state in epilepsy (interictal, preonset, onset, ictal) by comparing a vector of signal features for each data segment to four prototypical feature vectors obtained by Wendling model simulations. The classification result is validated against expert visual assessment. Model-driven brain state classification achieved a classification performance significantly above chance level (mean sensitivity 0.99 on model data, 0.77 on rat data, 0.56 on human data in a four-way classification task). Model-driven prototypes showed similarity with data-driven prototypes, which we obtained from real data for rats and humans. Our results indicate similar electrophysiological patterns of epileptic states in the human brain and the animal model that are well-reproduced by the computational model, and captured by a key set of signal features, enabling fully automated and unsupervised brain state classification in epilepsy.

The seizure severity score: a quantitative tool for comparing seizures and their response to therapy

Objective.Epilepsy is a neurological disorder characterized by recurrent seizures which vary widely in severity, from clinically silent to prolonged convulsions. Measuring severity is crucial for guiding therapy, particularly when complete control is not possible. Seizure diaries, the current standard for guiding therapy, are insensitive to the duration of events or the propagation of seizure activity across the brain. We present a quantitative seizure severity score that incorporates electroencephalography (EEG) and clinical data and demonstrate how it can guide epilepsy therapies.Approach.We collected intracranial EEG and clinical semiology data from 54 epilepsy patients who had 256 seizures during invasive, in-hospital presurgical evaluation. We applied an absolute slope algorithm to EEG recordings to identify seizing channels. From this data, we developed a seizure severity score that combines seizure duration, spread, and semiology using non-negative matrix factorization. For validation, we assessed its correlation with independent measures of epilepsy burden: seizure types, epilepsy duration, a pharmacokinetic model of medication load, and response to epilepsy surgery. We investigated the association between the seizure severity score and preictal network features.Main results.The seizure severity score augmented clinical classification by objectively delineating seizure duration and spread from recordings in available electrodes. Lower preictal medication loads were associated with higher seizure severity scores (p= 0.018, 97.5% confidence interval = [-1.242, -0.116]) and lower pre-surgical severity was associated with better surgical outcome (p= 0.042). In 85% of patients with multiple seizure types, greater preictal change from baseline was associated with higher severity.Significance.We present a quantitative measure of seizure severity that includes EEG and clinical features, validated on gold standard in-patient recordings. We provide a framework for extending our tool's utility to ambulatory EEG devices, for linking it to seizure semiology measured by wearable sensors, and as a tool to advance data-driven epilepsy care.

Temporal stability of intracranial electroencephalographic abnormality maps for localizing epileptogenic tissue

OBJECTIVE

Identifying abnormalities on interictal intracranial electroencephalogram (iEEG), by comparing patient data to a normative map, has shown promise for the localization of epileptogenic tissue and prediction of outcome. The approach typically uses short interictal segments of approximately 1 min. However, the temporal stability of findings has not been established.

METHODS

Here, we generated a normative map of iEEG in nonpathological brain tissue from 249 patients. We computed regional band power abnormalities in a separate cohort of 39 patients for the duration of their monitoring period (.92-8.62 days of iEEG data, mean = 4.58 days per patient, >4800 hours recording). To assess the localizing value of band power abnormality, we computedD RS -a measure of how different the surgically resected and spared tissue was in terms of band power abnormalities-over time.

RESULTS

In each patient, theD RS value was relatively consistent over time. The medianD RS of the entire recording period separated seizure-free (International League Against Epilepsy [ILAE] = 1) and not-seizure-free (ILAE> 1) patients well (area under the curve [AUC] = .69). This effect was similar interictally (AUC = .69) and peri-ictally (AUC = .71).

SIGNIFICANCE

Our results suggest that band power abnormality D_RS, as a predictor of outcomes from epilepsy surgery, is a relatively robust metric over time. These findings add further support for abnormality mapping of neurophysiology data during presurgical evaluation.

Chronic intracranial EEG recordings and interictal spike rate reveal multiscale temporal modulations in seizure states

Many biological processes are modulated by rhythms on circadian and multidien timescales. In focal epilepsy, various seizure features, such as spread and duration, can change from one seizure to the next within the same patient. However, the specific timescales of this variability, as well as the specific seizure characteristics that change over time, are unclear. Here, in a cross-sectional observational study, we analysed within-patient seizure variability in 10 patients with chronic intracranial EEG recordings (185-767 days of recording time, 57-452 analysed seizures/patient). We characterized the seizure evolutions as sequences of a finite number of patient-specific functional seizure network states. We then compared seizure network state occurrence and duration to (1) time since implantation and (2) patient-specific circadian and multidien cycles in interictal spike rate. In most patients, the occurrence or duration of at least one seizure network state was associated with the time since implantation. Some patients had one or more seizure network states that were associated with phases of circadian and/or multidien spike rate cycles. A given seizure network state's occurrence and duration were usually not associated with the same timescale. Our results suggest that different time-varying factors modulate within-patient seizure evolutions over multiple timescales, with separate processes modulating a seizure network state's occurrence and duration. These findings imply that the development of time-adaptive treatments in epilepsy must account for several separate properties of epileptic seizures and similar principles likely apply to other neurological conditions.

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.

A library of quantitative markers of seizure severity

OBJECTIVE

Understanding fluctuations in seizure severity within individuals is important for determining treatment outcomes and responses to therapy, as well as assessing novel treatments for epilepsy. Current methods for grading seizure severity rely on qualitative interpretations from patients and clinicians. Quantitative measures of seizure severity would complement existing approaches to electroencephalographic (EEG) monitoring, outcome monitoring, and seizure prediction. Therefore, we developed a library of quantitative EEG markers that assess the spread and intensity of abnormal electrical activity during and after seizures.

METHODS

We analyzed intracranial EEG (iEEG) recordings of 1009 seizures from 63 patients. For each seizure, we computed 16 markers of seizure severity that capture the signal magnitude, spread, duration, and postictal suppression of seizures.

RESULTS

Quantitative EEG markers of seizure severity distinguished focal versus subclinical seizures across patients. In individual patients, 53% had a moderate to large difference (rank sumr > .3 ,p < .05 ) between focal and subclinical seizures in three or more markers. Circadian and longer term changes in severity were found for the majority of patients.

SIGNIFICANCE

We demonstrate the feasibility of using quantitative iEEG markers to measure seizure severity. Our quantitative markers distinguish between seizure types and are therefore sensitive to established qualitative differences in seizure severity. Our results also suggest that seizure severity is modulated over different timescales. We envisage that our proposed seizure severity library will be expanded and updated in collaboration with the epilepsy research community to include more measures and modalities.

Multilayer brain networks can identify the epileptogenic zone and seizure dynamics

Seizure generation, propagation, and termination occur through spatiotemporal brain networks. In this paper, we demonstrate the significance of large-scale brain interactions in high-frequency (80-200Hz) for the identification of the epileptogenic zone (EZ) and seizure evolution. To incorporate the continuity of neural dynamics, here we have modeled brain connectivity constructed from stereoelectroencephalography (SEEG) data during seizures using multilayer networks. After introducing a new measure of brain connectivity for temporal networks, named multilayer eigenvector centrality (mlEVC), we applied a consensus hierarchical clustering on the developed model to identify the EZ as a cluster of nodes with distinctive brain connectivity in the ictal period. Our algorithm could successfully predict electrodes inside the resected volume as EZ for 88% of participants, who all were seizure-free for at least 12 months after surgery. Our findings illustrated significant and unique desynchronization between EZ and the rest of the brain in the early to mid-seizure. We showed that aging and the duration of epilepsy intensify this desynchronization, which can be the outcome of abnormal neuroplasticity. Additionally, we illustrated that seizures evolve with various network topologies, confirming the existence of different epileptogenic networks in each patient. Our findings suggest not only the importance of early intervention in epilepsy but possible factors that correlate with disease severity. Moreover, by analyzing the propagation patterns of different seizures, we demonstrate the necessity of collecting sufficient data for identifying epileptogenic networks.

Epileptic-network-based prediction and control of seizures in humans

Epilepsy is now conceptualized as a network disease. The epileptic brain network comprises structurally and functionally connected cortical and subcortical brain regions - spanning lobes and hemispheres -, whose connections and dynamics evolve in time. With this concept, focal and generalized seizures as well as other related pathophysiological phenomena are thought to emerge from, spread via, and be terminated by network vertices and edges that also generate and sustain normal, physiological brain dynamics. Research over the last years has advanced concepts and techniques to identify and characterize the evolving epileptic brain network and its constituents on various spatial and temporal scales. Network-based approaches further our understanding of how seizures emerge from the evolving epileptic brain network, and they provide both novel insights into pre-seizure dynamics and important clues for success or failure of measures for network-based seizure control and prevention. In this review, we summarize the current state of knowledge and address several important challenges that would need to be addressed to move network-based prediction and control of seizures closer to clinical translation.

Unsupervised EEG preictal interval identification in patients with drug-resistant epilepsy

Typical seizure prediction models aim at discriminating interictal brain activity from pre-seizure electrographic patterns. Given the lack of a preictal clinical definition, a fixed interval is widely used to develop these models. Recent studies reporting preictal interval selection among a range of fixed intervals show inter- and intra-patient preictal interval variability, possibly reflecting the heterogeneity of the seizure generation process. Obtaining accurate labels of the preictal interval can be used to train supervised prediction models and, hence, avoid setting a fixed preictal interval for all seizures within the same patient. Unsupervised learning methods hold great promise for exploring preictal alterations on a seizure-specific scale. Multivariate and univariate linear and nonlinear features were extracted from scalp electroencephalography (EEG) signals collected from 41 patients with drug-resistant epilepsy undergoing presurgical monitoring. Nonlinear dimensionality reduction was performed for each group of features and each of the 226 seizures. We applied different clustering methods in searching for preictal clusters located until 2 h before the seizure onset. We identified preictal patterns in 90% of patients and 51% of the visually inspected seizures. The preictal clusters manifested a seizure-specific profile with varying duration (22.9 ± 21.0 min) and starting time before seizure onset (47.6 ± 27.3 min). Searching for preictal patterns on the EEG trace using unsupervised methods showed that it is possible to identify seizure-specific preictal signatures for some patients and some seizures within the same patient.

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.

Multiple mechanisms shape the relationship between pathway and duration of focal seizures

A seizure's electrographic dynamics are characterized by its spatiotemporal evolution, also termed dynamical 'pathway', and the time it takes to complete that pathway, which results in the seizure's duration. Both seizure pathways and durations have been shown to vary within the same patient. However, it is unclear whether seizures following the same pathway will have the same duration or if these features can vary independently. We compared within-subject variability in these seizure features using (i) epilepsy monitoring unit intracranial EEG (iEEG) recordings of 31 patients (mean: 6.7 days, 16.5 seizures/subject), (ii) NeuroVista chronic iEEG recordings of 10 patients (mean: 521.2 days, 252.6 seizures/subject) and (iii) chronic iEEG recordings of three dogs with focal-onset seizures (mean: 324.4 days, 62.3 seizures/subject). While the strength of the relationship between seizure pathways and durations was highly subject-specific, in most subjects, changes in seizure pathways were only weakly to moderately associated with differences in seizure durations. The relationship between seizure pathways and durations was strengthened by seizures that were 'truncated' versions, both in pathway and duration, of other seizures. However, the relationship was weakened by seizures that had a common pathway, but different durations ('elasticity'), or had similar durations, but followed different pathways ('semblance'). Even in subjects with distinct populations of short and long seizures, seizure durations were not a reliable indicator of different seizure pathways. These findings suggest that seizure pathways and durations are modulated by multiple different mechanisms. Uncovering such mechanisms may reveal novel therapeutic targets for reducing seizure duration and severity.

Fluctuations in EEG band power at subject-specific timescales over minutes to days explain changes in seizure evolutions

Epilepsy is recognised as a dynamic disease, where both seizure susceptibility and seizure characteristics themselves change over time. Specifically, we recently quantified the variable electrographic spatio-temporal seizure evolutions that exist within individual patients. This variability appears to follow subject-specific circadian, or longer, timescale modulations. It is therefore important to know whether continuously recorded interictaliEEG features can capture signatures of these modulations over different timescales. In this study, we analyse continuous intracranial electroencephalographic (iEEG) recordings from video-telemetry units and find fluctuations in iEEG band power over timescales ranging from minutes up to 12 days. As expected and in agreement with previous studies, we find that all subjects show a circadian fluctuation in their iEEG band power. We additionally detect other fluctuations of similar magnitude on subject-specific timescales. Importantly, we find that a combination of these fluctuations on different timescales can explain changes in seizure evolutions in most subjects above chance level. These results suggest that subject-specific fluctuations in iEEG band power over timescales of minutes to days may serve as markers of seizure modulating processes. We hope that future study can link these detected fluctuations to their biological driver(s). There is a critical need to better understand seizure modulating processes, as this will enable the development of novel treatment strategies that could minimise the seizure spread, duration or severity and therefore the clinical impact of seizures.

Normative brain mapping of interictal intracranial EEG to localize epileptogenic tissue

The identification of abnormal electrographic activity is important in a wide range of neurological disorders, including epilepsy for localizing epileptogenic tissue. However, this identification may be challenging during non-seizure (interictal) periods, especially if abnormalities are subtle compared to the repertoire of possible healthy brain dynamics. Here, we investigate if such interictal abnormalities become more salient by quantitatively accounting for the range of healthy brain dynamics in a location-specific manner. To this end, we constructed a normative map of brain dynamics, in terms of relative band power, from interictal intracranial recordings from 234 participants (21 598 electrode contacts). We then compared interictal recordings from 62 patients with epilepsy to the normative map to identify abnormal regions. We proposed that if the most abnormal regions were spared by surgery, then patients would be more likely to experience continued seizures postoperatively. We first confirmed that the spatial variations of band power in the normative map across brain regions were consistent with healthy variations reported in the literature. Second, when accounting for the normative variations, regions that were spared by surgery were more abnormal than those resected only in patients with persistent postoperative seizures (t = -3.6, P = 0.0003), confirming our hypothesis. Third, we found that this effect discriminated patient outcomes (area under curve 0.75 P = 0.0003). Normative mapping is a well-established practice in neuroscientific research. Our study suggests that this approach is feasible to detect interictal abnormalities in intracranial EEG, and of potential clinical value to identify pathological tissue in epilepsy. Finally, we make our normative intracranial map publicly available to facilitate future investigations in epilepsy and beyond.

Perspectives on Understanding Aberrant Brain Networks in Epilepsy

Epilepsy is a neurological disorder affecting approximately 70 million people worldwide. It is characterized by seizures that are complex aberrant dynamical events typically treated with drugs and surgery. Unfortunately, not all patients become seizure-free, and there is an opportunity for novel approaches to treat epilepsy using a network view of the brain. The traditional seizure focus theory presumed that seizures originated within a discrete cortical area with subsequent recruitment of adjacent cortices with seizure progression. However, a more recent view challenges this concept, suggesting that epilepsy is a network disease, and both focal and generalized seizures arise from aberrant activity in a distributed network. Changes in the anatomical configuration or widespread neural activities spanning lobes and hemispheres could make the brain more susceptible to seizures. In this perspective paper, we summarize the current state of knowledge, address several important challenges that could further improve our understanding of the human brain in epilepsy, and invite novel studies addressing these challenges.

Interictal EEG and ECG for SUDEP Risk Assessment: A Retrospective Multicenter Cohort Study

Objective

Sudden unexpected death in epilepsy (SUDEP) is the leading cause of epilepsy-related mortality. Although lots of effort has been made in identifying clinical risk factors for SUDEP in the literature, there are few validated methods to predict individual SUDEP risk. Prolonged postictal EEG suppression (PGES) is a potential SUDEP biomarker, but its occurrence is infrequent and requires epilepsy monitoring unit admission. We use machine learning methods to examine SUDEP risk using interictal EEG and ECG recordings from SUDEP cases and matched living epilepsy controls.

Methods

This multicenter, retrospective, cohort study examined interictal EEG and ECG recordings from 30 SUDEP cases and 58 age-matched living epilepsy patient controls. We trained machine learning models with interictal EEG and ECG features to predict the retrospective SUDEP risk for each patient. We assessed cross-validated classification accuracy and the area under the receiver operating characteristic (AUC) curve.

Results

The logistic regression (LR) classifier produced the overall best performance, outperforming the support vector machine (SVM), random forest (RF), and convolutional neural network (CNN). Among the 30 patients with SUDEP [14 females; mean age (SD), 31 (8.47) years] and 58 living epilepsy controls [26 females (43%); mean age (SD) 31 (8.5) years], the LR model achieved the median AUC of 0.77 [interquartile range (IQR), 0.73–0.80] in five-fold cross-validation using interictal alpha and low gamma power ratio of the EEG and heart rate variability (HRV) features extracted from the ECG. The LR model achieved the mean AUC of 0.79 in leave-one-center-out prediction.

Conclusions

Our results support that machine learning-driven models may quantify SUDEP risk for epilepsy patients, future refinements in our model may help predict individualized SUDEP risk and help clinicians correlate predictive scores with the clinical data. Low-cost and noninvasive interictal biomarkers of SUDEP risk may help clinicians to identify high-risk patients and initiate preventive strategies.

Long-term changes in functional connectivity improve prediction of responses to intracranial stimulation of the human brain

OBJECTIVE

Targeted electrical stimulation of the brain perturbs neural networks and modulates their rhythmic activity both at the site of stimulation and at remote brain regions. Understanding, or even predicting, this neuromodulatory effect is crucial for any therapeutic use of brain stimulation. The objective of this study was to investigate if brain network properties prior to stimulation sessions hold associative and predictive value in understanding the neuromodulatory effect of electrical stimulation in a clinical context.

APPROACH

We analysed the stimulation responses in 131 stimulation sessions across 66 patients with focal epilepsy recorded through intracranial electroencephalogram (iEEG). We considered functional and structural connectivity features as predictors of the response at every iEEG contact. Taking advantage of multiple recordings over days, we also investigated how slow changes in interictal functional connectivity (FC) ahead of the stimulation, representing the long-term variability of FC, relate to stimulation responses.

MAIN RESULTS

The long-term variability of FC exhibits strong association with the stimulation-induced increases in delta and theta band power. Furthermore, we show through cross-validation that long-term variability of FC improves prediction of responses above the performance of spatial predictors alone.

SIGNIFICANCE

This study highlights the importance of the slow dynamics of functional connectivity in the prediction of brain stimulation responses. Furthermore, these findings can enhance the patient-specific design of effective neuromodulatory protocols for therapeutic interventions.

Epileptic seizure onset predicts its duration

BACKGROUND

Epilepsy is characterized by recurrent seizures that have a variety of manifestations. The severity of, and risks for patients associated with, seizures are largely linked to the duration of seizures. Methods that determine seizure duration based on seizure onsets could be used to help mitigate the risks associated with what might be extended seizures by guiding timely interventions.

METHODS

Using long-term intracranial electroencephalography (iEEG) recordings, this article presents a method for predicting whether a seizure is going to be long or short by analyzing the seizure onset. The definition of long and short depends on each patient's seizure distribution. By analyzing 2954 seizures from 10 patients, patient-specific classifiers were built to predict seizure duration given the first few seconds from the onset.

RESULTS

The proposed methodology achieved an average area under the receiver operating characteristic curve (AUC) performance of 0.7 for the 5 of 10 patients with above chance prediction performance (p value from 0.04 to 10^-9 ).

CONCLUSIONS

Our results imply that the duration of seizures can be predicted from the onset in some patients. This could form the basis of methods for predicting status epilepticus or optimizing the amount of electrical stimulation delivered by seizure control devices.

Medial septal GABAergic neurons reduce seizure duration upon optogenetic closed-loop stimulation

Seizures can emerge from multiple or large foci in temporal lobe epilepsy, complicating focally targeted strategies such as surgical resection or the modulation of the activity of specific hippocampal neuronal populations through genetic or optogenetic techniques. Here, we evaluate a strategy in which optogenetic activation of medial septal GABAergic neurons, which provide extensive projections throughout the hippocampus, is used to control seizures. We utilized the chronic intrahippocampal kainate mouse model of temporal lobe epilepsy, which results in spontaneous seizures and as is often the case in human patients, presents with hippocampal sclerosis. Medial septal GABAergic neuron populations were immunohistochemically labelled and were not reduced in epileptic conditions. Genetic labelling with mRuby of medial septal GABAergic neuron synaptic puncta and imaging across the rostral to caudal extent of the hippocampus, also indicated an unchanged number of putative synapses in epilepsy. Furthermore, optogenetic stimulation of medial septal GABAergic neurons consistently modulated oscillations across multiple hippocampal locations in control and epileptic conditions. Finally, wireless optogenetic stimulation of medial septal GABAergic neurons, upon electrographic detection of spontaneous hippocampal seizures, resulted in reduced seizure durations. We propose medial septal GABAergic neurons as a novel target for optogenetic control of seizures in temporal lobe epilepsy.

Extracting the transition network of epileptic seizure onset

In presurgical monitoring, focal seizure onset is visually assessed from intracranial electroencephalogram (EEG), typically based on the selection of channels that show the strongest changes in amplitude and frequency. As epileptic seizure dynamics is increasingly considered to reflect changes in potentially distributed neural networks, it becomes important to also assess the interrelationships between channels. We propose a workflow to quantitatively extract the nodes and edges contributing to the seizure onset using an across-seizure scoring. We propose a quantification of the consistency of EEG channel contributions to seizure onset within a patient. The workflow is exemplified using recordings from patients with different degrees of seizure-onset consistency.

Data-driven method to infer the seizure propagation patterns in an epileptic brain from intracranial electroencephalography

Surgical interventions in epileptic patients aimed at the removal of the epileptogenic zone have success rates at only 60-70%. This failure can be partly attributed to the insufficient spatial sampling by the implanted intracranial electrodes during the clinical evaluation, leading to an incomplete picture of spatio-temporal seizure organization in the regions that are not directly observed. Utilizing the partial observations of the seizure spreading through the brain network, complemented by the assumption that the epileptic seizures spread along the structural connections, we infer if and when are the unobserved regions recruited in the seizure. To this end we introduce a data-driven model of seizure recruitment and propagation across a weighted network, which we invert using the Bayesian inference framework. Using a leave-one-out cross-validation scheme on a cohort of 45 patients we demonstrate that the method can improve the predictions of the states of the unobserved regions compared to an empirical estimate that does not use the structural information, yet it is on the same level as the estimate that takes the structure into account. Furthermore, a comparison with the performed surgical resection and the surgery outcome indicates a link between the inferred excitable regions and the actual epileptogenic zone. The results emphasize the importance of the structural connectome in the large-scale spatio-temporal organization of epileptic seizures and introduce a novel way to integrate the patient-specific connectome and intracranial seizure recordings in a whole-brain computational model of seizure spread.

Focal to bilateral tonic-clonic seizures are associated with widespread network abnormality in temporal lobe epilepsy

OBJECTIVE

Our objective was to identify whether the whole-brain structural network alterations in patients with temporal lobe epilepsy (TLE) and focal to bilateral tonic-clonic seizures (FBTCS) differ from alterations in patients without FBTCS.

METHODS

We dichotomized a cohort of 83 drug-resistant patients with TLE into those with and without FBTCS and compared each group to 29 healthy controls. For each subject, we used diffusion-weighted magnetic resonance imaging to construct whole-brain structural networks. First, we measured the extent of alterations by performing FBTCS-negative (FBTCS-) versus control and FBTCS-positive (FBTCS+) versus control comparisons, thereby delineating altered subnetworks of the whole-brain structural network. Second, by standardizing each patient's networks using control networks, we measured the subject-specific abnormality at every brain region in the network, thereby quantifying the spatial localization and the amount of abnormality in every patient.

RESULTS

Both FBTCS+ and FBTCS- patient groups had altered subnetworks with reduced fractional anisotropy and increased mean diffusivity compared to controls. The altered subnetwork in FBTCS+ patients was more widespread than in FBTCS- patients (441 connections altered at t > 3, p < .001 in FBTCS+ compared to 21 connections altered at t > 3, p = .01 in FBTCS-). Significantly greater abnormalities-aggregated over the entire brain network as well as assessed at the resolution of individual brain areas-were present in FBTCS+ patients (p < .001, d = .82, 95% confidence interval = .32-1.3). In contrast, the fewer abnormalities present in FBTCS- patients were mainly localized to the temporal and frontal areas.

SIGNIFICANCE

The whole-brain structural network is altered to a greater and more widespread extent in patients with TLE and FBTCS. We suggest that these abnormal networks may serve as an underlying structural basis or consequence of the greater seizure spread observed in FBTCS.

Variability of seizure-like activity in an in vitro model of epilepsy depends on the electrical recording method

Background

Hippocampal and cortical slice-based models are widely used to study seizures and epilepsy. Seizure detection and quantification are essential components for studying mechanisms of epilepsy and assessing therapeutic interventions. To obtain meaningful signals and maximize experimental throughput, variability should be minimized. Some electrical recording methods require insertion of an electrode into neuronal tissue, change in slice chemical microenvironment, and transients in temperature and pH. These perturbations can cause acute and long-term alterations of the neuronal network which may be reflected in the variability of the recorded signal.

New method

In this study we investigated the effect of experimental perturbations in three local field potential (LFP) recording methods including substrate micro-wires (s-MWs), multiple electrode arrays (MEAs), and inserted micro wire electrodes (i-MW). These methods enabled us to isolate effects of different perturbations. We used organotypic hippocampal slices (OHCs) as an in-vitro model of posttraumatic epilepsy. To investigate the effect of the disturbances caused by the recording method on the paroxysmal events, we introduced jitter analysis, which is sensitive to small differences in the seizure spike timing.

Results

Medium replacement can introduce long-lasting perturbations. Electrode insertion increased variability on a shorter time scale. OHCs also underwent spontaneous state transitions characterized by transient increases in variability.

Comparison with existing methods

This new method of seizure waveform analysis allows for more sensitive assessment of variability of ictal events than simply measuring seizure frequency and duration.

Conclusion

We demonstrated that some of the variability in OHC recordings are due to experimental perturbations while some are spontaneous and independent of recording method.

Circadian/multidien Molecular Oscillations and Rhythmicity of Epilepsy (MORE)

The occurrence of seizures at specific times of the day has been consistently observed for centuries in individuals with epilepsy. Electrophysiological recordings provide evidence that seizures have a higher probability of occurring at a given time during the night and day cycle in individuals with epilepsy here referred to as the seizure rush hour. Which mechanisms underlie such circadian rhythmicity of seizures? Why don't they occur every day at the same time? Which mechanisms may underlie their occurrence outside the rush hour? In this commentary, I present a hypothesis: MORE - Molecular Oscillations and Rhythmicity of Epilepsy, a conceptual framework to study and understand the mechanisms underlying the circadian rhythmicity of seizures and their probabilistic nature. The core of the hypothesis is the existence of ~24-hour oscillations of gene and protein expression throughout the body in different cells and organs. The orchestrated molecular oscillations control the rhythmicity of numerous body events, such as feeding and sleep. The concept developed here is that molecular oscillations may favor seizure genesis at preferred times, generating the condition for a seizure rush hour. However, the condition is not sufficient, as other factors are necessary for a seizure to occur. Studying these molecular oscillations may help us understand seizure genesis mechanisms and find new therapeutic targets and predictive biomarkers. The MORE hypothesis can be generalized to comorbidities and the slower multidien (week/month period) rhythmicity of seizures, a phenomenon addressed in another article in this issue of Epilepsia.

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.

Seizure pathways change on circadian and slower timescales in individual patients with focal epilepsy

Personalized medicine requires that treatments adapt to not only the patient but also changing factors within each individual. Although epilepsy is a dynamic disorder characterized by pathological fluctuations in brain state, surprisingly little is known about whether and how seizures vary in the same patient. We quantitatively compared within-patient seizure network evolutions using intracranial electroencephalographic (iEEG) recordings of over 500 seizures from 31 patients with focal epilepsy (mean 16.5 seizures per patient). In all patients, we found variability in seizure paths through the space of possible network dynamics. Seizures with similar pathways tended to occur closer together in time, and a simple model suggested that seizure pathways change on circadian and/or slower timescales in the majority of patients. These temporal relationships occurred independent of whether the patient underwent antiepileptic medication reduction. Our results suggest that various modulatory processes, operating at different timescales, shape within-patient seizure evolutions, leading to variable seizure pathways that may require tailored treatment approaches.

The sensitivity of network statistics to incomplete electrode sampling on intracranial EEG

Network neuroscience applied to epilepsy holds promise to map pathological networks, localize seizure generators, and inform targeted interventions to control seizures. However, incomplete sampling of the epileptic brain because of sparse placement of intracranial electrodes may affect model results. In this study, we evaluate the sensitivity of several published network measures to incomplete spatial sampling and propose an algorithm using network subsampling to determine confidence in model results. We retrospectively evaluated intracranial EEG data from 28 patients implanted with grid, strip, and depth electrodes during evaluation for epilepsy surgery. We recalculated global and local network metrics after randomly and systematically removing subsets of intracranial EEG electrode contacts. We found that sensitivity to incomplete sampling varied significantly across network metrics. This sensitivity was largely independent of whether seizure onset zone contacts were targeted or spared from removal. We present an algorithm using random subsampling to compute patient-specific confidence intervals for network localizations. Our findings highlight the difference in robustness between commonly used network metrics and provide tools to assess confidence in intracranial network localization. We present these techniques as an important step toward translating personalized network models of seizures into rigorous, quantitative approaches to invasive therapy.