1
|
Ye C, Zhang Y, Ran C, Ma T. Recent Progress in Brain Network Models for Medical Applications: A Review. HEALTH DATA SCIENCE 2024; 4:0157. [PMID: 38979037 PMCID: PMC11227951 DOI: 10.34133/hds.0157] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/31/2023] [Accepted: 05/28/2024] [Indexed: 07/10/2024]
Abstract
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.
Collapse
Affiliation(s)
- Chenfei Ye
- International Research Institute for Artificial Intelligence,
Harbin Institute of Technology at Shenzhen, Shenzhen, China
| | - Yixuan Zhang
- Department of Electronic and Information Engineering,
Harbin Institute of Technology at Shenzhen, Shenzhen, China
| | - Chen Ran
- Department of Electronic and Information Engineering,
Harbin Institute of Technology at Shenzhen, Shenzhen, China
| | - Ting Ma
- International Research Institute for Artificial Intelligence,
Harbin Institute of Technology at Shenzhen, Shenzhen, China
- Department of Electronic and Information Engineering,
Harbin Institute of Technology at Shenzhen, Shenzhen, China
- Peng Cheng Laboratory, Shenzhen, China
- Guangdong Provincial Key Laboratory of Aerospace Communication and Networking Technology,
Harbin Institute of Technology at Shenzhen, China
| |
Collapse
|
2
|
Millán AP, van Straaten ECW, Stam CJ, Nissen IA, Idema S, Van Mieghem P, Hillebrand A. Individualized epidemic spreading models predict epilepsy surgery outcomes: A pseudo-prospective study. Netw Neurosci 2024; 8:437-465. [PMID: 38952815 PMCID: PMC11142635 DOI: 10.1162/netn_a_00361] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2023] [Accepted: 01/18/2024] [Indexed: 07/03/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Ana P. Millán
- Amsterdam UMC, Vrije Universiteit Amsterdam, Department of Clinical Neurophysiology and MEG Center, Amsterdam, The Netherlands
- Institute “Carlos I” for Theoretical and Computational Physics, and Electromagnetism and Matter Physics Department, University of Granada, Granada, Spain
| | - Elisabeth C. W. van Straaten
- Amsterdam UMC, Vrije Universiteit Amsterdam, Department of Clinical Neurophysiology and MEG Center, Amsterdam, The Netherlands
- Amsterdam Neuroscience, Systems and Network Neurosciences, Amsterdam, The Netherlands
- Amsterdam Neuroscience, Neurodegeneration, Amsterdam, The Netherlands
| | - Cornelis J. Stam
- Amsterdam UMC, Vrije Universiteit Amsterdam, Department of Clinical Neurophysiology and MEG Center, Amsterdam, The Netherlands
- Amsterdam Neuroscience, Brain Imaging, Amsterdam, The Netherlands
- Amsterdam Neuroscience, Neurodegeneration, Amsterdam, The Netherlands
| | - Ida A. Nissen
- Amsterdam UMC, Vrije Universiteit Amsterdam, Department of Clinical Neurophysiology and MEG Center, Amsterdam, The Netherlands
| | - Sander Idema
- Amsterdam UMC, Vrije Universiteit Amsterdam, Department of Neurosurgery, Amsterdam, The Netherlands
- Amsterdam Neuroscience, Systems and Network Neurosciences, Amsterdam, The Netherlands
- Amsterdam Neuroscience, Cancer Biology and Immonology, Amsterdam, The Netherlands
| | - Piet Van Mieghem
- Faculty of Electrical Engineering, Mathematics and Computer Science, Delft University of Technology, Delft, The Netherlands
| | - Arjan Hillebrand
- Amsterdam UMC, Vrije Universiteit Amsterdam, Department of Clinical Neurophysiology and MEG Center, Amsterdam, The Netherlands
- Amsterdam Neuroscience, Brain Imaging, Amsterdam, The Netherlands
- Amsterdam Neuroscience, Systems and Network Neurosciences, Amsterdam, The Netherlands
| |
Collapse
|
3
|
Saggio ML, Jirsa V. Bifurcations and bursting in the Epileptor. PLoS Comput Biol 2024; 20:e1011903. [PMID: 38446814 PMCID: PMC10947678 DOI: 10.1371/journal.pcbi.1011903] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2023] [Revised: 03/18/2024] [Accepted: 02/08/2024] [Indexed: 03/08/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Maria Luisa Saggio
- Institut de Neurosciences des Systemes INS UMR1106, AMU, INSERM, Marseille, France
| | - Viktor Jirsa
- Institut de Neurosciences des Systemes INS UMR1106, AMU, INSERM, Marseille, France
| |
Collapse
|
4
|
Brogin JAF, Faber J, Reyes-Garcia SZ, Cavalheiro EA, Bueno DD. Epileptic seizure suppression: A computational approach for identification and control using real data. PLoS One 2024; 19:e0298762. [PMID: 38416729 PMCID: PMC10901337 DOI: 10.1371/journal.pone.0298762] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2022] [Accepted: 01/31/2024] [Indexed: 03/01/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- João A. F. Brogin
- Department of Mechanical Engineering, São Paulo State University (UNESP), School of Engineering of Ilha Solteira, Ilha Solteira, São Paulo, Brazil
| | - Jean Faber
- Department of Neurology and Neurosurgery, Federal University of São Paulo (UNIFESP), São Paulo, São Paulo, Brazil
| | - Selvin Z. Reyes-Garcia
- Departamento de Ciencias Morfológicas, Facultad de Ciencias Médicas, Universidad Nacional Autónoma de Honduras, Tegucigalpa, Honduras
| | - Esper A. Cavalheiro
- Department of Neurology and Neurosurgery, Federal University of São Paulo (UNIFESP), São Paulo, São Paulo, Brazil
| | - Douglas D. Bueno
- Department of Mathematics, São Paulo State University (UNESP), School of Engineering of Ilha Solteira, Ilha Solteira, São Paulo, Brazil
| |
Collapse
|
5
|
Courson J, Quoy M, Timofeeva Y, Manos T. An exploratory computational analysis in mice brain networks of widespread epileptic seizure onset locations along with potential strategies for effective intervention and propagation control. Front Comput Neurosci 2024; 18:1360009. [PMID: 38468870 PMCID: PMC10925689 DOI: 10.3389/fncom.2024.1360009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2023] [Accepted: 02/08/2024] [Indexed: 03/13/2024] Open
Abstract
Mean-field models have been developed to replicate key features of epileptic seizure dynamics. However, the precise mechanisms and the role of the brain area responsible for seizure onset and propagation remain incompletely understood. In this study, we employ computational methods within The Virtual Brain framework and the Epileptor model to explore how the location and connectivity of an Epileptogenic Zone (EZ) in a mouse brain are related to focal seizures (seizures that start in one brain area and may or may not remain localized), with a specific focus on the hippocampal region known for its association with epileptic seizures. We then devise computational strategies to confine seizures (prevent widespread propagation), simulating medical-like treatments such as tissue resection and the application of an anti-seizure drugs or neurostimulation to suppress hyperexcitability. Through selectively removing (blocking) specific connections informed by the structural connectome and graph network measurements or by locally reducing outgoing connection weights of EZ areas, we demonstrate that seizures can be kept constrained around the EZ region. We successfully identified the minimal connections necessary to prevent widespread seizures, with a particular focus on minimizing surgical or medical intervention while simultaneously preserving the original structural connectivity and maximizing brain functionality.
Collapse
Affiliation(s)
- Juliette Courson
- ETIS Lab, ENSEA, CNRS, UMR8051, CY Cergy-Paris University, Cergy, France
- Laboratoire de Physique Théorique et Modélisation, UMR 8089, CY Cergy Paris Université, CNRS, Cergy-Pontoise, France
- Department of Computer Science, University of Warwick, Coventry, United Kingdom
| | - Mathias Quoy
- ETIS Lab, ENSEA, CNRS, UMR8051, CY Cergy-Paris University, Cergy, France
- IPAL CNRS Singapore, Singapore, Singapore
| | - Yulia Timofeeva
- Department of Computer Science, University of Warwick, Coventry, United Kingdom
- Department of Clinical and Experimental Epilepsy, UCL Queen Square Institute of Neurology, University College London, London, United Kingdom
| | - Thanos Manos
- ETIS Lab, ENSEA, CNRS, UMR8051, CY Cergy-Paris University, Cergy, France
| |
Collapse
|
6
|
Li Z, Zhang H, Niu S, Xing Y. Localizing epileptogenic zones with high-frequency oscillations and directed connectivity. Seizure 2023; 111:9-16. [PMID: 37487273 DOI: 10.1016/j.seizure.2023.07.013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2023] [Revised: 07/19/2023] [Accepted: 07/20/2023] [Indexed: 07/26/2023] Open
Abstract
PURPOSE Precise localization of the epileptogenic zone (EZ) is essential for epilepsy surgery. Existing methods often fail to detect slow onset patterns or similar neural activities presented in the recorded signals. To address this issue, we propose a new measure to quantify epileptogenicity, i.e., the connectivity high-frequency epileptogenicity index (cHFEI). METHODS The cHFEI method combines directed connectivity and high-frequency oscillations (HFOs) to measure the epileptogenicity of regions involved in a brain network. By applying this method to stereoelectroencephalography (SEEG) recordings of 49 seizures in 20 patients, we calculated the accuracy, sensitivity, and precision with a visually identified epileptogenic zone as a reference. The performance was evaluated by the confusion matrix and the area under the receiver operating characteristic (ROC) curve. RESULTS Epileptic network estimation based on cHFEI successfully distinguished brain regions involved in seizure onset from the propagation network. Moreover, cHFEI outperformed other existing detection methods in the estimation of EZs in all patients, with an average area under the ROC curve of 0.88 and an accuracy of 0.85. CONCLUSIONS cHFEI can characterize EZ in a robust manner despite various seizure onset patterns and has potential application in epilepsy therapy.
Collapse
Affiliation(s)
- Zhaohui Li
- School of Information Science and Engineering, Yanshan University, Qinhuangdao 066004, China; Hebei Key Laboratory of information transmission and signal processing, Yanshan University, Qinhuangdao 066004, China.
| | - Hao Zhang
- School of Information Science and Engineering, Yanshan University, Qinhuangdao 066004, China
| | - Shipeng Niu
- School of Information Science and Engineering, Yanshan University, Qinhuangdao 066004, China
| | - Yanyu Xing
- School of Information Science and Engineering, Yanshan University, Qinhuangdao 066004, China
| |
Collapse
|
7
|
Stöber TM, Batulin D, Triesch J, Narayanan R, Jedlicka P. Degeneracy in epilepsy: multiple routes to hyperexcitable brain circuits and their repair. Commun Biol 2023; 6:479. [PMID: 37137938 PMCID: PMC10156698 DOI: 10.1038/s42003-023-04823-0] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2022] [Accepted: 04/06/2023] [Indexed: 05/05/2023] Open
Abstract
Due to its complex and multifaceted nature, developing effective treatments for epilepsy is still a major challenge. To deal with this complexity we introduce the concept of degeneracy to the field of epilepsy research: the ability of disparate elements to cause an analogous function or malfunction. Here, we review examples of epilepsy-related degeneracy at multiple levels of brain organisation, ranging from the cellular to the network and systems level. Based on these insights, we outline new multiscale and population modelling approaches to disentangle the complex web of interactions underlying epilepsy and to design personalised multitarget therapies.
Collapse
Affiliation(s)
- Tristan Manfred Stöber
- Frankfurt Institute for Advanced Studies, 60438, Frankfurt am Main, Germany
- Institute for Neural Computation, Faculty of Computer Science, Ruhr University Bochum, 44801, Bochum, Germany
- Epilepsy Center Frankfurt Rhine-Main, Department of Neurology, Goethe University, 60590, Frankfurt, Germany
| | - Danylo Batulin
- Frankfurt Institute for Advanced Studies, 60438, Frankfurt am Main, Germany
- CePTER - Center for Personalized Translational Epilepsy Research, Goethe University, 60590, Frankfurt, Germany
- Faculty of Computer Science and Mathematics, Goethe University, 60486, Frankfurt, Germany
| | - Jochen Triesch
- Frankfurt Institute for Advanced Studies, 60438, Frankfurt am Main, Germany
| | - Rishikesh Narayanan
- Cellular Neurophysiology Laboratory, Molecular Biophysics Unit, Indian Institute of Science, Bangalore, 560012, India
| | - Peter Jedlicka
- ICAR3R - Interdisciplinary Centre for 3Rs in Animal Research, Faculty of Medicine, Justus Liebig University Giessen, 35390, Giessen, Germany.
- Institute of Clinical Neuroanatomy, Neuroscience Center, Goethe University, 60590, Frankfurt am Main, Germany.
| |
Collapse
|
8
|
Fan D, Wu H, Luan G, Wang Q. The distribution and heterogeneity of excitability in focal epileptic network potentially contribute to the seizure propagation. Front Psychiatry 2023; 14:1137704. [PMID: 36998622 PMCID: PMC10043226 DOI: 10.3389/fpsyt.2023.1137704] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 01/04/2023] [Accepted: 02/21/2023] [Indexed: 03/17/2023] Open
Abstract
IntroductionExisting dynamical models can explain the transmigration mechanisms involved in seizures but are limited to a single modality. Combining models with networks can reproduce scaled epileptic dynamics. And the structure and coupling interactions of the network, as well as the heterogeneity of both the node and network activities, may influence the final state of the network model.MethodsWe built a fully connected network with focal nodes prominently interacting and established a timescale separated epileptic network model. The factors affecting epileptic network seizure were explored by varying the connectivity patterns of focal network nodes and modulating the distribution of network excitability.ResultsThe whole brain network topology as the brain activity foundation affects the consistent delayed clustering seizure propagation. In addition, the network size and distribution heterogeneity of the focal excitatory nodes can influence seizure frequency. With the increasing of the network size and averaged excitability level of focal network, the seizure period decreases. In contrast, the larger heterogeneity of excitability for focal network nodes can lower the functional activity level (average degree) of focal network. There are also subtle effects of focal network topologies (connection patterns of excitatory nodes) that cannot be ignored along with non-focal nodes.DiscussionUnraveling the role of excitatory factors in seizure onset and propagation can be used to understand the dynamic mechanisms and neuromodulation of epilepsy, with profound implications for the treatment of epilepsy and even for the understanding of the brain.
Collapse
Affiliation(s)
- Denggui Fan
- School of Mathematics and Physics, University of Science and Technology Beijing, Beijing, China
| | - Hongyu Wu
- School of Mathematics and Physics, University of Science and Technology Beijing, Beijing, China
| | - Guoming Luan
- Epilepsy Center, Sanbo Brain Hospital, Capital Medical University, Beijing, China
- *Correspondence: Guoming Luan, ; Qingyun Wang,
| | - Qingyun Wang
- Department of Dynamics and Control, Beihang University, Beijing, China
- *Correspondence: Guoming Luan, ; Qingyun Wang,
| |
Collapse
|
9
|
Hashemi M, Vattikonda AN, Jha J, Sip V, Woodman MM, Bartolomei F, Jirsa VK. Amortized Bayesian inference on generative dynamical network models of epilepsy using deep neural density estimators. Neural Netw 2023; 163:178-194. [PMID: 37060871 DOI: 10.1016/j.neunet.2023.03.040] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2022] [Revised: 03/24/2023] [Accepted: 03/30/2023] [Indexed: 04/03/2023]
Abstract
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.
Collapse
|
10
|
Moosavi SA, Truccolo W. Criticality in probabilistic models of spreading dynamics in brain networks: Epileptic seizures. PLoS Comput Biol 2023; 19:e1010852. [PMID: 36749796 PMCID: PMC9904505 DOI: 10.1371/journal.pcbi.1010852] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2022] [Accepted: 01/05/2023] [Indexed: 02/08/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- S Amin Moosavi
- Department of Neuroscience, Brown University, Providence, Rhode Island, United States of America
| | - Wilson Truccolo
- Department of Neuroscience, Brown University, Providence, Rhode Island, United States of America
- Carney Institute for Brain Science, Brown University, Providence, Rhode Island, United States of America
- * E-mail:
| |
Collapse
|
11
|
Wang HE, Woodman M, Triebkorn P, Lemarechal JD, Jha J, Dollomaja B, Vattikonda AN, Sip V, Medina Villalon S, Hashemi M, Guye M, Makhalova J, Bartolomei F, Jirsa V. Delineating epileptogenic networks using brain imaging data and personalized modeling in drug-resistant epilepsy. Sci Transl Med 2023; 15:eabp8982. [PMID: 36696482 DOI: 10.1126/scitranslmed.abp8982] [Citation(s) in RCA: 13] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
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.
Collapse
Affiliation(s)
- Huifang E Wang
- Aix-Marseille Université, Institut National de la Santé et de la Recherche Médicale, Institut de Neurosciences des Systèmes (INS) UMR1106, Marseille 13005, France
| | - Marmaduke Woodman
- Aix-Marseille Université, Institut National de la Santé et de la Recherche Médicale, Institut de Neurosciences des Systèmes (INS) UMR1106, Marseille 13005, France
| | - Paul Triebkorn
- Aix-Marseille Université, Institut National de la Santé et de la Recherche Médicale, Institut de Neurosciences des Systèmes (INS) UMR1106, Marseille 13005, France
| | - Jean-Didier Lemarechal
- Aix-Marseille Université, Institut National de la Santé et de la Recherche Médicale, Institut de Neurosciences des Systèmes (INS) UMR1106, Marseille 13005, France.,Sorbonne Université, Institut du Cerveau - Paris Brain Institute - ICM, Inserm, CNRS, Centre MEG-EEG and Experimental Neurosurgery team, Paris F-75013, France
| | - Jayant Jha
- Aix-Marseille Université, Institut National de la Santé et de la Recherche Médicale, Institut de Neurosciences des Systèmes (INS) UMR1106, Marseille 13005, France
| | - Borana Dollomaja
- Aix-Marseille Université, Institut National de la Santé et de la Recherche Médicale, Institut de Neurosciences des Systèmes (INS) UMR1106, Marseille 13005, France
| | - Anirudh Nihalani Vattikonda
- Aix-Marseille Université, Institut National de la Santé et de la Recherche Médicale, Institut de Neurosciences des Systèmes (INS) UMR1106, Marseille 13005, France
| | - Viktor Sip
- Aix-Marseille Université, Institut National de la Santé et de la Recherche Médicale, Institut de Neurosciences des Systèmes (INS) UMR1106, Marseille 13005, France
| | - Samuel Medina Villalon
- Aix-Marseille Université, Institut National de la Santé et de la Recherche Médicale, Institut de Neurosciences des Systèmes (INS) UMR1106, Marseille 13005, France.,APHM, Epileptology and Clinical Neurophysiology Department, Timone Hospital, Marseille 13005, France
| | - Meysam Hashemi
- Aix-Marseille Université, Institut National de la Santé et de la Recherche Médicale, Institut de Neurosciences des Systèmes (INS) UMR1106, Marseille 13005, France
| | - Maxime Guye
- Aix-Marseille Université, CNRS, CRMBM, Marseille 13005, France.,APHM, Timone University Hospital, CEMEREM, Marseille 13005, France
| | - Julia Makhalova
- APHM, Epileptology and Clinical Neurophysiology Department, Timone Hospital, Marseille 13005, France.,Aix-Marseille Université, CNRS, CRMBM, Marseille 13005, France.,APHM, Timone University Hospital, CEMEREM, Marseille 13005, France
| | - Fabrice Bartolomei
- Aix-Marseille Université, Institut National de la Santé et de la Recherche Médicale, Institut de Neurosciences des Systèmes (INS) UMR1106, Marseille 13005, France.,APHM, Epileptology and Clinical Neurophysiology Department, Timone Hospital, Marseille 13005, France
| | - Viktor Jirsa
- Aix-Marseille Université, Institut National de la Santé et de la Recherche Médicale, Institut de Neurosciences des Systèmes (INS) UMR1106, Marseille 13005, France
| |
Collapse
|
12
|
Lagarde S, Bénar CG, Wendling F, Bartolomei F. Interictal Functional Connectivity in Focal Refractory Epilepsies Investigated by Intracranial EEG. Brain Connect 2022; 12:850-869. [PMID: 35972755 PMCID: PMC9807250 DOI: 10.1089/brain.2021.0190] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Stanislas Lagarde
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France.,Department of Epileptology and Cerebral Rythmology, APHM, Timone Hospital, Marseille, France.,Address correspondence to: Stanislas Lagarde, Department of Epileptology and Cerebral Rythmology, APHM, Timone Hospital, 264 Rue Saint-Pierre, 13005 Marseille, France
| | | | | | - Fabrice Bartolomei
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France.,Department of Epileptology and Cerebral Rythmology, APHM, Timone Hospital, Marseille, France
| |
Collapse
|
13
|
Brogin JAF, Faber J, Bueno DD. Estimating the Parameters of the Epileptor Model for Epileptic Seizure Suppression. Neuroinformatics 2022; 20:919-941. [PMID: 35303252 DOI: 10.1007/s12021-022-09583-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/10/2022] [Indexed: 12/31/2022]
Abstract
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.
Collapse
Affiliation(s)
- João Angelo Ferres Brogin
- Department of Mechanical Engineering, São Paulo State University (UNESP), 56 Brasil Avenue, Ilha Solteira, 15385-000, São Paulo, Brazil.
| | - Jean Faber
- Department of Neurology and Neurosurgery, Federal University of São Paulo, 667 Pedro de Toledo Street, São Paulo, 04039-032, São Paulo, Brazil
| | - Douglas D Bueno
- Department of Mathematics, São Paulo State University (UNESP), 56 Brasil Avenue, Ilha Solteira, 15385-000, São Paulo, Brazil
| |
Collapse
|
14
|
Jha J, Hashemi M, Vattikonda AN, Wang H, Jirsa V. Fully Bayesian estimation of virtual brain parameters with self-tuning Hamiltonian Monte Carlo. MACHINE LEARNING: SCIENCE AND TECHNOLOGY 2022. [DOI: 10.1088/2632-2153/ac9037] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022] Open
Abstract
Abstract
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.
Collapse
|
15
|
Stanslaski S, Farooqi H, Sanabria DE, Netoff TI. Fully Closed Loop Test Environment for Adaptive Implantable Neural Stimulators Using Computational Models. J Med Device 2022; 16:034501. [PMID: 35646224 PMCID: PMC9125865 DOI: 10.1115/1.4054083] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2021] [Revised: 02/26/2022] [Indexed: 07/22/2023] Open
Abstract
Implantable brain stimulation devices continue to be developed to treat and monitor brain conditions. As the complexity of these devices grows to include adaptive neuromodulation therapy, validating the operation and verifying the correctness of these systems becomes more complicated. The new complexities lie in the functioning of the device being dependent on the interaction with the patient and environmental factors such as noise and artifacts. Here, we present a hardware-in-the-loop (HIL) testing framework that employs computational models of pathological neural dynamics to test adaptive deep brain stimulation (DBS) devices prior to animal or human testing. A brain stimulation and recording electrode array is placed in the saline tank and connected to an adaptive neuromodulation system that measures and processes the synthetic signals and delivers stimulation back into the saline tank. A data acquisition system is used to detect the stimulation and provide feedback to the computational model in order to simulate the effects of stimulation on the neural dynamics. In this study, we used real-time computational models to emulate the dynamics of epileptic seizures observed in the anterior nucleus of the thalamus (ANT) in epilepsy patients and beta band (11-35 Hz) oscillations observed in the subthalamic nucleus (STN) of Parkinson's disease (PD) patients. These models simulated neuronal responses to electrical stimulation pulses and the saline tank tested hardware interactions between the detection algorithms and stimulation interference. We tested and validated the operation of adaptive DBS algorithms for seizure and beta band power suppression embedded in an implantable DBS system (Medtronic Summit RC+S). This study highlights the utility of the proposed hardware-in-the-loop framework to systematically test the adaptive DBS systems in the presence of system aggressors such as environmental noise and stimulation-induced electrical artifacts. This testing procedure can help ensure correctness and robustness of adaptive DBS devices prior to animal and human testing.
Collapse
Affiliation(s)
- Scott Stanslaski
- Department of Biomedical Engineering, University of Minnesota, Minneapolis, MN 55455; Neuromodulation Department, Medtronic PLC, Minneapolis, MN 55432
| | - Hafsa Farooqi
- Department of Neurology, University of Minnesota, Minneapolis, MN 55455
| | | | - Theoden Ivan Netoff
- Department of Biomedical Engineering, University of Minnesota, Minneapolis, MN 55455
| |
Collapse
|
16
|
Moosavi SA, Jirsa VK, Truccolo W. Critical dynamics in the spread of focal epileptic seizures: Network connectivity, neural excitability and phase transitions. PLoS One 2022; 17:e0272902. [PMID: 35998146 PMCID: PMC9397939 DOI: 10.1371/journal.pone.0272902] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2021] [Accepted: 07/29/2022] [Indexed: 11/24/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- S. Amin Moosavi
- Department of Neuroscience, Brown University, Providence, RI, United States of America
| | - Viktor K. Jirsa
- Aix Marseille University, INSERM, INS, Institut de Neurosciences de Système, Marseille, France
| | - Wilson Truccolo
- Department of Neuroscience, Brown University, Providence, RI, United States of America
- Carney Institute for Brain Science, Brown University, Providence, RI, United States of America
- * E-mail:
| |
Collapse
|
17
|
Tu T, Peng Z, Ren J, Zhang H. Cerebral Cavernous Malformation: Immune and Inflammatory Perspectives. Front Immunol 2022; 13:922281. [PMID: 35844490 PMCID: PMC9280619 DOI: 10.3389/fimmu.2022.922281] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2022] [Accepted: 05/18/2022] [Indexed: 12/03/2022] Open
Abstract
Cerebral cavernous malformation (CCM) is a type of vascular anomaly that arises due to the dyshomeostasis of brain capillary networks. In the past two decades, many advances have been made in this research field. Notably, as a more reasonable current view, the CCM lesions should be attributed to the results of a great number of additional events related to the homeostasis disorder of the endothelial cell. Indeed, one of the most fascinating concerns in the research field is the inflammatory perturbation in the immune microenvironment, which would affect the disease progression as well as the patients’ outcomes. In this work, we focused on this topic, and underlined the immune-related factors’ contribution to the CCM pathologic progression.
Collapse
Affiliation(s)
- Tianqi Tu
- Department of Neurosurgery, Xuanwu Hospital, Capital Medical University, Beijing, China
- China International Neuroscience Institute (China-INI), Xuanwu Hospital, Capital Medical University, Beijing, China
| | - Zhenghong Peng
- Health Management Department, The Affiliated Hospital of Southwest Medical University, Luzhou, China
| | - Jian Ren
- Department of Neurosurgery, Xuanwu Hospital, Capital Medical University, Beijing, China
- China International Neuroscience Institute (China-INI), Xuanwu Hospital, Capital Medical University, Beijing, China
| | - Hongqi Zhang
- Department of Neurosurgery, Xuanwu Hospital, Capital Medical University, Beijing, China
- China International Neuroscience Institute (China-INI), Xuanwu Hospital, Capital Medical University, Beijing, China
- *Correspondence: Hongqi Zhang,
| |
Collapse
|
18
|
Wang C, Chen S, Huang L, Yu L. Prediction and control of focal seizure spread: Random walk with restart on heterogeneous brain networks. Phys Rev E 2022; 105:064412. [PMID: 35854502 DOI: 10.1103/physreve.105.064412] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2022] [Accepted: 05/25/2022] [Indexed: 06/15/2023]
Abstract
Whole-brain models offer a promising method of predicting seizure spread, which is critical for successful surgical treatment of focal epilepsy. Existing methods are largely based on structural connectome, which ignores the effects of heterogeneity within the regional excitability of brains. In this study we used a whole-brain model to show that heterogeneity in nodal excitability had a significant impact on seizure propagation in the networks and compromised the prediction accuracy with structural connections. We then addressed this problem with an algorithm based on random walk with restart on graphs. We demonstrated that by establishing a relationship between the restarting probability and the excitability for each node, this algorithm could significantly improve the seizure spread prediction accuracy in heterogeneous networks and was more robust against the extent of heterogeneity. We also strategized surgical seizure control as a process to identify and remove the key nodes (connections) responsible for the early spread of seizures from the focal region. Compared to strategies based on structural connections, virtual surgery with a strategy based on a modified random walk with extended restart generated outcomes with a high success rate while maintaining low damage to the brain by removing fewer anatomical connections. These findings may have potential applications in developing personalized surgery strategies for epilepsy.
Collapse
Affiliation(s)
- Chen Wang
- School of Physical Science and Technology, Lanzhou University, Lanzhou 730000, China
| | - Sida Chen
- School of Physical Science and Technology, Lanzhou University, Lanzhou 730000, China
| | - Liang Huang
- School of Physical Science and Technology, Lanzhou University, Lanzhou 730000, China
- Lanzhou Center for Theoretical Physics, Key Laboratory of Theoretical Physics of Gansu Province, Lanzhou University, Lanzhou, Gansu 730000, China
| | - Lianchun Yu
- School of Physical Science and Technology, Lanzhou University, Lanzhou 730000, China
- Lanzhou Center for Theoretical Physics, Key Laboratory of Theoretical Physics of Gansu Province, Lanzhou University, Lanzhou, Gansu 730000, China
| |
Collapse
|
19
|
Armonaite K, Bertoli M, Paulon L, Gianni E, Balsi M, Conti L, Tecchio F. Neuronal Electrical Ongoing Activity as Cortical Areas Signature: An Insight from MNI Intracerebral Recording Atlas. Cereb Cortex 2021; 32:2895-2906. [PMID: 34727186 DOI: 10.1093/cercor/bhab389] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2021] [Revised: 09/28/2021] [Accepted: 09/29/2021] [Indexed: 11/14/2022] Open
Abstract
The time course of the neuronal activity in the brain network, the neurodynamics, reflects the structure and functionality of the generating neuronal pools. Here, using the intracranial stereo-electroencephalographic (sEEG) recordings of the public Montreal Neurological Institute (MNI) atlas, we investigated the neurodynamics of primary motor (M1), somatosensory (S1) and auditory (A1) cortices measuring power spectral densities (PSD) and Higuchi fractal dimension (HFD) in the same subject (M1 vs. S1 in 16 subjects, M1 vs. A1 in 9, S1 vs. A1 in 6). We observed specific spectral features in M1, which prevailed above beta band, S1 in the alpha band, and A1 in the delta band. M1 HFD was higher than S1, both higher than A1. A clear distinction of neurodynamics properties of specific primary cortices supports the efforts in cortical parceling based on this expression of the local cytoarchitecture and connectivity. In this perspective, we selected within the MNI intracortical database a first set of primary motor, somatosensory and auditory cortices' representatives to query in recognizing ongoing patterns of neuronal communication. Potential clinical impact stands primarily in exploiting such exchange patterns to enhance the efficacy of neuromodulation intervention to cure symptoms secondary to neuronal activity unbalances.
Collapse
Affiliation(s)
| | - Massimo Bertoli
- Laboratory of Electrophysiology for Translational NeuroScience (LET'S), Institute of Cognitive Sciences and Technologies - Consiglio Nazionale delle Ricerche, Rome 00185, Italy.,Department of Neuroscience, Imaging and Clinical Sciences, University 'Gabriele D'Annunzio' of Chieti-Pescara, Chieti 66100, Italy
| | - Luca Paulon
- Laboratory of Electrophysiology for Translational NeuroScience (LET'S), Institute of Cognitive Sciences and Technologies - Consiglio Nazionale delle Ricerche, Rome 00185, Italy
| | - Eugenia Gianni
- Laboratory of Electrophysiology for Translational NeuroScience (LET'S), Institute of Cognitive Sciences and Technologies - Consiglio Nazionale delle Ricerche, Rome 00185, Italy.,Unit of Neurology, Neurophysiology, Neurobiology, Department of Medicine, Università Campus Bio-Medico di Roma, Rome 00128, Italy
| | - Marco Balsi
- Department of Information Engineering, Electronics and Telecommunications, Sapienza University, Rome 00185, Italy
| | - Livio Conti
- Faculty of Engineering, Uninettuno University, Rome 00186, Italy.,INFN - Istituto Nazionale di Fisica Nucleare, Sezione Roma Tor Vergata, Rome 00133, Italy
| | - Franca Tecchio
- Laboratory of Electrophysiology for Translational NeuroScience (LET'S), Institute of Cognitive Sciences and Technologies - Consiglio Nazionale delle Ricerche, Rome 00185, Italy
| |
Collapse
|
20
|
Vattikonda AN, Hashemi M, Sip V, Woodman MM, Bartolomei F, Jirsa VK. Identifying spatio-temporal seizure propagation patterns in epilepsy using Bayesian inference. Commun Biol 2021; 4:1244. [PMID: 34725441 PMCID: PMC8560929 DOI: 10.1038/s42003-021-02751-5] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2020] [Accepted: 10/04/2021] [Indexed: 01/24/2023] Open
Abstract
Focal drug resistant epilepsy is a neurological disorder characterized by seizures caused by abnormal activity originating in one or more regions together called as epileptogenic zone. Treatment for such patients involves surgical resection of affected regions. Epileptogenic zone is typically identified using stereotactic EEG recordings from the electrodes implanted into the patient's brain. Identifying the epileptogenic zone is a challenging problem due to the spatial sparsity of electrode implantation. We propose a probabilistic hierarchical model of seizure propagation patterns, based on a phenomenological model of seizure dynamics called Epileptor. Using Bayesian inference, the Epileptor model is optimized to build patient specific virtual models that best fit to the log power of intracranial recordings. First, accuracy of the model predictions and identifiability of the model are investigated using synthetic data. Then, model predictions are evaluated against a retrospective patient cohort of 25 patients with varying surgical outcomes. In the patients who are seizure free after surgery, model predictions showed good match with the clinical hypothesis. In patients where surgery failed to achieve seizure freedom model predictions showed a strong mismatch. Our results demonstrate that proposed probabilistic model could be a valuable tool to aid the clinicians in identifying the seizure focus.
Collapse
Affiliation(s)
- Anirudh N Vattikonda
- Aix Marseille Univ, INSERM, INS, Institut de Neurosciences des Systèmes, Marseille, France
| | - Meysam Hashemi
- Aix Marseille Univ, INSERM, INS, Institut de Neurosciences des Systèmes, Marseille, France
| | - Viktor Sip
- Aix Marseille Univ, INSERM, INS, Institut de Neurosciences des Systèmes, Marseille, France
| | - Marmaduke M Woodman
- Aix Marseille Univ, INSERM, INS, Institut de Neurosciences des Systèmes, Marseille, France
| | - Fabrice Bartolomei
- Aix Marseille Univ, INSERM, INS, Institut de Neurosciences des Systèmes, Marseille, France
- Epileptology Department and Clinical Neurophysiology Department, Assistance publique des Hopitaux de Marseille, Marseille, France
| | - Viktor K Jirsa
- Aix Marseille Univ, INSERM, INS, Institut de Neurosciences des Systèmes, Marseille, France.
| |
Collapse
|
21
|
Gerster M, Taher H, Škoch A, Hlinka J, Guye M, Bartolomei F, Jirsa V, Zakharova A, Olmi S. Patient-Specific Network Connectivity Combined With a Next Generation Neural Mass Model to Test Clinical Hypothesis of Seizure Propagation. Front Syst Neurosci 2021; 15:675272. [PMID: 34539355 PMCID: PMC8440880 DOI: 10.3389/fnsys.2021.675272] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2021] [Accepted: 07/07/2021] [Indexed: 11/13/2022] Open
Abstract
Dynamics underlying epileptic seizures span multiple scales in space and time, therefore, understanding seizure mechanisms requires identifying the relations between seizure components within and across these scales, together with the analysis of their dynamical repertoire. In this view, mathematical models have been developed, ranging from single neuron to neural population. In this study, we consider a neural mass model able to exactly reproduce the dynamics of heterogeneous spiking neural networks. We combine mathematical modeling with structural information from non invasive brain imaging, thus building large-scale brain network models to explore emergent dynamics and test the clinical hypothesis. We provide a comprehensive study on the effect of external drives on neuronal networks exhibiting multistability, in order to investigate the role played by the neuroanatomical connectivity matrices in shaping the emergent dynamics. In particular, we systematically investigate the conditions under which the network displays a transition from a low activity regime to a high activity state, which we identify with a seizure-like event. This approach allows us to study the biophysical parameters and variables leading to multiple recruitment events at the network level. We further exploit topological network measures in order to explain the differences and the analogies among the subjects and their brain regions, in showing recruitment events at different parameter values. We demonstrate, along with the example of diffusion-weighted magnetic resonance imaging (dMRI) connectomes of 20 healthy subjects and 15 epileptic patients, that individual variations in structural connectivity, when linked with mathematical dynamic models, have the capacity to explain changes in spatiotemporal organization of brain dynamics, as observed in network-based brain disorders. In particular, for epileptic patients, by means of the integration of the clinical hypotheses on the epileptogenic zone (EZ), i.e., the local network where highly synchronous seizures originate, we have identified the sequence of recruitment events and discussed their links with the topological properties of the specific connectomes. The predictions made on the basis of the implemented set of exact mean-field equations turn out to be in line with the clinical pre-surgical evaluation on recruited secondary networks.
Collapse
Affiliation(s)
- Moritz Gerster
- Institut für Theoretische Physik, Technische Universität Berlin, Berlin, Germany
| | - Halgurd Taher
- Inria Sophia Antipolis Méditerranée Research Centre, MathNeuro Team, Valbonne, France
| | - Antonín Škoch
- National Institute of Mental Health, Klecany, Czechia
- MR Unit, Department of Diagnostic and Interventional Radiology, Institute for Clinical and Experimental Medicine, Prague, Czechia
| | - Jaroslav Hlinka
- National Institute of Mental Health, Klecany, Czechia
- Institute of Computer Science of the Czech Academy of Sciences, Prague, Czechia
| | - Maxime Guye
- Faculté de Médecine de la Timone, Centre de Résonance Magnétique et Biologique et Médicale (CRMBM, UMR CNRS-AMU 7339), Medical School of Marseille, Aix-Marseille Université, Marseille, France
- Assistance Publique -Hôpitaux de Marseille, Hôpital de la Timone, Pôle d'Imagerie, Marseille, France
| | - Fabrice Bartolomei
- Assistance Publique - Hôpitaux de Marseille, Hôpital de la Timone, Service de Neurophysiologie Clinique, Marseille, France
| | - Viktor Jirsa
- Aix Marseille Université, Inserm, Institut de Neurosciences des Systèmes, UMRS 1106, Marseille, France
| | - Anna Zakharova
- Institut für Theoretische Physik, Technische Universität Berlin, Berlin, Germany
| | - Simona Olmi
- Inria Sophia Antipolis Méditerranée Research Centre, MathNeuro Team, Valbonne, France
- Consiglio Nazionale delle Ricerche, Istituto dei Sistemi Complessi, Sesto Fiorentino, Italy
| |
Collapse
|
22
|
Hashemi M, Vattikonda AN, Sip V, Diaz-Pier S, Peyser A, Wang H, Guye M, Bartolomei F, Woodman MM, Jirsa VK. On the influence of prior information evaluated by fully Bayesian criteria in a personalized whole-brain model of epilepsy spread. PLoS Comput Biol 2021; 17:e1009129. [PMID: 34260596 PMCID: PMC8312957 DOI: 10.1371/journal.pcbi.1009129] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2020] [Revised: 07/26/2021] [Accepted: 05/29/2021] [Indexed: 11/18/2022] Open
Abstract
Individualized anatomical information has been used as prior knowledge in Bayesian inference paradigms of whole-brain network models. However, the actual sensitivity to such personalized information in priors is still unknown. In this study, we introduce the use of fully Bayesian information criteria and leave-one-out cross-validation technique on the subject-specific information to assess different epileptogenicity hypotheses regarding the location of pathological brain areas based on a priori knowledge from dynamical system properties. The Bayesian Virtual Epileptic Patient (BVEP) model, which relies on the fusion of structural data of individuals, a generative model of epileptiform discharges, and a self-tuning Monte Carlo sampling algorithm, is used to infer the spatial map of epileptogenicity across different brain areas. Our results indicate that measuring the out-of-sample prediction accuracy of the BVEP model with informative priors enables reliable and efficient evaluation of potential hypotheses regarding the degree of epileptogenicity across different brain regions. In contrast, while using uninformative priors, the information criteria are unable to provide strong evidence about the epileptogenicity of brain areas. We also show that the fully Bayesian criteria correctly assess different hypotheses about both structural and functional components of whole-brain models that differ across individuals. The fully Bayesian information-theory based approach used in this study suggests a patient-specific strategy for epileptogenicity hypothesis testing in generative brain network models of epilepsy to improve surgical outcomes.
Collapse
Affiliation(s)
- Meysam Hashemi
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France
| | | | - Viktor Sip
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France
| | - Sandra Diaz-Pier
- SimLab Neuroscience, Jülich Supercomputing Centre (JSC), Institute for Advanced Simulation, JARA, Forschungszentrum Jülich GmbH, Jülich, Germany
| | - Alexander Peyser
- SimLab Neuroscience, Jülich Supercomputing Centre (JSC), Institute for Advanced Simulation, JARA, Forschungszentrum Jülich GmbH, Jülich, Germany
- Google, München, Germany
| | - Huifang Wang
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France
| | - Maxime Guye
- Aix Marseille Univ, CNRS, CRMBM, Marseille, France
| | - Fabrice Bartolomei
- Epileptology Department, and Clinical Neurophysiology Department, Assistance Publique des Hôpitaux de Marseille, Marseille, France
| | | | - Viktor K. Jirsa
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France
| |
Collapse
|
23
|
Depannemaecker D, Destexhe A, Jirsa V, Bernard C. Modeling seizures: From single neurons to networks. Seizure 2021; 90:4-8. [PMID: 34219016 DOI: 10.1016/j.seizure.2021.06.015] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2021] [Revised: 06/11/2021] [Accepted: 06/11/2021] [Indexed: 11/26/2022] Open
Abstract
Dynamical system tools offer a complementary approach to detailed biophysical seizure modeling, with a high potential for clinical applications. This review describes the theoretical framework that provides a basis for theorizing certain properties of seizures and for their classification according to their dynamical properties at onset and offset. We describe various modeling approaches spanning different scales, from single neurons to large-scale networks. This narrative review provides an accessible overview of this field, including non-exhaustive examples of key recent works.
Collapse
Affiliation(s)
- Damien Depannemaecker
- Paris-Saclay University, French National Centre for Scientific Research (CNRS), Institute of Neuroscience (NeuroPSI), 91198 Gif sur Yvette, France.
| | - Alain Destexhe
- Paris-Saclay University, French National Centre for Scientific Research (CNRS), Institute of Neuroscience (NeuroPSI), 91198 Gif sur Yvette, France.
| | - Viktor Jirsa
- Aix Marseille Univ, INSERM, INS, Institut des Neurosciences des Systèmes, Marseille, France.
| | - Christophe Bernard
- Aix Marseille Univ, INSERM, INS, Institut des Neurosciences des Systèmes, Marseille, France.
| |
Collapse
|
24
|
Pérez-Cervera A, Hlinka J. Perturbations both trigger and delay seizures due to generic properties of slow-fast relaxation oscillators. PLoS Comput Biol 2021; 17:e1008521. [PMID: 33780437 PMCID: PMC8032201 DOI: 10.1371/journal.pcbi.1008521] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2020] [Revised: 04/08/2021] [Accepted: 02/22/2021] [Indexed: 01/24/2023] Open
Abstract
The mechanisms underlying the emergence of seizures are one of the most important unresolved issues in epilepsy research. In this paper, we study how perturbations, exogenous or endogenous, may promote or delay seizure emergence. To this aim, due to the increasingly adopted view of epileptic dynamics in terms of slow-fast systems, we perform a theoretical analysis of the phase response of a generic relaxation oscillator. As relaxation oscillators are effectively bistable systems at the fast time scale, it is intuitive that perturbations of the non-seizing state with a suitable direction and amplitude may cause an immediate transition to seizure. By contrast, and perhaps less intuitively, smaller amplitude perturbations have been found to delay the spontaneous seizure initiation. By studying the isochrons of relaxation oscillators, we show that this is a generic phenomenon, with the size of such delay depending on the slow flow component. Therefore, depending on perturbation amplitudes, frequency and timing, a train of perturbations causes an occurrence increase, decrease or complete suppression of seizures. This dependence lends itself to analysis and mechanistic understanding through methods outlined in this paper. We illustrate this methodology by computing the isochrons, phase response curves and the response to perturbations in several epileptic models possessing different slow vector fields. While our theoretical results are applicable to any planar relaxation oscillator, in the motivating context of epilepsy they elucidate mechanisms of triggering and abating seizures, thus suggesting stimulation strategies with effects ranging from mere delaying to full suppression of seizures. Despite its simplicity, the modelling of epileptic dynamics as a slow-fast transition between low and high activity states mediated by some slow feedback variable is a relatively novel albeit fruitful approach. This study is the first, to our knowledge, characterizing the response of such slow-fast models of epileptic brain to perturbations by computing its isochrons. Besides its numerical computation, we theoretically determine which factors shape the geometry of isochrons for planar slow-fast oscillators. As a consequence, we introduce a theoretical approach providing a clear understanding of the response of perturbations of slow-fast oscillators. Within the epilepsy context, this elucidates the origin of the contradictory role of interictal epileptiform discharges in the transition to seizure, manifested by both pro-convulsive and anti-convulsive effect depending on the amplitude, frequency and timing. More generally, this paper provides theoretical framework highlighting the role of the slow flow component on the response to perturbations in relaxation oscillators, pointing to the general phenomena in such slow-fast oscillators ubiquitous in biological systems.
Collapse
Affiliation(s)
- Alberto Pérez-Cervera
- Department of Complex Systems, Institute of Computer Science of the Czech Academy of Sciences, Prague, Czech Republic
- Center for Cognition and Decision Making, Institute for Cognitive Neuroscience, National Research University Higher School of Economics, Moscow, Russia
- * E-mail: (AP); (JH)
| | - Jaroslav Hlinka
- Department of Complex Systems, Institute of Computer Science of the Czech Academy of Sciences, Prague, Czech Republic
- National Institute of Mental Health, Klecany, Czech Republic
- * E-mail: (AP); (JH)
| |
Collapse
|
25
|
Spiegler A, Abadchi JK, Mohajerani M, Jirsa VK. In silico exploration of mouse brain dynamics by focal stimulation reflects the organization of functional networks and sensory processing. Netw Neurosci 2021; 4:807-851. [PMID: 33615092 PMCID: PMC7888484 DOI: 10.1162/netn_a_00152] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2019] [Accepted: 06/08/2020] [Indexed: 12/12/2022] Open
Abstract
Resting-state functional networks such as the default mode network (DMN) dominate spontaneous brain dynamics. To date, the mechanisms linking brain structure and brain dynamics and functions in cognition, perception, and action remain unknown, mainly due to the uncontrolled and erratic nature of the resting state. Here we used a stimulation paradigm to probe the brain’s resting behavior, providing insights on state-space stability and multiplicity of network trajectories after stimulation. We performed explorations on a mouse model to map spatiotemporal brain dynamics as a function of the stimulation site. We demonstrated the emergence of known functional networks in brain responses. Several responses heavily relied on the DMN and were suggestive of the DMN playing a mechanistic role between functional networks. We probed the simulated brain responses to the stimulation of regions along the information processing chains of sensory systems from periphery up to primary sensory cortices. Moreover, we compared simulated dynamics against in vivo brain responses to optogenetic stimulation. Our results underwrite the importance of anatomical connectivity in the functional organization of brain networks and demonstrate how functionally differentiated information processing chains arise from the same system. We demonstrate how functionally differentiated information processing chains arise from the same anatomical network. The main result of the in-silico mouse brain simulations is the emergence of specific functional networks based on structural data from the mouse brain. When the brain is stimulated, for example, by sensory inputs or direct electrical stimulation, the brain initially responds with activities in specific regions. The brain’s anatomical connectivity constrains the subsequent pattern formation. We built a high-resolution mouse brain network model. The model structure originated from experimental data. We systematically explored the mouse model and investigated the simulated brain dynamics after stimulation. Known functional networks emerged in the simulated brain responses. The default mode network occurred in almost all characteristic response patterns. Simulated brain response dynamics and in-vivo response dynamics of the mouse brain to optogenetic stimulation showed similarities even without parameter tuning. Anatomical connectivity and dynamics shape the functional organization of brain networks.
Collapse
Affiliation(s)
- Andreas Spiegler
- Department of Neurophysiology and Pathophysiology, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Javad Karimi Abadchi
- Canadian Center for Behavioural Neuroscience, University of Lethbridge, Alberta, Canada
| | - Majid Mohajerani
- Canadian Center for Behavioural Neuroscience, University of Lethbridge, Alberta, Canada
| | - Viktor K Jirsa
- Institut de Neurosciences des Systèmes, UMR Inserm 1106, Aix-Marseille Université, Faculté de Médecine, Marseille, France
| |
Collapse
|
26
|
Li Z, Li S, Yu T, Li X. Measuring the Coupling Direction between Neural Oscillations with Weighted Symbolic Transfer Entropy. ENTROPY (BASEL, SWITZERLAND) 2020; 22:e22121442. [PMID: 33371251 PMCID: PMC7767336 DOI: 10.3390/e22121442] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/09/2020] [Revised: 12/13/2020] [Accepted: 12/16/2020] [Indexed: 05/30/2023]
Abstract
Neural oscillations reflect rhythmic fluctuations in the synchronization of neuronal populations and play a significant role in neural processing. To further understand the dynamic interactions between different regions in the brain, it is necessary to estimate the coupling direction between neural oscillations. Here, we developed a novel method, termed weighted symbolic transfer entropy (WSTE), that combines symbolic transfer entropy (STE) and weighted probability distribution to measure the directionality between two neuronal populations. The traditional STE ignores the degree of difference between the amplitude values of a time series. In our proposed WSTE method, this information is picked up by utilizing a weighted probability distribution. The simulation analysis shows that the WSTE method can effectively estimate the coupling direction between two neural oscillations. In comparison with STE, the new method is more sensitive to the coupling strength and is more robust against noise. When applied to epileptic electrocorticography data, a significant coupling direction from the anterior nucleus of thalamus (ANT) to the seizure onset zone (SOZ) was detected during seizures. Considering the superiorities of the WSTE method, it is greatly advantageous to measure the coupling direction between neural oscillations and consequently characterize the information flow between different brain regions.
Collapse
Affiliation(s)
- Zhaohui Li
- School of Information Science and Engineering (School of Software), Yanshan University, Qinhuangdao 066004, China; (Z.L.); (S.L.)
- Hebei Key Laboratory of Information Transmission and Signal Processing, Yanshan University, Qinhuangdao 066004, China
| | - Shuaifei Li
- School of Information Science and Engineering (School of Software), Yanshan University, Qinhuangdao 066004, China; (Z.L.); (S.L.)
| | - Tao Yu
- Beijing Institute of Functional Neurosurgery, Capital Medical University, Beijing 100053, China;
| | - Xiaoli Li
- State Key Laboratory of Cognitive Neuroscience and Learning, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
27
|
Köksal Ersöz E, Modolo J, Bartolomei F, Wendling F. Neural mass modeling of slow-fast dynamics of seizure initiation and abortion. PLoS Comput Biol 2020; 16:e1008430. [PMID: 33166277 PMCID: PMC7676664 DOI: 10.1371/journal.pcbi.1008430] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2020] [Revised: 11/19/2020] [Accepted: 10/08/2020] [Indexed: 12/31/2022] Open
Abstract
Epilepsy is a dynamic and complex neurological disease affecting about 1% of the worldwide population, among which 30% of the patients are drug-resistant. Epilepsy is characterized by recurrent episodes of paroxysmal neural discharges (the so-called seizures), which manifest themselves through a large-amplitude rhythmic activity observed in depth-EEG recordings, in particular in local field potentials (LFPs). The signature characterizing the transition to seizures involves complex oscillatory patterns, which could serve as a marker to prevent seizure initiation by triggering appropriate therapeutic neurostimulation methods. To investigate such protocols, neurophysiological lumped-parameter models at the mesoscopic scale, namely neural mass models, are powerful tools that not only mimic the LFP signals but also give insights on the neural mechanisms related to different stages of seizures. Here, we analyze the multiple time-scale dynamics of a neural mass model and explain the underlying structure of the complex oscillations observed before seizure initiation. We investigate population-specific effects of the stimulation and the dependence of stimulation parameters on synaptic timescales. In particular, we show that intermediate stimulation frequencies (>20 Hz) can abort seizures if the timescale difference is pronounced. Those results have the potential in the design of therapeutic brain stimulation protocols based on the neurophysiological properties of tissue.
Collapse
Affiliation(s)
| | - Julien Modolo
- University of Rennes, Inserm-U1099, LTSI, Rennes, France
| | - Fabrice Bartolomei
- Aix Marseille University, Inserm, INS, Institut de Neurosciences des Systèmes, Marseille, France
- APHM, Timone Hospital, Clinical Neurophysiology, Marseille, France
| | | |
Collapse
|
28
|
Wang HE, Scholly J, Triebkorn P, Sip V, Medina Villalon S, Woodman MM, Le Troter A, Guye M, Bartolomei F, Jirsa V. VEP atlas: An anatomic and functional human brain atlas dedicated to epilepsy patients. J Neurosci Methods 2020; 348:108983. [PMID: 33121983 DOI: 10.1016/j.jneumeth.2020.108983] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2020] [Revised: 09/27/2020] [Accepted: 10/18/2020] [Indexed: 01/15/2023]
Abstract
BACKGROUND Several automated parcellation atlases of the human brain have been developed over the past decades, based on various criteria, and have been applied in basic and clinical research. NEW METHOD Here we present the Virtual Epileptic Patient (VEP) atlas that offers a new automated brain region parcellation and labeling, which has been developed for the specific use in the domains of epileptology and functional neurosurgery and is able to apply at individual patient's level. RESULTS It comprises 162 brain regions, including 73 cortical and 8 subcortical regions per hemisphere. We demonstrate the successful application of the VEP atlas in a cohort of 50 retrospective patients. The structural organization is complemented by the functional variation of stereotactic intracerebral EEG (SEEG) signal data features establishing brain region-specific 3d-maps. COMPARISON WITH EXISTING METHODS The VEP atlas integrates both anatomical and functional definitions in the same atlas, adapted to applications for epilepsy patients and individualizable. CONCLUSION The covariation of structural and functional organization is the basis for current efforts of patient-specific large-scale brain network modeling exploiting virtual brain technologies for the identification of the epileptogenic regions in an ongoing prospective clinical trial EPINOV.
Collapse
Affiliation(s)
- Huifang E Wang
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France.
| | - Julia Scholly
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France; Epileptology Department, and Clinical Neurophysiology Department, Assistance Publique des Hôpitaux de Marseille, Marseille, France
| | - Paul Triebkorn
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France
| | - Viktor Sip
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France
| | - Samuel Medina Villalon
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France; Epileptology Department, and Clinical Neurophysiology Department, Assistance Publique des Hôpitaux de Marseille, Marseille, France
| | | | | | - Maxime Guye
- Aix Marseille Univ, CNRS, CRMBM, Marseille, France
| | - Fabrice Bartolomei
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France; Epileptology Department, and Clinical Neurophysiology Department, Assistance Publique des Hôpitaux de Marseille, Marseille, France
| | - Viktor Jirsa
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France.
| |
Collapse
|
29
|
Creaser J, Lin C, Ridler T, Brown JT, D’Souza W, Seneviratne U, Cook M, Terry JR, Tsaneva-Atanasova K. Domino-like transient dynamics at seizure onset in epilepsy. PLoS Comput Biol 2020; 16:e1008206. [PMID: 32986695 PMCID: PMC7544071 DOI: 10.1371/journal.pcbi.1008206] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2019] [Revised: 10/08/2020] [Accepted: 07/29/2020] [Indexed: 12/20/2022] Open
Abstract
The International League Against Epilepsy (ILAE) groups seizures into "focal", "generalized" and "unknown" based on whether the seizure onset is confined to a brain region in one hemisphere, arises in several brain region simultaneously, or is not known, respectively. This separation fails to account for the rich diversity of clinically and experimentally observed spatiotemporal patterns of seizure onset and even less so for the properties of the brain networks generating them. We consider three different patterns of domino-like seizure onset in Idiopathic Generalized Epilepsy (IGE) and present a novel approach to classification of seizures. To understand how these patterns are generated on networks requires understanding of the relationship between intrinsic node dynamics and coupling between nodes in the presence of noise, which currently is unknown. We investigate this interplay here in the framework of domino-like recruitment across a network. In particular, we use a phenomenological model of seizure onset with heterogeneous coupling and node properties, and show that in combination they generate a range of domino-like onset patterns observed in the IGE seizures. We further explore the individual contribution of heterogeneous node dynamics and coupling by interpreting in-vitro experimental data in which the speed of onset can be chemically modulated. This work contributes to a better understanding of possible drivers for the spatiotemporal patterns observed at seizure onset and may ultimately contribute to a more personalized approach to classification of seizure types in clinical practice.
Collapse
Affiliation(s)
- Jennifer Creaser
- Department of Mathematics, University of Exeter, Exeter, EX4 4QF, UK
- EPSRC Centre for Predictive modeling in Healthcare, University of Exeter, Exeter, EX4 4QJ, UK
| | - Congping Lin
- Center for Mathematical Sciences, Huazhong University of Science and Technology, Wuhan, Hubei 430074, China
- Hubei Key Lab of Engineering Modeling and Scientific Computing, Huazhong University of Science and Technology, Wuhan, Hubei 430074, China
| | - Thomas Ridler
- Institute of Biomedical and Clinical Sciences, College of Medicine and Health, University of Exeter, EX4 4PS, UK
| | - Jonathan T. Brown
- Institute of Biomedical and Clinical Sciences, College of Medicine and Health, University of Exeter, EX4 4PS, UK
| | - Wendyl D’Souza
- Department of Medicine, St. Vincent’s Hospital, University of Melbourne, Melbourne, VIC 3065, Australia
| | - Udaya Seneviratne
- Department of Medicine, St. Vincent’s Hospital, University of Melbourne, Melbourne, VIC 3065, Australia
- Department of Neuroscience, Monash Medical Centre, Melbourne, VIC 3168, Australia
| | - Mark Cook
- Department of Medicine, St. Vincent’s Hospital, University of Melbourne, Melbourne, VIC 3065, Australia
- Graeme Clark Institute, University of Melbourne, Parkville, VIC 3010, Australia
| | - John R. Terry
- EPSRC Centre for Predictive modeling in Healthcare, University of Exeter, Exeter, EX4 4QJ, UK
- Centre for Systems Modelling and Quantitative Biomedicine, University of Birmingham, Birmingham, B15 2TT, UK
| | - Krasimira Tsaneva-Atanasova
- Department of Mathematics, University of Exeter, Exeter, EX4 4QF, UK
- EPSRC Centre for Predictive modeling in Healthcare, University of Exeter, Exeter, EX4 4QJ, UK
- Living System Institute, University of Exeter, Exeter, EX4 4QJ, UK
- Institute for Advanced Study, Technical University of Munich, Lichtenbergstrasse 2a, D-85748 Garching, Germany
- Department of Bioinformatics and Mathematical Modelling, Institute of Biophysics and Biomedical Engineering, Bulgarian Academy of Sciences, 105 Acad. G. Bonchev Str, 1113 Sofia, Bulgaria
| |
Collapse
|
30
|
Balatskaya A, Roehri N, Lagarde S, Pizzo F, Medina S, Wendling F, Bénar CG, Bartolomei F. The “Connectivity Epileptogenicity Index ” (cEI), a method for mapping the different seizure onset patterns in StereoElectroEncephalography recorded seizures. Clin Neurophysiol 2020; 131:1947-1955. [DOI: 10.1016/j.clinph.2020.05.029] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2020] [Revised: 04/15/2020] [Accepted: 05/15/2020] [Indexed: 10/24/2022]
|
31
|
Soriano J, Kubo T, Ikeda K. Feasibility of focal brain cooling for partial epilepsy with secondary generalization: A computational study . ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2020; 2020:2532-2535. [PMID: 33018522 DOI: 10.1109/embc44109.2020.9175157] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Experiments with animal models of epilepsy have consistently shown that focal cooling of epilepsy-induced brain region reversibly suppresses or terminates epileptic discharge activity. Recently, we formulated a physiologically plausible temperature dependence in a neural mass model that can reproduce the effect of focal cooling on epileptic discharge activity. This can be used to implement a temperature control in an implantable cooling device for thermal neuromodulation of the epileptogenic zone in patients with partial epilepsy when seizure activity is detected. However, there have been no experiments that looked into the effect of focal cooling in animal models of epilepsy with secondary generalization in which the seizure activity spreads from the pathologic region to other regions of the brain. Using the temperature-dependent neural mass model and a physiological coupling model, we show that focal cooling stops the propagation of low-frequency discharge activity; on the other hand, it increases the amount of coupling required to propagate high-frequency discharge activity. Moreover, discharge activities that are propagated with cooling are lower in both magnitude and frequency compared to those propagated without cooling. These results suggest the feasibility of focal cooling as an effective alternative therapeutic treatment for medically intractable partial epilepsy even with secondary generalization.Clinical Relevance- The computational study establishes focal cooling of the brain region with partial epilepsy not only suppresses epileptic discharges but can also prevent its generalization to other brain regions.
Collapse
|
32
|
An S, Kang C, Lee HW. Artificial Intelligence and Computational Approaches for Epilepsy. J Epilepsy Res 2020; 10:8-17. [PMID: 32983950 PMCID: PMC7494883 DOI: 10.14581/jer.20003] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2020] [Revised: 06/18/2020] [Accepted: 07/14/2020] [Indexed: 12/30/2022] Open
Abstract
Studies on treatment of epilepsy have been actively conducted in multiple avenues, but there are limitations in improving its efficacy due to between-subject variability in which treatment outcomes vary from patient to patient. Accordingly, there is a growing interest in precision medicine that provides accurate diagnosis for seizure types and optimal treatment for an individual epilepsy patient. Among these approaches, computational studies making this feasible are rapidly progressing in particular and have been widely applied in epilepsy. These computational studies are being conducted in two main streams: 1) artificial intelligence-based studies implementing computational machines with specific functions, such as automatic diagnosis and prognosis prediction for an individual patient, using machine learning techniques based on large amounts of data obtained from multiple patients and 2) patient-specific modeling-based studies implementing biophysical in-silico platforms to understand pathological mechanisms and derive the optimal treatment for each patient by reproducing the brain network dynamics of the particular patient per se based on individual patient's data. These computational approaches are important as it can integrate multiple types of data acquired from patients and analysis results into a single platform. If these kinds of methods are efficiently operated, it would suggest a novel paradigm for precision medicine.
Collapse
Affiliation(s)
- Sora An
- Department of Neurology, Ewha Womans University School of Medicine and Ewha Medical Research Institute, Seoul, Korea.,Department of Medical Science, Ewha Womans University School of Medicine and Ewha Medical Research Institute, Seoul, Korea
| | - Chaewon Kang
- Department of Neurology, Ewha Womans University School of Medicine and Ewha Medical Research Institute, Seoul, Korea.,Department of Computational Medicine, Ewha Womans University School of Medicine and Ewha Medical Research Institute, Seoul, Korea
| | - Hyang Woon Lee
- Department of Neurology, Ewha Womans University School of Medicine and Ewha Medical Research Institute, Seoul, Korea.,Department of Medical Science, Ewha Womans University School of Medicine and Ewha Medical Research Institute, Seoul, Korea.,Department of Computational Medicine, Ewha Womans University School of Medicine and Ewha Medical Research Institute, Seoul, Korea
| |
Collapse
|
33
|
Dynamical Mechanisms of Interictal Resting-State Functional Connectivity in Epilepsy. J Neurosci 2020; 40:5572-5588. [PMID: 32513827 PMCID: PMC7363471 DOI: 10.1523/jneurosci.0905-19.2020] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2019] [Revised: 05/31/2020] [Accepted: 06/02/2020] [Indexed: 12/18/2022] Open
Abstract
Drug-resistant focal epilepsy is a large-scale brain networks disorder characterized by altered spatiotemporal patterns of functional connectivity (FC), even during interictal resting state (RS). Although RS-FC-based metrics can detect these changes, results from RS functional magnetic resonance imaging (RS-fMRI) studies are unclear and difficult to interpret, and the underlying dynamical mechanisms are still largely unknown. To better capture the RS dynamics, we phenomenologically extended the neural mass model of partial seizures, the Epileptor, by including two neuron subpopulations of epileptogenic and nonepileptogenic type, making it capable of producing physiological oscillations in addition to the epileptiform activity. Using the neuroinformatics platform The Virtual Brain, we reconstructed 14 epileptic and 5 healthy human (of either sex) brain network models (BNMs), based on individual anatomical connectivity and clinically defined epileptogenic heatmaps. Through systematic parameter exploration and fitting to neuroimaging data, we demonstrated that epileptic brains during interictal RS are associated with lower global excitability induced by a shift in the working point of the model, indicating that epileptic brains operate closer to a stable equilibrium point than healthy brains. Moreover, we showed that functional networks are unaffected by interictal spikes, corroborating previous experimental findings; additionally, we observed higher excitability in epileptogenic regions, in agreement with the data. We shed light on new dynamical mechanisms responsible for altered RS-FC in epilepsy, involving the following two key factors: (1) a shift of excitability of the whole brain leading to increased stability; and (2) a locally increased excitability in the epileptogenic regions supporting the mixture of hyperconnectivity and hypoconnectivity in these areas. SIGNIFICANCE STATEMENT Advances in functional neuroimaging provide compelling evidence for epilepsy-related brain network alterations, even during the interictal resting state (RS). However, the dynamical mechanisms underlying these changes are still elusive. To identify local and network processes behind the RS-functional connectivity (FC) spatiotemporal patterns, we systematically manipulated the local excitability and the global coupling in the virtual human epileptic patient brain network models (BNMs), complemented by the analysis of the impact of interictal spikes and fitting to the neuroimaging data. Our results suggest that a global shift of the dynamic working point of the brain model, coupled with locally hyperexcitable node dynamics of the epileptogenic networks, provides a mechanistic explanation of the epileptic processes during the interictal RS period. These, in turn, are associated with the changes in FC.
Collapse
|
34
|
Hashemi M, Vattikonda AN, Sip V, Guye M, Bartolomei F, Woodman MM, Jirsa VK. The Bayesian Virtual Epileptic Patient: A probabilistic framework designed to infer the spatial map of epileptogenicity in a personalized large-scale brain model of epilepsy spread. Neuroimage 2020; 217:116839. [PMID: 32387625 DOI: 10.1016/j.neuroimage.2020.116839] [Citation(s) in RCA: 37] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2019] [Revised: 04/02/2020] [Accepted: 04/07/2020] [Indexed: 12/28/2022] Open
Abstract
Despite the importance and frequent use of Bayesian frameworks in brain network modeling for parameter inference and model prediction, the advanced sampling algorithms implemented in probabilistic programming languages to overcome the inference difficulties have received relatively little attention in this context. In this technical note, we propose a probabilistic framework, namely the Bayesian Virtual Epileptic Patient (BVEP), which relies on the fusion of structural data of individuals to infer the spatial map of epileptogenicity in a personalized large-scale brain model of epilepsy spread. To invert the individualized whole-brain model employed in this study, we use the recently developed algorithms known as No-U-Turn Sampler (NUTS) as well as Automatic Differentiation Variational Inference (ADVI). Our results indicate that NUTS and ADVI accurately estimate the degree of epileptogenicity of brain regions, therefore, the hypothetical brain areas responsible for the seizure initiation and propagation, while the convergence diagnostics and posterior behavior analysis validate the reliability of the estimations. Moreover, we illustrate the efficiency of the transformed non-centered parameters in comparison to centered form of parameterization. The Bayesian framework used in this work proposes an appropriate patient-specific strategy for estimating the epileptogenicity of the brain regions to improve outcome after epilepsy surgery.
Collapse
Affiliation(s)
- M Hashemi
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France.
| | - A N Vattikonda
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France
| | - V Sip
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France
| | - M Guye
- Aix Marseille Univ, CNRS, CRMBM, Marseille, France
| | - F Bartolomei
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France; Epileptology Department, and Clinical Neurophysiology Department, Assistance Publique des Hôpitaux de Marseille, Marseille, France
| | - M M Woodman
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France
| | - V K Jirsa
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France.
| |
Collapse
|
35
|
Deeba F, Sanz-Leon P, Robinson PA. Unified dynamics of interictal events and absence seizures. Phys Rev E 2019; 100:022407. [PMID: 31574631 DOI: 10.1103/physreve.100.022407] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2018] [Indexed: 01/09/2023]
Abstract
The dynamics of interictal events between absence seizures and their relationship to seizures themselves are investigated by employing a neural field model of the corticothalamic system. Interictal events are modeled as being due to transient parameter excursions beyond the seizure threshold, in the present case by sufficiently temporally varying the connection strength between the cerebral cortex and the thalamus. Increasing connection strength drives the system into ∼3-Hz seizure oscillations via a supercritical Hopf bifurcation once the linear instability threshold is passed. Depending on the time course of the excursion above threshold, different interictal activity event dynamics are seen in the time series of corticothalamic fields. These resemble experimental interictal time series observed via electroencephalography. It is found that the morphology of these events depends on the magnitude and duration of the excursion above threshold. For a large-amplitude excursion of short duration, events resemble interictal spikes, where one large spike is seen, followed by small damped oscillations. For a short excursion with long duration, events like observed interictal periodic sharp waves are seen. When both amplitude and duration above threshold are large, seizure oscillations are seen. Using these outcomes, proximity to seizure can be estimated and tracked.
Collapse
Affiliation(s)
- F Deeba
- Department of Physics, Dhaka University of Engineering and Technology, Gazipur 1700, Bangladesh; School of Physics, University of Sydney, New South Wales 2006, Australia; and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | - P Sanz-Leon
- School of Physics, University of Sydney, New South Wales 2006, Australia, and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| | - P A Robinson
- School of Physics, University of Sydney, New South Wales 2006, Australia, and Center for Integrative Brain Function, University of Sydney, New South Wales 2006, Australia
| |
Collapse
|
36
|
Jirsa VK, McIntosh AR, Huys R. Grand Unified Theories of the Brain Need Better Understanding of Behavior: The Two-Tiered Emergence of Function. ECOLOGICAL PSYCHOLOGY 2019. [DOI: 10.1080/10407413.2019.1615207] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Viktor K. Jirsa
- Institut de Neurosciences des Systèmes, UMR INSERM 1106, Aix-Marseille Université Faculté de Médecine
| | | | - Raoul Huys
- Université de Toulouse, UMR 5549 CERCO (Centre de Recherche Cerveau et Cognition), UPS, CNRS
| |
Collapse
|
37
|
Lagarde S, Roehri N, Lambert I, Trebuchon A, McGonigal A, Carron R, Scavarda D, Milh M, Pizzo F, Colombet B, Giusiano B, Medina Villalon S, Guye M, Bénar CG, Bartolomei F. Interictal stereotactic-EEG functional connectivity in refractory focal epilepsies. Brain 2019; 141:2966-2980. [PMID: 30107499 DOI: 10.1093/brain/awy214] [Citation(s) in RCA: 109] [Impact Index Per Article: 21.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2017] [Accepted: 06/25/2018] [Indexed: 12/28/2022] Open
Abstract
Drug-refractory focal epilepsies are network diseases associated with functional connectivity alterations both during ictal and interictal periods. A large majority of studies on the interictal/resting state have focused on functional MRI-based functional connectivity. Few studies have used electrophysiology, despite its high temporal capacities. In particular, stereotactic-EEG is highly suitable to study functional connectivity because it permits direct intracranial electrophysiological recordings with relative large-scale sampling. Most previous studies in stereotactic-EEG have been directed towards temporal lobe epilepsy, which does not represent the whole spectrum of drug-refractory epilepsies. The present study aims at filling this gap, investigating interictal functional connectivity alterations behind cortical epileptic organization and its association with post-surgical prognosis. To this purpose, we studied a large cohort of 59 patients with malformation of cortical development explored by stereotactic-EEG with a wide spatial sampling (76 distinct brain areas were recorded, median of 13.2 per patient). We computed functional connectivity using non-linear correlation. We focused on three zones defined by stereotactic-EEG ictal activity: the epileptogenic zone, the propagation zone and the non-involved zone. First, we compared within-zone and between-zones functional connectivity. Second, we analysed the directionality of functional connectivity between these zones. Third, we measured the associations between functional connectivity measures and clinical variables, especially post-surgical prognosis. Our study confirms that functional connectivity differs according to the zone under investigation. We found: (i) a gradual decrease of the within-zone functional connectivity with higher values for epileptogenic zone and propagation zone, and lower for non-involved zones; (ii) preferential coupling between structures of the epileptogenic zone; (iii) preferential coupling between epileptogenic zone and propagation zone; and (iv) poorer post-surgical outcome in patients with higher functional connectivity of non-involved zone (within- non-involved zone, between non-involved zone and propagation zone functional connectivity). Our work suggests that, even during the interictal state, functional connectivity is reinforced within epileptic cortices (epileptogenic zone and propagation zone) with a gradual organization. Moreover, larger functional connectivity alterations, suggesting more diffuse disease, are associated with poorer post-surgical prognosis. This is consistent with computational studies suggesting that connectivity is crucial in order to model the spatiotemporal dynamics of seizures.10.1093/brain/awy214_video1awy214media15833456182001.
Collapse
Affiliation(s)
- Stanislas Lagarde
- APHM, Timone Hospital, Clinical Neurophysiology, Marseille, France.,Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France
| | - Nicolas Roehri
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France
| | - Isabelle Lambert
- APHM, Timone Hospital, Clinical Neurophysiology, Marseille, France.,Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France
| | - Agnès Trebuchon
- APHM, Timone Hospital, Clinical Neurophysiology, Marseille, France.,Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France
| | - Aileen McGonigal
- APHM, Timone Hospital, Clinical Neurophysiology, Marseille, France.,Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France
| | - Romain Carron
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France.,APHM, Timone Hospital, Stereotactic and Functional Neurosurgery, Marseille, France
| | - Didier Scavarda
- APHM, Timone Hospital, Paediatric Neurosurgery, Marseille, France
| | - Mathieu Milh
- APHM, Timone Hospital, Paediatric Neurology, Marseille, France
| | - Francesca Pizzo
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France
| | - Bruno Colombet
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France
| | - Bernard Giusiano
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France
| | - Samuel Medina Villalon
- APHM, Timone Hospital, Clinical Neurophysiology, Marseille, France.,Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France
| | - Maxime Guye
- APHM, Timone Hospital, Clinical Neurophysiology, Marseille, France.,Aix Marseille Univ, CNRS, CRMBM, Marseille, France.,APHM, Timone Hospital, CEMEREM, Marseille, France
| | - Christian-G Bénar
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France
| | - Fabrice Bartolomei
- APHM, Timone Hospital, Clinical Neurophysiology, Marseille, France.,Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France
| |
Collapse
|
38
|
An S, Bartolomei F, Guye M, Jirsa V. Optimization of surgical intervention outside the epileptogenic zone in the Virtual Epileptic Patient (VEP). PLoS Comput Biol 2019; 15:e1007051. [PMID: 31242177 PMCID: PMC6594587 DOI: 10.1371/journal.pcbi.1007051] [Citation(s) in RCA: 42] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2018] [Accepted: 04/17/2019] [Indexed: 01/18/2023] Open
Abstract
Studies to improve the efficacy of epilepsy surgery have focused on better refining the localization of the epileptogenic zone (EZ) with the aim of effectively resecting it. However, in a considerable number of patients, EZs are distributed across multiple brain regions and may involve eloquent areas that cannot be removed due to the risk of neurological complications. There is a clear need for developing alternative approaches to induce seizure relief, but minimal impact on normal brain functions. Here, we develop a personalized in-silico network approach, that suggests effective and safe surgical interventions for each patient. Based on the clinically identified EZ, we employ modularity analysis to identify target brain regions and fiber tracts involved in seizure propagation. We then construct and simulate a patient-specific brain network model comprising phenomenological neural mass models at the nodes, and patient-specific structural brain connectivity using the neuroinformatics platform The Virtual Brain (TVB), in order to evaluate effectiveness and safety of the target zones (TZs). In particular, we assess safety via electrical stimulation for pre- and post-surgical condition to quantify the impact on the signal transmission properties of the network. We demonstrate the existence of a large repertoire of efficient surgical interventions resulting in reduction of degree of seizure spread, but only a small subset of them proves safe. The identification of novel surgical interventions through modularity analysis and brain network simulations may provide exciting solutions to the treatment of inoperable epilepsies. We propose a personalized in-silico surgical approach able to suggest effective and safe surgical options for each epilepsy patient. In particular, we focus on deriving effective alternative methods for those cases where EZs are inoperable because of issues related with neurological complications. Based on modularity analysis using structural brain connectivity from each patient, TZs that would be considered as surgical sites are obtained. The acquired TZs are evaluated by personalized brain network simulations in terms of effectiveness and safety. Through the feedback approach combining modularity analysis and brain network simulations, the optimized TZ options that minimize seizure propagation while not affecting normal brain functions are obtained. Our study has a great importance in that it demonstrates the possibility of computational neuroscience field being able to construct a paradigm for personalized medicine by deriving innovative surgical options suitable for each patient and predicting the surgical outcomes.
Collapse
Affiliation(s)
- Sora An
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France
| | | | - Maxime Guye
- Aix Marseille Univ, CNRS, CRMBM UMR 7339, Marseille, France
| | - Viktor Jirsa
- Aix Marseille Univ, INSERM, INS, Inst Neurosci Syst, Marseille, France
- * E-mail:
| |
Collapse
|
39
|
Grinenko O, Li J, Mosher JC, Wang IZ, Bulacio JC, Gonzalez-Martinez J, Nair D, Najm I, Leahy RM, Chauvel P. A fingerprint of the epileptogenic zone in human epilepsies. Brain 2019; 141:117-131. [PMID: 29253102 PMCID: PMC5837527 DOI: 10.1093/brain/awx306] [Citation(s) in RCA: 118] [Impact Index Per Article: 23.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2017] [Accepted: 09/27/2017] [Indexed: 11/14/2022] Open
Abstract
Defining a bio-electrical marker for the brain area responsible for initiating a seizure remains an unsolved problem. Fast gamma activity has been identified as the most specific marker for seizure onset, but conflicting results have been reported. In this study, we describe an alternative marker, based on an objective description of interictal to ictal transition, with the aim of identifying a time-frequency pattern or ‘fingerprint’ that can differentiate the epileptogenic zone from areas of propagation. Seventeen patients who underwent stereoelectroencephalography were included in the study. Each had seizure onset characterized by sustained gamma activity and were seizure-free after tailored resection or laser ablation. We postulated that the epileptogenic zone was always located inside the resection region based on seizure freedom following surgery. To characterize the ictal frequency pattern, we applied the Morlet wavelet transform to data from each pair of adjacent intracerebral electrode contacts. Based on a visual assessment of the time-frequency plots, we hypothesized that a specific time-frequency pattern in the epileptogenic zone should include a combination of (i) sharp transients or spikes; preceding (ii) multiband fast activity concurrent; with (iii) suppression of lower frequencies. To test this hypothesis, we developed software that automatically extracted each of these features from the time-frequency data. We then used a support vector machine to classify each contact-pair as being within epileptogenic zone or not, based on these features. Our machine learning system identified this pattern in 15 of 17 patients. The total number of identified contacts across all patients was 64, with 58 localized inside the resected area. Subsequent quantitative analysis showed strong correlation between maximum frequency of fast activity and suppression inside the resection but not outside. We did not observe significant discrimination power using only the maximum frequency or the timing of fast activity to differentiate contacts either between resected and non-resected regions or between contacts identified as epileptogenic versus non-epileptogenic. Instead of identifying a single frequency or a single timing trait, we observed the more complex pattern described above that distinguishes the epileptogenic zone. This pattern encompasses interictal to ictal transition and may extend until seizure end. Its time-frequency characteristics can be explained in light of recent models emphasizing the role of fast inhibitory interneurons acting on pyramidal cells as a prominent mechanism in seizure triggering. The pattern clearly differentiates the epileptogenic zone from areas of propagation and, as such, represents an epileptogenic zone ‘fingerprint’.
Collapse
Affiliation(s)
- Olesya Grinenko
- Epilepsy Center, Cleveland Clinic Neurological Institute, Cleveland OH, USA
| | - Jian Li
- Signal and Image Processing Institute, University of Southern California, Los Angeles CA, USA
| | - John C Mosher
- Epilepsy Center, Cleveland Clinic Neurological Institute, Cleveland OH, USA
| | - Irene Z Wang
- Epilepsy Center, Cleveland Clinic Neurological Institute, Cleveland OH, USA
| | - Juan C Bulacio
- Epilepsy Center, Cleveland Clinic Neurological Institute, Cleveland OH, USA
| | | | - Dileep Nair
- Epilepsy Center, Cleveland Clinic Neurological Institute, Cleveland OH, USA
| | - Imad Najm
- Epilepsy Center, Cleveland Clinic Neurological Institute, Cleveland OH, USA
| | - Richard M Leahy
- Signal and Image Processing Institute, University of Southern California, Los Angeles CA, USA
| | - Patrick Chauvel
- Epilepsy Center, Cleveland Clinic Neurological Institute, Cleveland OH, USA
| |
Collapse
|
40
|
Olmi S, Petkoski S, Guye M, Bartolomei F, Jirsa V. Controlling seizure propagation in large-scale brain networks. PLoS Comput Biol 2019; 15:e1006805. [PMID: 30802239 PMCID: PMC6405161 DOI: 10.1371/journal.pcbi.1006805] [Citation(s) in RCA: 66] [Impact Index Per Article: 13.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2018] [Revised: 03/07/2019] [Accepted: 01/18/2019] [Indexed: 01/26/2023] Open
Abstract
Information transmission in the human brain is a fundamentally dynamic network process. In partial epilepsy, this process is perturbed and highly synchronous seizures originate in a local network, the so-called epileptogenic zone (EZ), before recruiting other close or distant brain regions. We studied patient-specific brain network models of 15 drug-resistant epilepsy patients with implanted stereotactic electroencephalography (SEEG) electrodes. Each personalized brain model was derived from structural data of magnetic resonance imaging (MRI) and diffusion tensor weighted imaging (DTI), comprising 88 nodes equipped with region specific neural mass models capable of demonstrating a range of epileptiform discharges. Each patient's virtual brain was further personalized through the integration of the clinically hypothesized EZ. Subsequent simulations and connectivity modulations were performed and uncovered a finite repertoire of seizure propagation patterns. Across patients, we found that (i) patient-specific network connectivity is predictive for the subsequent seizure propagation pattern; (ii) seizure propagation is characterized by a systematic sequence of brain states; (iii) propagation can be controlled by an optimal intervention on the connectivity matrix; (iv) the degree of invasiveness can be significantly reduced via the proposed seizure control as compared to traditional resective surgery. To stop seizures, neurosurgeons typically resect the EZ completely. We showed that stability analysis of the network dynamics, employing structural and dynamical information, estimates reliably the spatiotemporal properties of seizure propagation. This suggests novel less invasive paradigms of surgical interventions to treat and manage partial epilepsy.
Collapse
Affiliation(s)
- Simona Olmi
- Inria Sophia Antipolis Méditerranée Research Centre, MathNeuro Team, 2004 route des Lucioles-Boîte Postale 93 06902 Sophia Antipolis, Cedex, France
- CNR - Consiglio Nazionale delle Ricerche - Istituto dei Sistemi Complessi, 50019, Sesto Fiorentino, Italy
| | - Spase Petkoski
- Aix Marseille Université, Inserm, Institut de Neurosciences des Systèmes, UMR_S 1106, 13005, Marseille, France
| | - Maxime Guye
- Faculté de Médecine de la Timone, centre de Résonance Magnétique et Biologique et Médicale (CRMBM, UMR CNRS-AMU 7339), Medical School of Marseille, Aix-Marseille Université, 13005, Marseille, France
- Assistance Publique - Hôpitaux de Marseille, Hôpital de la Timone, Pôle d’Imagerie, CHU, 13005, Marseille, France
| | - Fabrice Bartolomei
- Assistance Publique - Hôpitaux de Marseille, Hôpital de la Timone, Service de Neurophysiologie Clinique, CHU, 13005 Marseille, France
| | - Viktor Jirsa
- Aix Marseille Université, Inserm, Institut de Neurosciences des Systèmes, UMR_S 1106, 13005, Marseille, France
| |
Collapse
|
41
|
Buchin A, Kerr CC, Huberfeld G, Miles R, Gutkin B. Adaptation and Inhibition Control Pathological Synchronization in a Model of Focal Epileptic Seizure. eNeuro 2018; 5:ENEURO.0019-18.2018. [PMID: 30302390 PMCID: PMC6173584 DOI: 10.1523/eneuro.0019-18.2018] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2018] [Revised: 06/07/2018] [Accepted: 06/07/2018] [Indexed: 01/12/2023] Open
Abstract
Pharmacoresistant epilepsy is a common neurological disorder in which increased neuronal intrinsic excitability and synaptic excitation lead to pathologically synchronous behavior in the brain. In the majority of experimental and theoretical epilepsy models, epilepsy is associated with reduced inhibition in the pathological neural circuits, yet effects of intrinsic excitability are usually not explicitly analyzed. Here we present a novel neural mass model that includes intrinsic excitability in the form of spike-frequency adaptation in the excitatory population. We validated our model using local field potential (LFP) data recorded from human hippocampal/subicular slices. We found that synaptic conductances and slow adaptation in the excitatory population both play essential roles for generating seizures and pre-ictal oscillations. Using bifurcation analysis, we found that transitions towards seizure and back to the resting state take place via Andronov-Hopf bifurcations. These simulations therefore suggest that single neuron adaptation as well as synaptic inhibition are responsible for orchestrating seizure dynamics and transition towards the epileptic state.
Collapse
Affiliation(s)
- Anatoly Buchin
- University of Washington, Department of Physiology and Biophysics (United States, Seattle), 1959 NE Pacific St, 98195
| | - Cliff C. Kerr
- University of Sydney, School of Physics (Australia, Sydney), Physics Rd, NSW 2006
| | - Gilles Huberfeld
- Sorbonne Université-UPMC, Pitié-Salpêtrière Hô, Neurophysiology Department (France, Paris), 47-83 Boulevard de l’Hôpital, 75013
- Institut national de la santé et de la recherche médicale Unit 1129 “Infantile Epilepsies and Brain Plasticity”, Paris Descartes University, Sorbonne Paris Cité University group, (France, Paris), 149 rue de Sévres 75015
| | - Richard Miles
- Brain and Spine Institute, Cortex and Epilepsie Group (France, Paris), 47 Boulevard Hôpital, 75013
| | - Boris Gutkin
- Paris Sciences & Lettres Research University, Laboratoire des Neurosciences Cognitives, Group for Neural Theory (France, Paris), 29, rue d'Ulm, 75005 France
- National Research University Higher School of Economics, Center for Cognition and Decision Making (Russia, Moscow), 20 Myasnitskaya, 109316
| |
Collapse
|
42
|
Rosch RE, Hunter PR, Baldeweg T, Friston KJ, Meyer MP. Calcium imaging and dynamic causal modelling reveal brain-wide changes in effective connectivity and synaptic dynamics during epileptic seizures. PLoS Comput Biol 2018; 14:e1006375. [PMID: 30138336 PMCID: PMC6124808 DOI: 10.1371/journal.pcbi.1006375] [Citation(s) in RCA: 41] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2018] [Revised: 09/05/2018] [Accepted: 07/18/2018] [Indexed: 12/31/2022] Open
Abstract
Pathophysiological explanations of epilepsy typically focus on either the micro/mesoscale (e.g. excitation-inhibition imbalance), or on the macroscale (e.g. network architecture). Linking abnormalities across spatial scales remains difficult, partly because of technical limitations in measuring neuronal signatures concurrently at the scales involved. Here we use light sheet imaging of the larval zebrafish brain during acute epileptic seizure induced with pentylenetetrazole. Spectral changes of spontaneous neuronal activity during the seizure are then modelled using neural mass models, allowing Bayesian inference on changes in effective network connectivity and their underlying synaptic dynamics. This dynamic causal modelling of seizures in the zebrafish brain reveals concurrent changes in synaptic coupling at macro- and mesoscale. Fluctuations of both synaptic connection strength and their temporal dynamics are required to explain observed seizure patterns. These findings highlight distinct changes in local (intrinsic) and long-range (extrinsic) synaptic transmission dynamics as a possible seizure pathomechanism and illustrate how our Bayesian model inversion approach can be used to link existing neural mass models of seizure activity and novel experimental methods.
Collapse
Affiliation(s)
- Richard E. Rosch
- Wellcome Trust Centre for Neuroimaging, University College London, London, United Kingdom
- Developmental Neurosciences Programme, UCL Great Ormond Street Institute of Child Health, University College London, London, United Kingdom
| | - Paul R. Hunter
- Department of Developmental Neurobiology & MRC Center for Neurodevelopmental Disorders, Institute of Psychiatry, Psychology and Neuroscience, King’s College London, London, United Kingdom
| | - Torsten Baldeweg
- Developmental Neurosciences Programme, UCL Great Ormond Street Institute of Child Health, University College London, London, United Kingdom
| | - Karl J. Friston
- Wellcome Trust Centre for Neuroimaging, University College London, London, United Kingdom
| | - Martin P. Meyer
- Department of Developmental Neurobiology & MRC Center for Neurodevelopmental Disorders, Institute of Psychiatry, Psychology and Neuroscience, King’s College London, London, United Kingdom
| |
Collapse
|
43
|
Chizhov AV, Zefirov AV, Amakhin DV, Smirnova EY, Zaitsev AV. Minimal model of interictal and ictal discharges "Epileptor-2". PLoS Comput Biol 2018; 14:e1006186. [PMID: 29851959 PMCID: PMC6005638 DOI: 10.1371/journal.pcbi.1006186] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2017] [Revised: 06/18/2018] [Accepted: 05/09/2018] [Indexed: 12/01/2022] Open
Abstract
Seizures occur in a recurrent manner with intermittent states of interictal and ictal discharges (IIDs and IDs). The transitions to and from IDs are determined by a set of processes, including synaptic interaction and ionic dynamics. Although mathematical models of separate types of epileptic discharges have been developed, modeling the transitions between states remains a challenge. A simple generic mathematical model of seizure dynamics (Epileptor) has recently been proposed by Jirsa et al. (2014); however, it is formulated in terms of abstract variables. In this paper, a minimal population-type model of IIDs and IDs is proposed that is as simple to use as the Epileptor, but the suggested model attributes physical meaning to the variables. The model is expressed in ordinary differential equations for extracellular potassium and intracellular sodium concentrations, membrane potential, and short-term synaptic depression variables. A quadratic integrate-and-fire model driven by the population input current is used to reproduce spike trains in a representative neuron. In simulations, potassium accumulation governs the transition from the silent state to the state of an ID. Each ID is composed of clustered IID-like events. The sodium accumulates during discharge and activates the sodium-potassium pump, which terminates the ID by restoring the potassium gradient and thus polarizing the neuronal membranes. The whole-cell and cell-attached recordings of a 4-AP-based in vitro model of epilepsy confirmed the primary model assumptions and predictions. The mathematical analysis revealed that the IID-like events are large-amplitude stochastic oscillations, which in the case of ID generation are controlled by slow oscillations of ionic concentrations. The IDs originate in the conditions of elevated potassium concentrations in a bath solution via a saddle-node-on-invariant-circle-like bifurcation for a non-smooth dynamical system. By providing a minimal biophysical description of ionic dynamics and network interactions, the model may serve as a hierarchical base from a simple to more complex modeling of seizures. In pathological conditions of epilepsy, the functioning of the neural network crucially depends on the ionic concentrations inside and outside neurons. A number of factors that affect neuronal activity is large. That is why the development of a minimal model that reproduces typical seizures could structure further experimental and analytical studies of the pathological mechanisms. Here, on a base of known biophysical models, we present a simple population-type model that includes only four principal variables, the extracellular potassium concentration, the intracellular sodium concentration, the membrane potential and the synaptic resource diminishing due to short-term synaptic depression. A simple modeled neuron is used as an observer of the population activity. We validate the model assumptions with in vitro experiments. Our model reproduces ictal and interictal events, where the latter result in bursts of spikes in single neurons, and the former represent the cluster of spike bursts. Mathematical analysis reveals that the bursts are spontaneous large-amplitude oscillations, which may cluster after a saddle-node on invariant circle bifurcation in the pro-epileptic conditions. Our consideration has significant bearing in understanding pathological neuronal network dynamics.
Collapse
Affiliation(s)
- Anton V. Chizhov
- Laboratory of Molecular Mechanisms of Neural Interactions, Sechenov Institute of Evolutionary Physiology and Biochemistry of the Russian Academy of Sciences, Saint Petersburg, Russia
- Computational Physics Laboratory, Ioffe Institute, Saint Petersburg, Russia
- * E-mail:
| | - Artyom V. Zefirov
- Laboratory of Molecular Mechanisms of Neural Interactions, Sechenov Institute of Evolutionary Physiology and Biochemistry of the Russian Academy of Sciences, Saint Petersburg, Russia
- Computational Physics Laboratory, Ioffe Institute, Saint Petersburg, Russia
| | - Dmitry V. Amakhin
- Laboratory of Molecular Mechanisms of Neural Interactions, Sechenov Institute of Evolutionary Physiology and Biochemistry of the Russian Academy of Sciences, Saint Petersburg, Russia
| | - Elena Yu. Smirnova
- Laboratory of Molecular Mechanisms of Neural Interactions, Sechenov Institute of Evolutionary Physiology and Biochemistry of the Russian Academy of Sciences, Saint Petersburg, Russia
- Computational Physics Laboratory, Ioffe Institute, Saint Petersburg, Russia
| | - Aleksey V. Zaitsev
- Laboratory of Molecular Mechanisms of Neural Interactions, Sechenov Institute of Evolutionary Physiology and Biochemistry of the Russian Academy of Sciences, Saint Petersburg, Russia
- Institute of Experimental Medicine, Almazov National Medical Research Centre, Saint Petersburg, Russia
| |
Collapse
|
44
|
Deeba F, Sanz-Leon P, Robinson PA. Dependence of absence seizure dynamics on physiological parameter evolution. J Theor Biol 2018; 454:11-21. [PMID: 29807025 DOI: 10.1016/j.jtbi.2018.05.029] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2018] [Revised: 05/22/2018] [Accepted: 05/24/2018] [Indexed: 12/30/2022]
Abstract
A neural field model of the corticothalamic system is applied to investigate the temporal and spectral characteristics of absence seizures in the presence of a temporally varying connection strength between the cerebral cortex and thalamus. Increasing connection strength drives the system into an absence seizure-like state once a threshold is passed and a supercritical Hopf bifurcation occurs. The dynamics and spectral characteristics of the resulting model seizures are explored as functions of maximum connection strength, time above threshold, and the rate at which the connection strength increases (ramp rate). Our results enable spectral and temporal characteristics of seizures to be related to changes in the underlying physiological evolution of connections via nonlinear dynamics and neural field theory. Spectral analysis reveals that the power of the harmonics and the duration of the oscillations increase as the maximum connection strength and the time above threshold increase. It is also found that the time to reach the stable limit-cycle seizure oscillation from the instability threshold decreases with the square root of the ramp rate.
Collapse
Affiliation(s)
- F Deeba
- School of Physics, University of Sydney, NSW 2006, Australia; Center for Integrative Brain Function, University of Sydney, NSW 2006, Australia.
| | - Paula Sanz-Leon
- School of Physics, University of Sydney, NSW 2006, Australia; Center for Integrative Brain Function, University of Sydney, NSW 2006, Australia
| | - P A Robinson
- School of Physics, University of Sydney, NSW 2006, Australia; Center for Integrative Brain Function, University of Sydney, NSW 2006, Australia
| |
Collapse
|
45
|
Vaugier L, Lagarde S, McGonigal A, Trébuchon A, Milh M, Lépine A, Scavarda D, Carron R, Bartolomei F. The role of stereoelectroencephalography (SEEG) in reevaluation of epilepsy surgery failures. Epilepsy Behav 2018. [PMID: 29526579 DOI: 10.1016/j.yebeh.2018.02.005] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
Abstract
Management of patients after initial epilepsy surgical failure is challenging. In this study, we report our experience in using the stereoelectroencephalography (SEEG) method in the reevaluation of patients after initial epilepsy surgical failure. We selected 28 patients examined through SEEG in our department for drug-resistant focal epilepsy following initial epilepsy surgical failure. For each patient, the residual seizure onset zone (rSOZ) as defined by SEEG was classified as either contiguous if the seizure onset zone (SOZ) was focal and close to the surgical cavity (same lobe) or noncontiguous in cases where the SOZ included site(s) distant from the surgical cavity. The rSOZ was defined according to visual analysis of SEEG traces completed by an estimation of the epileptogenicity index (EI). A second surgical procedure was performed in 12 patients (45%). A favorable outcome (Engel class I or II) was obtained in 9/12 patients (6 in Engel class I, 50%). The proportion of patients that had reoperation was higher in the contiguous group (80%) than in the noncontiguous group (22%) (p=0.02). A rSOZ localized in close relation to the initial surgical resection zone (contiguous group) was found in 10 patients (35%). Among them, 8 have since undergone reoperation, and a good outcome (Engel class I) was achieved in 5/8 (63%). A rSOZ involving a distant region from the first surgery was observed in 18 patients (65%) (noncontiguous group). Among them, only 4 have undergone reoperation, leading to a failure in 2 (Engel class III or IV) and a good outcome in 2 (IA). Ten patients had a first standard temporal lobectomy, and in 50% of these cases, the insula was involved in the rSOZ. Stereoelectroencephalography offers a unique way to evaluate the rSOZ at the individual level and thus guide further surgical decision-making. The best results are observed in patients having a focal rSOZ close to the site of the surgical resection in the first surgery.
Collapse
Affiliation(s)
- Lisa Vaugier
- APHM, Timone Hospital, Clinical Neurophysiology and Epileptology Department, Marseille, France
| | - Stanislas Lagarde
- APHM, Timone Hospital, Clinical Neurophysiology and Epileptology Department, Marseille, France; Aix Marseille Univ, Inserm, INS, Institut de Neurosciences des Systèmes, Marseille, France
| | - Aileen McGonigal
- APHM, Timone Hospital, Clinical Neurophysiology and Epileptology Department, Marseille, France; Aix Marseille Univ, Inserm, INS, Institut de Neurosciences des Systèmes, Marseille, France
| | - Agnès Trébuchon
- APHM, Timone Hospital, Clinical Neurophysiology and Epileptology Department, Marseille, France; Aix Marseille Univ, Inserm, INS, Institut de Neurosciences des Systèmes, Marseille, France
| | - Mathieu Milh
- APHM, Timone Hospital, Clinical Neurophysiology and Epileptology Department, Marseille, France; Aix Marseille Univ, Inserm, INS, Institut de Neurosciences des Systèmes, Marseille, France
| | - Anne Lépine
- APHM, Timone Hospital, Clinical Neurophysiology and Epileptology Department, Marseille, France; Aix Marseille Univ, Inserm, INS, Institut de Neurosciences des Systèmes, Marseille, France
| | - Didier Scavarda
- APHM, Timone Hospital, Clinical Neurophysiology and Epileptology Department, Marseille, France; Aix Marseille Univ, Inserm, INS, Institut de Neurosciences des Systèmes, Marseille, France; APHM, Timone Hospital, Paediatric Neurosurgery Department, Marseille, France
| | - Romain Carron
- APHM, Timone Hospital, Clinical Neurophysiology and Epileptology Department, Marseille, France; Aix Marseille Univ, Inserm, INS, Institut de Neurosciences des Systèmes, Marseille, France; APHM, Timone Hospital, Functional and Stereotactical Neurosurgery Department, Marseille, France
| | - Fabrice Bartolomei
- APHM, Timone Hospital, Clinical Neurophysiology and Epileptology Department, Marseille, France; Aix Marseille Univ, Inserm, INS, Institut de Neurosciences des Systèmes, Marseille, France.
| |
Collapse
|
46
|
Proix T, Jirsa VK, Bartolomei F, Guye M, Truccolo W. Predicting the spatiotemporal diversity of seizure propagation and termination in human focal epilepsy. Nat Commun 2018. [PMID: 29540685 PMCID: PMC5852068 DOI: 10.1038/s41467-018-02973-y] [Citation(s) in RCA: 69] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Abstract
Recent studies have shown that seizures can spread and terminate across brain areas via a rich diversity of spatiotemporal patterns. In particular, while the location of the seizure onset area is usually invariant across seizures in an individual patient, the source of traveling (2–3 Hz) spike-and-wave discharges during seizures can either move with the slower propagating ictal wavefront or remain stationary at the seizure onset area. Furthermore, although many focal seizures terminate synchronously across brain areas, some evolve into distinct ictal clusters and terminate asynchronously. Here, we introduce a unifying perspective based on a new neural field model of epileptic seizure dynamics. Two main mechanisms, the co-existence of wave propagation in excitable media and coupled-oscillator dynamics, together with the interaction of multiple time scales, account for the reported diversity. We confirm our predictions in seizures and tractography data obtained from patients with pharmacologically resistant epilepsy. Our results contribute toward patient-specific seizure modeling. A major goal of epilepsy research is understanding the spatiotemporal dynamics of seizure. Here, the authors extend the Epileptor neural mass model into a neural field model, in order to provide a unified and patient-specific model of seizure initiation, propagation, and termination.
Collapse
Affiliation(s)
- Timothée Proix
- Department of Neuroscience, Brown University, Providence, RI, 02912, USA.,Institute for Brain Science, Brown University, Providence, RI, 02912, USA.,Center for Neurorestoration & Neurotechnology, U.S. Department of Veterans Affairs, Providence, RI, 02912, USA
| | - Viktor K Jirsa
- Institut de Neurosciences des Systèmes (INS), Inserm, Aix Marseille Univ, Marseille, 13005, France
| | - Fabrice Bartolomei
- Institut de Neurosciences des Systèmes (INS), Inserm, Aix Marseille Univ, Marseille, 13005, France
| | - Maxime Guye
- CNRS, CRMBM UMR 7339, Aix Marseille Univ, Marseille, 13005, France
| | - Wilson Truccolo
- Department of Neuroscience, Brown University, Providence, RI, 02912, USA. .,Institute for Brain Science, Brown University, Providence, RI, 02912, USA. .,Center for Neurorestoration & Neurotechnology, U.S. Department of Veterans Affairs, Providence, RI, 02912, USA.
| |
Collapse
|
47
|
Rosch R, Baldeweg T, Moeller F, Baier G. Network dynamics in the healthy and epileptic developing brain. Netw Neurosci 2018; 2:41-59. [PMID: 29911676 PMCID: PMC5989999 DOI: 10.1162/netn_a_00026] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2017] [Accepted: 09/09/2017] [Indexed: 12/29/2022] Open
Abstract
Electroencephalography (EEG) allows recording of cortical activity at high temporal resolution. EEG recordings can be summarized along different dimensions using network-level quantitative measures, such as channel-to-channel correlation, or band power distributions across channels. These reveal network patterns that unfold over a range of different timescales and can be tracked dynamically. Here we describe the dynamics of network state transitions in EEG recordings of spontaneous brain activity in normally developing infants and infants with severe early infantile epileptic encephalopathies (n = 8, age: 1–8 months). We describe differences in measures of EEG dynamics derived from band power, and correlation-based summaries of network-wide brain activity. We further show that EEGs from different patient groups and controls may be distinguishable on a small set of the novel quantitative measures introduced here, which describe dynamic network state switching. Quantitative measures related to the sharpness of switching from one correlation pattern to another show the largest differences between groups. These findings reveal that the early epileptic encephalopathies are associated with characteristic dynamic features at the network level. Quantitative network-based analyses like the one presented here may in the future inform the clinical use of quantitative EEG for diagnosis.
Collapse
Affiliation(s)
- Richard Rosch
- Wellcome Trust Centre for Neuroimaging, University College London, United Kingdom.,Developmental Neurosciences Programme, UCL Great Ormond Street Institute of Child Health, University College London, United Kingdom
| | - Torsten Baldeweg
- Developmental Neurosciences Programme, UCL Great Ormond Street Institute of Child Health, University College London, United Kingdom
| | - Friederike Moeller
- Department of Clinical Neurophysiology, Great Ormond Street Hospital, London, United Kingdom
| | - Gerold Baier
- Cell and Developmental Biology, University College London, United Kingdom
| |
Collapse
|
48
|
Gu S, Cieslak M, Baird B, Muldoon SF, Grafton ST, Pasqualetti F, Bassett DS. The Energy Landscape of Neurophysiological Activity Implicit in Brain Network Structure. Sci Rep 2018; 8:2507. [PMID: 29410486 PMCID: PMC5802783 DOI: 10.1038/s41598-018-20123-8] [Citation(s) in RCA: 54] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2017] [Accepted: 01/08/2018] [Indexed: 01/03/2023] Open
Abstract
A critical mystery in neuroscience lies in determining how anatomical structure impacts the complex functional dynamics of the brain. How does large-scale brain circuitry constrain states of neuronal activity and transitions between those states? We address these questions using a maximum entropy model of brain dynamics informed by white matter tractography. We demonstrate that the most probable brain states - characterized by minimal energy - display common activation profiles across brain areas: local spatially-contiguous sets of brain regions reminiscent of cognitive systems are co-activated frequently. The predicted activation rate of these systems is highly correlated with the observed activation rate measured in a separate resting state fMRI data set, validating the utility of the maximum entropy model in describing neurophysiological dynamics. This approach also offers a formal notion of the energy of activity within a system, and the energy of activity shared between systems. We observe that within- and between-system energies cleanly separate cognitive systems into distinct categories, optimized for differential contributions to integrated versus segregated function. These results support the notion that energetic and structural constraints circumscribe brain dynamics, offering insights into the roles that cognitive systems play in driving whole-brain activation patterns.
Collapse
Affiliation(s)
- Shi Gu
- Department of Computer Science and Engineering, University of Electronic Science and Technology of China, Chengdu, Sichuan, 611731, China
- Department of Psychiatry, University of Pennsylvania, Philadelphia, PA, 19104, USA
- Department of Bioengineering, University of Pennsylvania, Philadelphia, PA, 19104, USA
| | - Matthew Cieslak
- Department of Psychological and Brain Sciences, University of California, Santa Barbara, CA, 93106, USA
| | - Benjamin Baird
- Center for Sleep and Consciousness, University of Wisconsin - Madison, Madison, WI, 53706, USA
| | - Sarah F Muldoon
- Department of Mathematics and CDSE Program, University at Buffalo, SUNY, Buffalo, NY, 14260, USA
| | - Scott T Grafton
- Department of Psychological and Brain Sciences, University of California, Santa Barbara, CA, 93106, USA
| | - Fabio Pasqualetti
- Department of Mechanical Engineering, University of California, Riverside, CA, 92521, USA
| | - Danielle S Bassett
- Department of Bioengineering, University of Pennsylvania, Philadelphia, PA, 19104, USA.
- Department of Electrical & Systems Engineering, University of Pennsylvania, Philadelphia, PA, 19104, USA.
- Department of Neurology, University of Pennsylvania, Philadelphia, PA, 19104, USA.
| |
Collapse
|
49
|
Klein P, Dingledine R, Aronica E, Bernard C, Blümcke I, Boison D, Brodie MJ, Brooks-Kayal AR, Engel J, Forcelli PA, Hirsch LJ, Kaminski RM, Klitgaard H, Kobow K, Lowenstein DH, Pearl PL, Pitkänen A, Puhakka N, Rogawski MA, Schmidt D, Sillanpää M, Sloviter RS, Steinhäuser C, Vezzani A, Walker MC, Löscher W. Commonalities in epileptogenic processes from different acute brain insults: Do they translate? Epilepsia 2018; 59:37-66. [PMID: 29247482 PMCID: PMC5993212 DOI: 10.1111/epi.13965] [Citation(s) in RCA: 189] [Impact Index Per Article: 31.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 11/01/2017] [Indexed: 12/12/2022]
Abstract
The most common forms of acquired epilepsies arise following acute brain insults such as traumatic brain injury, stroke, or central nervous system infections. Treatment is effective for only 60%-70% of patients and remains symptomatic despite decades of effort to develop epilepsy prevention therapies. Recent preclinical efforts are focused on likely primary drivers of epileptogenesis, namely inflammation, neuron loss, plasticity, and circuit reorganization. This review suggests a path to identify neuronal and molecular targets for clinical testing of specific hypotheses about epileptogenesis and its prevention or modification. Acquired human epilepsies with different etiologies share some features with animal models. We identify these commonalities and discuss their relevance to the development of successful epilepsy prevention or disease modification strategies. Risk factors for developing epilepsy that appear common to multiple acute injury etiologies include intracranial bleeding, disruption of the blood-brain barrier, more severe injury, and early seizures within 1 week of injury. In diverse human epilepsies and animal models, seizures appear to propagate within a limbic or thalamocortical/corticocortical network. Common histopathologic features of epilepsy of diverse and mostly focal origin are microglial activation and astrogliosis, heterotopic neurons in the white matter, loss of neurons, and the presence of inflammatory cellular infiltrates. Astrocytes exhibit smaller K+ conductances and lose gap junction coupling in many animal models as well as in sclerotic hippocampi from temporal lobe epilepsy patients. There is increasing evidence that epilepsy can be prevented or aborted in preclinical animal models of acquired epilepsy by interfering with processes that appear common to multiple acute injury etiologies, for example, in post-status epilepticus models of focal epilepsy by transient treatment with a trkB/PLCγ1 inhibitor, isoflurane, or HMGB1 antibodies and by topical administration of adenosine, in the cortical fluid percussion injury model by focal cooling, and in the albumin posttraumatic epilepsy model by losartan. Preclinical studies further highlight the roles of mTOR1 pathways, JAK-STAT3, IL-1R/TLR4 signaling, and other inflammatory pathways in the genesis or modulation of epilepsy after brain injury. The wealth of commonalities, diversity of molecular targets identified preclinically, and likely multidimensional nature of epileptogenesis argue for a combinatorial strategy in prevention therapy. Going forward, the identification of impending epilepsy biomarkers to allow better patient selection, together with better alignment with multisite preclinical trials in animal models, should guide the clinical testing of new hypotheses for epileptogenesis and its prevention.
Collapse
Affiliation(s)
- Pavel Klein
- Mid-Atlantic Epilepsy and Sleep Center, Bethesda, MD, USA
| | | | - Eleonora Aronica
- Department of (Neuro) Pathology, Academic Medical Center and Swammerdam Institute for Life Sciences, Center for Neuroscience, University of Amsterdam, Amsterdam, The Netherlands
- Stichting Epilepsie Instellingen Nederland (SEIN), Heemstede, The Netherlands
| | - Christophe Bernard
- Aix Marseille Univ, Inserm, INS, Instit Neurosci Syst, Marseille, 13005, France
| | - Ingmar Blümcke
- Department of Neuropathology, University Hospital Erlangen, Erlangen, Germany
| | - Detlev Boison
- Robert Stone Dow Neurobiology Laboratories, Legacy Research Institute, Portland, OR, USA
| | - Martin J Brodie
- Epilepsy Unit, West Glasgow Ambulatory Care Hospital-Yorkhill, Glasgow, UK
| | - Amy R Brooks-Kayal
- Division of Neurology, Departments of Pediatrics and Neurology, University of Colorado School of Medicine, Aurora, CO, USA
- Children's Hospital Colorado, Aurora, CO, USA
- Neuroscience Program, University of Colorado Anschutz Medical Campus, Aurora, CO, USA
| | - Jerome Engel
- Departments of Neurology, Neurobiology, and Psychiatry and Biobehavioral Sciences, David Geffen School of Medicine, Brain Research Institute, University of California, Los Angeles, CA, USA
| | | | | | | | | | - Katja Kobow
- Department of Neuropathology, University Hospital Erlangen, Erlangen, Germany
| | | | - Phillip L Pearl
- Department of Neurology, Boston Children's Hospital, Harvard Medical School, Boston, MA, USA
| | - Asla Pitkänen
- Department of Neurobiology, A. I. Virtanen Institute for Molecular Sciences, University of Eastern Finland, Kuopio, Finland
| | - Noora Puhakka
- Department of Neurobiology, A. I. Virtanen Institute for Molecular Sciences, University of Eastern Finland, Kuopio, Finland
| | - Michael A Rogawski
- Department of Neurology, University of California, Davis, Sacramento, CA, USA
| | | | - Matti Sillanpää
- Departments of Child Neurology and General Practice, University of Turku and Turku University Hospital, Turku, Finland
| | - Robert S Sloviter
- Department of Neurobiology, Morehouse School of Medicine, Atlanta, GA, USA
| | - Christian Steinhäuser
- Institute of Cellular Neurosciences, Medical Faculty, University of Bonn, Bonn, Germany
| | - Annamaria Vezzani
- Department of Neuroscience, IRCCS-Istituto di Ricerche Farmacologiche Mario Negri, Institute for Pharmacological Research, Milan,, Italy
| | - Matthew C Walker
- Department of Clinical and Experimental Epilepsy, UCL Institute of Neurology, London, UK
| | - Wolfgang Löscher
- Department of Pharmacology, Toxicology, and Pharmacy, University of Veterinary Medicine, Hannover, Germany
- Center for Systems Neuroscience, Hannover, Germany
| |
Collapse
|
50
|
The Virtual Mouse Brain: A Computational Neuroinformatics Platform to Study Whole Mouse Brain Dynamics. eNeuro 2017; 4:eN-MNT-0111-17. [PMID: 28664183 PMCID: PMC5489253 DOI: 10.1523/eneuro.0111-17.2017] [Citation(s) in RCA: 38] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2017] [Revised: 05/29/2017] [Accepted: 05/31/2017] [Indexed: 12/17/2022] Open
Abstract
Connectome-based modeling of large-scale brain network dynamics enables causal in silico interrogation of the brain’s structure-function relationship, necessitating the close integration of diverse neuroinformatics fields. Here we extend the open-source simulation software The Virtual Brain (TVB) to whole mouse brain network modeling based on individual diffusion magnetic resonance imaging (dMRI)-based or tracer-based detailed mouse connectomes. We provide practical examples on how to use The Virtual Mouse Brain (TVMB) to simulate brain activity, such as seizure propagation and the switching behavior of the resting state dynamics in health and disease. TVMB enables theoretically driven experimental planning and ways to test predictions in the numerous strains of mice available to study brain function in normal and pathological conditions.
Collapse
|