Dynamic Network Drivers of Seizure Generation, Propagation and Termination in Human Neocortical Epilepsy

Brain connectivity analysis in preictal phases of seizure induced by pentylenetetrazol in rats

Abnormal patterns of brain connectivity characterize epilepsy. However, little is known about these patterns during the stages preceding a seizure induced by pentylenetetrazol (PTZ). To investigate brain connectivity in male Wistar rats during the preictal phase of PTZ-induced seizures (60 mg/kg), we recorded local field potentials in the primary motor (M1) cortex, the ventral anterior (VA) nucleus of the thalamus, the hippocampal CA1 area, and the dentate gyrus (DG) during the baseline period and after PTZ administration. While there were no changes in power density between the baseline and preictal periods, we observed an increase in directional functional connectivity in theta from the hippocampal formation to M1 and VA, as well as in middle gamma from DG to CA1 and from CA1 to M1, and also in slow gamma from M1 to CA1. These findings are supported by increased phase coherence between DG-M1 in theta and CA1-M1 in middle gamma, as well as enhanced phase-amplitude coupling of delta-middle gamma in M1 and delta-fast gamma in CA1. Interestingly, we also noted a slight decrease in phase synchrony between CA1 and VA in slow gamma. Together, these results demonstrate increased functional connectivity between brain regions during the PTZ-induced preictal period, with this increase being particularly driven by the hippocampal formation.

Subject-based assessment of large-scale integration dynamics in epileptic brain networks: insights from the intrinsic ignition framework

This study examined the dynamic properties of brain regions involved in the genesis and spread of seizures in 10 individuals diagnosed with pharmacoresistant focal epilepsy. The patients and 30 healthy controls underwent resting-state functional magnetic resonance imaging scans and the brain's functional network dynamics were analyzed using the intrinsic ignition framework. Comparative statistical analyses examined the differences in the integration and metastability measures in both groups in the whole brain and specific local brain regions. Invasive electroencephalography evaluations validated the findings of significant global and regional changes in the patient's brain network dynamics. There was a marked increase in global integration and metastability across the brain, reflecting substantial alterations in the overall connectivity and flexibility of the functional networks. Specific brain regions exhibited paradoxical dynamics within the seizure onset zone, with decreased intrinsic ignition and increased metastability. Increased intrinsic ignition was observed in remote brain regions, suggesting a reorganization of the brain network hubs and potential pathways for seizure propagation. Using the intrinsic ignition framework provided insights into dynamic alterations in the brain networks of patients with epilepsy. These have increased our understanding of the mechanisms underlying epileptic seizures and may guide the development of diagnostic biomarkers and targeted therapeutic interventions.

Explaining slow seizure propagation with white matter tractography

Epileptic seizures recorded with stereo-EEG can take a fraction of a second or several seconds to propagate from one region to another. What explains such propagation patterns? We combine tractography and stereo-EEG to determine the relationship between seizure propagation and the white matter architecture and to describe seizure propagation mechanisms. Patient-specific spatiotemporal seizure propagation maps were combined with tractography from diffusion imaging of matched subjects from the Human Connectome Project. The onset of seizure activity was marked on a channel-by-channel basis by two board-certified neurologists for all channels involved in the seizure. We measured the tract connectivity (number of tracts) between regions-of-interest pairs among the seizure onset zone, regions of seizure spread and non-involved regions. We also investigated how tract-connected the seizure onset zone is to regions of early seizure spread compared with regions of late spread. Comparisons were made after correcting for differences in distance. Sixty-nine seizures were marked across 26 patients with drug-resistant epilepsy; 11 were seizure free after surgery (Engel IA) and 15 were not (Engel IB-Engel IV). The seizure onset zone was more tract-connected to regions of seizure spread than to non-involved regions (P < 0.0001); however, regions of seizure spread were not differentially tract-connected to other regions of seizure spread compared with non-involved regions. In seizure-free patients only, regions of seizure spread were more tract-connected to the seizure onset zone than to other regions of spread (P < 0.0001). Over the temporal evolution of a seizure, the seizure onset zone was significantly more tract-connected to regions of early spread compared with regions of late spread in seizure-free patients only (P < 0.0001). By integrating information on structure, we demonstrate that seizure propagation is likely to be mediated by white matter tracts. The pattern of connectivity between seizure onset zone, regions of spread and non-involved regions demonstrates that the onset zone might be largely responsible for seizures propagating throughout the brain, rather than seizures propagating to intermediate points, from which further propagation takes place. Our findings also suggest that seizure propagation over seconds might be the result of a continuous bombardment of action potentials from the seizure onset zone to regions of spread. In non-seizure-free patients, the paucity of tracts from the presumed seizure onset zone to regions of spread suggests that the onset zone was missed. Fully understanding the structure-propagation relationship might eventually provide insight into selecting the correct targets for epilepsy surgery.

A data augmentation procedure to improve detection of spike ripples in brain voltage recordings

Epilepsy is a major neurological disorder characterized by recurrent, spontaneous seizures. For patients with drug-resistant epilepsy, treatments include neurostimulation or surgical removal of the epileptogenic zone (EZ), the brain region responsible for seizure generation. Precise targeting of the EZ requires reliable biomarkers. Spike ripples - high-frequency oscillations that co-occur with large amplitude epileptic discharges - have gained prominence as a candidate biomarker. However, spike ripple detection remains a challenge. The gold-standard approach requires an expert manually visualize and interpret brain voltage recordings, which limits reproducibility and high-throughput analysis. Addressing these limitations requires more objective, efficient, and automated methods for spike ripple detection, including approaches that utilize deep neural networks. Despite advancements, dataset heterogeneity and scarcity severely limit machine learning performance. Our study explores long-short term memory (LSTM) neural network architectures for spike ripple detection, leveraging data augmentation to improve classifier performance. We highlight the potential of combining training on augmented and in vivo data for enhanced spike ripple detection and ultimately improving diagnostic accuracy in epilepsy treatment.

Connectome reorganization associated with temporal lobe pathology and its surgical resection

Network neuroscience offers a unique framework to understand the organizational principles of the human brain. Despite recent progress, our understanding of how the brain is modulated by focal lesions remains incomplete. Resection of the temporal lobe is the most effective treatment to control seizures in pharmaco-resistant temporal lobe epilepsy (TLE), making this syndrome a powerful model to study lesional effects on network organization in young and middle-aged adults. Here, we assessed the downstream consequences of a focal lesion and its surgical resection on the brain's structural connectome, and explored how this reorganization relates to clinical variables at the individual patient level. We included adults with pharmaco-resistant TLE (n = 37) who underwent anterior temporal lobectomy between two imaging time points, as well as age- and sex-matched healthy controls who underwent comparable imaging (n = 31). Core to our analysis was the projection of high-dimensional structural connectome data-derived from diffusion MRI tractography from each subject-into lower-dimensional gradients. We then compared connectome gradients in patients relative to controls before surgery, tracked surgically-induced connectome reconfiguration from pre- to postoperative time points, and examined associations to patient-specific clinical and imaging phenotypes. Before surgery, individuals with TLE presented with marked connectome changes in bilateral temporo-parietal regions, reflecting an increased segregation of the ipsilateral anterior temporal lobe from the rest of the brain. Surgery-induced connectome reorganization was localized to this temporo-parietal subnetwork, but primarily involved postoperative integration of contralateral regions with the rest of the brain. Using a partial least-squares analysis, we uncovered a latent clinical imaging signature underlying this pre- to postoperative connectome reorganization, showing that patients who displayed postoperative integration in bilateral fronto-occipital cortices also had greater preoperative ipsilateral hippocampal atrophy, lower seizure frequency and secondarily generalized seizures. Our results bridge the effects of focal brain lesions and their surgical resections with large-scale network reorganization and interindividual clinical variability, thus offering new avenues to examine the fundamental malleability of the human brain.

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.

Hub overload and failure as a final common pathway in neurological brain network disorders

Understanding the concept of network hubs and their role in brain disease is now rapidly becoming important for clinical neurology. Hub nodes in brain networks are areas highly connected to the rest of the brain, which handle a large part of all the network traffic. They also show high levels of neural activity and metabolism, which makes them vulnerable to many different types of pathology. The present review examines recent evidence for the prevalence and nature of hub involvement in a variety of neurological disorders, emphasizing common themes across different types of pathology. In focal epilepsy, pathological hubs may play a role in spreading of seizure activity, and removal of such hub nodes is associated with improved outcome. In stroke, damage to hubs is associated with impaired cognitive recovery. Breakdown of optimal brain network organization in multiple sclerosis is accompanied by cognitive dysfunction. In Alzheimer's disease, hyperactive hub nodes are directly associated with amyloid-beta and tau pathology. Early and reliable detection of hub pathology and disturbed connectivity in Alzheimer's disease with imaging and neurophysiological techniques opens up opportunities to detect patients with a network hyperexcitability profile, who could benefit from treatment with anti-epileptic drugs.

Community detection in the human connectome: Method types, differences and their impact on inference

Community structure is a fundamental topological characteristic of optimally organized brain networks. Currently, there is no clear standard or systematic approach for selecting the most appropriate community detection method. Furthermore, the impact of method choice on the accuracy and robustness of estimated communities (and network modularity), as well as method-dependent relationships between network communities and cognitive and other individual measures, are not well understood. This study analyzed large datasets of real brain networks (estimated from resting-state fMRI fromn $$ n $$ = 5251 pre/early adolescents in the adolescent brain cognitive development [ABCD] study), andn $$ n $$ = 5338 synthetic networks with heterogeneous, data-inspired topologies, with the goal to investigate and compare three classes of community detection methods: (i) modularity maximization-based (Newman and Louvain), (ii) probabilistic (Bayesian inference within the framework of stochastic block modeling (SBM)), and (iii) geometric (based on graph Ricci flow). Extensive comparisons between methods and their individual accuracy (relative to the ground truth in synthetic networks), and reliability (when applied to multiple fMRI runs from the same brains) suggest that the underlying brain network topology plays a critical role in the accuracy, reliability and agreement of community detection methods. Consistent method (dis)similarities, and their correlations with topological properties, were estimated across fMRI runs. Based on synthetic graphs, most methods performed similarly and had comparable high accuracy only in some topological regimes, specifically those corresponding to developed connectomes with at least quasi-optimal community organization. In contrast, in densely and/or weakly connected networks with difficult to detect communities, the methods yielded highly dissimilar results, with Bayesian inference within SBM having significantly higher accuracy compared to all others. Associations between method-specific modularity and demographic, anthropometric, physiological and cognitive parameters showed mostly method invariance but some method dependence as well. Although method sensitivity to different levels of community structure may in part explain method-dependent associations between modularity estimates and parameters of interest, method dependence also highlights potential issues of reliability and reproducibility. These findings suggest that a probabilistic approach, such as Bayesian inference in the framework of SBM, may provide consistently reliable estimates of community structure across network topologies. In addition, to maximize robustness of biological inferences, identified network communities and their cognitive, behavioral and other correlates should be confirmed with multiple reliable detection methods.

High frequency oscillation network dynamics predict outcome in non-palliative epilepsy surgery

High frequency oscillations are a promising biomarker of outcome in intractable epilepsy. Prior high frequency oscillation work focused on counting high frequency oscillations on individual channels, and it is still unclear how to translate those results into clinical care. We show that high frequency oscillations arise as network discharges that have valuable properties as predictive biomarkers. Here, we develop a tool to predict patient outcome before surgical resection is performed, based on only prospective information. In addition to determining high frequency oscillation rate on every channel, we performed a correlational analysis to evaluate the functional connectivity of high frequency oscillations in 28 patients with intracranial electrodes. We found that high frequency oscillations were often not solitary events on a single channel, but part of a local network discharge. Eigenvector and outcloseness centrality were used to rank channel importance within the connectivity network, then used to compare patient outcome by comparison with the seizure onset zone or a proportion within the proposed resected channels (critical resection percentage). Combining the knowledge of each patient's seizure onset zone resection plan along with our computed high frequency oscillation network centralities and high frequency oscillation rate, we develop a Naïve Bayes model that predicts outcome (positive predictive value: 100%) better than predicting based upon fully resecting the seizure onset zone (positive predictive value: 71%). Surgical margins had a large effect on outcomes: non-palliative patients in whom most of the seizure onset zone was resected ('definitive surgery', ≥ 80% resected) had predictable outcomes, whereas palliative surgeries (<80% resected) were not predictable. These results suggest that the addition of network properties of high frequency oscillations is more accurate in predicting patient outcome than seizure onset zone alone in patients with most of the seizure onset zone removed and offer great promise for informing clinical decisions in surgery for refractory epilepsy.

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.

Local and Remote Chemogenetic Suppression of Hippocampal Seizures in Rats

BACKGROUND

Innovative treatments of refractory epilepsy are widely desired, for which chemogenetic technology can provide region- and cell-type-specific modulation with relative noninvasiveness.

OBJECTIVES

We aimed to explore the specific applications of chemogenetics for locally and remotely networks controlling hippocampal seizures.

METHODS

A virus coding for a modified human Gi-coupled M4 muscarinic receptor (hM4Di) on pyramidal cells was injected into either the right hippocampal CA3 or the bilateral anterior nucleus of the thalamus (ANT) in rats. After one month, seizures were induced by 4-aminopyridine (4-AP) injection into the right CA3. Simultaneously, clozapine-N-oxide (CNO) (2.5 mg/kg) or clozapine (0.1 mg/kg), the specific ligands acting on hM4Di, were injected intraperitoneally. We also set up hM4Di control and clozapine control groups to eliminate the influence of viral transfection and the ligand alone on the experimental results.

RESULTS

For both local and remote controls, the mean seizure duration was significantly reduced upon ligand application in the experimental groups. Seizure frequency, on the other hand, only showed a significant decrease in local control, with a lower frequency in the clozapine group than in the CNO group. Both the effects of CNO and clozapine were time-dependent, and clozapine was faster than CNO in local seizure control.

CONCLUSION

This study shows the potency of chemogenetics to attenuate hippocampal seizures locally or remotely by activating the transfected hM4Di receptor with CNO or clozapine. ANT is suggested as a potentially safe chemogenetic application target in the epileptic network for focal hippocampal seizures.

SEEG-based epileptic seizure network modeling and analysis for pre-surgery evaluation

Stereo-electroencephalography is a minimally invasive technique for patients with refractory epilepsy pursuing surgery to reduce or control seizures. Electrodes are implanted based on pre-surgery evaluations and can collect deep brain activities for surgery decisions. This paper presents a methodology to analyze stereo-electroencephalography and assist clinicians by recommending the optimal surgical option and target areas for focal epilepsy patients. A seizure network (graph) model is proposed to characterize the spatial distribution and temporal changes of ictal events. The network nodes and edges correspond to specific epileptogenic regions and propagation/impact pathways (weighted by directed transfer function), respectively. We then employ a K-means clustering strategy to group nodes into a few clusters, from which the target surgical areas can be identified. Ten patients with different types of focal seizures were thoroughly analyzed. Promising consistency between results of our method's recommendations, clinical decisions and surgery outcomes were observed.

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.

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.

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.

Dynamic flexibility and controllability of network communities in juvenile myoclonic epilepsy

Juvenile myoclonic epilepsy (JME) is the most common syndrome within the idiopathic generalized epilepsy spectrum, manifested by myoclonic and generalized tonic-clonic seizures and spike-and-wave discharges (SWDs) on electroencephalography (EEG). Currently, the pathophysiological concepts addressing SWD generation in JME are still incomplete. In this work, we characterize the temporal and spatial organization of functional networks and their dynamic properties as derived from high-density EEG (hdEEG) recordings and MRI in 40 JME patients (25.4 ± 7.6 years, 25 females). The adopted approach allows for the construction of a precise dynamic model of ictal transformation in JME at the cortical and deep brain nuclei source levels. We implement Louvain algorithm to attribute brain regions with similar topological properties to modules during separate time windows before and during SWD generation. Afterwards, we quantify how modular assignments evolve and steer through different states towards the ictal state by measuring characteristics of flexibility and controllability. We find antagonistic dynamics of flexibility and controllability within network modules as they evolve towards and undergo ictal transformation. Prior to SWD generation, we observe concomitantly increasing flexibility (F(1,39) = 25.3, corrected p < 0.001) and decreasing controllability (F(1,39) = 55.3, p < 0.001) within the fronto-parietal module in γ-band. On a step further, during interictal SWDs as compared to preceding time windows, we notice decreasing flexibility (F(1,39) = 11.9, p < 0.001) and increasing controllability (F(1,39) = 10.1, p < 0.001) within the fronto-temporal module in γ-band. During ictal SWDs as compared to prior time windows, we demonstrate significantly decreasing flexibility (F(1,14) = 31.6; p < 0.001) and increasing controllability (F(1,14) = 44.7, p < 0.001) within the basal ganglia module. Furthermore, we show that flexibility and controllability within the fronto-temporal module of the interictal SWDs relate to seizure frequency and cognitive performance in JME patients. Our results demonstrate that detection of network modules and quantification of their dynamic properties is relevant to track the generation of SWDs. The observed flexibility and controllability dynamics reflect the reorganization of de-/synchronized connections and the ability of evolving network modules to reach a seizure-free state, respectively. These findings may advance the elaboration of network-based biomarkers and more targeted therapeutic neuromodulatory approaches in JME.

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.

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.

Progressive alterations in electrophysiological and epileptic network properties during the development of temporal lobe epilepsy in rats

OBJECTIVE

Refractory temporal lobe epilepsy (TLE) with recurring seizures causing continuing pathological changes in neural reorganization. There is an incomplete understanding of how spatiotemporal electrophysiological characteristics changes during the development of TLE. Long-term multi-site epilepsy patients' data is hard to obtain. Thus, our study relied on animal models to reveal the changes in electrophysiological and epileptic network characteristics systematically.

METHODS

Long-term local field potentials (LFPs) were recorded over a period of 1 to 4 months from 6 pilocarpine-treated TLE rats. We compared variations of seizure onset zone (SOZ), seizure onset pattern (SOP), the latency of seizure onsets, and functional connectivity network from 10-channel LFPs between the early and late stages. Moreover, three machine learning classifiers trained by early-stage data were used to test seizure detection performance in the late stage.

RESULTS

Compared to the early stage, the earliest seizure onset was more frequently detected in hippocampus areas in the late stage. The latency of seizure onsets between electrodes became shorter. Low-voltage fast activity (LVFA) was the most common SOP and the proportion of it increased in the late stage. Different brain states were observed during seizures using Granger causality (GC). Moreover, seizure detection classifiers trained by early-stage data were less accurate when tested in late-stage data.

SIGNIFICANCE

Neuromodulation especially closed-loop deep brain stimulation (DBS) is effective in the treatment of refractory TLE. Although the frequency or amplitude of the stimulation is generally adjusted in existing closed-loop DBS devices in clinical usage, the adjustment rarely considers the pathological progression of chronic TLE. This suggests that an important factor affecting the therapeutic effect of neuromodulation may have been overlooked. The present study reveals time-varying electrophysiological and epileptic network properties in chronic TLE rats and indicates that classifiers of seizure detection and neuromodulation parameters might be designed to adapt to the current state dynamically with the progression of epilepsy.

Unruptured brain arteriovenous malformations causing seizures localize to one common brain network

Seizures are a frequent symptom of unruptured brain arteriovenous malformations (bAVMs). However, the brain regions responsible for these seizures remain unclear. To identify the brain regions causally involved in bAVM-related seizures, we retrospectively reviewed 220 patients with unruptured bAVMs. Using voxel-based lesion-symptom mapping (VLSM) analyses, we tested whether individual brain regions were associated with unruptured bAVM-related seizures. The result revealed that unruptured bAVMs causing seizures are anatomically heterogeneous at the voxel level. Subsequently, lesion network mapping (LNM) analyses was performed to determine whether bAVMs causing seizures belonged to a distributed brain network. LNM analyses indicated that these lesions were located in a functional network characterized by connectivity to the left caudate and precuneus. Moreover, the discrimination performance of the identified seizure network was evaluated in discovery set by calculating the individualized network damage score and was tested in validation set. Based on the calculated network damage scores, patients were divided into low-, medium-, and high-risk groups. The prevalence of seizures significantly differed among the three risk categories in both discovery (p = .003) and validation set (p = .004). Finally, we calculated the percentage of voxels in the canonical resting-state networks that overlapped with the seizure-susceptible brain regions to investigate the involvement of resting-state networks. With an involvement percentage over 50%, the frontoparietal control (82.9%), limbic function (76.7%), and default mode network (69.3%) were considered to be impacted in bAVM-related seizures. Our study identified the seizure-susceptible brain regions for unruptured bAVMs, which could be a plausible neuroimaging biomarker in predicting possible seizures.

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.

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.

A time-reversed model selection approach to time series forecasting

In this paper, we introduce a novel model selection approach to time series forecasting. For linear stationary processes, such as AR processes, the direction of time is independent of the model parameters. By combining theoretical principles of time-reversibility in time series with conventional modeling approaches such as information criteria, we construct a criterion that employs the backwards prediction (backcast) as a proxy for the forecast. Hereby, we aim to adopt a theoretically grounded, data-driven approach to model selection. The novel criterion is named the backwards validated information criterion (BVIC). The BVIC identifies suitable models by trading off a measure of goodness-of-fit and a models ability to predict backwards. We test the performance of the BVIC by conducting experiments on synthetic and real data. In each experiment, the BVIC is examined in contrast to conventionally employed criteria. Our experimental results suggest that the BVIC has comparable performance as conventional information criteria. Specifically, in most of the experiments performed, we did not find statistically significant differences between the forecast error of the BVIC under certain parameterizations and that of the different information criteria. Nonetheless, it is worth emphasizing that the BVIC guarantees are established by design where the model order penalization term depends on strong mathematical properties of time-reversible time series forecasting properties and a finite data assessment. In particular, the penalization term is replaced by a weighted trade-off between functional dimensions pertaining to forecasting.That said, we observed that the BVIC recovered more accurately the real order of the underlying process than the other criteria, which rely on a static penalization of the model order. Lastly, leveraging the latter property we perform the assessment of the order model (or, memory) of time series pertaining to epileptic seizures recorded using electrocorticographic data. Our results provide converging evidence that the order of the model increases during the epileptic events.

Interictal dynamic network transitions in mesial temporal lobe epilepsy

OBJECTIVE

To reveal the possible routine of brain network dynamic alterations in patients with mesial temporal lobe epilepsy (mTLE) and to establish a predicted model of seizure recurrence during interictal periods.

METHODS

Seventy-nine unilateral mTLE patients with hippocampal sclerosis and 97 healthy controls from two centers were retrospectively enrolled. Dynamic brain configuration analyses were performed with resting-state functional magnetic resonance imaging (MRI) data to quantify the functional stability over time and the dynamic interactions between brain regions. Relationships between seizure frequency and ipsilateral hippocampal module allegiance were evaluated using a machine learning predictive model.

RESULTS

Compared to the healthy controls, patients with mTLE displayed an overall higher dynamic network, switching mainly in the epileptogenic regions (false discovery rate [FDR] corrected p-FDR < .05). Moreover, the dynamic network configuration in mTLE was characterized by decreased recruitment (intra-network communication), and increased integration (inter-network communication) among hippocampal systems and large-scale higher-order brain networks (p-FDR < .05). We further found that the dynamic interactions between the hippocampal system and the default-mode network (DMN) or control networks exhibited an opposite distribution pattern (p-FDR < .05). Strikingly, we showed that there was a robust association between predicted seizure frequency based on the ipsilateral hippocampal-DMN dynamics model and actual seizure frequency (p-perm < .001).

SIGNIFICANCE

These findings suggest that the interictal brain of mTLE is characterized by dynamical shifts toward unstable state. Our study provides novel insights into the brain dynamic network alterations and supports the potential use of DMN dynamic parameters as candidate neuroimaging markers in monitoring the seizure frequency clinically during interictal periods.

Graph theoretical measures of fast ripples support the epileptic network hypothesis

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

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

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

Network connectivity predicts effectiveness of responsive neurostimulation in focal epilepsy

Responsive neurostimulation is a promising treatment for drug-resistant focal epilepsy; however, clinical outcomes are highly variable across individuals. The therapeutic mechanism of responsive neurostimulation likely involves modulatory effects on brain networks; however, with no known biomarkers that predict clinical response, patient selection remains empiric. This study aimed to determine whether functional brain connectivity measured non-invasively prior to device implantation predicts clinical response to responsive neurostimulation therapy. Resting-state magnetoencephalography was obtained in 31 participants with subsequent responsive neurostimulation device implantation between 15 August 2014 and 1 October 2020. Functional connectivity was computed across multiple spatial scales (global, hemispheric, and lobar) using pre-implantation magnetoencephalography and normalized to maps of healthy controls. Normalized functional connectivity was investigated as a predictor of clinical response, defined as percent change in self-reported seizure frequency in the most recent year of clinic visits relative to pre-responsive neurostimulation baseline. Area under the receiver operating characteristic curve quantified the performance of functional connectivity in predicting responders (≥50% reduction in seizure frequency) and non-responders (<50%). Leave-one-out cross-validation was furthermore performed to characterize model performance. The relationship between seizure frequency reduction and frequency-specific functional connectivity was further assessed as a continuous measure. Across participants, stimulation was enabled for a median duration of 52.2 (interquartile range, 27.0-62.3) months. Demographics, seizure characteristics, and responsive neurostimulation lead configurations were matched across 22 responders and 9 non-responders. Global functional connectivity in the alpha and beta bands were lower in non-responders as compared with responders (alpha, p_fdr < 0.001; beta, p_fdr < 0.001). The classification of responsive neurostimulation outcome was improved by combining feature inputs; the best model incorporated four features (i.e. mean and dispersion of alpha and beta bands) and yielded an area under the receiver operating characteristic curve of 0.970 (0.919-1.00). The leave-one-out cross-validation analysis of this four-feature model yielded a sensitivity of 86.3%, specificity of 77.8%, positive predictive value of 90.5%, and negative predictive value of 70%. Global functional connectivity in alpha band correlated with seizure frequency reduction (alpha, P = 0.010). Global functional connectivity predicted responder status more strongly, as compared with hemispheric predictors. Lobar functional connectivity was not a predictor. These findings suggest that non-invasive functional connectivity may be a candidate personalized biomarker that has the potential to predict responsive neurostimulation effectiveness and to identify patients most likely to benefit from responsive neurostimulation therapy. Follow-up large-cohort, prospective studies are required to validate this biomarker. These findings furthermore support an emerging view that the therapeutic mechanism of responsive neurostimulation involves network-level effects in the brain.

Effects of anterior thalamic nuclei stimulation on hippocampal activity: Chronic recording in a patient with drug-resistant focal epilepsy

•

Devices for RNS and thalamic DBS were implanted in a single person with epilepsy.

•

RNS electrocorticography enabled characterization of acute and chronic DBS effects.

•

DBS caused acute, phasic, frequency-dependent responses in hippocampus and cortex.

•

DBS modulated functional connectivity and suppressed epileptiform activity over time.

•

Chronic electrocorticography elucidates progressive effects of thalamic stimulation.

Implanted neurostimulation devices are gaining traction as palliative treatment options for certain forms of drug-resistant epilepsy, but clinical utility of these devices is hindered by incomplete mechanistic understanding of their therapeutic effects. Approved devices for anterior thalamic nuclei deep brain stimulation (ANT DBS) are thought to work at a network level, but limited sensing capability precludes characterization of neurophysiological effects outside the thalamus. Here, we describe a patient with drug-resistant temporal lobe epilepsy who was implanted with a responsive neurostimulation device (RNS System), involving hippocampal and ipsilateral temporal neocortical leads, and subsequently received ANT DBS. Over 1.5 years, RNS System electrocorticography enabled multiscale characterization of neurophysiological effects of thalamic stimulation. In brain regions sampled by the RNS System, ANT DBS produced acute, phasic, frequency-dependent responses, including suppression of hippocampal low frequency local field potentials. ANT DBS modulated functional connectivity between hippocampus and neocortex. Finally, ANT DBS progressively suppressed hippocampal epileptiform activity in relation to the extent of hippocampal theta suppression, which informs stimulation parameter selection for ANT DBS. Taken together, this unique clinical scenario, involving hippocampal recordings of unprecedented chronicity alongside ANT DBS, sheds light on the therapeutic mechanism of thalamic stimulation and highlights capabilities needed in next-generation devices.

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.

Ictal gamma-band interactions localize ictogenic nodes of the epileptic network in focal cortical dysplasia

OBJECTIVE

Epilepsy surgery fails in > 30% of patients with focal cortical dysplasia (FCD). The seizure persistence after surgery can be attributed to the inability to precisely localize the tissue with an endogenous potential to generate seizures. In this study, we aimed to identify the critical components of the epileptic network that were actively involved in seizure genesis.

METHODS

The directed transfer function was applied to intracranial EEG recordings and the effective connectivity was determined with a high temporal and frequency resolution. Pre-ictal network properties were compared with ictal epochs to identify regions actively generating ictal activity and discriminate them from the areas of propagation.

RESULTS

Analysis of 276 seizures from 30 patients revealed the existence of a seizure-related network reconfiguration in the gamma-band (25-170 Hz; p < 0.005) - ictogenic nodes. Unlike seizure onset zone, resecting the majority of ictogenic nodes correlated with favorable outcomes (p < 0.012).

CONCLUSION

The prerequisite to successful epilepsy surgery is the accurate identification of brain areas from which seizures arise. We show that in FCD-related epilepsy, gamma-band network markers can reliably identify and distinguish ictogenic areas in macroelectrode recordings, improve intracranial EEG interpretation and better delineate the epileptogenic zone.

SIGNIFICANCE

Ictogenic nodes localize the critical parts of the epileptogenic tissue and increase the diagnostic yield of intracranial evaluation.

Prediction of Seizure Recurrence. A Note of Caution

Great strides have been made recently in documenting that machine-learning programs can predict seizure occurrence in people who have epilepsy. Along with this progress have come claims that appear to us to be a bit premature. We anticipate that many people will benefit from seizure prediction. We also doubt that all will benefit. Although machine learning is a useful tool for aiding discovery, we believe that the greatest progress will come from deeper understanding of seizures, epilepsy, and the EEG features that enable seizure prediction. In this essay, we lay out reasons for optimism and skepticism.

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.

Recursive dynamic functional connectivity reveals a characteristic correlation structure in human scalp EEG

Time-varying neurophysiological activity has been classically explored using correlation based sliding window analysis. However, this method employs only lower order statistics to track dynamic functional connectivity of the brain. We introduce recursive dynamic functional connectivity (rdFC) that incorporates higher order statistics to generate a multi-order connectivity pattern by analyzing neurophysiological data at multiple time scales. The technique builds a hierarchical graph between various temporal scales as opposed to traditional approaches that analyze each scale independently. We examined more than a million rdFC patterns obtained from morphologically diverse EEGs of 2378 subjects of varied age and neurological health. Spatiotemporal evaluation of these patterns revealed three dominant connectivity patterns that represent a universal underlying correlation structure seen across subjects and scalp locations. The three patterns are both mathematically equivalent and observed with equal prevalence in the data. The patterns were observed across a range of distances on the scalp indicating that they represent a spatially scale-invariant correlation structure. Moreover, the number of patterns representing the correlation structure has been shown to be linked with the number of nodes used to generate them. We also show evidence that temporal changes in the rdFC patterns are linked with seizure dynamics.

Static and dynamic functional connectivity supports the configuration of brain networks associated with creative cognition

Creative cognition is recognized to involve the integration of multiple spontaneous cognitive processes and is manifested as complex networks within and between the distributed brain regions. We propose that the processing of creative cognition involves the static and dynamic re-configuration of brain networks associated with complex cognitive processes. We applied the sliding-window approach followed by a community detection algorithm and novel measures of network flexibility on the blood-oxygen level dependent (BOLD) signal of 8 major functional brain networks to reveal static and dynamic alterations in the network reconfiguration during creative cognition using functional magnetic resonance imaging (fMRI). Our results demonstrate the temporal connectivity of the dynamic large-scale creative networks between default mode network (DMN), salience network, and cerebellar network during creative cognition, and advance our understanding of the network neuroscience of creative cognition.

Transcranial direct current stimulation modulates brain functional connectivity in autism

Autism spectrum disorder (ASD) is characterized by deficits in social interactions, impairments in language and communication, and highly restricted behavioral interests. Transcranial direct current stimulation (tDCS) is a widely used form of noninvasive stimulation and may have therapeutic potential for ASD. So far, despite the widespread use of this technique in the neuroscience field, its effects on network-level neural activity and the underlying mechanisms of any effects are still unclear. In the present study, we used electroencephalography (EEG) to investigate tDCS induced brain network changes in children with ASD before and after active and sham stimulation. We recorded 5 min of resting state EEG before and after a single session of tDCS (of approximately 20 min) over dorsolateral prefrontal cortex (DLPFC). Two network-based methods were applied to investigate tDCS modulation on brain networks: 1) temporal network dynamics were analyzed by comparing "flexibility" changes before vs after stimulation, and 2) frequency specific network changes were identified using non-negative matrix factorization (NMF). We found 1) an increase in network flexibility following tDCS (rapid network configuration of dynamic network communities), 2) specific increase in interhemispheric connectivity within the alpha frequency band following tDCS. Together, these results demonstrate that tDCS could help modify both local and global brain network dynamics, and highlight stimulation-induced differences in the manifestation of network reconfiguration. Meanwhile, frequency-specific subnetworks, as a way to index local and global information processing, highlight the core modulatory effects of tDCS on the modular architecture of the functional connectivity patterns within higher frequency bands.

Transfer Function Models for the Localization of Seizure Onset Zone From Cortico-Cortical Evoked Potentials

Surgical resection of the seizure onset zone (SOZ) could potentially lead to seizure-freedom in medically refractory epilepsy patients. However, localizing the SOZ can be a time consuming and tedious process involving visual inspection of intracranial electroencephalographic (iEEG) recordings captured during passive patient monitoring. Cortical stimulation is currently performed on patients undergoing invasive EEG monitoring for the main purpose of mapping functional brain networks such as language and motor networks. We hypothesized that evoked responses from single pulse electrical stimulation (SPES) can also be used to localize the SOZ as they may express the natural frequencies and connectivity of the iEEG network. To test our hypothesis, we constructed patient specific transfer function models from the evoked responses recorded from 22 epilepsy patients that underwent SPES evaluation and iEEG monitoring. We then computed the frequency and connectivity dependent “peak gain” of the system as measured by the H∞  norm from systems theory. We found that in cases for which clinicians had high confidence in localizing the SOZ, the highest peak gain transfer functions with the smallest “floor gain” (gain at which the dipped H∞   3dB below DC gain) corresponded to when the clinically annotated SOZ and early spread regions were stimulated. In more complex cases, there was a large spread of the peak-to-floor (PF) ratios when the clinically annotated SOZ was stimulated. Interestingly for patients who had successful surgeries, our ratio of gains, agreed with clinical localization, no matter the complexity of the case. For patients with failed surgeries, the PF ratio did not match clinical annotations. Our findings suggest that transfer function gains and their corresponding frequency responses computed from SPES evoked responses may improve SOZ localization and thus surgical outcomes.

Semi-automatic Extraction of Functional Dynamic Networks Describing Patient's Epileptic Seizures

Intracranial electroencephalography (EEG) studies using stereotactic EEG (SEEG) have shown that during seizures, epileptic activity spreads across several anatomical regions from the seizure onset zone toward remote brain areas. A full and objective characterization of this patient-specific time-varying network is crucial for optimal surgical treatment. Functional connectivity (FC) analysis of SEEG signals recorded during seizures enables to describe the statistical relations between all pairs of recorded signals. However, extracting meaningful information from those large datasets is time consuming and requires high expertise. In the present study, we first propose a novel method named Brain-wide Time-varying Network Decomposition (BTND) to characterize the dynamic epileptogenic networks activated during seizures in individual patients recorded with SEEG electrodes. The method provides a number of pathological FC subgraphs with their temporal course of activation. The method can be applied to several seizures of the patient to extract reproducible subgraphs. Second, we compare the activated subgraphs obtained by the BTND method with visual interpretation of SEEG signals recorded in 27 seizures from nine different patients. As a whole, we found that activated subgraphs corresponded to brain regions involved during the course of the seizures and their time course was highly consistent with classical visual interpretation. We believe that the proposed method can complement the visual analysis of SEEG signals recorded during seizures by highlighting and characterizing the most significant parts of epileptic networks with their activation dynamics.

Large-Scale Desynchronization During Interictal Epileptic Discharges Recorded With Intracranial EEG

It is increasingly recognized that deep understanding of epileptic seizures requires both localizing and characterizing the functional network of the region where they are initiated, i. e., the epileptic focus. Previous investigations of the epileptogenic focus' functional connectivity have yielded contrasting results, reporting both pathological increases and decreases during resting periods and seizures. In this study, we shifted paradigm to investigate the time course of connectivity in relation to interictal epileptiform discharges. We recruited 35 epileptic patients undergoing intracranial EEG (iEEG) investigation as part of their presurgical evaluation. For each patient, 50 interictal epileptic discharges (IEDs) were marked and iEEG signals were epoched around those markers. Signals were narrow-band filtered and time resolved phase-locking values were computed to track the dynamics of functional connectivity during IEDs. Results show that IEDs are associated with a transient decrease in global functional connectivity, time-locked to the peak of the discharge and specific to the high range of the gamma frequency band. Disruption of the long-range connectivity between the epileptic focus and other brain areas might be an important process for the generation of epileptic activity. Transient desynchronization could be a potential biomarker of the epileptogenic focus since 1) the functional connectivity involving the focus decreases significantly more than the connectivity outside the focus and 2) patients with good surgical outcome appear to have a significantly more disconnected focus than patients with bad outcomes.

Reconfiguration of human evolving large-scale epileptic brain networks prior to seizures: an evaluation with node centralities

Previous research has indicated that temporal changes of centrality of specific nodes in human evolving large-scale epileptic brain networks carry information predictive of impending seizures. Centrality is a fundamental network-theoretical concept that allows one to assess the role a node plays in a network. This concept allows for various interpretations, which is reflected in a number of centrality indices. Here we aim to achieve a more general understanding of local and global network reconfigurations during the pre-seizure period as indicated by changes of different node centrality indices. To this end, we investigate—in a time-resolved manner—evolving large-scale epileptic brain networks that we derived from multi-day, multi-electrode intracranial electroencephalograpic recordings from a large but inhomogeneous group of subjects with pharmacoresistant epilepsies with different anatomical origins. We estimate multiple centrality indices to assess the various roles the nodes play while the networks transit from the seizure-free to the pre-seizure period. Our findings allow us to formulate several major scenarios for the reconfiguration of an evolving epileptic brain network prior to seizures, which indicate that there is likely not a single network mechanism underlying seizure generation. Rather, local and global aspects of the pre-seizure network reconfiguration affect virtually all network constituents, from the various brain regions to the functional connections between them.

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

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

Temporal Lobe Epilepsy Surgical Outcomes Can Be Inferred Based on Structural Connectome Hubs: A Machine Learning Study

OBJECTIVE

Medial temporal lobe epilepsy (TLE) is the most common form of medication-resistant focal epilepsy in adults. Despite removal of medial temporal structures, more than one-third of patients continue to have disabling seizures postoperatively. Seizure refractoriness implies that extramedial regions are capable of influencing the brain network and generating seizures. We tested whether abnormalities of structural network integration could be associated with surgical outcomes.

METHODS

Presurgical magnetic resonance images from 121 patients with drug-resistant TLE across 3 independent epilepsy centers were used to train feed-forward neural network models based on tissue volume or graph-theory measures from whole-brain diffusion tensor imaging structural connectomes. An independent dataset of 47 patients with TLE from 3 other epilepsy centers was used to assess the predictive values of each model and regional anatomical contributions toward surgical treatment results.

RESULTS

The receiver operating characteristic area under the curve based on regional betweenness centrality was 0.88, significantly higher than a random model or models based on gray matter volumes, degree, strength, and clustering coefficient. Nodes most strongly contributing to the predictive models involved the bilateral parahippocampal gyri, as well as the superior temporal gyri.

INTERPRETATION

Network integration in the medial and lateral temporal regions was related to surgical outcomes. Patients with abnormally integrated structural network nodes were less likely to achieve seizure freedom. These findings are in line with previous observations related to network abnormalities in TLE and expand on the notion of underlying aberrant plasticity. Our findings provide additional information on the mechanisms of surgical refractoriness. ANN NEUROL 2020;88:970-983.

Electrical stimulation of the endopiriform nucleus attenuates epilepsy in rats by network modulation

OBJECTIVE

Neuromodulatory anterior thalamic deep brain stimulation (DBS) is an effective therapy for intractable epilepsy, but few patients achieve complete seizure control with thalamic DBS. Other stimulation sites may be considered for anti-seizure DBS. We investigated bilateral low-frequency stimulation of the endopiriform nuclei (LFS-EPN) to control seizures induced by intracortically implanted cobalt wire in rats.

METHODS

Chronic epilepsy was induced by cobalt wire implantation in the motor cortex unilaterally. Bipolar-stimulating electrodes were implanted into the EPN bilaterally. Continuous electroencephalography (EEG) was recorded using electrodes placed into bilateral motor cortex and hippocampus CA1 areas. Spontaneous seizures were monitored by long-term video-EEG, and behavioral seizures were classified based on the Racine scale. Continuous 1-Hz LFS-EPN began on the third day after electrode implantation and was controlled by a multi-channel stimulator. Stimulation continued until the rats had no seizures for three consecutive days.

RESULTS

Compared with the control and sham stimulation groups, the LFS-EPN group experienced significantly fewer seizures per day and the mean Racine score of seizures was lower due to fewer generalized seizures. Ictal discharges at the epileptogenic site had significantly reduced theta band power in the LFS-EPN group compared to the other groups.

INTERPRETATION

Bilateral LFS-EPN attenuates cobalt wire-induced seizures in rats by modulating epileptic networks. Reduced ictal theta power of the EEG broadband spectrum at the lesion site may be associated with the anti-epileptogenic mechanism of LFS-EPN. Bilateral EPN DBS may have therapeutic applications in human partial epilepsies.

Edge-centric functional network representations of human cerebral cortex reveal overlapping system-level architecture

Network neuroscience has relied on a node-centric network model in which cells, populations and regions are linked to one another via anatomical or functional connections. This model cannot account for interactions of edges with one another. In this study, we developed an edge-centric network model that generates constructs 'edge time series' and 'edge functional connectivity' (eFC). Using network analysis, we show that, at rest, eFC is consistent across datasets and reproducible within the same individual over multiple scan sessions. We demonstrate that clustering eFC yields communities of edges that naturally divide the brain into overlapping clusters, with regions in sensorimotor and attentional networks exhibiting the greatest levels of overlap. We show that eFC is systematically modulated by variation in sensory input. In future work, the edge-centric approach could be useful for identifying novel biomarkers of disease, characterizing individual variation and mapping the architecture of highly resolved neural circuits.

Signal diffusion along connectome gradients and inter-hub routing differentially contribute to dynamic human brain function

Human cognition is dynamic, alternating over time between externally-focused states and more abstract, often self-generated, patterns of thought. Although cognitive neuroscience has documented how networks anchor particular modes of brain function, mechanisms that describe transitions between distinct functional states remain poorly understood. Here, we examined how time-varying changes in brain function emerge within the constraints imposed by macroscale structural network organization. Studying a large cohort of healthy adults (n = 326), we capitalized on manifold learning techniques that identify low dimensional representations of structural connectome organization and we decomposed neurophysiological activity into distinct functional states and their transition patterns using Hidden Markov Models. Structural connectome organization predicted dynamic transitions anchored in sensorimotor systems and those between sensorimotor and transmodal states. Connectome topology analyses revealed that transitions involving sensorimotor states traversed short and intermediary distances and adhered strongly to communication mechanisms of network diffusion. Conversely, transitions between transmodal states involved spatially distributed hubs and increasingly engaged long-range routing. These findings establish that the structure of the cortex is optimized to allow neural states the freedom to vary between distinct modes of processing, and so provides a key insight into the neural mechanisms that give rise to the flexibility of human cognition.

Early seizure spread and epilepsy surgery: A systematic review

OBJECTIVE

A fundamental question in epilepsy surgery is how to delineate the margins of cortex that must be resected to result in seizure freedom. Whether and which areas showing seizure activity early in ictus must be removed to avoid postoperative recurrence of seizures is an area of ongoing research. Seizure spread dynamics in the initial seconds of ictus are often correlated with postoperative outcome; there is neither a consensus definition of early spread nor a concise summary of the existing literature linking seizure spread to postsurgical seizure outcomes. The present study is intended to summarize the literature that links seizure spread to postoperative seizure outcome and to provide a framework for quantitative assessment of early seizure spread.

METHODS

A systematic review was carried out according to Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines. A Medline search identified clinical studies reporting data on seizure spread measured by intracranial electrodes, having at least 10 subjects and reporting at least 1-year postoperative outcome in the English literature from 1990 to 2019. Studies were evaluated regarding support for a primary hypothesis: Areas of early seizure spread represent cortex with seizure-generating potential.

RESULTS

The search yielded 4562 studies: 15 studies met inclusion criteria and 7 studies supported the primary hypothesis. The methods and metrics used to describe seizure spread were heterogenous. The timeframe of seizure spread associated with seizure outcome ranged from 1-14 seconds, with large, well-designed, retrospective studies pointing to 3-10 seconds as most likely to provide meaningful correlates of postoperative seizure freedom.

SIGNIFICANCE

The complex correlation between electrophysiologic seizure spread and the potential for seizure generation needs further elucidation. Prospective cohort studies or trials are needed to evaluate epilepsy surgery targeting cortex involved in the first 3-10 seconds of ictus.

Recurrence quantification analysis of dynamic brain networks

Evidence suggests that brain network dynamics are a key determinant of brain function and dysfunction. Here we propose a new framework to assess the dynamics of brain networks based on recurrence analysis. Our framework uses recurrence plots and recurrence quantification analysis to characterize dynamic networks. For resting-state magnetoencephalographic dynamic functional networks (dFNs), we have found that functional networks recur more quickly in people with epilepsy than in healthy controls. This suggests that recurrence of dFNs may be used as a biomarker of epilepsy. For stereo electroencephalography data, we have found that dFNs involved in epileptic seizures emerge before seizure onset, and recurrence analysis allows us to detect seizures. We further observe distinct dFNs before and after seizures, which may inform neurostimulation strategies to prevent seizures. Our framework can also be used for understanding dFNs in healthy brain function and in other neurological disorders besides epilepsy.

Identification of neural networks preferentially engaged by epileptogenic mass lesions through lesion network mapping analysis

Lesion network mapping (LNM) has been applied to true lesions (e.g., cerebrovascular lesions in stroke) to identify functionally connected brain networks. No previous studies have utilized LNM for analysis of intra-axial mass lesions. Here, we implemented LNM for identification of potentially vulnerable epileptogenic networks in mass lesions causing medically-refractory epilepsy (MRE). Intra-axial brain lesions were manually segmented in patients with MRE seen at our institution (EL_INST). These lesions were then normalized to standard space and used as seeds in a high-resolution normative resting state functional magnetic resonance imaging template. The resulting connectivity maps were first thresholded (p_{Bonferroni_cor} < 0.05) and binarized; the thresholded binarized connectivity maps were subsequently summed to produce overall group connectivity maps, which were compared with established resting-state networks to identify potential networks prone to epileptogenicity. To validate our data, this approach was also applied to an external dataset of epileptogenic lesions identified from the literature (EL_LIT). As an additional exploratory analysis, we also segmented and computed the connectivity of institutional non-epileptogenic lesions (NEL_INST), calculating voxel-wise odds ratios (VORs) to identify voxels more likely to be functionally-connected with EL_INST versus NEL_INST. To ensure connectivity results were not driven by anatomical overlap, the extent of lesion overlap between EL_INST, and EL_LIT and NEL_INST was assessed using the Dice Similarity Coefficient (DSC, lower index ~ less overlap). Twenty-eight patients from our institution were included (EL_INST: 17 patients, 17 lesions, 10 low-grade glioma, 3 cavernoma, 4 focal cortical dysplasia; NEL_INST: 11 patients, 33 lesions, all brain metastases). An additional 23 cases (25 lesions) with similar characteristics to the EL_INST data were identified from the literature (EL_LIT). Despite minimal anatomical overlap of lesions, both EL_INST and EL_LIT showed greatest functional connectivity overlap with structures in the Default Mode Network, Frontoparietal Network, Ventral Attention Network, and the Limbic Network-with percentage volume overlap of 19.5%, 19.1%, 19.1%, and 12.5%, respectively-suggesting them as networks consistently engaged by epileptogenic mass lesions. Our exploratory analysis moreover showed that the mesial frontal lobes, parahippocampal gyrus, and lateral temporal neocortex were at least twice as likely to be functionally connected with the EL_INST compared to the NEL_INST group (i.e. Peak VOR > 2.0); canonical resting-state networks preferentially engaged by EL_INSTs were the Limbic and the Frontoparietal Networks (Mean VOR > 1.5). In this proof of concept study, we demonstrate the feasibility of LNM for intra-axial mass lesions by showing that ELs have discrete functional connections and may preferentially engage in discrete resting-state networks. Thus, the underlying normative neural circuitry may, in part, explain the propensity of particular lesions toward the development of MRE. If prospectively validated, this has ramifications for patient counseling along with both approach and timing of surgery for lesions in locations prone to development of MRE.

Time-variant Epileptic Brain Functional Connectivity of Focal and Generalized Seizure in Chronic Temporal Lobe Epilepsy Rat. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2020;2020:2833-2836. [PMID: 33018596 DOI: 10.1109/embc44109.2020.9175924]

Seizure types and characteristics may vary with time in a patient with distinct mechanisms underlying the propagation of ictal activity. Similarly, we found that both focal and generalized seizures coexist in some pilocarpine-induced chronic temporal lobe epilepsy (TLE) rats. In different seizure patterns, mapping complex networks and analyzing epileptic characteristics involved in seizure propagation are likely to reflect seizure propagation mechanisms, and indicate the establishment of stimulation strategy for epilepsy treatment, especially on the selection of stimulation targets. In our study, we used Granger causality method to track the time-variant epileptic brain functional connectivity in focal and generalized seizures from multi-site local field potentials (LFPs). Results showed that these two major types of seizures had different propagation patterns during ictal period. When comparing them, generalized seizures involved in a network with more complex relationships and spread to more extensive brain regions than in local seizures at mid-ictal stage. Moreover, we observed that focal seizures had a focused causal hub with strong interactions, while generalized seizures had relative distributed causal hubs to drive the development of seizure during seizure-onset stage. These findings suggest that stimulation strategy might need to be adapted to different seizure types thus allowing for retuning abnormal epileptic brain network and obtaining better treatment effect on seizure suppression.