Individual brain structure and modelling predict seizure propagation

Virtual epilepsy patient cohort: Generation and evaluation

Epilepsy is a prevalent brain disorder, characterized by sudden, abnormal brain activity, making it difficult to live with. One-third of people with epilepsy do not respond to anti-epileptic drugs. Drug-resistant epilepsy is treated with brain surgery. Successful surgical treatment relies on identifying brain regions responsible for seizure onset, known as epileptogenic zones (EZ). Despite various methods for EZ estimation, evaluating their efficacy remains challenging due to a lack of ground truth for empirical data. To address this, we generated and evaluated a cohort of 30 virtual epilepsy patients, using patient-specific anatomical and functional data from 30 real drug-resistant epilepsy patients. This personalized modeling approach, based on each patient's brain data, is called a virtual brain twin. For each virtual patient, we provided data that included anatomically parcellated brain regions, structural connectivity, reconstructed intracranial electrodes, simulated brain activity at both the brain region and electrode levels, and key parameters of the virtual brain twin. These key parameters, which include the EZ hypothesis, serve as the ground truth for simulated brain activity. For each virtual brain twin, we generated synthetic spontaneous seizures, stimulation-induced seizures and interictal activity. We systematically evaluated these simulated brain signals by quantitatively comparing them against their corresponding empirical intracranial recordings. Simulated signals based on patient-specific EZ captured spatio-temporal seizure generation and propagation. Through in-silico exploration of stimulation parameters, we also demonstrated the role of patient-specific stimulation location and amplitude in reproducing empirically stimulated seizures. The virtual epileptic cohort is openly available, and can be used to systematically evaluate methods for the estimation of EZ or source localization using ground truth EZ parameters and source signals.

Pioneering new paths: the role of generative modelling in neurological disease research

Recently, deep generative modelling has become an increasingly powerful tool with seminal work in a myriad of disciplines. This powerful modelling approach is supposed to not only have the potential to solve current problems in the medical field but also to enable personalised precision medicine and revolutionise healthcare through applications such as digital twins of patients. Here, the core concepts of generative modelling and popular modelling approaches are first introduced to consider the potential based on methodological concepts for the generation of synthetic data and the ability to learn a representation of observed data. These potentials will be reviewed using current applications in neuroimaging for data synthesis and disease decomposition in Alzheimer's disease and multiple sclerosis. Finally, challenges for further research and applications will be discussed, including computational and data requirements, model evaluation, and potential privacy risks.

Voxel-Based Lesion Symptom Mapping to Predict Poststroke Epilepsy After Mechanical Thrombectomy

INTRODUCTION

Poststroke epilepsy (PSE) is an important long-term complication after stroke. Data regarding predictors of PSE in patients with large-vessel occlusion stroke receiving mechanical thrombectomy (MT) are scarce. Voxel-based lesion symptom mapping on brain MRI might be a valuable tool in the risk prediction of PSE. This study aims to assess PSE risk after acute stroke treated with MT via voxel- and volumetric-based analyses.

METHODS

In this bi-center study from two tertiary-care stroke centers, we included consecutive acute ischemic stroke patients who had received MT between 2011 and 2017, and had postinterventional brain MRI as well as long-term follow-up data available. Infarct volume and location were assessed on MRI. Following semiautomated lesion outlining and generation of binarized lesion masks, lesion symptom mapping was applied to identify relevant topographical lesion patterns in PSE.

RESULTS

Of 348 analyzed patients, 97 cases had to be excluded due to insufficient image quality and inaccurate registration results. Finally, lesion maps from 251 patients (median age: 66, 45.4% women) were considered for lesion symptom mapping, including maps from 26 patients with PSE (10.4%). Mean infarct volume was higher in PSE patients (119.2 cm3 vs. 43.9 cm3, p < 0.0001). Lesion symptom mapping identified the orbitofrontal region, the operculum, and the temporal pole as brain regions associated with PSE.

CONCLUSION

Apart from infarct volume, lesion symptom mapping on postinterventional brain MRI identified specific brain regions associated with PSE after large vessel occlusion stroke. This information might be helpful for PSE risk stratification and follow-up care in this specific population.

The different subtypes of temporal lobe seizures networks

Temporal lobe epilepsies (TLEs) are among the forms of epilepsy most frequently encountered in surgical evaluations, characterized by a wide range of anatomical, functional, and electroclinical subtypes. Traditional classifications, such as mesial and lateral TLE, have been broadened by advances in stereoelectroencephalography (SEEG), revealing more complex forms such as mesio-lateral and temporal-plus seizures. These findings support the concept of epileptogenic networks, emphasizing interconnected regions rather than isolated focal areas in the genesis of seizures. Quantitative tools, such as the epileptogenicity index (EI), are improving the accuracy of SEEG interpretation, which is closely correlated with surgical results. Temporal-plus epilepsies, in particular, require full SEEG exploration due to their broader involvement in the network, necessitating tailored surgical approaches. A better understanding of TLEs subtypes and epileptogenic networks is an essential basis for advancing minimally invasive surgical techniques, including laser interstitial thermal therapy (LITT) and neuromodulation. These methods rely on the precise localization of epileptogenic networks. This network-based framework represents an important step towards optimizing surgical outcomes and advancing personalized epilepsy care.

Localization of the epileptogenic network from scalp EEG using a patient-specific whole-brain model

Computational modeling is a key tool for elucidating the neuronal mechanisms underlying epileptic activity. Despite considerable progress, existing models often lack realistic accuracy in representing electrophysiological epileptic activity. In this study, we used a comprehensive human brain model based on a neural mass model, which is tailored to the layered structure of the neocortex and incorporates patient-specific imaging data. This approach allowed the simulation of scalp EEGs in an epileptic patient suffering from type 2 focal cortical dysplasia (FCD). The simulation specifically addressed epileptic activity induced by FCD, faithfully reproducing intracranial interictal epileptiform discharges (IEDs) recorded with electrocorticography. For constructing the patient-specific scalp EEG, we carefully defined a clear delineation of the epileptogenic zone by numerical simulations to ensure fidelity to the topography, polarity, and diffusion characteristics of IEDs. This nuanced approach improves the accuracy of the simulated EEG signal, provides a more accurate representation of epileptic activity, and enhances our understanding of the mechanism behind the epileptogenic networks. The accuracy of the model was confirmed by a postoperative reevaluation with a secondary EEG simulation that was consistent with the lesion's removal. Ultimately, this personalized approach may prove instrumental in optimizing and tailoring epilepsy treatment strategies.

EEG signature orchestrating expression of ictal behavior in mesial temporal lobe epilepsy

OBJECTIVES

We investigated EEG features differentiating clinical seizures (CSs) from subclinical seizures (SCSs) to explore the mechanisms underlying the generation of ictal behavior in mesial temporal lobe epilepsy (mTLE).

METHODS

Peri-ictal state of power spectral density (PSD) within seizure onset zone (SOZ) and propagation zone (PZ) were compared between SCSs and CSs. Functional connectivity was analyzed using the nonlinear correlation coefficient h2, outgoing links (OUT) and ingoing links (IN). The EEG epochs of CSs-early part and SCSs were equally divided into four segments to reveal dynamic EEG changes.

RESULTS

During pre-ictal state, PSD at 30-80 Hz in SOZ was higher in CSs than in SCSs. The preictal OUT and IN values in SOZ at 30-80 Hz were greater in CSs than in SCSs. During CSs-early part, PSD displayed an initial increase in SOZ but a late increase in PZ, with enhanced high-frequency activity in temporal regions and increased low-frequency activity in insula.

CONCLUSION

The enhanced pre-ictal gamma activity within the epileptic network was able to distinguish CSs from SCSs. The unique temporospatial alterations within the epileptic network drive the expression of ictal behavior in mTLE.

SIGNIFICANCE

The distinct EEG features between SCSs and CSs offer transformative insights into the mechanisms driving ictal behavior.

Precision neurology

Over the past several decades, high-resolution brain imaging, blood and cerebrospinal fluid analyses, and other advanced technologies have changed diagnosis from an exercise depending primarily on the history and physical examination to a computer- and online resource-aided process that relies on larger and larger quantities of data. In addition, randomized controlled trials (RCT) at a population level have led to many new drugs and devices to treat neurological disease, including disease-modifying therapies. We are now at a crossroads. Combinatorially profound increases in data about individuals has led to an alternative to population-based RCTs. Genotyping and comprehensive "deep" phenotyping can sort individuals into smaller groups, enabling precise medical decisions at a personal level. In neurology, precision medicine that includes prediction, prevention and personalization requires that genomic and phenomic information further incorporate imaging and behavioral data. In this article, we review the genomic, phenomic, and computational aspects of precision medicine for neurology. After defining biological markers, we discuss some applications of these "-omic" and neuroimaging measures, and then outline the role of computation and ultimately brain simulation. We conclude the article with a discussion of the relation between precision medicine and value-based care.

Mechanistic insight of curcumin: a potential pharmacological candidate for epilepsy

Recurrent spontaneous seizures with an extended epileptic discharge are the hallmarks of epilepsy. At present, there are several available anti-epileptic drugs (AEDs) in the market. Still no adequate treatment for epilepsy treatment is available. The main disadvantages of AEDs are their associated adverse effects. It is a challenge to develop new therapies that can reduce seizures by modulating the underlying mechanisms with no adverse effects. In the last decade, the neuromodulatory potential of phytoconstituents has sparked their usage in the treatment of central nervous system disorders. Curcumin is an active polyphenolic component that interacts at cellular and molecular levels. Curcumin's neuroprotective properties have been discovered in recent preclinical and clinical studies due to its immunomodulatory effects. Curcumin has the propensity to modulate signaling pathways involved in cell survival and manage oxidative stress, apoptosis, and inflammatory mechanisms. Further, curcumin can persuade epigenetic alterations, including histone modifications (acetylation/deacetylation), which are the changes responsible for the altered expression of genes facilitating the process of epileptogenesis. The bioavailability of curcumin in the brain is a concern that needs to be tackled. Therefore, nanonization has emerged as a novel drug delivery system to enhance the pharmacokinetics of curcumin. In the present review, we reviewed curcumin's modulatory effects on potential biomarkers involved in epileptogenesis including dendritic cells, T cell subsets, cytokines, chemokines, apoptosis mediators, antioxidant mechanisms, and cognition impairment. Also, we have discussed the nanocarrier systems for encapsulating curcumin, offering a promising approach to enhance bioavailability of curcumin.

Advances in physiological and clinical relevance of hiPSC-derived brain models for precision medicine pipelines

Precision, or personalized, medicine aims to stratify patients based on variable pathogenic signatures to optimize the effectiveness of disease prevention and treatment. This approach is favorable in the context of brain disorders, which are often heterogeneous in their pathophysiological features, patterns of disease progression and treatment response, resulting in limited therapeutic standard-of-care. Here we highlight the transformative role that human induced pluripotent stem cell (hiPSC)-derived neural models are poised to play in advancing precision medicine for brain disorders, particularly emerging innovations that improve the relevance of hiPSC models to human physiology. hiPSCs derived from accessible patient somatic cells can produce various neural cell types and tissues; current efforts to increase the complexity of these models, incorporating region-specific neural tissues and non-neural cell types of the brain microenvironment, are providing increasingly relevant insights into human-specific neurobiology. Continued advances in tissue engineering combined with innovations in genomics, high-throughput screening and imaging strengthen the physiological relevance of hiPSC models and thus their ability to uncover disease mechanisms, therapeutic vulnerabilities, and tissue and fluid-based biomarkers that will have real impact on neurological disease treatment. True physiological understanding, however, necessitates integration of hiPSC-neural models with patient biophysical data, including quantitative neuroimaging representations. We discuss recent innovations in cellular neuroscience that can provide these direct connections through generative AI modeling. Our focus is to highlight the great potential of synergy between these emerging innovations to pave the way for personalized medicine becoming a viable option for patients suffering from neuropathologies, particularly rare epileptic and neurodegenerative disorders.

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.

Computational neurosciences and quantitative systems pharmacology: a powerful combination for supporting drug development in neurodegenerative diseases

Successful clinical development of new therapeutic interventions is notoriously difficult, especially in neurodegenerative diseases, where predictive biomarkers are scarce and functional improvement is often based on patient's perception, captured by structured interviews. As a consequence, mechanistic modeling of the processes relevant to therapeutic interventions in CNS disorders has been lagging behind other disease indications, probably because of the perceived complexity of the brain. However in this report, we develop the argument that a combination of Computational Neurosciences and Quantitative Systems Pharmacology (QSP) modeling of molecular pathways is a powerful simulation tool to enhance the probability of successful drug development for neurodegenerative diseases. Computational Neurosciences aims to predict action potential dynamics and neuronal circuit activation that are ultimately linked to behavioral changes and clinically relevant functional outcomes. These processes can not only be affected by the disease state, but also by common genotype variants on neurotransmitter-related proteins and the psycho-active medications often prescribed in these patient populations. Quantitative Systems Pharmacology (QSP) modeling of molecular pathways allows to simulate key pathological drivers of dementia, such as protein aggregation and neuroinflammatory responses. They often impact neurotransmitter homeostasis and voltage-gated ion-channels or lead to mitochondrial dysfunction, ultimately leading to changes in action potential dynamics and clinical readouts. Combining these two modeling approaches can lead to better actionable understanding of the many non-linear pharmacodynamic processes active in the human diseased brain. Practical applications include a rational selection of the optimal doses in combination therapies, identification of subjects more likely to respond to treatment, a more balanced stratification of treatment arms in terms of comedications, disease status and common genotype variants and re-analysis of small clinical trials to uncover a possible clinical signal. Ultimately this will lead to a higher success rate of bringing new therapeutics to the right patient populations.

Transient alteration of Awareness triggered by direct electrical stimulation of the brain

BACKGROUND

Awareness is a state of consciousness that enables a subject to interact with the environment. Transient alteration of awareness (AA) is a disabling sign of many types of epileptic seizures. The brain mechanisms of awareness and its alteration are not well known.

OBJECTIVE/HYPOTHESIS

Transient and isolated AA induced by electrical brain stimulation during a stereoelectroencephalography (SEEG) recording represents an ideal model for studying the associated modifications of functional connectivity and locating the hubs of awareness networks.

METHODS

We investigated the SEEG signals-based brain functional connectivity (FC) changes vs background occurring during AA triggered by three thalamic and two insular stimulations in three patients explored by SEEG in the frame of presurgical evaluation for focal drug-resistant epilepsy. The results were compared to the stimulations of the same sites that did not induce clinical changes (negative stimulations).

RESULTS

We observed decreased node strength in the pulvinar, insula, and parietal associative cortices during the thalamic and insular stimulations that induced AA. The link strengths characterizing functional coupling between the thalamus and the insular, prefrontal, temporal, or parietal associative cortices were also decreased. In contrast, there was an increased synchronization between the precuneus and the temporal lateral cortex. These FC changes were absent during the negative stimulations.

CONCLUSION

Our study highlights the role of the pulvinar, insular, and parietal hubs in maintaining the awareness networks and paves the way for invasive or non-invasive neuromodulation protocols to reduce AA manifestations during epileptic seizures.

Predicting radiofrequency thermocoagulation surgical outcomes in refractory focal epilepsy patients using functional coupled neural mass model

Objective

The success rate of achieving seizure freedom after radiofrequency thermocoagulation surgery for patients with refractory focal epilepsy is about 20-40%. This study aims to enhance the prediction of surgical outcomes based on preoperative decisions through network model simulation, providing a reference for clinicians to validate and optimize surgical plans.

Methods

Twelve patients with epilepsy who underwent radiofrequency thermocoagulation were retrospectively reviewed in this study. A coupled model based on model subsets of the neural mass model was constructed by calculating partial directed coherence as the coupling matrix from stereoelectroencephalography (SEEG) signals. Multi-channel time-varying model parameters of excitation and inhibitions were identified by fitting the real SEEG signals with the coupled model. Further incorporating these model parameters, the coupled model virtually removed contacts destroyed in radiofrequency thermocoagulation or selected randomly. Subsequently, the coupled model after virtual surgery was simulated.

Results

The identified excitatory and inhibitory parameters showed significant difference before and after seizure onset (p < 0.05), and the trends of parameter changes aligned with the seizure process. Additionally, excitatory parameters of epileptogenic contacts were higher than that of non-epileptogenic contacts, and opposite findings were noticed for inhibitory parameters. The simulated signals of postoperative models to predict surgical outcomes yielded an area under the curve (AUC) of 83.33% and an accuracy of 91.67%.

Conclusion

The multi-channel coupled model proposed in this study with physiological characteristics showed a desirable performance for preoperatively predicting patients' prognoses.

Recent Progress in Brain Network Models for Medical Applications: A Review

Importance: Pathological perturbations of the brain often spread via connectome to fundamentally alter functional consequences. By integrating multimodal neuroimaging data with mathematical neural mass modeling, brain network models (BNMs) enable to quantitatively characterize aberrant network dynamics underlying multiple neurological and psychiatric disorders. We delved into the advancements of BNM-based medical applications, discussed the prevalent challenges within this field, and provided possible solutions and future directions. Highlights: This paper reviewed the theoretical foundations and current medical applications of computational BNMs. Composed of neural mass models, the BNM framework allows to investigate large-scale brain dynamics behind brain diseases by linking the simulated functional signals to the empirical neurophysiological data, and has shown promise in exploring neuropathological mechanisms, elucidating therapeutic effects, and predicting disease outcome. Despite that several limitations existed, one promising trend of this research field is to precisely guide clinical neuromodulation treatment based on individual BNM simulation. Conclusion: BNM carries the potential to help understand the mechanism underlying how neuropathology affects brain network dynamics, further contributing to decision-making in clinical diagnosis and treatment. Several constraints must be addressed and surmounted to pave the way for its utilization in the clinic.

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.

Parameter estimation in a whole-brain network model of epilepsy: Comparison of parallel global optimization solvers

The Virtual Epileptic Patient (VEP) refers to a computer-based representation of a patient with epilepsy that combines personalized anatomical data with dynamical models of abnormal brain activities. It is capable of generating spatio-temporal seizure patterns that resemble those recorded with invasive methods such as stereoelectro EEG data, allowing for the evaluation of clinical hypotheses before planning surgery. This study highlights the effectiveness of calibrating VEP models using a global optimization approach. The approach utilizes SaCeSS, a cooperative metaheuristic algorithm capable of parallel computation, to yield high-quality solutions without requiring excessive computational time. Through extensive benchmarking on synthetic data, our proposal successfully solved a set of different configurations of VEP models, demonstrating better scalability and superior performance against other parallel solvers. These results were further enhanced using a Bayesian optimization framework for hyperparameter tuning, with significant gains in terms of both accuracy and computational cost. Additionally, we added a scalable uncertainty quantification phase after model calibration, and used it to assess the variability in estimated parameters across different problems. Overall, this study has the potential to improve the estimation of pathological brain areas in drug-resistant epilepsy, thereby to inform the clinical decision-making process.

Individualized epidemic spreading models predict epilepsy surgery outcomes: A pseudo-prospective study

Epilepsy surgery is the treatment of choice for drug-resistant epilepsy patients, but up to 50% of patients continue to have seizures one year after the resection. In order to aid presurgical planning and predict postsurgical outcome on a patient-by-patient basis, we developed a framework of individualized computational models that combines epidemic spreading with patient-specific connectivity and epileptogeneity maps: the Epidemic Spreading Seizure and Epilepsy Surgery framework (ESSES). ESSES parameters were fitted in a retrospective study (N = 15) to reproduce invasive electroencephalography (iEEG)-recorded seizures. ESSES reproduced the iEEG-recorded seizures, and significantly better so for patients with good (seizure-free, SF) than bad (nonseizure-free, NSF) outcome. We illustrate here the clinical applicability of ESSES with a pseudo-prospective study (N = 34) with a blind setting (to the resection strategy and surgical outcome) that emulated presurgical conditions. By setting the model parameters in the retrospective study, ESSES could be applied also to patients without iEEG data. ESSES could predict the chances of good outcome after any resection by finding patient-specific model-based optimal resection strategies, which we found to be smaller for SF than NSF patients, suggesting an intrinsic difference in the network organization or presurgical evaluation results of NSF patients. The actual surgical plan overlapped more with the model-based optimal resection, and had a larger effect in decreasing modeled seizure propagation, for SF patients than for NSF patients. Overall, ESSES could correctly predict 75% of NSF and 80.8% of SF cases pseudo-prospectively. Our results show that individualised computational models may inform surgical planning by suggesting alternative resections and providing information on the likelihood of a good outcome after a proposed resection. This is the first time that such a model is validated with a fully independent cohort and without the need for iEEG recordings.

Treatment effects in epilepsy: a mathematical framework for understanding response over time

Epilepsy is a neurological disorder characterized by recurrent seizures, affecting over 65 million people worldwide. Treatment typically commences with the use of anti-seizure medications, including both mono- and poly-therapy. Should these fail, more invasive therapies such as surgery, electrical stimulation and focal drug delivery are often considered in an attempt to render the person seizure free. Although a significant portion ultimately benefit from these treatment options, treatment responses often fluctuate over time. The physiological mechanisms underlying these temporal variations are poorly understood, making prognosis a significant challenge when treating epilepsy. Here we use a dynamic network model of seizure transition to understand how seizure propensity may vary over time as a consequence of changes in excitability. Through computer simulations, we explore the relationship between the impact of treatment on dynamic network properties and their vulnerability over time that permit a return to states of high seizure propensity. For small networks we show vulnerability can be fully characterised by the size of the first transitive component (FTC). For larger networks, we find measures of network efficiency, incoherence and heterogeneity (degree variance) correlate with robustness of networks to increasing excitability. These results provide a set of potential prognostic markers for therapeutic interventions in epilepsy. Such markers could be used to support the development of personalized treatment strategies, ultimately contributing to understanding of long-term seizure freedom.

Biophysical models applied to dementia patients reveal links between geographical origin, gender, disease duration, and loss of neural inhibition

BACKGROUND

The hypothesis of decreased neural inhibition in dementia has been sparsely studied in functional magnetic resonance imaging (fMRI) data across patients with different dementia subtypes, and the role of social and demographic heterogeneities on this hypothesis remains to be addressed.

METHODS

We inferred regional inhibition by fitting a biophysical whole-brain model (dynamic mean field model with realistic inter-areal connectivity) to fMRI data from 414 participants, including patients with Alzheimer's disease, behavioral variant frontotemporal dementia, and controls. We then investigated the effect of disease condition, and demographic and clinical variables on the local inhibitory feedback, a variable related to the maintenance of balanced neural excitation/inhibition.

RESULTS

Decreased local inhibitory feedback was inferred from the biophysical modeling results in dementia patients, specific to brain areas presenting neurodegeneration. This loss of local inhibition correlated positively with years with disease, and showed differences regarding the gender and geographical origin of the patients. The model correctly reproduced known disease-related changes in functional connectivity.

CONCLUSIONS

Results suggest a critical link between abnormal neural and circuit-level excitability levels, the loss of grey matter observed in dementia, and the reorganization of functional connectivity, while highlighting the sensitivity of the underlying biophysical mechanism to demographic and clinical heterogeneities in the patient population.

Bifurcations and bursting in the Epileptor

The Epileptor is a phenomenological model for seizure activity that is used in a personalized large-scale brain modeling framework, the Virtual Epileptic Patient, with the aim of improving surgery outcomes for drug-resistant epileptic patients. Transitions between interictal and ictal states are modeled as bifurcations, enabling the definition of seizure classes in terms of onset/offset bifurcations. This establishes a taxonomy of seizures grounded in their essential underlying dynamics and the Epileptor replicates the activity of the most common class, as observed in patients with focal epilepsy, which is characterized by square-wave bursting properties. The Epileptor also encodes an additional mechanism to account for interictal spikes and spike and wave discharges. Here we use insights from a more generic model for square-wave bursting, based on the Unfolding Theory approach, to guide the bifurcation analysis of the Epileptor and gain a deeper understanding of the model and the role of its parameters. We show how the Epileptor's parameters can be modified to produce activities for other seizures classes of the taxonomy, as observed in patients, so that the large-scale brain models could be further personalized. Some of these classes have already been described in the literature in the Epileptor, others, predicted by the generic model, are new. Finally, we unveil how the interaction with the additional mechanism for spike and wave discharges alters the bifurcation structure of the main burster.

Epileptic seizure suppression: A computational approach for identification and control using real data

Epilepsy affects millions of people worldwide every year and remains an open subject for research. Current development on this field has focused on obtaining computational models to better understand its triggering mechanisms, attain realistic descriptions and study seizure suppression. Controllers have been successfully applied to mitigate epileptiform activity in dynamic models written in state-space notation, whose applicability is, however, restricted to signatures that are accurately described by them. Alternatively, autoregressive modeling (AR), a typical data-driven tool related to system identification (SI), can be directly applied to signals to generate more realistic models, and since it is inherently convertible into state-space representation, it can thus be used for the artificial reconstruction and attenuation of seizures as well. Considering this, the first objective of this work is to propose an SI approach using AR models to describe real epileptiform activity. The second objective is to provide a strategy for reconstructing and mitigating such activity artificially, considering non-hybrid and hybrid controllers - designed from ictal and interictal events, respectively. The results show that AR models of relatively low order represent epileptiform activities fairly well and both controllers are effective in attenuating the undesired activity while simultaneously driving the signal to an interictal condition. These findings may lead to customized models based on each signal, brain region or patient, from which it is possible to better define shape, frequency and duration of external stimuli that are necessary to attenuate seizures.

Identifying sources of human interictal discharges with travelling wave and white matter propagation

Interictal epileptiform discharges have been shown to propagate from focal epileptogenic sources as travelling waves or through more rapid white matter conduction. We hypothesize that both modes of propagation are necessary to explain interictal discharge timing delays. We propose a method that, for the first time, incorporates both propagation modes to identify unique potential sources of interictal activity. We retrospectively analysed 38 focal epilepsy patients who underwent intracranial EEG recordings and diffusion-weighted imaging for epilepsy surgery evaluation. Interictal discharges were detected and localized to the most likely source based on relative delays in time of arrival across electrodes, incorporating travelling waves and white matter propagation. We assessed the influence of white matter propagation on distance of spread, timing and clinical interpretation of interictal activity. To evaluate accuracy, we compared our source localization results to earliest spiking regions to predict seizure outcomes. White matter propagation helps to explain the timing delays observed in interictal discharge sequences, underlying rapid and distant propagation. Sources identified based on differences in time of receipt of interictal discharges are often distinct from the leading electrode location. Receipt of activity propagating rapidly via white matter can occur earlier than more local activity propagating via slower cortical travelling waves. In our cohort, our source localization approach was more accurate in predicting seizure outcomes than the leading electrode location. Inclusion of white matter in addition to travelling wave propagation in our model of discharge spread did not improve overall accuracy but allowed for identification of unique and at times distant potential sources of activity, particularly in patients with persistent postoperative seizures. Since distant white matter propagation can occur more rapidly than local travelling wave propagation, combined modes of propagation within an interictal discharge sequence can decouple the commonly assumed relationship between spike timing and distance from the source. Our findings thus highlight the clinical importance of recognizing the presence of dual modes of propagation during interictal discharges, as this may be a cause of clinical mislocalization.

Critical-like Brain Dynamics in a Continuum from Second- to First-Order Phase Transition

The classic brain criticality hypothesis postulates that the brain benefits from operating near a continuous second-order phase transition. Slow feedback regulation of neuronal activity could, however, lead to a discontinuous first-order transition and thereby bistable activity. Observations of bistability in awake brain activity have nonetheless remained scarce and its functional significance unclear. Moreover, there is no empirical evidence to support the hypothesis that the human brain could flexibly operate near either a first- or second-order phase transition despite such a continuum being common in models. Here, using computational modeling, we found bistable synchronization dynamics to emerge through elevated positive feedback and occur exclusively in a regimen of critical-like dynamics. We then assessed bistability in vivo with resting-state MEG in healthy adults (7 females, 11 males) and stereo-electroencephalography in epilepsy patients (28 females, 36 males). This analysis revealed that a large fraction of the neocortices exhibited varying degrees of bistability in neuronal oscillations from 3 to 200 Hz. In line with our modeling results, the neuronal bistability was positively correlated with classic assessment of brain criticality across narrow-band frequencies. Excessive bistability was predictive of epileptic pathophysiology in the patients, whereas moderate bistability was positively correlated with task performance in the healthy subjects. These empirical findings thus reveal the human brain as a one-of-a-kind complex system that exhibits critical-like dynamics in a continuum between continuous and discontinuous phase transitions.SIGNIFICANCE STATEMENT In the model, while synchrony per se was controlled by connectivity, increasing positive local feedback led to gradually emerging bistable synchrony with scale-free dynamics, suggesting a continuum between second- and first-order phase transitions in synchrony dynamics inside a critical-like regimen. In resting-state MEG and SEEG, bistability of ongoing neuronal oscillations was pervasive across brain areas and frequency bands and was observed only with concurring critical-like dynamics as the modeling predicted. As evidence for functional relevance, moderate bistability was positively correlated with executive functioning in the healthy subjects, and excessive bistability was associated with epileptic pathophysiology. These findings show that critical-like neuronal dynamics in vivo involves both continuous and discontinuous phase transitions in a frequency-, neuroanatomy-, and state-dependent manner.

Complementary structural and functional abnormalities to localise epileptogenic tissue

BACKGROUND

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

METHODS

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

FINDINGS

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

INTERPRETATION

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

FUNDING

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

Complementary structural and functional abnormalities to localise epileptogenic tissue

Background

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

Methods

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

Findings

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

Interpretation

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

Funding

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

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

BACKGROUND AND OBJECTIVES

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

METHODS

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

RESULTS

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

DISCUSSION

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

CLASSIFICATION OF EVIDENCE

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

Personalized dynamic network models of the human brain as a future tool for planning and optimizing epilepsy therapy

Epilepsy is a common neurological disorder, with one third of patients not responding to currently available antiepileptic drugs. The proportion of pharmacoresistant epilepsies has remained unchanged for many decades. To cure epilepsy and control seizures requires a paradigm shift in the development of new approaches to epilepsy diagnosis and treatment. Contemporary medicine has benefited from the exponential growth of computational modeling, and the application of network dynamics theory to understanding and treating human brain disorders. In epilepsy, the introduction of these approaches has led to personalized epileptic network modeling that can explore the patient's seizure genesis and predict the functional impact of resection on its individual network's propensity to seize. The application of the dynamic systems approach to neurostimulation therapy of epilepsy allows designing stimulation strategies that consider the patient's seizure dynamics and long-term fluctuations in the stability of their epileptic networks. In this article, we review, in a nontechnical fashion suitable for a broad neuroscientific audience, recent progress in personalized dynamic brain network modeling that is shaping the future approach to the diagnosis and treatment of epilepsy.

Does seizure propagate within or across intrinsic brain networks? An intracranial EEG study

BACKGROUND

Understanding the spatiotemporal propagation profiles of seizures is crucial for the preoperative assessment of epilepsy patients. The present study aimed to investigate whether seizures exhibit propagation patterns that align with intrinsic networks (INs).

METHODS

A quantitative analysis was conducted to examine ictal fast activity (IFA). The Epileptogenicity Index (EI) was employed to assess the epileptogenicity, spectral features, and temporal characteristics of IFA. Intra-network and inter-network comparisons were made regarding the IFA-related metrics. Additionally, the metrics were correlated with Euclidean distance. Network connection maps were generated to visualize seizures originating from different INs, allowing for comparisons between distinct groups.

RESULTS

Data for 81 seizures in 43 subjects were captured using stereoelectroencephalography implantation. Three metrics were compared: EI, time involvement (TI), and energy ratio index (ERI). Intra-network channels exhibited higher EI, earlier involvement of IFA, and stronger high-frequency energy. These findings were further validated through subgroup analyses stratified by neuropathology, seizure type, and seizure origination lobe. Correlation analyses revealed a negative association between distance and both EI and ERI, while distance exhibited a positive correlation with TI. Seizures originating from different INs exhibited varying propagation characteristics.

CONCLUSIONS

The study findings highlight the dominant role of intra-network dynamics over inter-network during seizure propagation. These results contribute to our understanding of seizure dynamics and their relationship with INs.

The Interictal Suppression Hypothesis in focal epilepsy: network-level supporting evidence

Why are people with focal epilepsy not continuously having seizures? Previous neuronal signalling work has implicated gamma-aminobutyric acid balance as integral to seizure generation and termination, but is a high-level distributed brain network involved in suppressing seizures? Recent intracranial electrographic evidence has suggested that seizure-onset zones have increased inward connectivity that could be associated with interictal suppression of seizure activity. Accordingly, we hypothesize that seizure-onset zones are actively suppressed by the rest of the brain network during interictal states. Full testing of this hypothesis would require collaboration across multiple domains of neuroscience. We focused on partially testing this hypothesis at the electrographic network level within 81 individuals with drug-resistant focal epilepsy undergoing presurgical evaluation. We used intracranial electrographic resting-state and neurostimulation recordings to evaluate the network connectivity of seizure onset, early propagation and non-involved zones. We then used diffusion imaging to acquire estimates of white-matter connectivity to evaluate structure-function coupling effects on connectivity findings. Finally, we generated a resting-state classification model to assist clinicians in detecting seizure-onset and propagation zones without the need for multiple ictal recordings. Our findings indicate that seizure onset and early propagation zones demonstrate markedly increased inwards connectivity and decreased outwards connectivity using both resting-state (one-way ANOVA, P-value = 3.13 × 10-13) and neurostimulation analyses to evaluate evoked responses (one-way ANOVA, P-value = 2.5 × 10-3). When controlling for the distance between regions, the difference between inwards and outwards connectivity remained stable up to 80 mm between brain connections (two-way repeated measures ANOVA, group effect P-value of 2.6 × 10-12). Structure-function coupling analyses revealed that seizure-onset zones exhibit abnormally enhanced coupling (hypercoupling) of surrounding regions compared to presumably healthy tissue (two-way repeated measures ANOVA, interaction effect P-value of 9.76 × 10-21). Using these observations, our support vector classification models achieved a maximum held-out testing set accuracy of 92.0 ± 2.2% to classify early propagation and seizure-onset zones. These results suggest that seizure-onset zones are actively segregated and suppressed by a widespread brain network. Furthermore, this electrographically observed functional suppression is disproportionate to any observed structural connectivity alterations of the seizure-onset zones. These findings have implications for the identification of seizure-onset zones using only brief electrographic recordings to reduce patient morbidity and augment the presurgical evaluation of drug-resistant epilepsy. Further testing of the interictal suppression hypothesis can provide insight into potential new resective, ablative and neuromodulation approaches to improve surgical success rates in those suffering from drug-resistant focal epilepsy.

Observers for Phenomenological Models of Epileptic Seizures

Progress towards effective treatment of epileptic seizures has seen much improvement in the past decade. In particular, the emergence of phenomenological models of epileptic seizures specifically designed to capture the electrical seizure dynamics in the Epileptor model is inspiring new approaches to predicting and controlling seizures. These new models present in various forms and contain important but unmeasurable variables that control the occurrence of seizures. These models have been used mostly as nodes in large networks to study the complex brain behaviour of seizures. In order to use this model for the purposes of seizure forecasting or to control seizures through deep brain stimulation, the states of the model will need to be estimated. Although devices such as EEG electrodes can be related to some of the states of the model, most remain unmeasured and would require an observer (as defined in control theory) for their estimation. Additionally, we would like to consider the case for large nodes of systems where the number of electrodes is far smaller than the number of nodes being estimated. In this paper, we provide methods towards obtaining the full states of these phenomenological models using nonlinear observers. In particular, we explore the effectiveness of the Extended Kalman Filter for small networks of nodes of a smoothed sixth order Epileptor model. We show that observer design is possible for this family of systems and identify the difficulties in doing so.Clinical relevance-The methods presented here can be applied with an individual epileptic patient's EEG to reveal previously hidden biomarkers of epilepsy for seizure forecasting.

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

SIGNIFICANCE

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

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.

Linear feedback control of spreading dynamics in stochastic nonlinear network models: epileptic seizures

The development of models and approaches for controlling the spreading dynamics of epileptic seizures is an essential step towards new therapies for people with pharmacologically resistant epilepsy. Beyond resective neurosurgery, in which epileptogenic zones (EZs) are the target of surgery, closed-loop control based on intracranial electrical stimulation, applied at the very early stage of seizure evolution, has been the main alternative, e.g. the RNS system from NeuroPace (Mountain View, CA). In this approach the electrical stimulation is delivered to target brain areas after detection of seizure initiation in the EZ. Here, we examined, on model simulations, some of the closed-loop control aspects of the problem. Seizure dynamics and spread are typically modeled with highly nonlinear dynamics on complex brain networks. Despite the nonlinearity and complexity, currently available optimal feedback control approaches are mostly based on linear approximations. Alternative machine learning control approaches might require amounts of data beyond what is commonly available in the intended application. We thus examined how standard linear feedback control approaches perform when applied to nonlinear models of neural dynamics of seizure generation and spread. In particular, we considered patient-specific epileptor network models for seizure onset and spread. The models incorporate network connectivity derived from (diffusion MRI) white-matter tractography, have been shown to capture the qualitative dynamics of epileptic seizures and can be fit to patient data. For control, we considered simple linear quadratic Gaussian (LQG) regulators. The LQG control was based on a discrete-time state-space model derived from the linearization of the patient-specific epileptor network model around a stable fixed point in the regime of preictal dynamics. We show in simulations that LQG regulators acting on the EZ node during the initial seizure period tend to be unstable. The LQG solution for the control law in this case leads to global feedback to the EZ-node actuator. However, if the LQG solution is constrained to depend on only local feedback originating from the EZ node itself, the controller is stable. In this case, we demonstrate that localized LQG can easily terminate the seizure at the early stage and prevent spread. In the context of optimal feedback control based on linear approximations, our results point to the need for investigating in more detail feedback localization and additional relevant control targets beyond epileptogenic zones.

Personalized whole brain modeling of status epilepticus

How status epilepticus (SE) is generated and propagates in the brain is not known. As for seizures, a patient-specific approach is necessary, and the analysis should be performed at the whole brain level. Personalized brain models can be used to study seizure genesis and propagation at the whole brain scale in The Virtual Brain (TVB), using the Epileptor mathematical construct. Building on the fact that SE is part of the repertoire of activities that the Epileptor can generate, we present the first attempt to model SE at the whole brain scale in TVB, using data from a patient who experienced SE during presurgical evaluation. Simulations reproduced the patterns found with SEEG recordings. We find that if, as expected, the pattern of SE propagation correlates with the properties of the patient's structural connectome, SE propagation also depends upon the global state of the network, i.e., that SE propagation is an emergent property. We conclude that individual brain virtualization can be used to study SE genesis and propagation. This type of theoretical approach may be used to design novel interventional approaches to stop SE. This paper was presented at the 8th London-Innsbruck Colloquium on Status Epilepticus and Acute Seizures held in September 2022.

Personalised virtual brain models in epilepsy

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

Amortized Bayesian inference on generative dynamical network models of epilepsy using deep neural density estimators

Whole-brain modeling of epilepsy combines personalized anatomical data with dynamical models of abnormal activities to generate spatio-temporal seizure patterns as observed in brain imaging data. Such a parametric simulator is equipped with a stochastic generative process, which itself provides the basis for inference and prediction of the local and global brain dynamics affected by disorders. However, the calculation of likelihood function at whole-brain scale is often intractable. Thus, likelihood-free algorithms are required to efficiently estimate the parameters pertaining to the hypothetical areas, ideally including the uncertainty. In this study, we introduce the simulation-based inference for the virtual epileptic patient model (SBI-VEP), enabling us to amortize the approximate posterior of the generative process from a low-dimensional representation of whole-brain epileptic patterns. The state-of-the-art deep learning algorithms for conditional density estimation are used to readily retrieve the statistical relationships between parameters and observations through a sequence of invertible transformations. We show that the SBI-VEP is able to efficiently estimate the posterior distribution of parameters linked to the extent of the epileptogenic and propagation zones from sparse intracranial electroencephalography recordings. The presented Bayesian methodology can deal with non-linear latent dynamics and parameter degeneracy, paving the way for fast and reliable inference on brain disorders from neuroimaging modalities.

Criticality in probabilistic models of spreading dynamics in brain networks: Epileptic seizures

The spread of seizures across brain networks is the main impairing factor, often leading to loss-of-consciousness, in people with epilepsy. Despite advances in recording and modeling brain activity, uncovering the nature of seizure spreading dynamics remains an important challenge to understanding and treating pharmacologically resistant epilepsy. To address this challenge, we introduce a new probabilistic model that captures the spreading dynamics in patient-specific complex networks. Network connectivity and interaction time delays between brain areas were estimated from white-matter tractography. The model's computational tractability allows it to play an important complementary role to more detailed models of seizure dynamics. We illustrate model fitting and predictive performance in the context of patient-specific Epileptor networks. We derive the phase diagram of spread size (order parameter) as a function of brain excitability and global connectivity strength, for different patient-specific networks. Phase diagrams allow the prediction of whether a seizure will spread depending on excitability and connectivity strength. In addition, model simulations predict the temporal order of seizure spread across network nodes. Furthermore, we show that the order parameter can exhibit both discontinuous and continuous (critical) phase transitions as neural excitability and connectivity strength are varied. Existence of a critical point, where response functions and fluctuations in spread size show power-law divergence with respect to control parameters, is supported by mean-field approximations and finite-size scaling analyses. Notably, the critical point separates two distinct regimes of spreading dynamics characterized by unimodal and bimodal spread-size distributions. Our study sheds new light on the nature of phase transitions and fluctuations in seizure spreading dynamics. We expect it to play an important role in the development of closed-loop stimulation approaches for preventing seizure spread in pharmacologically resistant epilepsy. Our findings may also be of interest to related models of spreading dynamics in epidemiology, biology, finance, and statistical physics.

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

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

Interictal Functional Connectivity in Focal Refractory Epilepsies Investigated by Intracranial EEG

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

Advancing brain network models to reconcile functional neuroimaging and clinical research

Functional magnetic resonance imaging (fMRI) captures information on brain function beyond the anatomical alterations that are traditionally visually examined by neuroradiologists. However, the fMRI signals are complex in addition to being noisy, so fMRI still faces limitations for clinical applications. Here we review methods that have been proposed as potential solutions so far, namely statistical, biophysical and decoding models, with their strengths and weaknesses. We especially evaluate the ability of these models to directly predict clinical variables from their parameters (predictability) and to extract clinically relevant information regarding biological mechanisms and relevant features for classification and prediction (interpretability). We then provide guidelines for useful applications and pitfalls of such fMRI-based models in a clinical research context, looking beyond the current state of the art. In particular, we argue that the clinical relevance of fMRI calls for a new generation of models for fMRI data, which combine the strengths of both biophysical and decoding models. This leads to reliable and biologically meaningful model parameters, which thus fulfills the need for simultaneous interpretability and predictability. In our view, this synergy is fundamental for the discovery of new pharmacological and interventional targets, as well as the use of models as biomarkers in neurology and psychiatry.

Tracking Multisite Seizure Propagation Using Ictal High-Gamma Activity

PURPOSE

Spatial patterns of long-range seizure propagation in epileptic networks have not been well characterized. Here, we use ictal high-gamma activity (HGA) as a proxy of intense neuronal population firing to map the spatial evolution of seizure recruitment.

METHODS

Ictal HGA (80-150 Hz) was analyzed in 13 patients with 72 seizures recorded by stereotactic depth electrodes, using previously validated methods. Distinct spatial clusters of channels with the ictal high-gamma signature were identified, and seizure hubs were defined as stereotypically recruited nonoverlapping clusters. Clusters correlated with asynchronous seizure terminations to provide supportive evidence for independent seizure activity at these sites. The spatial overlap between seizure hubs and interictal ripples was compared.

RESULTS

Ictal HGA was detected in 71% of seizures and 10% of implanted contacts, enabling tracking of contiguous and noncontiguous seizure recruitment. Multiple seizure hubs were identified in 54% of cases, including 43% of patients thought preoperatively to have unifocal epilepsy. Noncontiguous recruitment was associated with asynchronous seizure termination (odds ratio = 19.7; p = 0.029). Interictal ripples demonstrated greater spatial overlap with ictal HGA in cases with single seizure hubs compared with those with multiple hubs (100% vs. 66% per patient; p = 0.03).

CONCLUSIONS

Ictal HGA may serve as a useful adjunctive biomarker to distinguish contiguous seizure spread from propagation to remote seizure sites. High-gamma sites were found to cluster in stereotyped seizure hubs rather than being broadly distributed. Multiple hubs were common even in cases that were considered unifocal.

A perturbative approach to study information communication in brain networks

How communication among neuronal ensembles shapes functional brain dynamics is a question of fundamental importance to neuroscience. Communication in the brain can be viewed as a product of the interaction of node activities with the structural network over which these activities flow. The study of these interactions is, however, restricted by the difficulties in describing the complex dynamics of the brain. There is thus a need to develop methods to study these network-dynamical interactions and how they impact information flow, without having to ascertain dynamics a priori or resort to restrictive analytical approaches. Here, we adapt a recently established network analysis method based on perturbations, it to a neuroscientific setting to study how information flow in the brain can raise from properties of underlying structure. For proof-of-concept, we apply the approach on in silico whole-brain models. We expound on the functional implications of the distributions of metrics that capture network-dynamical interactions, termed net influence and flow. We also study the network-dynamical interactions at the level of resting-state networks. An attractive feature of this method is its simplicity, which allows a direct translation to an experimental or clinical setting, such as for identifying targets for stimulation studies or therapeutic interventions.

Estimating the Parameters of the Epileptor Model for Epileptic Seizure Suppression

Epilepsy is one of the most common brain disorders worldwide, affecting millions of people every year. Given the partially successful existing treatments for epileptiform activity suppression, dynamic mathematical models have been proposed with the purpose of better understanding the factors that might trigger an epileptic seizure and how to mitigate it, among which Epileptor stands out, due to its relative simplicity and consistency with experimental observations. Recent studies using this model have provided evidence that establishing a feedback-based control approach is possible. However, for this strategy to work properly, Epileptor's parameters, which describe the dynamic characteristics of a seizure, must be known beforehand. Therefore, this work proposes a methodology for estimating such parameters based on a successive optimization technique. The results show that it is feasible to approximate their values as they converge to reference values based on different initial conditions, which are modeled by an uncertainty factor or noise addition. Also, interictal (healthy) and ictal (ongoing seizure) conditions, as well as time resolution, must be taken into account for an appropriate estimation. At last, integrating such a parameter estimation approach with observers and controllers for purposes of seizure suppression is carried out, which might provide an interesting alternative for seizure suppression in practice in the future.

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

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

Fully Bayesian estimation of virtual brain parameters with self-tuning Hamiltonian Monte Carlo

Virtual brain models are data-driven patient-specific brain models integrating individual brain imaging data with neural mass modeling in a single computational framework, capable of autonomously generating brain activity and its associated brain imaging signals. Along the example of epilepsy, we develop an efficient and accurate Bayesian methodology estimating the parameters linked to the extent of the epileptogenic zone. State-of-the-art advances in Bayesian inference using Hamiltonian Monte Carlo (HMC) algorithms have remained elusive for large-scale differential-equations based models due to their slow convergence. We propose appropriate priors and a novel reparameterization to facilitate efficient exploration of the posterior distribution in terms of computational time and convergence diagnostics. The methodology is illustrated for in-silico dataset and then, applied to infer the personalized model parameters based on the empirical stereotactic electroencephalography (SEEG) recordings of retrospective patients. This improved methodology may pave the way to render HMC methods sufficiently easy and efficient to use, thus applicable in personalized medicine.

Connectivity Alterations in Emotional and Cognitive Networks During a Manic State Induced by Direct Electrical Stimulation

Mania is characterized by affective and cognitive alterations, with heightened external and self-awareness that are opposite to the alteration of awareness during epileptic seizures. Electrical stimulations carried out routinely during stereotactic intracerebral EEG (SEEG) recordings for presurgical evaluation of epilepsy may represent a unique opportunity to study the pathophysiology of such complex emotional-behavioral phenomenon, particularly difficult to reproduce in experimental setting. We investigated SEEG signals-based functional connectivity between different brain regions involved in emotions and in consciousness processing during a manic state induced by electrical stimulation in a patient with drug-resistant focal epilepsy. The stimulation inducing manic state and an asymptomatic stimulation of the same site, as well as a seizure with alteration of awareness (AOA) were analyzed. Functional connectivity analysis was performed by measuring interdependencies (nonlinear regression analysis based on the h2 coefficient) between broadband SEEG signals and within typical sub-bands, before and after stimulation, or before and during the seizure with AOA, respectively. Stimulation of the right lateral prefrontal cortex induced a manic state lasting several hours. Its onset was associated with significant increase of broadband-signal functional coupling between the right hemispheric limbic nodes, the temporal pole and the claustrum, whereas significant decorrelation between the right lateral prefrontal and the anterior cingulate cortex was observed in theta-band. In contrast, ictal alteration of awareness was associated with increased broadband and sub-bands synchronization within and between the internal and external awareness networks, including the anterior and middle cingulate, the mesial and lateral prefrontal, the inferior parietal and the temporopolar cortex. Our data suggest the existence of network- and frequency-specific functional connectivity patterns during manic state. A transient desynchronization of theta activity between the external and internal awareness network hubs is likely to increase awareness, with potential therapeutic effect.

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

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

Volumetric and structural connectivity abnormalities co-localise in TLE

Patients with temporal lobe epilepsy (TLE) exhibit both volumetric and structural connectivity abnormalities relative to healthy controls. How these abnormalities inter-relate and their mechanisms are unclear. We computed grey matter volumetric changes and white matter structural connectivity abnormalities in 144 patients with unilateral TLE and 96 healthy controls. Regional volumes were calculated using T1-weighted MRI, while structural connectivity was derived using white matter fibre tractography from diffusion-weighted MRI. For each regional volume and each connection strength, we calculated the effect size between patient and control groups in a group-level analysis. We then applied hierarchical regression to investigate the relationship between volumetric and structural connectivity abnormalities in individuals. Additionally, we quantified whether abnormalities co-localised within individual patients by computing Dice similarity scores. In TLE, white matter connectivity abnormalities were greater when joining two grey matter regions with abnormal volumes. Similarly, grey matter volumetric abnormalities were greater when joined by abnormal white matter connections. The extent of volumetric and connectivity abnormalities related to epilepsy duration, but co-localisation did not. Co-localisation was primarily driven by neighbouring abnormalities in the ipsilateral hemisphere. Overall, volumetric and structural connectivity abnormalities were related in TLE. Our results suggest that shared mechanisms may underlie changes in both volume and connectivity alterations in patients with TLE.

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.

Normalizing the brain connectome for communication through synchronization

Networks in neuroscience determine how brain function unfolds, and their perturbations lead to psychiatric disorders and brain disease. Brain networks are characterized by their connectomes, which comprise the totality of all connections, and are commonly described by graph theory. This approach is deeply rooted in a particle view of information processing, based on the quantification of informational bits such as firing rates. Oscillations and brain rhythms demand, however, a wave perspective of information processing based on synchronization. We extend traditional graph theory to a dual, particle-wave, perspective, integrate time delays due to finite transmission speeds, and derive a normalization of the connectome. When applied to the database of the Human Connectome Project, it explains the emergence of frequency-specific network cores including the visual and default mode networks. These findings are robust across human subjects (N = 100) and are a fundamental network property within the wave picture. The normalized connectome comprises the particle view in the limit of infinite transmission speeds and opens the applicability of graph theory to a wide range of novel network phenomena, including physiological and pathological brain rhythms. These two perspectives are orthogonal, but not incommensurable, when understood within the novel, here-proposed, generalized framework of structural connectivity.