Virtual Cortical Resection Reveals Push-Pull Network Control Preceding Seizure Evolution

Spontaneous HFO Sequences Reveal Propagation Pathways for Precise Delineation of Epileptogenic Networks

Epilepsy, a neurological disorder affecting millions worldwide, poses great challenges in precisely delineating the epileptogenic zone - the brain region generating seizures - for effective treatment. High-frequency oscillations (HFOs) are emerging as promising biomarkers; however, the clinical utility is hindered by the difficulties in distinguishing pathological HFOs from non- epileptiform activities at single electrode and single patient resolution and understanding their dynamic role in epileptic networks. Here, we introduce an HFO-sequencing approach to analyze spontaneous HFOs traversing cortical regions in 40 drug-resistant epilepsy patients. This data- driven method automatically detected over 8.9 million HFOs, pinpointing pathological HFO- networks, and unveiled intricate millisecond-scale spatiotemporal dynamics, stability, and functional connectivity of HFOs in prolonged intracranial EEG recordings. These HFO sequences demonstrated a significant improvement in localization of epileptic tissue, with an 818.47% increase in concordance with seizure-onset zone (mean error: 2.92 mm), compared to conventional benchmarks. They also accurately predicted seizure outcomes for 90% AUC based on pre-surgical information using generalized linear models. Importantly, this mapping remained reliable even with short recordings (mean standard deviation: 3.23 mm for 30-minute segments). Furthermore, HFO sequences exhibited distinct yet highly repetitive spatiotemporal patterns, characterized by pronounced synchrony and predominant inward information flow from periphery towards areas involved in propagation, suggesting a crucial role for excitation-inhibition balance in HFO initiation and progression. Together, these findings shed light on the intricate organization of epileptic network and highlight the potential of HFO-sequencing as a translational tool for improved diagnosis, surgical targeting, and ultimately, better outcomes for vulnerable patients with drug-resistant epilepsy.

One Sentence Summary

Pathological fast brain oscillations travel like traffic along varied routes, outlining recurrently visited neural sites emerging as critical hotspots in epilepsy network.

Measurable transitions during seizures in intracranial EEG: A stereoelectroencephalography and SPECT study

OBJECTIVE

Ictal Single Photon Emission Computed Tomography (SPECT) and stereo-electroencephalography (SEEG) are diagnostic techniques used for the management of patients with drug-resistant focal epilepsies. While hyperperfusion patterns in ictal SPECT studies reveal seizure onset and propagation pathways, the role of ictal hypoperfusion remains poorly understood. The goal of this study was to systematically characterize the spatio-temporal information flow dynamics between differently perfused brain regions using stereo-EEG recordings.

METHODS

We identified seizure-free patients after resective epilepsy surgery who had prior ictal SPECT and SEEG investigations. We estimated directional connectivity between the epileptogenic-zone (EZ), non-resected areas of hyperperfusion, hypoperfusion, and baseline perfusion during the interictal, preictal, ictal, and postictal periods.

RESULTS

Compared to the background, we noted significant information flow (1) during the preictal period from the EZ to the baseline and hyperperfused regions, (2) during the ictal onset from the EZ to all three regions, and (3) during the period of seizure evolution from the area of hypoperfusion to all three regions.

CONCLUSIONS

Hypoperfused brain regions were found to indirectly interact with the EZ during the ictal period.

SIGNIFICANCE

Our unique study, combining intracranial electrophysiology and perfusion imaging, presents compelling evidence of dynamic changes in directional connectivity between brain regions during the transition from interictal to ictal states.

Epileptic tissue localization using graph-based networks in the high frequency oscillation range of intracranial electroencephalography

PURPOSE

High frequency oscillations (HFOs) are an emerging biomarker of epilepsy. However, very few studies have investigated the functional connectivity of interictal iEEG signals in the frequency range of HFOs. Here, we study the corresponding functional networks using graph theory, and we assess their predictive value for automatic electrode classification in a cohort of 20 drug resistant patients.

METHODS

Coherence-based connectivity analysis was performed on the iEEG recordings, and six different local graph measures were computed in both sub-bands of the HFO frequency range (80-250 Hz and 250-500 Hz). Correlation analysis was implemented between the local graph measures and the ripple and fast ripple rates. Finally, the WEKA software was employed for training and testing different predictive models on the aforementioned local graph measures.

RESULTS

The ripple rate was significantly correlated with five out of six local graph measures in the functional network. For fast ripples, their rate was also significantly (but negatively) correlated with most of the local metrics. The results from WEKA showed that the Logistic Regression algorithm was able to classify highly HFO-contaminated electrodes with an accuracy of 82.5 % for ripples and 75.4 % for fast ripples.

CONCLUSION

Functional connectivity networks in the HFO band could represent an alternative to the direct use of distinct HFO events, while also providing important insights about hub epileptic areas that can represent possible surgical targets. Automatic electrode classification through FC-based classifiers can help bypass the burden of manual HFO annotation, providing at the same time similar amount of information about the epileptic tissue.

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.

Brain functional connectivity-based prediction of vagus nerve stimulation efficacy in pediatric pharmacoresistant epilepsy

OBJECTIVE

Although vagus nerve stimulation (VNS) is a common and widely used therapy for pharmacoresistant epilepsy, the reported efficacy of VNS in pediatric patients varies, so it is unclear which children will respond to VNS therapy. This study aimed to identify functional brain network features associated with VNS action to distinguish VNS responders from nonresponders using scalp electroencephalogram (EEG) data.

METHODS

Twenty-three children were included in this study, 16 in the discovery cohort and 7 in the test cohort. Using partial correlation value as a measure of whole-brain functional connectivity, we identified the differential edges between responders and nonresponders. Results derived from this were used as input to generate a support vector machine-learning classifier to predict VNS outcomes.

RESULTS

The postcentral gyrus in the left and right parietal lobe regions was identified as the most significant differential brain region between VNS responders and nonresponders (p < 0.001). The resultant classifier demonstrated a mean AUC value of 0.88, a mean sensitivity rate of 91.4%, and a mean specificity rate of 84.3% on fivefold cross-validation in the discovery cohort. In the testing cohort, our study demonstrated an AUC value of 0.91, a sensitivity rate of 86.6%, and a specificity rate of 79.3%. Furthermore, for prediction accuracy, our model can achieve 81.4% accuracy at the epoch level and 100% accuracy at the patient level.

SIGNIFICANCE

This study provides the first treatment response prediction model for VNS using scalp EEG data with ictal recordings and offers new insights into its mechanism of action. Our results suggest that brain functional connectivity features can help predict therapeutic response to VNS therapy. With further validation, our model could facilitate the selection of targeted pediatric patients and help avoid risky and costly procedures for patients who are unlikely to benefit from VNS therapy.

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.

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.

Functional connectivity discriminates epileptogenic states and predicts surgical outcome in children with drug resistant epilepsy

Normal brain functioning emerges from a complex interplay among regions forming networks. In epilepsy, these networks are disrupted causing seizures. Highly connected nodes in these networks are epilepsy surgery targets. Here, we assess whether functional connectivity (FC) using intracranial electroencephalography can quantify brain regions epileptogenicity and predict surgical outcome in children with drug resistant epilepsy (DRE). We computed FC between electrodes on different states (i.e. interictal without spikes, interictal with spikes, pre-ictal, ictal, and post-ictal) and frequency bands. We then estimated the electrodes' nodal strength. We compared nodal strength between states, inside and outside resection for good- (n = 22, Engel I) and poor-outcome (n = 9, Engel II-IV) patients, respectively, and tested their utility to predict the epileptogenic zone and outcome. We observed a hierarchical epileptogenic organization among states for nodal strength: lower FC during interictal and pre-ictal states followed by higher FC during ictal and post-ictal states (p < 0.05). We further observed higher FC inside resection (p < 0.05) for good-outcome patients on different states and bands, and no differences for poor-outcome patients. Resection of nodes with high FC was predictive of outcome (positive and negative predictive values: 47-100%). Our findings suggest that FC can discriminate epileptogenic states and predict outcome in patients with DRE.

Evolution of brain network dynamics in early Parkinson's disease with mild cognitive impairment

How mild cognitive impairment (MCI) is instantiated in dynamically interacting and spatially distributed functional brain networks remains an unexplored mystery in early Parkinson's disease (PD). We applied a machine-learning technology based on personalized sliding-window algorithm to track continuously time-varying and overlapping subnetworks under the functional brain networks calculated form resting state electroencephalogram data within a sample of 33 early PD patients (13 early PD patients with MCI and 20 early PD patients without MCI). We decoded a set of subnetworks that captured surprisingly dynamically varying and integrated interactions among certain brain lobes. We observed that the master expressed subnetworks were particularly transient, and flexibly switching between high and low expression during integration into a dynamic brain network. This transience was particularly salient in a subnetwork predominantly linking temporal-parietal-occipital lobes, which decreases in both expression and flexibility in early PD patients with MCI and expresses their degree of cognitive impairment. Moreover, MCI induced a regularly interrupted, slow evolution of subnetworks in functional brain network dynamics in early PD at the individual level, and the dynamic expression characteristics of subnetworks also reflected the degree of cognitive impairment in patients with early PD. Collectively, these results provide novel and deeper insights regarding MCI-induced abnormal dynamical interaction and large-scale changes in functional brain network of early PD.

Remote effects of temporal lobe epilepsy surgery: Long-term morphological changes after surgical resection

OBJECTIVE

Epilepsy surgery is an effective treatment for drug-resistant patients. However, how different surgical approaches affect long-term brain structure remains poorly characterized. Here, we present a semiautomated method for quantifying structural changes after epilepsy surgery and compare the remote structural effects of two approaches, anterior temporal lobectomy (ATL), and selective amygdalohippocampectomy (SAH).

METHODS

We studied 36 temporal lobe epilepsy patients who underwent resective surgery (ATL = 22, SAH = 14). All patients received same-scanner MR imaging preoperatively and postoperatively (mean 2 years). To analyze postoperative structural changes, we segmented the resection zone and modified the Advanced Normalization Tools (ANTs) longitudinal cortical pipeline to account for resections. We compared global and regional annualized cortical thinning between surgical treatments.

RESULTS

Across procedures, there was significant cortical thinning in the ipsilateral insula, fusiform, pericalcarine, and several temporal lobe regions outside the resection zone as well as the contralateral hippocampus. Additionally, increased postoperative cortical thickness was seen in the supramarginal gyrus. Patients treated with ATL exhibited greater annualized cortical thinning compared with SAH cases (ATL: -0.08 ± 0.11 mm per year, SAH: -0.01 ± 0.02 mm per year, t = 2.99, P = 0.006). There were focal postoperative differences between the two treatment groups in the ipsilateral insula (P = 0.039, corrected). Annualized cortical thinning rates correlated with preoperative cortical thickness (r = 0.60, P < 0.001) and had weaker associations with age at surgery (r = -0.33, P = 0.051) and disease duration (r = -0.42, P = 0.058).

SIGNIFICANCE

Our evidence suggests that selective procedures are associated with less cortical thinning and that earlier surgical intervention may reduce long-term impacts on brain structure.

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.

An individual data-driven virtual resection model based on epileptic network dynamics in children with intractable epilepsy: a magnetoencephalography interictal activity application

Epilepsy surgery continues to be a recommended treatment for intractable (medication-resistant) epilepsy; however, 30-70% of epilepsy surgery patients can continue to have seizures. Surgical failures are often associated with incomplete resection or inaccurate localization of the epileptogenic zone. This retrospective study aims to improve surgical outcome through in silico testing of surgical hypotheses through a personalized computational neurosurgery model created from individualized patient's magnetoencephalography recording and MRI. The framework assesses the extent of the epileptic network and evaluates underlying spike dynamics, resulting in identification of one single brain volume as a candidate for resection. Dynamic-locked networks were utilized for virtual cortical resection. This in silico protocol was tested in a cohort of 24 paediatric patients with focal drug-resistant epilepsy who underwent epilepsy surgery. Of 24 patients who were included in the analysis, 79% (19 of 24) of the models agreed with the patient's clinical surgery outcome and 21% (5 of 24) were considered as model failures (accuracy 0.79, sensitivity 0.77, specificity 0.82). Patients with unsuccessful surgery outcome typically showed a model cluster outside of the resected cavity, while those with successful surgery showed the cluster model within the cavity. Two of the model failures showed the cluster in the vicinity of the resected tissue and either a functional disconnection or lack of precision of the magnetoencephalography-MRI overlapping could explain the results. Two other cases were seizure free for 1 year but developed late recurrence. This is the first study that provides in silico personalized protocol for epilepsy surgery planning using magnetoencephalography spike network analysis. This model could provide complementary information to the traditional pre-surgical assessment methods and increase the proportion of patients achieving seizure-free outcome from surgery.

Atypical connectome topography and signal flow in temporal lobe epilepsy

Temporal lobe epilepsy (TLE) is one of the most common pharmaco-resistant epilepsies in adults. While hippocampal pathology is the hallmark of this condition, emerging evidence indicates that brain alterations extend beyond the mesiotemporal epicenter and affect macroscale brain function and cognition. We studied macroscale functional reorganization in TLE, explored structural substrates, and examined cognitive associations. We investigated a multisite cohort of 95 patients with pharmaco-resistant TLE and 95 healthy controls using state-of-the-art multimodal 3T magnetic resonance imaging (MRI). We quantified macroscale functional topographic organization using connectome dimensionality reduction techniques and estimated directional functional flow using generative models of effective connectivity. We observed atypical functional topographies in patients with TLE relative to controls, manifesting as reduced functional differentiation between sensory/motor networks and transmodal systems such as the default mode network, with peak alterations in bilateral temporal and ventromedial prefrontal cortices. TLE-related topographic changes were consistent in all three included sites and reflected reductions in hierarchical flow patterns between cortical systems. Integration of parallel multimodal MRI data indicated that these findings were independent of TLE-related cortical grey matter atrophy, but mediated by microstructural alterations in the superficial white matter immediately beneath the cortex. The magnitude of functional perturbations was robustly associated with behavioral markers of memory function. Overall, this work provides converging evidence for macroscale functional imbalances, contributing microstructural alterations, and their associations with cognitive dysfunction in TLE.

Pre-surgical features of intrinsic brain networks predict single and joint epilepsy surgery outcomes

Despite the effectiveness of surgical interventions for the treatment of intractable focal temporal lobe epilepsy (TLE), the substrates that support good outcomes are poorly understood. While algorithms have been developed for the prediction of either seizure or cognitive/psychiatric outcomes alone, no study has reported on the functional and structural architecture that supports joint outcomes. We measured key aspects of pre-surgical whole brain functional/structural network architecture and evaluated their ability to predict post-operative seizure control in combination with cognitive/psychiatric outcomes. Pre-surgically, we identified the intrinsic connectivity networks (ICNs) unique to each person through independent component analysis (ICA), and computed: (1) the spatial-temporal match between each person's ICA components and established, canonical ICNs, (2) the connectivity strength within each identified person-specific ICN, (3) the gray matter (GM) volume underlying the person-specific ICNs, and (4) the amount of variance not explained by the canonical ICNs for each person. Post-surgical seizure control and reliable change indices of change (for language [naming, phonemic fluency], verbal episodic memory, and depression) served as binary outcome responses in random forest (RF) models. The above functional and structural measures served as input predictors. Our empirically derived ICN-based measures customized to the individual showed that good joint seizure and cognitive/psychiatric outcomes depended upon higher levels of brain reserve (GM volume) in specific networks. In contrast, singular outcomes relied on systematic, idiosyncratic variance in the case of seizure control, and the weakened pre-surgical presence of functional ICNs that encompassed the ictal temporal lobe in the case of cognitive/psychiatric outcomes. Our data made clear that the ICNs differed in their propensity to provide reserve for adaptive outcomes, with some providing structural (brain), and others functional (cognitive) reserve. Our customized methodology demonstrated that when substantial unique, patient-specific ICNs are present prior to surgery there is a reliable association with poor post-surgical seizure control. These ICNs are idiosyncratic in that they did not match the canonical, normative ICNs and, therefore, could not be defined functionally, with their location likely varying by patient. This important finding suggested the level of highly individualized ICN's in the epileptic brain may signal the emergence of epileptogenic activity after surgery.

Highly synchronized cortical circuit dynamics mediate spontaneous pain in mice

Cortical neural dynamics mediate information processing for the cerebral cortex, which is implicated in fundamental biological processes such as vision and olfaction, in addition to neurological and psychiatric diseases. Spontaneous pain is a key feature of human neuropathic pain. Whether spontaneous pain pushes the cortical network into an aberrant state and, if so, whether it can be brought back to a "normal" operating range to ameliorate pain are unknown. Using a clinically relevant mouse model of neuropathic pain with spontaneous pain-like behavior, we report that orofacial spontaneous pain activated a specific area within the primary somatosensory cortex (S1), displaying synchronized neural dynamics revealed by intravital two-photon calcium imaging. This synchronization was underpinned by local GABAergic interneuron hypoactivity. Pain-induced cortical synchronization could be attenuated by manipulating local S1 networks or clinically effective pain therapies. Specifically, both chemogenetic inhibition of pain-related c-Fos-expressing neurons and selective activation of GABAergic interneurons significantly attenuated S1 synchronization. Clinically effective pain therapies including carbamazepine and nerve root decompression could also dampen S1 synchronization. More important, restoring a "normal" range of neural dynamics through attenuation of pain-induced S1 synchronization alleviated pain-like behavior. These results suggest that spontaneous pain pushed the S1 regional network into a synchronized state, whereas reversal of this synchronization alleviated pain.

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.

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.

Revisiting the involvement of tau in complex neural network remodeling: analysis of the extracellular neuronal activity in organotypic brain slice co-cultures

Objective.Tau ablation has a protective effect in epilepsy due to inhibition of the hyperexcitability/hypersynchrony. Protection may also occur in transgenic models of Alzheimer's disease by reducing the epileptic activity and normalizing the excitation/inhibition imbalance. However, it is difficult to determine the exact functions of tau, because tau knockout (tau^KO) brain networks exhibit elusive phenotypes. In this study, we aimed to further explore the physiological role of tau using brain network remodeling.Approach.The effect of tau ablation was investigated in hippocampal-entorhinal slice co-cultures during network remodeling. We recorded the spontaneous extracellular neuronal activity over 2 weeks in single-slice cultures and co-cultures from control andtau^KOmice. We compared the burst activity and applied concepts and analytical tools intended for the analysis of the network synchrony and connectivity.Main results.Comparison of the control andtau^KOco-cultures revealed that tau ablation had an anti-synchrony effect on the hippocampal-entorhinal two-slice networks at late stages of culture, in line with the literature. Differences were also found between the single-slice and co-culture conditions, which indicated that tau ablation had differential effects at the sub-network scale. For instance, tau ablation was found to have an anti-synchrony effect on the co-cultured hippocampal slices throughout the culture, possibly due to a reduction in the excitation/inhibition ratio. Conversely, tau ablation led to increased synchrony in the entorhinal slices at early stages of the co-culture, possibly due to homogenization of the connectivity distribution.Significance.The new methodology presented here proved useful for investigating the role of tau in the remodeling of complex brain-derived neural networks. The results confirm previous findings and hypotheses concerning the effects of tau ablation on neural networks. Moreover, the results suggest, for the first time, that tau has multifaceted roles that vary in different brain sub-networks.

Interictal Functional Connectivity in Focal Refractory Epilepsies Investigated by Intracranial EEG

Introduction: Focal epilepsies are diseases of neuronal excitability affecting macroscopic networks of cortical and subcortical neural structures. These networks ("epileptogenic networks") can generate pathological electrophysiological activities during seizures, and also between seizures (interictal period). Many works attempt to describe these networks by using quantification methods, particularly based on the estimation of statistical relationships between signals produced by brain regions, namely functional connectivity (FC). Results: FC has been shown to be greatly altered during seizures and in the immediate peri-ictal period. An increasing number of studies have shown that FC is also altered during the interictal period depending on the degree of epileptogenicity of the structures. Furthermore, connectivity values could be correlated with other clinical variables including surgical outcome. Significance: This leads to a conceptual change and to consider epileptic areas as both hyperexcitable and abnormally connected. These data open the door to the use of interictal FC as a marker of epileptogenicity and as a complementary tool for predicting the effect of surgery. Aim: In this article, we review the available data concerning interictal FC estimated from intracranial electroencephalograhy (EEG) in focal epilepsies and discuss it in the light of data obtained from other modalities (EEG imaging) and modeling studies.

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.

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.

Synthetic Epileptic Brain Activities with TripleGAN

Epilepsy is a chronic noninfectious disease caused by sudden abnormal discharge of brain neurons, which leads to intermittent brain dysfunction. It is also one of the most common neurological diseases in the world. The automatic detection of epilepsy based on electroencephalogram through machine learning, correlation analysis, and temporal-frequency analysis plays an important role in epilepsy early warning and automatic recognition. In this study, we propose a method to realize EEG epilepsy recognition by means of triple genetic antagonism network (GAN). TripleGAN is used for EEG temporal domain, frequency domain, and temporal-frequency domain, respectively. The experiment was conducted through CHB-MIT datasets, which operated at the latest level in the same industry in the world. In the CHB-MIT dataset, the classification accuracy, sensitivity, and specificity exceeded 1.19%, 1.36%, and 0.27%, respectively. The crossobject ratio exceeded 0.53%, 2.2%, and 0.37%, respectively. It shows that the established deep learning model of TripleGAN has a good effect on EEG epilepsy classification through simulation and classification optimization of real signals.

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.

Normative intracranial EEG maps epileptogenic tissues in focal epilepsy

Planning surgery for patients with medically refractory epilepsy often requires recording seizures using intracranial EEG. Quantitative measures derived from interictal intracranial EEG yield potentially appealing biomarkers to guide these surgical procedures; however, their utility is limited by the sparsity of electrode implantation as well as the normal confounds of spatiotemporally varying neural activity and connectivity. We propose that comparing intracranial EEG recordings to a normative atlas of intracranial EEG activity and connectivity can reliably map abnormal regions, identify targets for invasive treatment and increase our understanding of human epilepsy. Merging data from the Penn Epilepsy Center and a public database from the Montreal Neurological Institute, we aggregated interictal intracranial EEG retrospectively across 166 subjects comprising >5000 channels. For each channel, we calculated the normalized spectral power and coherence in each canonical frequency band. We constructed an intracranial EEG atlas by mapping the distribution of each feature across the brain and tested the atlas against data from novel patients by generating a z-score for each channel. We demonstrate that for seizure onset zones within the mesial temporal lobe, measures of connectivity abnormality provide greater distinguishing value than univariate measures of abnormal neural activity. We also find that patients with a longer diagnosis of epilepsy have greater abnormalities in connectivity. By integrating measures of both single-channel activity and inter-regional functional connectivity, we find a better accuracy in predicting the seizure onset zones versus normal brain (area under the curve = 0.77) compared with either group of features alone. We propose that aggregating normative intracranial EEG data across epilepsy centres into a normative atlas provides a rigorous, quantitative method to map epileptic networks and guide invasive therapy. We publicly share our data, infrastructure and methods, and propose an international framework for leveraging big data in surgical planning for refractory epilepsy.

Interictal SEEG Resting-State Connectivity Localizes the Seizure Onset Zone and Predicts Seizure Outcome

Localization of epileptogenic zone currently requires prolonged intracranial recordings to capture seizure, which may take days to weeks. The authors developed a novel method to identify the seizure onset zone (SOZ) and predict seizure outcome using short-time resting-state stereotacticelectroencephalography (SEEG) data. In a cohort of 27 drug-resistant epilepsy patients, the authors estimated the information flow via directional connectivity and inferred the excitation-inhibition ratio from the 1/f power slope. They hypothesized that the antagonism of information flow at multiple frequencies between SOZ and non-SOZ underlying the relatively stable epilepsy resting state could be related to the disrupted excitation-inhibition balance. They found flatter 1/f power slope in non-SOZ regions compared to the SOZ, with dominant information flow from non-SOZ to SOZ regions. Greater differences in resting-state information flow between SOZ and non-SOZ regions are associated with favorable seizure outcome. By integrating a balanced random forest model with resting-state connectivity, their method localized the SOZ with an accuracy of 88% and predicted the seizure outcome with an accuracy of 92% using clinically determined SOZ. Overall, this study suggests that brief resting-state SEEG data can significantly facilitate the identification of SOZ and may eventually predict seizure outcomes without requiring long-term ictal recordings.

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.

Drug-resistant focal epilepsy in children is associated with increased modal controllability of the whole brain and epileptogenic regions

Network control theory provides a framework by which neurophysiological dynamics of the brain can be modelled as a function of the structural connectome constructed from diffusion MRI. Average controllability describes the ability of a region to drive the brain to easy-to-reach neurophysiological states whilst modal controllability describes the ability of a region to drive the brain to difficult-to-reach states. In this study, we identify increases in mean average and modal controllability in children with drug-resistant epilepsy compared to healthy controls. Using simulations, we purport that these changes may be a result of increased thalamocortical connectivity. At the node level, we demonstrate decreased modal controllability in the thalamus and posterior cingulate regions. In those undergoing resective surgery, we also demonstrate increased modal controllability of the resected parcels, a finding specific to patients who were rendered seizure free following surgery. Changes in controllability are a manifestation of brain network dysfunction in epilepsy and may be a useful construct to understand the pathophysiology of this archetypical network disease. Understanding the mechanisms underlying these controllability changes may also facilitate the design of network-focussed interventions that seek to normalise network structure and function.

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.

Quantitative Analysis of Visually Reviewed Normal Scalp EEG Predicts Seizure Freedom Following Anterior Temporal Lobectomy

Objective

Anterior temporal lobectomy (ATL) is a widely performed and successful intervention for drug‐resistant temporal lobe epilepsy (TLE). However, up to one third of patients experience seizure recurrence within 1 year after ATL. Despite the extensive literature on presurgical electroencephalography (EEG) and magnetic resonance imaging (MRI) abnormalities to prognosticate seizure freedom following ATL, the value of quantitative analysis of visually reviewed normal interictal EEG in such prognostication remains unclear. In this retrospective multicenter study, we investigate whether machine learning analysis of normal interictal scalp EEG studies can inform the prediction of postoperative seizure freedom outcomes in patients who have undergone ATL.

Methods

We analyzed normal presurgical scalp EEG recordings from 41 Mayo Clinic (MC) and 23 Cleveland Clinic (CC) patients. We used an unbiased automated algorithm to extract eyes closed awake epochs from scalp EEG studies that were free of any epileptiform activity and then extracted spectral EEG features representing (a) spectral power and (b) interhemispheric spectral coherence in frequencies between 1 and 25 Hz across several brain regions. We analyzed the differences between the seizure‐free and non–seizure‐free patients and employed a Naïve Bayes classifier using multiple spectral features to predict surgery outcomes. We trained the classifier using a leave‐one‐patient‐out cross‐validation scheme within the MC data set and then tested using the out‐of‐sample CC data set. Finally, we compared the predictive performance of normal scalp EEG‐derived features against MRI abnormalities.

Results

We found that several spectral power and coherence features showed significant differences correlated with surgical outcomes and that they were most pronounced in the 10–25 Hz range. The Naïve Bayes classification based on those features predicted 1‐year seizure freedom following ATL with area under the curve (AUC) values of 0.78 and 0.76 for the MC and CC data sets, respectively. Subsequent analyses revealed that (a) interhemispheric spectral coherence features in the 10–25 Hz range provided better predictability than other combinations and (b) normal scalp EEG‐derived features provided superior and potentially distinct predictive value when compared with MRI abnormalities (>10% higher F1 score).

Significance

These results support that quantitative analysis of even a normal presurgical scalp EEG may help prognosticate seizure freedom following ATL in patients with drug‐resistant TLE. Although the mechanism for this result is not known, the scalp EEG spectral and coherence properties predicting seizure freedom may represent activity arising from the neocortex or the networks responsible for temporal lobe seizure generation within vs outside the margins of an ATL.

Chaotic and stochastic dynamics of epileptiform-like activities in sclerotic hippocampus resected from patients with pharmacoresistant epilepsy

The types of epileptiform activity occurring in the sclerotic hippocampus with highest incidence are interictal-like events (II) and periodic ictal spiking (PIS). These activities are classified according to their event rates, but it is still unclear if these rate differences are consequences of underlying physiological mechanisms. Identifying new and more specific information related to these two activities may bring insights to a better understanding about the epileptogenic process and new diagnosis. We applied Poincaré map analysis and Recurrence Quantification Analysis (RQA) onto 35 in vitro electrophysiological signals recorded from slices of 12 hippocampal tissues surgically resected from patients with pharmacoresistant temporal lobe epilepsy. These analyzes showed that the II activity is related to chaotic dynamics, whereas the PIS activity is related to deterministic periodic dynamics. Additionally, it indicates that their different rates are consequence of different endogenous dynamics. Finally, by using two computational models we were able to simulate the transition between II and PIS activities. The RQA was applied to different periods of these simulations to compare the recurrences between artificial and real signals, showing that different ranges of regularity-chaoticity can be directly associated with the generation of PIS and II activities.

Temporal lobe epilepsy (TLE) is the most prevalent type of epilepsy in adults and hippocampal sclerosis is the major pathophysiological substrate of pharmaco-refractory TLE. Different patterns of epileptiform-like activity have been described in human hippocampal sclerosis, but the standard analysis applied to characterize the activities usually do not consider the nonlinear features that epileptiform patterns exhibit. Here, using Poincaré map and Recurrence Quantitative Analysis we characterized the most prevalent type of epileptiform-like activities—interictal-like events (II) and periodic ictal spiking (PIS), recorded in vitro from resected hippocampi of pharmacoresistant patients with TLE—according to their levels of stochasticity, chaoticity and determinism. The II activities showed to be more chaotic with complex rhythmicity than PIS activities. The nonlinear dynamic differences between II and PIS leads us to conjecture that they are expressions of different seizure susceptibility. We also identified that each hippocampal subfield expresses II and PIS activities in a specific and different way. Finally, from the modulation of internal parameters of two computational models, we show the conversion of one type of activity into the other, showing how specific neuron networks synchronize over time, leading to II and PIS activities and then into a generalized seizure.

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.

A framework For brain atlases: Lessons from seizure dynamics

Brain maps, or atlases, are essential tools for studying brain function and organization. The abundance of available atlases used across the neuroscience literature, however, creates an implicit challenge that may alter the hypotheses and predictions we make about neurological function and pathophysiology. Here, we demonstrate how parcellation scale, shape, anatomical coverage, and other atlas features may impact our prediction of the brain’s function from its underlying structure. We show how network topology, structure–function correlation (SFC), and the power to test specific hypotheses about epilepsy pathophysiology may change as a result of atlas choice and atlas features. Through the lens of our disease system, we propose a general framework and algorithm for atlas selection. This framework aims to maximize the descriptive, explanatory, and predictive validity of an atlas. Broadly, our framework strives to provide empirical guidance to neuroscience research utilizing the various atlases published over the last century.

Intracranial electroencephalographic biomarker predicts effective responsive neurostimulation for epilepsy prior to treatment

OBJECTIVE

Despite the overall success of responsive neurostimulation (RNS) therapy for drug-resistant focal epilepsy, clinical outcomes in individuals vary significantly and are hard to predict. Biomarkers that indicate the clinical efficacy of RNS-ideally before device implantation-are critically needed, but challenges include the intrinsic heterogeneity of the RNS patient population and variability in clinical management across epilepsy centers. The aim of this study is to use a multicenter dataset to evaluate a candidate biomarker from intracranial electroencephalographic (iEEG) recordings that predicts clinical outcome with subsequent RNS therapy.

METHODS

We assembled a federated dataset of iEEG recordings, collected prior to RNS implantation, from a retrospective cohort of 30 patients across three major epilepsy centers. Using ictal iEEG recordings, each center independently calculated network synchronizability, a candidate biomarker indicating the susceptibility of epileptic brain networks to RNS therapy.

RESULTS

Ictal measures of synchronizability in the high-γ band (95-105 Hz) significantly distinguish between good and poor RNS responders after at least 3 years of therapy under the current RNS therapy guidelines (area under the curve = .83). Additionally, ictal high-γ synchronizability is inversely associated with the degree of therapeutic response.

SIGNIFICANCE

This study provides a proof-of-concept roadmap for collaborative biomarker evaluation in federated data, where practical considerations impede full data sharing across centers. Our results suggest that network synchronizability can help predict therapeutic response to RNS therapy. With further validation, this biomarker could facilitate patient selection and help avert a costly, invasive intervention in patients who are unlikely to benefit.

Virtual intracranial EEG signals reconstructed from MEG with potential for epilepsy surgery

Modelling the interactions that arise from neural dynamics in seizure genesis is challenging but important in the effort to improve the success of epilepsy surgery. Dynamical network models developed from physiological evidence offer insights into rapidly evolving brain networks in the epileptic seizure. A limitation of previous studies in this field is the dependence on invasive cortical recordings with constrained spatial sampling of brain regions that might be involved in seizure dynamics. Here, we propose virtual intracranial electroencephalography (ViEEG), which combines non-invasive ictal magnetoencephalographic imaging (MEG), dynamical network models and a virtual resection technique. In this proof-of-concept study, we show that ViEEG signals reconstructed from MEG alone preserve critical temporospatial characteristics for dynamical approaches to identify brain areas involved in seizure generation. We show the non-invasive ViEEG approach may have some advantage over intracranial electroencephalography (iEEG). Future work may be designed to test the potential of the virtual iEEG approach for use in surgical management of epilepsy.

Reproducibility of graph measures derived from resting-state MEG functional connectivity metrics in sensor and source spaces

Prior studies have used graph analysis of resting‐state magnetoencephalography (MEG) to characterize abnormal brain networks in neurological disorders. However, a present challenge for researchers is the lack of guidance on which network construction strategies to employ. The reproducibility of graph measures is important for their use as clinical biomarkers. Furthermore, global graph measures should ideally not depend on whether the analysis was performed in the sensor or source space. Therefore, MEG data of the 89 healthy subjects of the Human Connectome Project were used to investigate test–retest reliability and sensor versus source association of global graph measures. Atlas‐based beamforming was used for source reconstruction, and functional connectivity (FC) was estimated for both sensor and source signals in six frequency bands using the debiased weighted phase lag index (dwPLI), amplitude envelope correlation (AEC), and leakage‐corrected AEC. Reliability was examined over multiple network density levels achieved with proportional weight and orthogonal minimum spanning tree thresholding. At a 100% density, graph measures for most FC metrics and frequency bands had fair to excellent reliability and significant sensor versus source association. The greatest reliability and sensor versus source association was obtained when using amplitude metrics. Reliability was similar between sensor and source spaces when using amplitude metrics but greater for the source than the sensor space in higher frequency bands when using the dwPLI. These results suggest that graph measures are useful biomarkers, particularly for investigating functional networks based on amplitude synchrony.

Transfer Entropy between Intracranial EEG Nodes Highlights Network Dynamics that Cause and Stop Epileptic Seizures

Transfer entropy (TE) is used to examine the connectivity between nodes and the roles of nodes in epileptic neural networks during rest, moments before seizure, during seizure, and moments after seizure. There is a set of nodes that dominate information flow to epileptogenic zone (EZ) nodes, regions that trigger seizure, and non-EZ nodes during rest. The TE from the dominant to the EZ nodes decreases shortly before a seizure event and reaches a minimum during seizure. During the seizure, the dominant nodes cease or only weakly interact with the EZ nodes. This supports the hypothesis that seizure occurs when some nodes stop inhibiting the EZ nodes. The TE from the dominant to the EZ nodes peaks immediately after seizure, suggesting that seizure may stop when the brain exerts the highest level of information flow/activation/communication to the EZ nodes. The information flow from the dominant to EZ nodes is different from that to non-EZ nodes. This TE dynamics entering and exiting seizures may identify more accurately the EZ nodes, which may improve surgical planning.

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.

Leveraging electrophysiologic correlates of word encoding to map seizure onset zone in focal epilepsy: Task-dependent changes in epileptiform activity, spectral features, and functional connectivity

OBJECTIVE

Verbal memory dysfunction is common in focal, drug-resistant epilepsy (DRE). Unfortunately, surgical removal of seizure-generating brain tissue can be associated with further memory decline. Therefore, localization of both the circuits generating seizures and those underlying cognitive functions is critical in presurgical evaluations for patients who may be candidates for resective surgery. We used intracranial electroencephalographic (iEEG) recordings during a verbal memory task to investigate word encoding in focal epilepsy. We hypothesized that engagement in a memory task would exaggerate local iEEG feature differences between the seizure onset zone (SOZ) and neighboring tissue as compared to wakeful rest ("nontask").

METHODS

Ten participants undergoing presurgical iEEG evaluation for DRE performed a free recall verbal memory task. We evaluated three iEEG features in SOZ and non-SOZ electrodes during successful word encoding and compared them with nontask recordings: interictal epileptiform spike (IES) rates, power in band (PIB), and relative entropy (REN; a functional connectivity measure).

RESULTS

We found a complex pattern of PIB and REN changes in SOZ and non-SOZ electrodes during successful word encoding compared to nontask. Successful word encoding was associated with a reduction in local electrographic functional connectivity (increased REN), which was most exaggerated in temporal lobe SOZ. The IES rates were reduced during task, but only in the non-SOZ electrodes. Compared with nontask, REN features during task yielded marginal improvements in SOZ classification.

SIGNIFICANCE

Previous studies have supported REN as a biomarker for epileptic brain. We show that REN differences between SOZ and non-SOZ are enhanced during a verbal memory task. We also show that IESs are reduced during task in non-SOZ, but not in SOZ. These findings support the hypothesis that SOZ and non-SOZ respond differently to task and warrant further exploration into the use of cognitive tasks to identify functioning memory circuits and localize SOZ.

Patient-Specific Network Connectivity Combined With a Next Generation Neural Mass Model to Test Clinical Hypothesis of Seizure Propagation

Dynamics underlying epileptic seizures span multiple scales in space and time, therefore, understanding seizure mechanisms requires identifying the relations between seizure components within and across these scales, together with the analysis of their dynamical repertoire. In this view, mathematical models have been developed, ranging from single neuron to neural population. In this study, we consider a neural mass model able to exactly reproduce the dynamics of heterogeneous spiking neural networks. We combine mathematical modeling with structural information from non invasive brain imaging, thus building large-scale brain network models to explore emergent dynamics and test the clinical hypothesis. We provide a comprehensive study on the effect of external drives on neuronal networks exhibiting multistability, in order to investigate the role played by the neuroanatomical connectivity matrices in shaping the emergent dynamics. In particular, we systematically investigate the conditions under which the network displays a transition from a low activity regime to a high activity state, which we identify with a seizure-like event. This approach allows us to study the biophysical parameters and variables leading to multiple recruitment events at the network level. We further exploit topological network measures in order to explain the differences and the analogies among the subjects and their brain regions, in showing recruitment events at different parameter values. We demonstrate, along with the example of diffusion-weighted magnetic resonance imaging (dMRI) connectomes of 20 healthy subjects and 15 epileptic patients, that individual variations in structural connectivity, when linked with mathematical dynamic models, have the capacity to explain changes in spatiotemporal organization of brain dynamics, as observed in network-based brain disorders. In particular, for epileptic patients, by means of the integration of the clinical hypotheses on the epileptogenic zone (EZ), i.e., the local network where highly synchronous seizures originate, we have identified the sequence of recruitment events and discussed their links with the topological properties of the specific connectomes. The predictions made on the basis of the implemented set of exact mean-field equations turn out to be in line with the clinical pre-surgical evaluation on recruited secondary networks.

Long-term brain network reorganization predicts responsive neurostimulation outcomes for focal epilepsy

Responsive neurostimulation (RNS) devices, able to detect imminent seizures and to rapidly deliver electrical stimulation to the brain, are effective in reducing seizures in some patients with focal epilepsy. However, therapeutic response to RNS is often slow, is highly variable, and defies prognostication based on clinical factors. A prevailing view holds that RNS efficacy is primarily mediated by acute seizure termination; yet, stimulations greatly outnumber seizures and occur mostly in the interictal state, suggesting chronic modulation of brain networks that generate seizures. Here, using years-long intracranial neural recordings collected during RNS therapy, we found that patients with the greatest therapeutic benefit undergo progressive, frequency-dependent reorganization of interictal functional connectivity. The extent of this reorganization scales directly with seizure reduction and emerges within the first year of RNS treatment, enabling potential early prediction of therapeutic response. Our findings reveal a mechanism for RNS that involves network plasticity and may inform development of next-generation devices for epilepsy.

Electrocorticography and stereo EEG provide distinct measures of brain connectivity: implications for network models

Brain network models derived from graph theory have the potential to guide functional neurosurgery, and to improve rates of post-operative seizure freedom for patients with epilepsy. A barrier to applying these models clinically is that intracranial EEG electrode implantation strategies vary by centre, region and country, from cortical grid & strip electrodes (Electrocorticography), to purely stereotactic depth electrodes (Stereo EEG), to a mixture of both. To determine whether models derived from one type of study are broadly applicable to others, we investigate the differences in brain networks mapped by electrocorticography and stereo EEG in a cohort of patients who underwent surgery for temporal lobe epilepsy and achieved a favourable outcome. We show that networks derived from electrocorticography and stereo EEG define distinct relationships between resected and spared tissue, which may be driven by sampling bias of temporal depth electrodes in patients with predominantly cortical grids. We propose a method of correcting for the effect of internodal distance that is specific to electrode type and explore how additional methods for spatially correcting for sampling bias affect network models. Ultimately, we find that smaller surgical targets tend to have lower connectivity with respect to the surrounding network, challenging notions that abnormal connectivity in the epileptogenic zone is typically high. Our findings suggest that effectively applying computational models to localize epileptic networks requires accounting for the effects of spatial sampling, particularly when analysing both electrocorticography and stereo EEG recordings in the same cohort, and that future network studies of epilepsy surgery should also account for differences in focality between resection and ablation. We propose that these findings are broadly relevant to intracranial EEG network modelling in epilepsy and an important step in translating them clinically into patient care.

To Stop the Seizure, You Must Become the Seizure: Closed-Loop Stimulation Phase Locked to Seizure Waves Disrupts Kindled Seizures While Open-Loop Stimulation Fails

Closed-Loop Stimulation of the Medial Septum Terminates Epileptic Seizures Takeuchi Y, Harangozó M, Pedraza L, Földi T, Kozák G, Li Q, Berényi A. Brain . 2021;144(3):885-908. doi:10.1093/brain/awaa450 Temporal lobe epilepsy with distributed hippocampal seizure foci is often intractable and its secondary generalization might lead to sudden death. Early termination through spatially extensive hippocampal intervention is not feasible directly, because of the large size and irregular shape of the hippocampus. In contrast, the medial septum is a promising target to govern hippocampal oscillations through its divergent connections to both hippocampi. Combining this “proxy intervention” concept and precisely timed stimulation, we report here that closed-loop medial septum electrical stimulation can quickly terminate intrahippocampal seizures and suppress secondary generalization in a rat kindling model. Precise stimulus timing governed by internal seizure rhythms was essential. Cell type-specific stimulation revealed that the precisely timed activation of medial septum GABAergic neurons underlaid the effects. Our concept of time-targeted proxy stimulation for intervening pathological oscillations can be extrapolated to other neurological and psychiatric disorders, and has potential for clinical translation.

The use of computational models in the management and prognosis of refractory epilepsy: A critical evaluation

PURPOSE

Drug resistant epilepsy (DRE) affects approximately 30 percent of individuals with epilepsy worldwide. Surgery remains the most effective treatment for individuals with DRE, but referral to surgery is low and only about 60 percent of individuals who undergo surgery experience seizure control postoperatively. The present paper evaluates the evidence for using computational models in the prediction of surgical resection sites and surgical outcomes for patients with DRE.

METHODS

We conducted a search in the Medline data base using the terms "refractory epilepsy", "drug-resistant epilepsy", "surgery", "computational model", and "artificial intelligence". Inclusion: original articles in English and case reports from 2000 to 2020. Reviews were excluded.

RESULTS

Clinical applications of computational models may lead to increased utilisation of surgical services through improving our ability to predict outcomes and by improving surgical outcomes outright. The identification and optimisation of nodes that are crucial for the genesis and propagation of epileptiform activity offers the most promising clinical applications of computational models discussed herein.

CONCLUSION

Advances in computational models may in the future significantly increase the application and efficacy of surgery for patients with DRE by optimising the site and amount of cortex to resect, but more research is needed before it achieves therapeutic utility.

The Effect of Anti-seizure Medications on the Propagation of Epileptic Activity: A Review

The propagation of epileptiform events is a highly interesting phenomenon from the pathophysiological point of view, as it involves several mechanisms of recruitment of neural networks. Extensive in vivo and in vitro research has been performed, suggesting that multiple networks as well as cellular candidate mechanisms govern this process, including the co-existence of wave propagation, coupled oscillator dynamics, and more. The clinical importance of seizure propagation stems mainly from the fact that the epileptic manifestations cannot be attributed solely to the activity in the seizure focus itself, but rather to the propagation of epileptic activity to other brain structures. Propagation, especially when causing secondary generalizations, poses a risk to patients due to recurrent falls, traumatic injuries, and poor neurological outcome. Anti-seizure medications (ASMs) affect propagation in diverse ways and with different potencies. Importantly, for drug-resistant patients, targeting seizure propagation may improve the quality of life even without a major reduction in simple focal events. Motivated by the extensive impact of this phenomenon, we sought to review the literature regarding the propagation of epileptic activity and specifically the effect of commonly used ASMs on it. Based on this body of knowledge, we propose a novel classification of ASMs into three main categories: major, minor, and intermediate efficacy in reducing the propagation of epileptiform activity.

Pairwise maximum entropy model explains the role of white matter structure in shaping emergent co-activation states

A major challenge in neuroscience is determining a quantitative relationship between the brain's white matter structural connectivity and emergent activity. We seek to uncover the intrinsic relationship among brain regions fundamental to their functional activity by constructing a pairwise maximum entropy model (MEM) of the inter-ictal activation patterns of five patients with medically refractory epilepsy over an average of ~14 hours of band-passed intracranial EEG (iEEG) recordings per patient. We find that the pairwise MEM accurately predicts iEEG electrodes' activation patterns' probability and their pairwise correlations. We demonstrate that the estimated pairwise MEM's interaction weights predict structural connectivity and its strength over several frequencies significantly beyond what is expected based solely on sampled regions' distance in most patients. Together, the pairwise MEM offers a framework for explaining iEEG functional connectivity and provides insight into how the brain's structural connectome gives rise to large-scale activation patterns by promoting co-activation between connected structures.

Time-evolving controllability of effective connectivity networks during seizure progression

Over one third of the estimated 3 million people with epilepsy in the United States are medication resistant. Responsive neurostimulation from chronically implanted electrodes provides a promising treatment alternative to resective surgery. However, determining optimal personalized stimulation parameters, including when and where to intervene to guarantee a positive patient outcome, is a major open challenge. Network neuroscience and control theory offer useful tools that may guide improvements in parameter selection for control of anomalous neural activity. Here we use a method to characterize dynamic controllability across consecutive effective connectivity (EC) networks based on regularized partial correlations between implanted electrodes during the onset, propagation, and termination regimes of 34 seizures. We estimate regularized partial correlation adjacency matrices from 1-s time windows of intracranial electrocorticography recordings using the Graphical Least Absolute Shrinkage and Selection Operator (GLASSO). Average and modal controllability metrics calculated from each resulting EC network track the time-varying controllability of the brain on an evolving landscape of conditionally dependent network interactions. We show that average controllability increases throughout a seizure and is negatively correlated with modal controllability throughout. Our results support the hypothesis that the energy required to drive the brain to a seizure-free state from an ictal state is smallest during seizure onset, yet we find that applying control energy at electrodes in the seizure onset zone may not always be energetically favorable. Our work suggests that a low-complexity model of time-evolving controllability may offer insights for developing and improving control strategies targeting seizure suppression.