Individual brain structure and modelling predict seizure propagation

Artificial Intelligence and Computational Approaches for Epilepsy

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.

Dynamical Mechanisms of Interictal Resting-State Functional Connectivity in Epilepsy

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.

Control of brain network dynamics across diverse scales of space and time

The human brain is composed of distinct regions that are each associated with particular functions and distinct propensities for the control of neural dynamics. However, the relation between these functions and control profiles is poorly understood, as is the variation in this relation across diverse scales of space and time. Here we probe the relation between control and dynamics in brain networks constructed from diffusion tensor imaging data in a large community sample of young adults. Specifically, we probe the control properties of each brain region and investigate their relationship with dynamics across various spatial scales using the Laplacian eigenspectrum. In addition, through analysis of regional modal controllability and partitioning of modes, we determine whether the associated dynamics are fast or slow, as well as whether they are alternating or monotone. We find that brain regions that facilitate the control of energetically easy transitions are associated with activity on short length scales and slow timescales. Conversely, brain regions that facilitate control of difficult transitions are associated with activity on long length scales and fast timescales. Built on linear dynamical models, our results offer parsimonious explanations for the activity propagation and network control profiles supported by regions of differing neuroanatomical structure.

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

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.

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

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

Dynamic network properties of the interictal brain determine whether seizures appear focal or generalised

Current explanatory concepts suggest seizures emerge from ongoing dynamics of brain networks. It is unclear how brain network properties determine focal or generalised seizure onset, or how network properties can be described in a clinically-useful manner. Understanding network properties would cast light on seizure-generating mechanisms and allow to quantify to which extent a seizure is focal or generalised. Functional brain networks were estimated in segments of scalp-EEG without interictal discharges (68 people with epilepsy, 38 controls). Simplified brain dynamics were simulated using a computer model. We introduce: Critical Coupling (C_c), the ability of a network to generate seizures; Onset Index (OI), the tendency of a region to generate seizures; and Participation Index (PI), the tendency of a region to become involved in seizures. C_c was lower in both patient groups compared with controls. OI and PI were more variable in focal-onset than generalised-onset cases. In focal cases, the regions with highest OI and PI corresponded to the side of seizure onset. Properties of interictal functional networks from scalp EEG can be estimated using a computer model and used to predict seizure likelihood and onset patterns. This may offer potential to enhance diagnosis through quantification of seizure type using inter-ictal recordings.

Controversies on the network theory of epilepsy: Debates held during the ICTALS 2019 conference

Debates on six controversial topics on the network theory of epilepsy were held during two debate sessions, as part of the International Conference for Technology and Analysis of Seizures, 2019 (ICTALS 2019) convened at the University of Exeter, UK, September 2-5 2019. The debate topics were (1) From pathologic to physiologic: is the epileptic network part of an existing large-scale brain network? (2) Are micro scale recordings pertinent for defining the epileptic network? (3) From seconds to years: do we need all temporal scales to define an epileptic network? (4) Is it necessary to fully define the epileptic network to control it? (5) Is controlling seizures sufficient to control the epileptic network? (6) Does the epileptic network want to be controlled? This article, written by the organizing committee for the debate sessions and the debaters, summarizes the arguments presented during the debates on these six topics.

Linking Structure and Function in Macroscale Brain Networks

Structure-function relationships are a fundamental principle of many naturally occurring systems. However, network neuroscience research suggests that there is an imperfect link between structural connectivity and functional connectivity in the brain. Here, we synthesize the current state of knowledge linking structure and function in macroscale brain networks and discuss the different types of models used to assess this relationship. We argue that current models do not include the requisite biological detail to completely predict function. Structural network reconstructions enriched with local molecular and cellular metadata, in concert with more nuanced representations of functions and properties, hold great potential for a truly multiscale understanding of the structure-function relationship.

Changing concepts in presurgical assessment for epilepsy surgery

Candidates for epilepsy surgery must undergo presurgical evaluation to establish whether and how surgical treatment can stop seizures without causing neurological deficits. Various techniques, including MRI, PET, single-photon emission CT, video-EEG, magnetoencephalography and invasive EEG, aim to identify the diseased brain tissue and the involved network. Recent technical and methodological developments, encompassing both advances in existing techniques and new combinations of technologies, are enhancing the ability to define the optimal resection strategy. Multimodal interpretation and predictive computer models are expected to aid surgical planning and patient counselling, and multimodal intraoperative guidance is likely to increase surgical precision. In this Review, we discuss how the knowledge derived from these new approaches is challenging our way of thinking about surgery to stop focal seizures. In particular, we highlight the importance of looking beyond the EEG seizure onset zone and considering focal epilepsy as a brain network disease in which long-range connections need to be taken into account. We also explore how new diagnostic techniques are revealing essential information in the brain that was previously hidden from view.

Interictal epileptiform discharges shape large-scale intercortical communication

Dynamic interactions between remote but functionally specialized brain regions enable complex information processing. This intercortical communication is disrupted in the neural networks of patients with focal epilepsy, and epileptic activity can exert widespread effects within the brain. Using large-scale human intracranial electroencephalography recordings, we show that interictal epileptiform discharges (IEDs) are significantly coupled with spindles in discrete, individualized brain regions outside of the epileptic network. We found that a substantial proportion of these localized spindles travel across the cortical surface. Brain regions that participate in this IED-driven oscillatory coupling express spindles that have a broader spatial extent and higher tendency to propagate than spindles occurring in uncoupled regions. These altered spatiotemporal oscillatory properties identify areas that are shaped by epileptic activity independent of IED or seizure detection. Our findings suggest that IED-spindle coupling may be an important mechanism of interictal global network dysfunction that could be targeted to prevent disruption of normal neural activity.

Quantification and Selection of Ictogenic Zones in Epilepsy Surgery

Network models of brain dynamics provide valuable insight into the healthy functioning of the brain and how this breaks down in disease. A pertinent example is the use of network models to understand seizure generation (ictogenesis) in epilepsy. Recently, computational models have emerged to aid our understanding of seizures and to predict the outcome of surgical perturbations to brain networks. Such approaches provide the opportunity to quantify the effect of removing regions of tissue from brain networks and thereby search for the optimal resection strategy. Here, we use computational models to elucidate how sets of nodes contribute to the ictogenicity of networks. In small networks we fully elucidate the ictogenicity of all possible sets of nodes and demonstrate that the distribution of ictogenicity across sets depends on network topology. However, the full elucidation is a combinatorial problem that becomes intractable for large networks. Therefore, we combine computational models with a genetic algorithm to search for minimal sets of nodes that contribute significantly to ictogenesis. We demonstrate the potential applicability of these methods in practice by identifying optimal sets of nodes to resect in networks derived from 20 individuals who underwent resective surgery for epilepsy. We show that they have the potential to aid epilepsy surgery by suggesting alternative resection sites as well as facilitating the avoidance of brain regions that should not be resected.

Transmission time delays organize the brain network synchronization

The timing of activity across brain regions can be described by its phases for oscillatory processes, and is of crucial importance for brain functioning. The structure of the brain constrains its dynamics through the delays due to propagation and the strengths of the white matter tracts. We use self-sustained delay-coupled, non-isochronous, nonlinearly damped and chaotic oscillators to study how spatio-temporal organization of the brain governs phase lags between the coherent activity of its regions. In silico results for the brain network model demonstrate a robust switching from in- to anti-phase synchronization by increasing the frequency, with a consistent lagging of the stronger connected regions. Relative phases are well predicted by an earlier analysis of Kuramoto oscillators, confirming the spatial heterogeneity of time delays as a crucial mechanism in shaping the functional brain architecture. Increased frequency and coupling are also shown to distort the oscillators by decreasing their amplitude, and stronger regions have lower, but more synchronized activity. These results indicate specific features in the phase relationships within the brain that need to hold for a wide range of local oscillatory dynamics, given that the time delays of the connectome are proportional to the lengths of the structural pathways. This article is part of the theme issue 'Nonlinear dynamics of delay systems'.

The road ahead in clinical network neuroscience

Clinical network neuroscience, the study of brain network topology in neurological and psychiatric diseases, has become a mainstay field within clinical neuroscience. Being a multidisciplinary group of clinical network neuroscience experts based in The Netherlands, we often discuss the current state of the art and possible avenues for future investigations. These discussions revolve around questions like "How do dynamic processes alter the underlying structural network?" and "Can we use network neuroscience for disease classification?" This opinion paper is an incomplete overview of these discussions and expands on ten questions that may potentially advance the field. By no means intended as a review of the current state of the field, it is instead meant as a conversation starter and source of inspiration to others.

Linking Molecular Pathways and Large-Scale Computational Modeling to Assess Candidate Disease Mechanisms and Pharmacodynamics in Alzheimer's Disease

Introduction: While the prevalence of neurodegenerative diseases associated with dementia such as Alzheimer's disease (AD) increases, our knowledge on the underlying mechanisms, outcome predictors, or therapeutic targets is limited. In this work, we demonstrate how computational multi-scale brain modeling links phenomena of different scales and therefore identifies potential disease mechanisms leading the way to improved diagnostics and treatment. Methods: The Virtual Brain (TVB; thevirtualbrain.org) neuroinformatics platform allows standardized large-scale structural connectivity-based simulations of whole brain dynamics. We provide proof of concept for a novel approach that quantitatively links the effects of altered molecular pathways onto neuronal population dynamics. As a novelty, we connect chemical compounds measured with positron emission tomography (PET) with neural function in TVB addressing the phenomenon of hyperexcitability in AD related to the protein amyloid beta (Abeta). We construct personalized virtual brains based on an averaged healthy connectome and individual PET derived distributions of Abeta in patients with mild cognitive impairment (MCI, N = 8) and Alzheimer's Disease (AD, N = 10) and in age-matched healthy controls (HC, N = 15) using data from ADNI-3 data base (http://adni.loni.usc.edu). In the personalized virtual brains, individual Abeta burden modulates regional Excitation-Inhibition balance, leading to local hyperexcitation with high Abeta loads. We analyze simulated regional neural activity and electroencephalograms (EEG). Results: Known empirical alterations of EEG in patients with AD compared to HCs were reproduced by simulations. The virtual AD group showed slower frequencies in simulated local field potentials and EEG compared to MCI and HC groups. The heterogeneity of the Abeta load is crucial for the virtual EEG slowing which is absent for control models with homogeneous Abeta distributions. Slowing phenomena primarily affect the network hubs, independent of the spatial distribution of Abeta. Modeling the N-methyl-D-aspartate (NMDA) receptor antagonism of memantine in local population models, reveals potential functional reversibility of the observed large-scale alterations (reflected by EEG slowing) in virtual AD brains. Discussion: We demonstrate how TVB enables the simulation of systems effects caused by pathogenetic molecular candidate mechanisms in human virtual brains.

Background EEG Connectivity Captures the Time-Course of Epileptogenesis in a Mouse Model of Epilepsy

Large-scale brain networks are increasingly recognized as important for the generation of seizures in epilepsy. However, how a network evolves from a healthy state through the process of epileptogenesis remains unclear. To address this question, here, we study longitudinal epicranial background EEG recordings (30 electrodes, EEG free from epileptiform activity) of a mouse model of mesial temporal lobe epilepsy. We analyze functional connectivity networks and observe that over the time course of epileptogenesis the networks become increasingly asymmetric. Furthermore, computational modelling reveals that a set of nodes, located outside of the region of initial insult, emerges as particularly important for the network dynamics. These findings are consistent with experimental observations, thus demonstrating that ictogenic mechanisms can be revealed on the EEG, that computational models can be used to monitor unfolding epileptogenesis and that both the primary focus and epileptic network play a role in epileptogenesis.

The path from scientific discovery to cures for epilepsy

The epilepsies are a complex group of disorders that can be caused by a myriad of genetic and acquired factors. As such, identifying interventions that will prevent development of epilepsy, as well as cure the disorder once established, will require a multifaceted approach. Here we discuss the progress in scientific discovery propelling us towards this goal, including identification of genetic risk factors and big data approaches that integrate clinical and molecular 'omics' datasets to identify common pathophysiological signatures and biomarkers. We discuss the many animal and cellular models of epilepsy, what they have taught us about pathophysiology, and the cutting edge cellular, optogenetic, chemogenetic and anti-seizure drug screening approaches that are being used to find new cures in these models. Finally, we reflect on the work that still needs to be done towards identify at-risk individuals early, targeting and stopping epileptogenesis, and optimizing promising treatment approaches. Ultimately, developing and implementing cures for epilepsy will require a coordinated and immense effort from clinicians and basic scientists, as well as industry, and should always be guided by the needs of individuals affected by epilepsy and their families. This article is part of the special issue entitled 'New Epilepsy Therapies for the 21st Century - From Antiseizure Drugs to Prevention, Modification and Cure of Epilepsy'.

Interictal stereotactic-EEG functional connectivity in refractory focal epilepsies

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.

The Human Brain Project-Synergy between neuroscience, computing, informatics, and brain-inspired technologies

The Human Brain Project (HBP) is a European flagship project with a 10-year horizon aiming to understand the human brain and to translate neuroscience knowledge into medicine and technology. To achieve such aims, the HBP explores the multilevel complexity of the brain in space and time; transfers the acquired knowledge to brain-derived applications in health, computing, and technology; and provides shared and open computing tools and data through the HBP European brain research infrastructure. We discuss how the HBP creates a transdisciplinary community of researchers united by the quest to understand the brain, with fascinating perspectives on societal benefits.

This Community Page article presents the Human Brain Project; a European Flagship project with a ten-year horizon aiming to understand the human brain and translate neuroscience knowledge into medicine and technology.

Characterizing the role of the structural connectome in seizure dynamics

How does the human brain's structural scaffold give rise to its intricate functional dynamics? This is a central question in translational neuroscience that is particularly relevant to epilepsy, a disorder affecting over 50 million subjects worldwide. Treatment for medication-resistant focal epilepsy is often structural-through surgery or laser ablation-but structural targets, particularly in patients without clear lesions, are largely based on functional mapping via intracranial EEG. Unfortunately, the relationship between structural and functional connectivity in the seizing brain is poorly understood. In this study, we quantify structure-function coupling, specifically between white matter connections and intracranial EEG, across pre-ictal and ictal periods in 45 seizures from nine patients with unilateral drug-resistant focal epilepsy. We use high angular resolution diffusion imaging (HARDI) tractography to construct structural connectivity networks and correlate these networks with time-varying broadband and frequency-specific functional networks derived from coregistered intracranial EEG. Across all frequency bands, we find significant increases in structure-function coupling from pre-ictal to ictal periods. We demonstrate that short-range structural connections are primarily responsible for this increase in coupling. Finally, we find that spatiotemporal patterns of structure-function coupling are highly stereotyped for each patient. These results suggest that seizures harness the underlying structural connectome as they propagate. Mapping the relationship between structural and functional connectivity in epilepsy may inform new therapies to halt seizure spread, and pave the way for targeted patient-specific interventions.

Optimization of surgical intervention outside the epileptogenic zone in the Virtual Epileptic Patient (VEP)

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.

Robust closed-loop control of spike-and-wave discharges in a thalamocortical computational model of absence epilepsy

In this paper, we investigate the abatement of spike-and-wave discharges in a thalamocortical model using a closed-loop brain stimulation method. We first explore the complex states and various transitions in the thalamocortical computational model of absence epilepsy by using bifurcation analysis. We demonstrate that the Hopf and double cycle bifurcations are the key dynamical mechanisms of the experimental observed bidirectional communications during absence seizures through top-down cortical excitation and thalamic feedforward inhibition. Then, we formulate the abatement of epileptic seizures to a closed-loop tracking control problem. Finally, we propose a neural network based sliding mode feedback control system to drive the dynamics of pathological cortical area to track the desired normal background activities. The control system is robust to uncertainties and disturbances, and its stability is guaranteed by Lyapunov stability theorem. Our results suggest that the seizure abatement can be modeled as a tracking control problem and solved by a robust closed-loop control method, which provides a promising brain stimulation strategy.

WONOEP appraisal: Network concept from an imaging perspective

Neuroimaging techniques applied to a variety of organisms-from zebrafish, to rodents to humans-can offer valuable insights into neuronal network properties and their dysfunction in epilepsy. A wide range of imaging methods used to monitor neuronal circuits and networks during evoked seizures in animal models and advances in functional magnetic resonance imaging (fMRI) applied to patients with epilepsy were discussed during the XIV Workshop on Neurobiology of Epilepsy (XIV WONOEP) organized in 2017 by the Neurobiology Commission of the International League Against Epilepsy (ILAE). We review the growing number of technological approaches developed, as well as the current state of knowledge gained from studies applying these advanced imaging approaches to epilepsy research.

The role that choice of model plays in predictions for epilepsy surgery

Mathematical modelling has been widely used to predict the effects of perturbations to brain networks. An important example is epilepsy surgery, where the perturbation in question is the removal of brain tissue in order to render the patient free of seizures. Different dynamical models have been proposed to represent transitions to ictal states in this context. However, our choice of which mathematical model to use to address this question relies on making assumptions regarding the mechanism that defines the transition from background to the seizure state. Since these mechanisms are unknown, it is important to understand how predictions from alternative dynamical descriptions compare. Herein we evaluate to what extent three different dynamical models provide consistent predictions for the effect of removing nodes from networks. We show that for small, directed, connected networks the three considered models provide consistent predictions. For larger networks, predictions are shown to be less consistent. However consistency is higher in networks that have sufficiently large differences in ictogenicity between nodes. We further demonstrate that heterogeneity in ictogenicity across nodes correlates with variability in the number of connections for each node.

Big data analytics for personalized medicine

Big Data are radically changing biomedical research. The unprecedented advances in automated collection of large-scale molecular and clinical data pose major challenges to data analysis and interpretation, calling for the development of new computational approaches. The creation of powerful systems for the effective use of biomedical Big Data in Personalized Medicine (a.k.a. Precision Medicine) will require significant scientific and technical developments, including infrastructure, engineering, project and financial management. We review here how the evolution of data-driven methods offers the possibility to address many of these problems, guiding the formulation of hypotheses on systems functioning and the generation of mechanistic models, and facilitating the design of clinical procedures in Personalized Medicine.

Quantitative System Pharmacology (QSP): An Integrative Framework for paradigm change in the treatment of the first-episode schizophrenia

Despite the lack of progress in the curative treatment of mental illness, especially schizophrenia, the accumulation of neuroscience data over the past decade suggests the re-conceptualization of schizophrenia. With the advent of new biomarkers and cognitive tools, new neuroscience technologies such as functional dynamic connectivity and the identification of subtle clinical features; it is now possible to detect early stages at risk or prodromes of a first psychotic episode. Current concepts reconceptualizes schizophrenia as a neurodevelopmental disorder at early onset, with polygenic risk and only symptomatic treatment for positive symptoms at this time. The use of such technologies in the future suggests new diagnostic and therapeutic options. Next steps include new pharmacological perspectives and potential contributions of new technologies such as quantitative system pharmacology brain computational modeling approach.

Neuroimaging and connectomics of drug-resistant epilepsy at multiple scales: From focal lesions to macroscale networks

Epilepsy is among the most common chronic neurologic disorders, with 30%-40% of patients having seizures despite antiepileptic drug treatment. The advent of brain imaging and network analyses has greatly improved the understanding of this condition. In particular, developments in magnetic resonance imaging (MRI) have provided measures for the noninvasive characterization and detection of lesions causing epilepsy. MRI techniques can probe structural and functional connectivity, and network analyses have shaped our understanding of whole-brain anomalies associated with focal epilepsies. This review considers the progress made by neuroimaging and connectomics in the study of drug-resistant epilepsies due to focal substrates, particularly temporal lobe epilepsy related to mesiotemporal sclerosis and extratemporal lobe epilepsies associated with malformations of cortical development. In these disorders, there is evidence of widespread disturbances of structural and functional connectivity that may contribute to the clinical and cognitive prognosis of individual patients. It is hoped that studying the interplay between macroscale network anomalies and lesional profiles will improve our understanding of focal epilepsies and assist treatment choices.

Controlling seizure propagation in large-scale brain networks

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.

Conflicting emergences. Weak vs. strong emergence for the modelling of brain function

The concept of "emergence" has become commonplace in the modelling of complex systems, both natural and man-made; a functional property" emerges" from a system when it cannot be readily explained by the properties of the system's sub-units. A bewildering array of adaptive and sophisticated behaviours can be observed from large ensembles of elementary agents such as ant colonies, bird flocks or by the interactions of elementary material units such as molecules or weather elements. Ultimately, emergence has been adopted as the ontological support of a number of attempts to model brain function. This manuscript aims to clarify the ontology of emergence and delve into its many facets, particularly into its "strong" and "weak" versions that underpin two different approaches to the modelling of behaviour. The first group of models is here represented by the "free energy" principle of brain function and the "integrated information theory" of consciousness. The second group is instead represented by computational models such as oscillatory networks that use mathematical scalable representations to generate emergent behaviours and are then able to bridge neurobiology with higher mental functions. Drawing on the epistemological literature, we observe that due to their loose mechanistic links with the underlying biology, models based on strong forms of emergence are at risk of metaphysical implausibility. This, in practical terms, translates into the over determination that occurs when the proposed model becomes only one of a large set of possible explanations for the observable phenomena. On the other hand, computational models that start from biologically plausible elementary units, hence are weakly emergent, are not limited by ontological faults and, if scalable and able to realistically simulate the hierarchies of brain output, represent a powerful vehicle for future neuroscientific research programmes.

Computer modelling of connectivity change suggests epileptogenesis mechanisms in idiopathic generalised epilepsy

Patients with idiopathic generalised epilepsy (IGE) typically have normal conventional magnetic resonance imaging (MRI), hence diagnosis based on MRI is challenging. Anatomical abnormalities underlying brain dysfunctions in IGE are unclear and their relation to the pathomechanisms of epileptogenesis is poorly understood.

In this study, we applied connectometry, an advanced quantitative neuroimaging technique for investigating localised changes in white-matter tissues in vivo. Analysing white matter structures of 32 subjects we incorporated our in vivo findings in a computational model of seizure dynamics to suggest a plausible mechanism of epileptogenesis.

Patients with IGE have significant bilateral alterations in major white-matter fascicles. In the cingulum, fornix, and superior longitudinal fasciculus, tract integrity is compromised, whereas in specific parts of tracts between thalamus and the precentral gyrus, tract integrity is enhanced in patients. Combining these alterations in a logistic regression model, we computed the decision boundary that discriminated patients and controls. The computational model, informed with the findings on the tract abnormalities, specifically highlighted the importance of enhanced cortico-reticular connections along with impaired cortico-cortical connections in inducing pathological seizure-like dynamics.

We emphasise taking directionality of brain connectivity into consideration towards understanding the pathological mechanisms; this is possible by combining neuroimaging and computational modelling. Our imaging evidence of structural alterations suggest the loss of cortico-cortical and enhancement of cortico-thalamic fibre integrity in IGE. We further suggest that impaired connectivity from cortical regions to the thalamic reticular nucleus offers a therapeutic target for selectively modifying the brain circuit for reversing the mechanisms leading to epileptogenesis.

•

Significant focal alterations along major white-matter fascicles in IGE patients are characterised.

•

Increased white matter integrity found in thalamo-cortical connections.

•

Decreased white matter integrity found in cortico-cortical connections.

•

Disease mechanism is investigated by combining the neuroimaging findings with a dynamical model of seizure activity.

•

Model implicates cortical projections to the thalamic reticular nucleus in IGE.

What Can Computational Models Contribute to Neuroimaging Data Analytics?

Over the past years, nonlinear dynamical models have significantly contributed to the general understanding of brain activity as well as brain disorders. Appropriately validated and optimized mathematical models can be used to mechanistically explain properties of brain structure and neuronal dynamics observed from neuroimaging data. A thorough exploration of the model parameter space and hypothesis testing with the methods of nonlinear dynamical systems and statistical physics can assist in classification and prediction of brain states. On the one hand, such a detailed investigation and systematic parameter variation are hardly feasible in experiments and data analysis. On the other hand, the model-based approach can establish a link between empirically discovered phenomena and more abstract concepts of attractors, multistability, bifurcations, synchronization, noise-induced dynamics, etc. Such a mathematical description allows to compare and differentiate brain structure and dynamics in health and disease, such that model parameters and dynamical regimes may serve as additional biomarkers of brain states and behavioral modes. In this perspective paper we first provide very brief overview of the recent progress and some open problems in neuroimaging data analytics with emphasis on the resting state brain activity. We then focus on a few recent contributions of mathematical modeling to our understanding of the brain dynamics and model-based approaches in medicine. Finally, we discuss the question stated in the title. We conclude that incorporating computational models in neuroimaging data analytics as well as in translational medicine could significantly contribute to the progress in these fields.

Seizure pathways: A model-based investigation

We present the results of a model inversion algorithm for electrocorticography (ECoG) data recorded during epileptic seizures. The states and parameters of neural mass models were tracked during a total of over 3000 seizures from twelve patients with focal epilepsy. These models provide an estimate of the effective connectivity within intracortical circuits over the time course of seizures. Observing the dynamics of effective connectivity provides insight into mechanisms of seizures. Estimation of patients seizure dynamics revealed: 1) a highly stereotyped pattern of evolution for each patient, 2) distinct sub-groups of onset mechanisms amongst patients, and 3) different offset mechanisms for long and short seizures. Stereotypical dynamics suggest that, once initiated, seizures follow a deterministic path through the parameter space of a neural model. Furthermore, distinct sub-populations of patients were identified based on characteristic motifs in the dynamics at seizure onset. There were also distinct patterns between long and short duration seizures that were related to seizure offset. Understanding how these different patterns of seizure evolution arise may provide new insights into brain function and guide treatment for epilepsy, since specific therapies may have preferential effects on the various parameters that could potentially be individualized. Methods that unite computational models with data provide a powerful means to generate testable hypotheses for further experimental research. This work provides a demonstration that the hidden connectivity parameters of a neural mass model can be dynamically inferred from data. Our results underscore the power of theoretical models to inform epilepsy management. It is our hope that this work guides further efforts to apply computational models to clinical data.

On the nature and use of models in network neuroscience

Network theory provides an intuitively appealing framework for studying relationships among interconnected brain mechanisms and their relevance to behaviour. As the space of its applications grows, so does the diversity of meanings of the term network model. This diversity can cause confusion, complicate efforts to assess model validity and efficacy, and hamper interdisciplinary collaboration. In this Review, we examine the field of network neuroscience, focusing on organizing principles that can help overcome these challenges. First, we describe the fundamental goals in constructing network models. Second, we review the most common forms of network models, which can be described parsimoniously along the following three primary dimensions: from data representations to first-principles theory; from biophysical realism to functional phenomenology; and from elementary descriptions to coarse-grained approximations. Third, we draw on biology, philosophy and other disciplines to establish validation principles for these models. We close with a discussion of opportunities to bridge model types and point to exciting frontiers for future pursuits.

Decision Support System for Seizure Onset Zone Localization Based on Channel Ranking and High-Frequency EEG Activity

Interictal high-frequency oscillations (HFO) are a promising biomarker that can help define the seizure onset zone (SOZ) and predict the surgical outcome after the epilepsy surgery. The utility of HFO in planning the surgery, though, is unclear. Reasons include the variability of the HFO across patients and brain regions and the influence of the sleep-wake cycle, which causes large fluctuations in the ratio between the HFO observed in SOZ and non-SOZ regions. To cope with these limitations, a rank-based solution is proposed to identify the SOZ by using the HFO in multichannel intracranial EEG. A time-varying index of the epileptic susceptibility of the different brain areas is derived from the HFO rate and a support vector machine is applied on this index to identify the SOZ. The solution is trained and tested on separate groups of patients to avoid the use of patient-specific information and provides optimal SOZ prediction using as little as 30 min of recordings per channel (window). Tested on 14 patients with various combinations of seizure type, epilepsy etiology, and SOZ arrangement (172.7 ± 90.1 h/channel per patient and 75.6 ± 23.5 channels/patient, mean ± S.D.), our solution identified the SOZ with 0.92 ± 0.03 accuracy and 0.91 ± 0.03 area under the ROC curve (mean ± S.D.) across patients. For each patient, the window onset time was varied over 72 continuous hours and the prediction of the SOZ remained insensitive to the onset time, thus showing potential for surgery planning.

Modeling Brain Dynamics in Brain Tumor Patients Using the Virtual Brain

Presurgical planning for brain tumor resection aims at delineating eloquent tissue in the vicinity of the lesion to spare during surgery. To this end, noninvasive neuroimaging techniques such as functional MRI and diffusion-weighted imaging fiber tracking are currently employed. However, taking into account this information is often still insufficient, as the complex nonlinear dynamics of the brain impede straightforward prediction of functional outcome after surgical intervention. Large-scale brain network modeling carries the potential to bridge this gap by integrating neuroimaging data with biophysically based models to predict collective brain dynamics. As a first step in this direction, an appropriate computational model has to be selected, after which suitable model parameter values have to be determined. To this end, we simulated large-scale brain dynamics in 25 human brain tumor patients and 11 human control participants using The Virtual Brain, an open-source neuroinformatics platform. Local and global model parameters of the Reduced Wong–Wang model were individually optimized and compared between brain tumor patients and control subjects. In addition, the relationship between model parameters and structural network topology and cognitive performance was assessed. Results showed (1) significantly improved prediction accuracy of individual functional connectivity when using individually optimized model parameters; (2) local model parameters that can differentiate between regions directly affected by a tumor, regions distant from a tumor, and regions in a healthy brain; and (3) interesting associations between individually optimized model parameters and structural network topology and cognitive performance.

How stimulation frequency and intensity impact on the long-lasting effects of coordinated reset stimulation

Several brain diseases are characterized by abnormally strong neuronal synchrony. Coordinated Reset (CR) stimulation was computationally designed to specifically counteract abnormal neuronal synchronization processes by desynchronization. In the presence of spike-timing-dependent plasticity (STDP) this may lead to a decrease of synaptic excitatory weights and ultimately to an anti-kindling, i.e. unlearning of abnormal synaptic connectivity and abnormal neuronal synchrony. The long-lasting desynchronizing impact of CR stimulation has been verified in pre-clinical and clinical proof of concept studies. However, as yet it is unclear how to optimally choose the CR stimulation frequency, i.e. the repetition rate at which the CR stimuli are delivered. This work presents the first computational study on the dependence of the acute and long-term outcome on the CR stimulation frequency in neuronal networks with STDP. For this purpose, CR stimulation was applied with Rapidly Varying Sequences (RVS) as well as with Slowly Varying Sequences (SVS) in a wide range of stimulation frequencies and intensities. Our findings demonstrate that acute desynchronization, achieved during stimulation, does not necessarily lead to long-term desynchronization after cessation of stimulation. By comparing the long-term effects of the two different CR protocols, the RVS CR stimulation turned out to be more robust against variations of the stimulation frequency. However, SVS CR stimulation can obtain stronger anti-kindling effects. We revealed specific parameter ranges that are favorable for long-term desynchronization. For instance, RVS CR stimulation at weak intensities and with stimulation frequencies in the range of the neuronal firing rates turned out to be effective and robust, in particular, if no closed loop adaptation of stimulation parameters is (technically) available. From a clinical standpoint, this may be relevant in the context of both invasive as well as non-invasive CR stimulation.

Personalized brain network models for assessing structure-function relationships

Many recent efforts in computational modeling of macro-scale brain dynamics have begun to take a data-driven approach by incorporating structural and/or functional information derived from subject data. Here, we discuss recent work using personalized brain network models to study structure-function relationships in human brains. We describe the steps necessary to build such models and show how this computational approach can provide previously unobtainable information through the ability to perform virtual experiments. Finally, we present examples of how personalized brain network models can be used to gain insight into the effects of local stimulation and improve surgical outcomes in epilepsy.

Credibility, Replicability, and Reproducibility in Simulation for Biomedicine and Clinical Applications in Neuroscience

Modeling and simulation in computational neuroscience is currently a research enterprise to better understand neural systems. It is not yet directly applicable to the problems of patients with brain disease. To be used for clinical applications, there must not only be considerable progress in the field but also a concerted effort to use best practices in order to demonstrate model credibility to regulatory bodies, to clinics and hospitals, to doctors, and to patients. In doing this for neuroscience, we can learn lessons from long-standing practices in other areas of simulation (aircraft, computer chips), from software engineering, and from other biomedical disciplines. In this manuscript, we introduce some basic concepts that will be important in the development of credible clinical neuroscience models: reproducibility and replicability; verification and validation; model configuration; and procedures and processes for credible mechanistic multiscale modeling. We also discuss how garnering strong community involvement can promote model credibility. Finally, in addition to direct usage with patients, we note the potential for simulation usage in the area of Simulation-Based Medical Education, an area which to date has been primarily reliant on physical models (mannequins) and scenario-based simulations rather than on numerical simulations.

Interpreting and Utilising Intersubject Variability in Brain Function

We consider between-subject variance in brain function as data rather than noise. We describe variability as a natural output of a noisy plastic system (the brain) where each subject embodies a particular parameterisation of that system. In this context, variability becomes an opportunity to: (i) better characterise typical versus atypical brain functions; (ii) reveal the different cognitive strategies and processing networks that can sustain similar tasks; and (iii) predict recovery capacity after brain damage by taking into account both damaged and spared processing pathways. This has many ramifications for understanding individual learning preferences and explaining the wide differences in human abilities and disabilities. Understanding variability boosts the translational potential of neuroimaging findings, in particular in clinical and educational neuroscience.

Predicting the spatiotemporal diversity of seizure propagation and termination in human focal epilepsy

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.

Seizure development in the acute intrahippocampal epileptic focus

Currently, an epileptic seizure is considered to involve a temporary network that exists for a finite period of time. Formation of this network evolves through spread of epileptiform activity from a seizure onset zone (SOZ). Propagation of seizures evoked by kainic acid injection in hippocampus to different brain areas was analyzed at macro- and micro-intervals. The mean latency of seizure occurrence in different brain areas varied between 0.5 sec and 85 sec (mean 14.9 ± 14.5 (SD)), and it increased after each consecutive seizure in areas located contralateral to the area of injection, but not in the ipsilateral sites. We have shown that only 41% of epileptic individual events in target brain areas were driven by epileptic events generated in the SOZ once the seizure began. Fifty-nine percent of epileptiform events in target areas occurred one millisecond before or after events in the SOZ. These data illustrate that during seizure maintenance, only some individual epileptiform events in areas outside of SOZ could be consistently triggered by the SOZ; and the majority must be triggered by a driver located outside the SOZ or brain areas involved in ictal activity could be coupled to each other via an unknown mechanism such as stochastic resonance.

From Maps to Multi-dimensional Network Mechanisms of Mental Disorders

The development of advanced neuroimaging techniques and their deployment in large cohorts has enabled an assessment of functional and structural brain network architecture at an unprecedented level of detail. Across many temporal and spatial scales, network neuroscience has emerged as a central focus of intellectual efforts, seeking meaningful descriptions of brain networks and explanatory sets of network features that underlie circuit function in health and dysfunction in disease. However, the tools of network science commonly deployed provide insight into brain function at a fundamentally descriptive level, often failing to identify (patho-)physiological mechanisms that link system-level phenomena to the multiple hierarchies of brain function. Here we describe recently developed techniques stemming from advances in complex systems and network science that have the potential to overcome this limitation, thereby contributing mechanistic insights into neuroanatomy, functional dynamics, and pathology. Finally, we build on the Research Domain Criteria framework, highlighting the notion that mental illnesses can be conceptualized as dysfunctions of neural circuitry present across conventional diagnostic boundaries, to sketch how network-based methods can be combined with pharmacological, intermediate phenotype, genetic, and magnetic stimulation studies to probe mechanisms of psychopathology.

Commonalities in epileptogenic processes from different acute brain insults: Do they translate?

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.