Epilepsies as dynamical diseases of brain systems: basic models of the transition between normal and epileptic activity

Optimal input reverberation and homeostatic self-organization toward the edge of synchronization

Transient or partial synchronization can be used to do computations, although a fully synchronized network is sometimes related to the onset of epileptic seizures. Here, we propose a homeostatic mechanism that is capable of maintaining a neuronal network at the edge of a synchronization transition, thereby avoiding the harmful consequences of a fully synchronized network. We model neurons by maps since they are dynamically richer than integrate-and-fire models and more computationally efficient than conductance-based approaches. We first describe the synchronization phase transition of a dense network of neurons with different tonic spiking frequencies coupled by gap junctions. We show that at the transition critical point, inputs optimally reverberate through the network activity through transient synchronization. Then, we introduce a local homeostatic dynamic in the synaptic coupling and show that it produces a robust self-organization toward the edge of this phase transition. We discuss the potential biological consequences of this self-organization process, such as its relation to the Brain Criticality hypothesis, its input processing capacity, and how its malfunction could lead to pathological synchronization and the onset of seizure-like activity.

Multi-scale modelling of the epileptic brain: advantages of computational therapy exploration

Objective: Epilepsy is a complex disease spanning across multiple scales, from ion channels in neurons to neuronal circuits across the entire brain. Over the past decades, computational models have been used to describe the pathophysiological activity of the epileptic brain from different aspects. Traditionally, each computational model can aid in optimizing therapeutic interventions, therefore, providing a particular view to design strategies for treating epilepsy. As a result, most studies are concerned with generating specific models of the epileptic brain that can help us understand the certain machinery of the pathological state. Those specific models vary in complexity and biological accuracy, with system-level models often lacking biological details.Approach: Here, we review various types of computational model of epilepsy and discuss their potential for different therapeutic approaches and scenarios, including drug discovery, surgical strategies, brain stimulation, and seizure prediction. We propose that we need to consider an integrated approach with a unified modelling framework across multiple scales to understand the epileptic brain. Our proposal is based on the recent increase in computational power, which has opened up the possibility of unifying those specific epileptic models into simulations with an unprecedented level of detail.Main results: A multi-scale epilepsy model can bridge the gap between biologically detailed models, used to address molecular and cellular questions, and brain-wide models based on abstract models which can account for complex neurological and behavioural observations.Significance: With these efforts, we move toward the next generation of epileptic brain models capable of connecting cellular features, such as ion channel properties, with standard clinical measures such as seizure severity.

Response to photic stimulation as a measure of cortical excitability in epilepsy patients

Studying states and state transitions in the brain is challenging due to nonlinear, complex dynamics. In this research, we analyze the brain's response to non-invasive perturbations. Perturbation techniques offer a powerful method for studying complex dynamics, though their translation to human brain data is under-explored. This method involves applying small inputs, in this case via photic stimulation, to a system and measuring its response. Sensitivity to perturbations can forewarn a state transition. Therefore, biomarkers of the brain's perturbation response or "cortical excitability" could be used to indicate seizure transitions. However, perturbing the brain often involves invasive intracranial surgeries or expensive equipment such as transcranial magnetic stimulation (TMS) which is only accessible to a minority of patient groups, or animal model studies. Photic stimulation is a widely used diagnostic technique in epilepsy that can be used as a non-invasive perturbation paradigm to probe brain dynamics during routine electroencephalography (EEG) studies in humans. This involves changing the frequency of strobing light, sometimes triggering a photo-paroxysmal response (PPR), which is an electrographic event that can be studied as a state transition to a seizure state. We investigate alterations in the response to these perturbations in patients with genetic generalized epilepsy (GGE), with (n = 10) and without (n = 10) PPR, and patients with psychogenic non-epileptic seizures (PNES; n = 10), compared to resting controls (n = 10). Metrics of EEG time-series data were evaluated as biomarkers of the perturbation response including variance, autocorrelation, and phase-based synchrony measures. We observed considerable differences in all group biomarker distributions during stimulation compared to controls. In particular, variance and autocorrelation demonstrated greater changes in epochs close to PPR transitions compared to earlier stimulation epochs. Comparison of PPR and spontaneous seizure morphology found them indistinguishable, suggesting PPR is a valid proxy for seizure dynamics. Also, as expected, posterior channels demonstrated the greatest change in synchrony measures, possibly reflecting underlying PPR pathophysiologic mechanisms. We clearly demonstrate observable changes at a group level in cortical excitability in epilepsy patients as a response to perturbation in EEG data. Our work re-frames photic stimulation as a non-invasive perturbation paradigm capable of inducing measurable changes to brain dynamics.

Heterogeneous correlate and potential diagnostic biomarker of tinnitus based on nonlinear dynamics of resting-state EEG recordings

Tinnitus is a heterogeneous condition of hearing a rattling sound when there is no auditory stimulus. This rattling sound is associated with abnormal synchronous oscillations in auditory and non-auditory cortical areas. Since tinnitus is a highly heterogeneous condition with no objective detection criteria, it is necessary to search for indicators that can be compared between and within participants for diagnostic purposes. This study introduces heterogeneous though comparable indicators of tinnitus through investigation of spontaneous fluctuations in resting-state brain dynamics. The proposed approach uses nonlinear measures of chaos theory, to detect tinnitus and cross correlation patterns to reflect many of the previously reported neural correlates of tinnitus. These indicators may serve as effective measures of tinnitus risk even at early ages before any symptom is reported. The approach quantifies differences in oscillatory brain dynamics of tinnitus and normal subjects. It demonstrates that the left temporal areas of subjects with tinnitus exhibit larger lyapunov exponent indicating irregularity of brain dynamics in these regions. More complex dynamics is further recognized in tinnitus cases through entropy. We use this evidence to distinguish tinnitus patients from normal participants. Besides, we illustrate that certain anticorrelation patterns appear in these nonlinear measures across temporal and frontal areas in the brain perhaps corresponding to increased/decreased connectivity in certain brain networks and a shift in the balance of excitation and inhibition in tinnitus. Additionally, the main correlations are lost in tinnitus participants compared to control group suggesting involvement of distinct neural mechanisms in generation and persistence of tinnitus.

Alterations in brain network functional connectivity and topological properties in DRE patients

Objective

The study aimed to find the difference in functional network topology on interictal electroencephalographic (EEG) between patients with drug-resistant epilepsy (DRE) and healthy people.

Methods

We retrospectively analyzed the medical records as well as EEG data of ten patients with DRE and recruited five sex-age-matched healthy controls (HC group). Each participant remained awake while undergoing video-electroencephalography (vEEG) monitoring. After excluding data that contained abnormal discharges, we screened EEG segments that were free of artifacts and put them together into 20-min segments. The screened data was bandpass filtered to different frequency bands (delta, theta, alpha, beta, and gamma). The weighted phase lag index (wPLI) and the network properties were calculated to evaluate changes in the topology of the functional network. Finally, the results were statistically analyzed, and the false discovery rate (FDR) was used to correct for differences after multiple comparisons.

Results

In the full frequency band (0.5-45 Hz), the functional connectivity in the DRE group during the interictal period was significantly lower than that in the HC group (p < 0.05). Compared to the HC group, in the full frequency band, the DRE group exhibited significantly decreased clustering coefficient (CC), node degree (D), and global efficiency (GE), while the characteristic path length (CPL) significantly increased (p < 0.05). In the sub-frequency bands, the functional connectivity of the DRE group was significantly lower than that of the HC group in the delta band but higher in the alpha, beta, and gamma bands (p < 0.05). The statistical results of network properties revealed that in the delta band, the DRE group had significantly decreased values for D, CC, and GE, but in the alpha, beta, and gamma bands, these values were significantly increased (p < 0.05). Additionally, the CPL of the DRE group significantly increased in the delta and theta bands but significantly decreased in the alpha, beta, and gamma bands (p < 0.05).

Conclusion

The topology structure of the functional network in DRE patients was significantly changed compared with healthy people, which was reflected in different frequency bands. It provided a theoretical basis for understanding the pathological network alterations of DRE.

Kindling in humans: Does secondary epileptogenesis occur?

The relevance of secondary epileptogenesis for human epilepsy remains a controversial subject decades after it was first described in animal models. Whether or not a previously normal brain region can become independently epileptogenic through a kindling-like process has not, and cannot, be definitely proven in humans. Rather than reliance on direct experimental evidence, attempts to answering this question must depend on observational data. In this review, observations based largely upon contemporary surgical series will advance the case for secondary epileptogenesis in humans. As will be argued, hypothalamic hamartoma-related epilepsy provides the strongest case for this process; all the stages of secondary epileptogenesis can be observed. Hippocampal sclerosis (HS) is another pathology where the question of secondary epileptogenesis frequently arises, and observations from bitemporal and dual pathology series are explored. The verdict here is far more difficult to reach, in large part because of the scarcity of longitudinal cohorts; moreover, recent experimental data have challenged the claim that HS is acquired consequent to recurrent seizures. Synaptic plasticity more than seizure-induced neuronal injury is the likely mechanism of secondary epileptogenesis. Postoperative running-down phenomenon provides the best evidence that a kindling-like process occurs in some patients, evidenced by its reversal. Finally, a network perspective of secondary epileptogenesis is considered, as well as the possible role for subcortical surgical interventions.

EEG microstate quantifiers and state space descriptors during anaesthesia in patients with postoperative delirium: a descriptive analysis

Postoperative delirium is a serious sequela of surgery and surgery-related anaesthesia. One recommended method to prevent postoperative delirium is using bi-frontal EEG recording. The single, processed index of depth of anaesthesia allows the anaesthetist to avoid episodes of suppression EEG and excessively deep anaesthesia. The study data presented here were based on multichannel (19 channels) EEG recordings during anaesthesia. This enabled the analysis of various parameters of global electrical brain activity. These parameters were used to compare microstate topographies under anaesthesia with those in healthy volunteers and to analyse changes in microstate quantifiers and EEG global state space descriptors with increasing exposure to anaesthesia. Seventy-three patients from the Surgery Depth of Anaesthesia and Cognitive Outcome study (SRCTN 36437985) received intraoperative multichannel EEG recordings. Altogether, 720 min of artefact-free EEG data, including 210 min (29.2%) of suppression EEG, were analysed. EEG microstate topographies, microstate quantifiers (duration, frequency of occurrence and global field power) and the state space descriptors sigma (overall EEG power), phi (generalized frequency) and omega (number of uncorrelated brain processes) were evaluated as a function of duration of exposure to anaesthesia, suppression EEG and subsequent development of postoperative delirium. The major analyses involved covariate-adjusted linear mixed-effects models. The older (71 ± 7 years), predominantly male (60%) patients received a median exposure of 210 (range: 75-675) min of anaesthesia. During seven postoperative days, 21 patients (29%) developed postoperative delirium. Microstate topographies under anaesthesia resembled topographies from healthy and much younger awake persons. With increasing duration of exposure to anaesthesia, single microstate quantifiers progressed differently in suppression or non-suppression EEG and in patients with or without subsequent postoperative delirium. The most pronounced changes occurred during enduring suppression EEG in patients with subsequent postoperative delirium: duration and frequency of occurrence of microstates C and D progressed in opposite directions, and the state space descriptors showed a pattern of declining uncorrelated brain processes (omega) combined with increasing EEG variance (sigma). With increasing exposure to general anaesthesia, multiple changes in the dynamics of microstates and global EEG parameters occurred. These changes varied partly between suppression and non-suppression EEG and between patients with or without subsequent postoperative delirium. Ongoing suppression EEG in patients with subsequent postoperative delirium was associated with reduced network complexity in combination with increased overall EEG power. Additionally, marked changes in quantifiers in microstate C and in microstate D occurred. These putatively adverse intraoperative trajectories in global electrical brain activity may be seen as preceding and ultimately predicting postoperative delirium.

Localizing epileptogenic zones with high-frequency oscillations and directed connectivity

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.

The Roles of Potassium and Calcium Currents in the Bistable Firing Transition

Healthy brains display a wide range of firing patterns, from synchronized oscillations during slow-wave sleep to desynchronized firing during movement. These physiological activities coexist with periods of pathological hyperactivity in the epileptic brain, where neurons can fire in synchronized bursts. Most cortical neurons are pyramidal regular spiking (RS) cells with frequency adaptation and do not exhibit bursts in current-clamp experiments (in vitro). In this work, we investigate the transition mechanism of spike-to-burst patterns due to slow potassium and calcium currents, considering a conductance-based model of a cortical RS cell. The joint influence of potassium and calcium ion channels on high synchronous patterns is investigated for different synaptic couplings (gsyn) and external current inputs (I). Our results suggest that slow potassium currents play an important role in the emergence of high-synchronous activities, as well as in the spike-to-burst firing pattern transitions. This transition is related to the bistable dynamics of the neuronal network, where physiological asynchronous states coexist with pathological burst synchronization. The hysteresis curve of the coefficient of variation of the inter-spike interval demonstrates that a burst can be initiated by firing states with neuronal synchronization. Furthermore, we notice that high-threshold (IL) and low-threshold (IT) ion channels play a role in increasing and decreasing the parameter conditions (gsyn and I) in which bistable dynamics occur, respectively. For high values of IL conductance, a synchronous burst appears when neurons are weakly coupled and receive more external input. On the other hand, when the conductance IT increases, higher coupling and lower I are necessary to produce burst synchronization. In light of our results, we suggest that channel subtype-specific pharmacological interactions can be useful to induce transitions from pathological high bursting states to healthy states.

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

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

The Role of Potassium and Calcium Currents in the Bistable Firing Transition

Healthy brains display a wide range of firing patterns, from synchronized oscillations during slowwave sleep to desynchronized firing during movement. These physiological activities coexist with periods of pathological hyperactivity in the epileptic brain, where neurons can fire in synchronized bursts. Most cortical neurons are pyramidal regular spiking cells (RS) with frequency adaptation and do not exhibit bursts in current-clamp experiments ( in vitro ). In this work, we investigate the transition mechanism of spike-to-burst patterns due to slow potassium and calcium currents, considering a conductance-based model of a cortical RS cell. The joint influence of potassium and calcium ion channels on high synchronous patterns is investigated for different synaptic couplings ( g syn ) and external current inputs ( I ). Our results suggest that slow potassium currents play an important role in the emergence of high-synchronous activities, as well as in the spike-to-burst firing pattern transitions. This transition is related to bistable dynamics of the neuronal network, where physiological asynchronous states coexist with pathological burst synchronization. The hysteresis curve of the coefficient of variation of the inter-spike interval demonstrates that a burst can be initiated by firing states with neuronal synchronization. Furthermore, we notice that high-threshold ( I L ) and low-threshold ( I T ) ion channels play a role in increasing and decreasing the parameter conditions ( g syn and I ) in which bistable dynamics occur, respectively. For high values of I L conductance, a synchronous burst appears when neurons are weakly coupled and receive more external input. On the other hand, when the conductance I T increases, higher coupling and lower I are necessary to produce burst synchronization. In light of our results, we suggest that channel subtype-specific pharmacological interactions can be useful to induce transitions from pathological high bursting states to healthy states.

Switching state-space modeling of neural signal dynamics

Linear parametric state-space models are a ubiquitous tool for analyzing neural time series data, providing a way to characterize the underlying brain dynamics with much greater statistical efficiency than non-parametric data analysis approaches. However, neural time series data are frequently time-varying, exhibiting rapid changes in dynamics, with transient activity that is often the key feature of interest in the data. Stationary methods can be adapted to time-varying scenarios by employing fixed-duration windows under an assumption of quasi-stationarity. But time-varying dynamics can be explicitly modeled by switching state-space models, i.e., by using a pool of state-space models with different dynamics selected by a probabilistic switching process. Unfortunately, exact solutions for state inference and parameter learning with switching state-space models are intractable. Here we revisit a switching state-space model inference approach first proposed by Ghahramani and Hinton. We provide explicit derivations for solving the inference problem iteratively after applying a variational approximation on the joint posterior of the hidden states and the switching process. We introduce a novel initialization procedure using an efficient leave-one-out strategy to compare among candidate models, which significantly improves performance compared to the existing method that relies on deterministic annealing. We then utilize this state inference solution within a generalized expectation-maximization algorithm to estimate model parameters of the switching process and the linear state-space models with dynamics potentially shared among candidate models. We perform extensive simulations under different settings to benchmark performance against existing switching inference methods and further validate the robustness of our switching inference solution outside the generative switching model class. Finally, we demonstrate the utility of our method for sleep spindle detection in real recordings, showing how switching state-space models can be used to detect and extract transient spindles from human sleep electroencephalograms in an unsupervised manner.

On temporal scale-free non-periodic stimulation and its mechanisms as an infinite improbability drive of the brain's functional connectogram

Rationalized development of electrical stimulation (ES) therapy is of paramount importance. Not only it will foster new techniques and technologies with increased levels of safety, efficacy, and efficiency, but it will also facilitate the translation from basic research to clinical practice. For such endeavor, design of new technologies must dialogue with state-of-the-art neuroscientific knowledge. By its turn, neuroscience is transitioning-a movement started a couple of decades earlier-into adopting a new conceptual framework for brain architecture, in which time and thus temporal patterns plays a central role in the neuronal representation of sampled data from the world. This article discusses how neuroscience has evolved to understand the importance of brain rhythms in the overall functional architecture of the nervous system and, consequently, that neuromodulation research should embrace this new conceptual framework. Based on such support, we revisit the literature on standard (fixed-frequency pulsatile stimuli) and mostly non-standard patterns of ES to put forward our own rationale on how temporally complex stimulation schemes may impact neuromodulation strategies. We then proceed to present a low frequency, on average (thus low energy), scale-free temporally randomized ES pattern for the treatment of experimental epilepsy, devised by our group and termed NPS (Non-periodic Stimulation). The approach has been shown to have robust anticonvulsant effects in different animal models of acute and chronic seizures (displaying dysfunctional hyperexcitable tissue), while also preserving neural function. In our understanding, accumulated mechanistic evidence suggests such a beneficial mechanism of action may be due to the natural-like characteristic of a scale-free temporal pattern that may robustly compete with aberrant epileptiform activity for the recruitment of neural circuits. Delivering temporally patterned or random stimuli within specific phases of the underlying oscillations (i.e., those involved in the communication within and across brain regions) could both potentiate and disrupt the formation of neuronal assemblies with random probability. The usage of infinite improbability drive here is obviously a reference to the "The Hitchhiker's Guide to the Galaxy" comedy science fiction classic, written by Douglas Adams. The parallel is that dynamically driving brain functional connectogram, through neuromodulation, in a manner that would not favor any specific neuronal assembly and/or circuit, could re-stabilize a system that is transitioning to fall under the control of a single attractor. We conclude by discussing future avenues of investigation and their potentially disruptive impact on neurotechnology, with a particular interest in NPS implications in neural plasticity, motor rehabilitation, and its potential for clinical translation.

Glycolytic Wave Patterns in a Simple Reaction-diffusion System with Inhomogeneous Influx: Dynamic Transitions

An inhomogeneous profile of chemostatted species generates a rich variety of patterns in glycolytic waves depicted in a Selkov reaction-diffusion framework here. A key role played by diffusion amplitude and symmetry in the chemostatted species profile in dictating the fate of local spatial dynamics involving periodic, quasiperiodic, and chaotic patterns and transitions among them are investigated systematically. More importantly, various dynamic transitions, including wave propagation direction changes, are illustrated in interesting situations. Besides numerical results, our analytical formulation of the amplitude equation connecting complex Ginzburg-Landau and Lambda-omega representation shed light on the phase dynamics of the system. This systematic study of the glycolytic reaction-diffusion wave is in line with previous experimental results in open spatial reactor and will provide a knowledge about the dynamics that shape and control biological information processing and related phenomena.

Combined effect of stimulation and electromagnetic induction on absence seizure inhibition in coupled thalamocortical circuits

Deep brain stimulation (DBS) and electromagnetic induction are new techniques that are increasingly used in modern epilepsy treatments; however, the mechanism of action remains unclear. In this study, we constructed a bidirectional-coupled cortico-thalamic model, based on which we proposed three regulation schemes: isolated regulation of DBS, isolated regulation of electromagnetic induction and combined regulation of the previous two. In particular, we introduced DBS with a lower amplitude and considered the influence of electromagnetic induction caused by the transmembrane current on the membrane potential. The most striking finding of this study is that the three therapeutic schemes could effectively control abnormal discharge, and combined regulation could reduce the occurrence of epileptic seizures more effectively. The present study bridges the gap between the bidirectional coupling model and combined control. In this way, the damage induced by electrical stimulation of the patient's brain tissue could be reduced, and the abnormal physiological discharge pattern of the cerebral cortex was simultaneously regulated by different techniques. This work opens new avenues for improving brain dysfunction in patients with epilepsy, expands ideas for promoting the development of neuroscience and is meaningful for improving the health of modern society and developing the field of science.

Pro-Ictal State in Human Temporal Lobe Epilepsy

BACKGROUND: Studies of continuous electroencephalography (EEG) suggest that seizures in individuals with focal-onset epilepsies preferentially occur during periods of heightened risk, typified by pathologic brain activities, termed pro-ictal states; however, the presence of (pathologic) pro-ictal states among a plethora of otherwise physiologic (e.g., sleep–wake cycle) states has not been established. METHODS: We studied a prospective, consecutive series of 15 patients with temporal lobe epilepsy who underwent limbic thalamic recordings in addition to routine (cortical) intracranial EEG for seizure localization. For each participant, pro-ictal (45 minutes before seizure onset) and interictal (4 hours removed from all seizures) EEG segments were divided into 10-minute, nonoverlapping windows, which were randomly distributed into training and validation cohorts in a 1:1 ratio. A deep neural classifier was applied to distinguish pro-ictal from interictal brain activities in a patient-specific fashion. RESULTS: We analyzed 1800 patient-hours of continuous thalamocortical EEG. Distinct pro-ictal states were detected in each participant. The median area under the receiver-operating characteristic curve of the classifier was 0.92 (interquartile range, 0.90–0.96). Pro-ictal states were distinguished at least 45 minutes before seizure onset in 13 of 15 participants; in 2 of 15 participants, they were distinguished up to 35 minutes prior. CONCLUSIONS: On the basis of thalamocortical EEG, pro-ictal states — pathologic brain activities during periods of heightened seizure risk — could be identified in patients with temporal lobe epilepsy and were detected, in our small sample, more than one half hour before seizure onset.

Diverse nature of interictal oscillations: EEG-based biomarkers in epilepsy

Interictal electroencephalogram (EEG) patterns, including high-frequency oscillations (HFOs), interictal spikes (ISs), and slow wave activities (SWAs), are defined as specific oscillations between seizure events. These interictal oscillations reflect specific dynamic changes in network excitability and play various roles in epilepsy. In this review, we briefly describe the electrographic characteristics of HFOs, ISs, and SWAs in the interictal state, and discuss the underlying cellular and network mechanisms. We also summarize representative evidence from experimental and clinical epilepsy to address their critical roles in ictogenesis and epileptogenesis, indicating their potential as electrophysiological biomarkers of epilepsy. Importantly, we put forwards some perspectives for further research in the field.

Unsupervised EEG preictal interval identification in patients with drug-resistant epilepsy

Typical seizure prediction models aim at discriminating interictal brain activity from pre-seizure electrographic patterns. Given the lack of a preictal clinical definition, a fixed interval is widely used to develop these models. Recent studies reporting preictal interval selection among a range of fixed intervals show inter- and intra-patient preictal interval variability, possibly reflecting the heterogeneity of the seizure generation process. Obtaining accurate labels of the preictal interval can be used to train supervised prediction models and, hence, avoid setting a fixed preictal interval for all seizures within the same patient. Unsupervised learning methods hold great promise for exploring preictal alterations on a seizure-specific scale. Multivariate and univariate linear and nonlinear features were extracted from scalp electroencephalography (EEG) signals collected from 41 patients with drug-resistant epilepsy undergoing presurgical monitoring. Nonlinear dimensionality reduction was performed for each group of features and each of the 226 seizures. We applied different clustering methods in searching for preictal clusters located until 2 h before the seizure onset. We identified preictal patterns in 90% of patients and 51% of the visually inspected seizures. The preictal clusters manifested a seizure-specific profile with varying duration (22.9 ± 21.0 min) and starting time before seizure onset (47.6 ± 27.3 min). Searching for preictal patterns on the EEG trace using unsupervised methods showed that it is possible to identify seizure-specific preictal signatures for some patients and some seizures within the same patient.

The adjustment mechanism of the spike and wave discharges in thalamic neurons: a simulation analysis

Different from many previous theoretical studies, this paper explores the regulatory mechanism of the spike and wave discharges (SWDs) in the reticular thalamic nucleus (TRN) by a dynamic computational model. We observe that the SWDs appears in the TRN by changing the coupling weights and delays in the thalamocortical circuit. The abundant poly-spikes wave discharges is also induced when the delay increases to large enough. These discharges can be inhibited by tuning the inhibitory output from the basal ganglia to the thalamus. The mechanisms of these waves can be explained in this model together with simulation results, which are different from the mechanisms in the cortex. The TRN is an important target in treating epilepsy, and the results may be a theoretical evidence for experimental study in the future.

Temporal Evolution of Multiday, Epileptic Functional Networks Prior to Seizure Occurrence

Epilepsy is one of the most common neurological disorders, characterized by the occurrence of repeated seizures. Given that epilepsy is considered a network disorder, tools derived from network neuroscience may confer the valuable ability to quantify the properties of epileptic brain networks. In this study, we use well-established brain network metrics (i.e., mean strength, variance of strength, eigenvector centrality, betweenness centrality) to characterize the temporal evolution of epileptic functional networks over several days prior to seizure occurrence. We infer the networks using long-term electroencephalographic recordings from 12 people with epilepsy. We found that brain network metrics are variable across days and show a circadian periodicity. In addition, we found that in 9 out of 12 patients the distribution of the variance of strength in the day (or even two last days) prior to seizure occurrence is significantly different compared to the corresponding distributions on all previous days. Our results suggest that brain network metrics computed fromelectroencephalographic recordings could potentially be used to characterize brain network changes that occur prior to seizures, and ultimately contribute to seizure warning systems.

Localizing targets for neuromodulation in drug-resistant epilepsy using intracranial EEG and computational model

Neuromodulation has emerged as a promising technique for the treatment of epilepsy. The target for neuromodulation is critical for the effectiveness of seizure control. About 30% of patients with drug-resistant epilepsy (DRE) fail to achieve seizure freedom after surgical intervention. It is difficult to find effective brain targets for neuromodulation in these patients because brain regions are damaged during surgery. In this study, we propose a novel approach for localizing neuromodulatory targets, which uses intracranial EEG and multi-unit computational models to simulate the dynamic behavior of epileptic networks through external stimulation. First, we validate our method on a multivariate autoregressive model and compare nine different methods of constructing brain networks. Our results show that the directed transfer function with surrogate analysis achieves the best performance. Intracranial EEGs of 11 DRE patients are further analyzed. These patients all underwent surgery. In three seizure-free patients, the localized targets are concordant with the resected regions. For the eight patients without seizure-free outcome, the localized targets in three of them are outside the resected regions. Finally, we provide candidate targets for neuromodulation in these patients without seizure-free outcome based on virtual resected epileptic network. We demonstrate the ability of our approach to locate optimal targets for neuromodulation. We hope that our approach can provide a new tool for localizing patient-specific targets for neuromodulation therapy in DRE.

Regulating absence seizures by tri-phase delay stimulation applied to globus pallidus internal

In this paper, a reduced globus pallidus internal (GPI)-corticothalamic (GCT) model is developed, and a tri-phase delay stimulation (TPDS) with sequentially applying three pulses on the GPI representing the inputs from the striatal D ₁ neurons, subthalamic nucleus (STN), and globus pallidus external (GPE), respectively, is proposed. The GPI is evidenced to control absence seizures characterized by 2 Hz-4 Hz spike and wave discharge (SWD). Hence, based on the basal ganglia-thalamocortical (BGCT) model, we firstly explore the triple effects of D _l-GPI, GPE-GPI, and STN-GPI pathways on seizure patterns. Then, using the GCT model, we apply the TPDS on the GPI to potentially investigate the alternative and improved approach if these pathways to the GPI are blocked. The results show that the striatum D ₁, GPE, and STN can indeed jointly and significantly affect seizure patterns. In particular, the TPDS can effectively reproduce the seizure pattern if the D ₁-GPI, GPE-GPI, and STN-GPI pathways are cut off. In addition, the seizure abatement can be obtained by well tuning the TPDS stimulation parameters. This implies that the TPDS can play the surrogate role similar to the modulation of basal ganglia, which hopefully can be helpful for the development of the brain-computer interface in the clinical application of epilepsy.

Tuning Neural Synchronization: The Role of Variable Oscillation Frequencies in Neural Circuits

Brain oscillations emerge during sensory and cognitive processes and have been classified into different frequency bands. Yet, even within the same frequency band and between nearby brain locations, the exact frequencies of brain oscillations can differ. These frequency differences (detuning) have been largely ignored and play little role in current functional theories of brain oscillations. This contrasts with the crucial role that detuning plays in synchronization theory, as originally derived in physical systems. Here, we propose that detuning is equally important to understand synchronization in biological systems. Detuning is a critical control parameter in synchronization, which is not only important in shaping phase-locking, but also in establishing preferred phase relations between oscillators. We review recent evidence that frequency differences between brain locations are ubiquitous and essential in shaping temporal neural coordination. With the rise of powerful experimental techniques to probe brain oscillations, the contributions of exact frequency and detuning across neural circuits will become increasingly clear and will play a key part in developing a new understanding of the role of oscillations in brain function.

Interictal dynamic network transitions in mesial temporal lobe epilepsy

OBJECTIVE

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

METHODS

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

RESULTS

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

SIGNIFICANCE

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

Epidemic models characterize seizure propagation and the effects of epilepsy surgery in individualized brain networks based on MEG and invasive EEG recordings

Epilepsy surgery is the treatment of choice for drug-resistant epilepsy patients. However, seizure-freedom is currently achieved in only 2/3 of the patients after surgery. In this study we have developed an individualized computational model based on MEG brain networks to explore seizure propagation and the efficacy of different virtual resections. Eventually, the goal is to obtain individualized models to optimize resection strategy and outcome. We have modelled seizure propagation as an epidemic process using the susceptible-infected (SI) model on individual brain networks derived from presurgical MEG. We included 10 patients who had received epilepsy surgery and for whom the surgery outcome at least one year after surgery was known. The model parameters were tuned in in order to reproduce the patient-specific seizure propagation patterns as recorded with invasive EEG. We defined a personalized search algorithm that combined structural and dynamical information to find resections that maximally decreased seizure propagation for a given resection size. The optimal resection for each patient was defined as the smallest resection leading to at least a 90% reduction in seizure propagation. The individualized model reproduced the basic aspects of seizure propagation for 9 out of 10 patients when using the resection area as the origin of epidemic spreading, and for 10 out of 10 patients with an alternative definition of the seed region. We found that, for 7 patients, the optimal resection was smaller than the resection area, and for 4 patients we also found that a resection smaller than the resection area could lead to a 100% decrease in propagation. Moreover, for two cases these alternative resections included nodes outside the resection area. Epidemic spreading models fitted with patient specific data can capture the fundamental aspects of clinically observed seizure propagation, and can be used to test virtual resections in silico. Combined with optimization algorithms, smaller or alternative resection strategies, that are individually targeted for each patient, can be determined with the ultimate goal to improve surgery outcome. MEG-based networks can provide a good approximation of structural connectivity for computational models of seizure propagation, and facilitate their clinical use.

Abnormal Static and Dynamic Functional Connectivity in Left and Right Temporal Lobe Epilepsy

Objective

Temporal lobe epilepsy (TLE) can be conceptualized as a network disease. However, the network characteristics in lateralization remain controversial.

Methods

In this study, resting-state functional MRI scans were acquired from 53 TLE patients [22 with left-side TLE (LTLE) and 31 with right-side TLE (RTLE)] and 37 matched healthy controls. We focused on the characteristics of static and dynamic functional connectivity, including static connectivity patterns and topological properties, as well as temporal properties of the dynamic connectivity state and the variability of the dynamic connectivity and network topological organization. Correlation analyses were conducted between abnormal static and dynamic properties and cognitive performances.

Results

The static functional connectivity analysis presented a significantly decreased cortical-cortical connectivity pattern and increased subcortical-cortical connectivity pattern in RTLE. The global-level network in RTLE showed a significant decrease in global efficiency. The dynamic functional connectivity analysis revealed that RTLE patients exhibited aberrant connectivity states, as well as increased variability in the subcortical-cortical connectivity. The global-level network in RTLE revealed increased variance in global efficiency and local efficiency. The static or dynamic functional connectivity in LTLE did not show any significant abnormalities. The altered dynamic properties were associated with worsening cognitive performance in language and conceptual thinking by the TLE patients.

Conclusion

Our findings demonstrated the presence of abnormalities in the static and dynamic functional connectivity of TLE patients. RTLE patients exhibited more pronounced aberrant connectivity patterns and topological properties, which might represent a mechanism for reconfiguration of brain networks in RTLE patients. These observations extended our understanding of the pathophysiological network mechanisms of TLE.

Temporally Targeted Interactions With Pathologic Oscillations as Therapeutical Targets in Epilepsy and Beyond

Self-organized neuronal oscillations rely on precisely orchestrated ensemble activity in reverberating neuronal networks. Chronic, non-malignant disorders of the brain are often coupled to pathological neuronal activity patterns. In addition to the characteristic behavioral symptoms, these disturbances are giving rise to both transient and persistent changes of various brain rhythms. Increasing evidence support the causal role of these "oscillopathies" in the phenotypic emergence of the disease symptoms, identifying neuronal network oscillations as potential therapeutic targets. While the kinetics of pharmacological therapy is not suitable to compensate the disease related fine-scale disturbances of network oscillations, external biophysical modalities (e.g., electrical stimulation) can alter spike timing in a temporally precise manner. These perturbations can warp rhythmic oscillatory patterns via resonance or entrainment. Properly timed phasic stimuli can even switch between the stable states of networks acting as multistable oscillators, substantially changing the emergent oscillatory patterns. Novel transcranial electric stimulation (TES) approaches offer more reliable neuronal control by allowing higher intensities with tolerable side-effect profiles. This precise temporal steerability combined with the non- or minimally invasive nature of these novel TES interventions make them promising therapeutic candidates for functional disorders of the brain. Here we review the key experimental findings and theoretical background concerning various pathological aspects of neuronal network activity leading to the generation of epileptic seizures. The conceptual and practical state of the art of temporally targeted brain stimulation is discussed focusing on the prevention and early termination of epileptic seizures.

Proteins related to ictogenesis and seizure clustering in chronic epilepsy

Seizure clustering is a common phenomenon in epilepsy. Protein expression profiles during a seizure cluster might reflect the pathomechanism underlying ictogenesis. We performed proteomic analyses to identify proteins with a specific temporal expression pattern in cluster phases and to demonstrate their potential pathomechanistic role. Pilocarpine epilepsy model mice with confirmed cluster pattern of spontaneous recurrent seizures by long-term video-electroencpehalography were sacrificed at the onset, peak, or end of a seizure cluster or in the seizure-free period. Proteomic analysis was performed in the hippocampus and the cortex. Differentially expressed proteins (DEPs) were identified and classified according to their temporal expression pattern. Among the five hippocampal (HC)-DEP classes, HC-class 1 (66 DEPs) represented disrupted cell homeostasis due to clustered seizures, HC-class 2 (63 DEPs) cluster-onset downregulated processes, HC-class 3 (42 DEPs) cluster-onset upregulated processes, and HC-class 4 (103 DEPs) consequences of clustered seizures. Especially, DEPs in HC-class 3 were hippocampus-specific and involved in axonogenesis, synaptic vesicle assembly, and neuronal projection, indicating their pathomechanistic roles in ictogenesis. Key proteins in HC-class 3 were highly interconnected and abundantly involved in those biological processes. This study described the seizure cluster-associated spatiotemporal regulation of protein expression. HC-class 3 provides insights regarding ictogenesis-related processes.

Impact of Antiepileptic Drugs on Cognition and Neuromagnetic Activity in Childhood Epilepsy With Centrotemporal Spikes: A Magnetoencephalography Study

Objective: Childhood epilepsy with centrotemporal spikes (CECTS), the most common childhood epilepsy, still lacks longitudinal imaging studies involving antiepileptic drugs (AEDs). In order to examine the effect of AEDs on cognition and brain activity. We investigated the neuromagnetic activities and cognitive profile in children with CECTS before and after 1 year of treatment.

Methods: Fifteen children with CECTS aged 6–12 years underwent high-sampling magnetoencephalography (MEG) recordings before treatment and at 1 year after treatment, and 12 completed the cognitive assessment (The Wechsler Intelligence Scale for Children). Next, magnetic source location and functional connectivity (FC) were investigated in order to characterize interictal neuromagnetic activity in the seven frequency sub-bands, including: delta (1–4 Hz), theta (4–8 Hz), alpha (8–12 Hz), beta (12–30 Hz), gamma (30–80 Hz), ripple (80–250 Hz), and fast ripple (250–500 Hz).

Results: After 1 year of treatment, children with CECTS had increased scores on full-scale intelligence quotient, verbal comprehension index (VCI) and perceptual reasoning index (PRI). Alterations of neural activity occurred in specific frequency bands. Source location, in the 30–80 Hz frequency band, was significantly increased in the posterior cingulate cortex (PCC) after treatment. Moreover, FC analysis demonstrated that after treatment, the connectivity between the PCC and the medial frontal cortex (MFC) was enhanced in the 8–12 Hz frequency band. Additionally, the whole-brain network distribution was more dispersed in the 80–250 Hz frequency band.

Conclusion: Intrinsic neural activity has frequency-dependent characteristic. AEDs have impact on regional activity and FC of the default mode network (DMN). Normalization of aberrant DMN in children with CECTS after treatment is likely the reason for improvement of cognitive function.

Canard solutions in neural mass models: consequences on critical regimes

Mathematical models at multiple temporal and spatial scales can unveil the fundamental mechanisms of critical transitions in brain activities. Neural mass models (NMMs) consider the average temporal dynamics of interconnected neuronal subpopulations without explicitly representing the underlying cellular activity. The mesoscopic level offered by the neural mass formulation has been used to model electroencephalographic (EEG) recordings and to investigate various cerebral mechanisms, such as the generation of physiological and pathological brain activities. In this work, we consider a NMM widely accepted in the context of epilepsy, which includes four interacting neuronal subpopulations with different synaptic kinetics. Due to the resulting three-time-scale structure, the model yields complex oscillations of relaxation and bursting types. By applying the principles of geometric singular perturbation theory, we unveil the existence of the canard solutions and detail how they organize the complex oscillations and excitability properties of the model. In particular, we show that boundaries between pathological epileptic discharges and physiological background activity are determined by the canard solutions. Finally we report the existence of canard-mediated small-amplitude frequency-specific oscillations in simulated local field potentials for decreased inhibition conditions. Interestingly, such oscillations are actually observed in intracerebral EEG signals recorded in epileptic patients during pre-ictal periods, close to seizure onsets.

Control of epileptic seizures by electrical stimulation: a model-based study

High frequency electrical stimulation of brain is commonly used in research experiments and clinical trials as a modern tool for control of epileptic seizures. However, the mechanistic basis by which periodic external stimuli alter the brain state is not well understood. This study provides a computational insight into the mechanism of seizure suppression by high frequency stimulation (HFS). In particular, a modified version of the Jansen-Rit neural mass model is employed, in which EEG signals can be considered as the input. The proposed model reproduces seizure-like activity in the output during the ictal period of the input signal. By applying a control signal to the model, a wide range of stimulation amplitudes and frequencies are systematically explored. Simulation results reveal that HFS can effectively suppress the seizure-like activity. Our results suggest that HFS has the ability of shifting the operating state of neural populations away from a critical condition. Furthermore, a closed-loop control strategy is proposed in this paper. The main objective has been to considerably reduce the control effort needed for blocking abnormal activity of the brain. Such an energy reduction could be of practical importance, to reduce possible side effects and increase battery life for implanted neurostimulators.

Epilepsy as a dynamical system, a most needed paradigm shift in epileptology

The idea of the epileptic brain being highly excitable and facilitated to synchronic activity has guided pharmacological treatment since the early twentieth century. Although tackling epilepsy's seizure-prone feature, by tonically modifying overall circuit excitability and/or connectivity, the last 50 years of drug development has not seen a substantial improvement in seizure suppression of refractory epilepsies. This review presents a new conceptual framework for epilepsy in which the temporal dynamics of the disease plays a more critical role in both its understanding and therapeutic strategies. The repetitive epileptiform pattern (characteristic during ictal activity) and other well-defined electrographic signatures (i.e., present during the interictal period) are discussed in terms of the sequential activation of the circuit motifs. Lessons learned from the physiological activation of neural circuitry are used to further corroborate the argument and explore the transition from proper function to a state of instability. Furthermore, the review explores how interfering in the temporally dependent abnormal connectivity between circuits may work as a therapeutic approach. We also review the use of probing stimulation to access network connectivity and evaluate its power to determine transitional states of the dynamical system as it moves towards regions of instability, especially when conventional electrographic monitoring is proven inefficient. Unorthodox cases, with little or no scalp electrographic correlate, in which ictogenic circuitry and/or seizure spread is temporally restricted to neurovegetative, cognitive, and motivational areas are shown as possible explanations for sudden death in epilepsy (SUDEP) and other psychiatric comorbidities. In short, this review presents a paradigm shift in the way that we address the disease and is aimed to encourage debate rather than narrow the rationale epilepsy is currently engaged in. This article is part of the Special Issue "NEWroscience 2018".

Nonperiodic stimulation for the treatment of refractory epilepsy: Applications, mechanisms, and novel insights

Electrical stimulation of the central nervous system is a promising alternative for the treatment of pharmacoresistant epilepsy. Successful clinical and experimental stimulation is most usually carried out as continuous trains of current or voltage pulses fired at rates of 100 Hz or above, since lower frequencies yield controversial results. On the other hand, stimulation frequency should be as low as possible, in order to maximize implant safety and battery efficiency. Moreover, the development of stimulation approaches has been largely empirical in general, while they should be engineered with the neurobiology of epilepsy in mind if a more robust, efficient, efficacious, and safe application is intended. In an attempt to reconcile evidence of therapeutic effect with the understanding of the underpinnings of epilepsy, our group has developed a nonstandard form of low-frequency stimulation with randomized interpulse intervals termed nonperiodic stimulation (NPS). The rationale was that an irregular temporal pattern would impair neural hypersynchronization, which is a hallmark of epilepsy. In this review, we start by briefly revisiting the literature on the molecular, cellular, and network level mechanisms of epileptic phenomena in order to highlight this often-overlooked emergent property of cardinal importance in the pathophysiology of the disease. We then review our own studies on the efficacy of NPS against acute and chronic experimental seizures and also on the anatomical and physiological mechanism of the method, paying special attention to the hypothesis that the lack of temporal regularity induces desynchronization. We also put forward a novel insight regarding the temporal structure of NPS that may better encompass the set of findings published by the group: the fact that intervals between stimulation pulses have a distribution that follows a power law and thus may induce natural-like activity that would compete with epileptiform discharge for the recruitment of networks. We end our discussion by mentioning ongoing research and future projects of our lab.

The role of interictal discharges in ictogenesis - A dynamical perspective

Interictal epileptiform discharge (IED) is a traditional hallmark of epileptic tissue that is generated by the synchronous activity of a population of neurons. Interictal epileptiform discharges represent a heterogeneous group of pathological activities that differ in shape, duration, spatiotemporal distribution, underlying cellular and network mechanisms, and their relationship to seizure genesis. The exact role of IEDs in epilepsy is still not well understood, and there remains a persistent dichotomy about the impact on IEDs on seizures. Proseizure, antiseizure, and no impact on ictogenesis have all been described in previous studies. In this article, we review the existing knowledge on the role of interictal discharges in seizure genesis, and we discuss how dynamical approaches to ictogenesis can explain the existing dichotomy about the multifaceted role of IEDs in ictogenesis. This article is part of the Special Issue "NEWroscience 2018".

Measuring the effects of sleep on epileptogenicity with multifrequency entropy

OBJECTIVE

We demonstrate that multifrequency entropy gives insight into the relationship between epileptogenicity and sleep, and forms the basis for an improved measure of medical assessment of sleep impairment in epilepsy patients.

METHODS

Multifrequency entropy was computed from electroencephalography measurements taken from 31 children with Benign Epilepsy with Centrotemporal Spikes and 31 non-epileptic controls while awake and during sleep. Values were compared in the epileptic zone and away from the epileptic zone in various sleep stages.

RESULTS

We find that (I) in lower frequencies, multifrequency entropy decreases during non-rapid eye movement sleep stages when compared with wakefulness in a general population of pediatric patients, (II) patients with Benign Epilepsy with Centrotemporal Spikes had lower multifrequency entropy across stages of sleep and wakefulness, and (III) the epileptic regions of the brain exhibit lower multifrequency entropy patterns than the rest of the brain in epilepsy patients.

CONCLUSIONS

Our results show that multifrequency entropy decreases during sleep, particularly sleep stage 2, confirming, in a pediatric population, an association between sleep, lower multifrequency entropy, and increased likelihood of seizure.

SIGNIFICANCE

We observed a correlation between lowered multifrequency entropy and increased epileptogenicity that lays preliminary groundwork for the detection of a digital biomarker for epileptogenicity.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Exploring the Inhibitory Effect of Low-frequency Magnetic Fields on Epileptiform Discharges in Juvenile Rat Hippocampus

Stimulation with a low frequency electromagnetic field (LF-EMF) has proven to represent a powerful method for the suppression of seizures, as demonstrated in select clinical and laboratory studies. However, the mechanism by which LF-EMF suppresses seizures remains unclear. The purpose of the present study was to explore the modulatory effect of LF-EMF on epileptiform discharges (EDs) using rat hippocampal slices and investigate the underlying mechanisms that mediate these effects. EDs in hippocampal slices was induced by magnesium-free (zero-Mg²⁺) artificial cerebrospinal fluid (ACSF) and recorded using an in vitro micro-electrode array (MEA). A small sub-decimeter coil was designed and incorporated in a flexible magnetic stimulation device that allowed electromagnetic fields with different parameters to be delivered to slices. After a stable ED event was recorded, magnetic fields of 0.5 Hz (30 min) with a magnetic intensity of 0.13 mT (5 V_pp voltage input) and 0.25 mT (20 V_pp voltage input) were applied. The results indicated that a high-amplitude 0.5 Hz magnetic field could lead to persistent suppression of ictal discharges (IDs), while low-amplitude magnetic fields did not influence IDs. The persistent suppression of complex ED was prevented if the magnetic fields were applied in the presence of 10 μmol/L bicuculline (BIC), a γ-aminobutyric acid type A (GABA_A) receptor antagonist, while the application of BIC subsequent to a magnetic field application led to the reappearance of ID. The addition of BIC resulted in EDs that had previously been inhibited by magnetic fields, reappearing. Low-frequency magnetic stimulation was able to inhibit the conversion from interictal discharges (IIDs) or preictal discharges (PIDs) to IDs. This suppression was attributed to the modulation of GABA_A receptor activity.

Prediction of Seizure Recurrence. A Note of Caution

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

Cross-Scale Causality and Information Transfer in Simulated Epileptic Seizures

An information-theoretic approach for detecting causality and information transfer was applied to phases and amplitudes of oscillatory components related to different time scales and obtained using the wavelet transform from a time series generated by the Epileptor model. Three main time scales and their causal interactions were identified in the simulated epileptic seizures, in agreement with the interactions of the model variables. An approach consisting of wavelet transform, conditional mutual information estimation, and surrogate data testing applied to a single time series generated by the model was demonstrated to be successful in the identification of all directional (causal) interactions between the three different time scales described in the model. Thus, the methodology was prepared for the identification of causal cross-frequency phase–phase and phase–amplitude interactions in experimental and clinical neural data.

Machine learning for detection of interictal epileptiform discharges

The electroencephalogram (EEG) is a fundamental tool in the diagnosis and classification of epilepsy. In particular, Interictal Epileptiform Discharges (IEDs) reflect an increased likelihood of seizures and are routinely assessed by visual analysis of the EEG. Visual assessment is, however, time consuming and prone to subjectivity, leading to a high misdiagnosis rate and motivating the development of automated approaches. Research towards automating IED detection started 45 years ago. Approaches range from mimetic methods to deep learning techniques. We review different approaches to IED detection, discussing their performance and limitations. Traditional machine learning and deep learning methods have yielded the best results so far and their application in the field is still growing. Standardization of datasets and outcome measures is necessary to compare models more objectively and decide which should be implemented in a clinical setting.

Computational Models in Electroencephalography

Computational models lie at the intersection of basic neuroscience and healthcare applications because they allow researchers to test hypotheses in silico and predict the outcome of experiments and interactions that are very hard to test in reality. Yet, what is meant by “computational model” is understood in many different ways by researchers in different fields of neuroscience and psychology, hindering communication and collaboration. In this review, we point out the state of the art of computational modeling in Electroencephalography (EEG) and outline how these models can be used to integrate findings from electrophysiology, network-level models, and behavior. On the one hand, computational models serve to investigate the mechanisms that generate brain activity, for example measured with EEG, such as the transient emergence of oscillations at different frequency bands and/or with different spatial topographies. On the other hand, computational models serve to design experiments and test hypotheses in silico. The final purpose of computational models of EEG is to obtain a comprehensive understanding of the mechanisms that underlie the EEG signal. This is crucial for an accurate interpretation of EEG measurements that may ultimately serve in the development of novel clinical applications.

Perturbations both trigger and delay seizures due to generic properties of slow-fast relaxation oscillators

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.

Bad Timing for Epileptic Networks: Role of Temporal Dynamics in Seizures and Cognitive Deficits

The precise coordination of neuronal activity is critical for optimal brain function. When such coordination fails, this can lead to dire consequences. In this review, I will present evidence that in epilepsy, failed coordination leads not only to seizures but also to alterations of the rhythmical patterns observed in the electroencephalogram and cognitive deficits. Restoring the dynamic coordination of epileptic networks could therefore both improve seizures and cognitive outcomes.

Assessing criticality in pre-seizure single-neuron activity of human epileptic cortex

Epileptic seizures are characterized by abnormal and excessive neural activity, where cortical network dynamics seem to become unstable. However, most of the time, during seizure-free periods, cortex of epilepsy patients shows perfectly stable dynamics. This raises the question of how recurring instability can arise in the light of this stable default state. In this work, we examine two potential scenarios of seizure generation: (i) epileptic cortical areas might generally operate closer to instability, which would make epilepsy patients generally more susceptible to seizures, or (ii) epileptic cortical areas might drift systematically towards instability before seizure onset. We analyzed single-unit spike recordings from both the epileptogenic (focal) and the nonfocal cortical hemispheres of 20 epilepsy patients. We quantified the distance to instability in the framework of criticality, using a novel estimator, which enables an unbiased inference from a small set of recorded neurons. Surprisingly, we found no evidence for either scenario: Neither did focal areas generally operate closer to instability, nor were seizures preceded by a drift towards instability. In fact, our results from both pre-seizure and seizure-free intervals suggest that despite epilepsy, human cortex operates in the stable, slightly subcritical regime, just like cortex of other healthy mammalians.

Autonomic Characteristics of Sudden Unexpected Death in Epilepsy in Children-A Systematic Review of Studies and Their Relevance to the Management of Epilepsy in Rett Syndrome

Aim: To systematically identify and critically appraise studies that investigate the autonomic characteristics of Sudden Unexpected Death in Epilepsy (SUDEP) in the pediatric population. We also wanted to explore how this information would be relevant to the management of epilepsy in patients with Rett Syndrome.

Method: Using PRISMA guidelines, a systematic review of PubMed, Scopus, Cochrane, PsycINFO, Embase, and Web of Science databases was performed to identify eligible studies. After extracting data from the included studies, a thematic analysis was undertaken to identify emerging themes. A quality appraisal was also done to assess the quality of the included studies.

Results: The systematic search revealed 41 records, and 15 full-text articles on the autonomic characteristics of SUDEP in children were included in the final analysis. Following thematic analysis, three themes were identified (I) modulation in sympathovagal tone, (II) pre- and post-ictal autonomic changes, and (III) other markers of autonomic dysregulation in children with epilepsy. Modulation in sympathovagal tone emerged as the theme with the highest frequency followed by pre- and post-ictal autonomic changes. While the themes provide additional insight into the management of epilepsy in the Rett Syndrome population, the quality of evidence concerning the autonomic characteristics of SUDEP in the pediatric population was low and underscores the importance of much needed research in this area.

Conclusion: The mechanism of SUDEP in the pediatric population is complex and involves an interplay between several components of the autonomic nervous system. While direct clinical inferences regarding pediatric SUDEP could not be made, the thematic analysis does suggest that in vulnerable populations such as Rett Syndrome, where there is already a pervasive autonomic dysregulation, pro-active surveillance of the autonomic profile in this patient group would be useful to better manage epilepsy and reduce the SUDEP risk.

Falling off a limit cycle using phase-agnostic stimuli: Applications to clinical oscillopathies

For over a century, physiological studies have shown that precisely timed pulses can switch off a biological oscillator. This empiric finding has shaped our mechanistic understanding of how perturbations start, stop, and reset biological oscillators and has led to treatments that suppress pathological oscillations using electrical pulses given within specified therapeutic phase windows. Here, we present evidence, using numerical simulations of models of epileptic seizures and reentrant tachycardia, that the phase window can be opened to the entire cycle using novel complex stimulus waveforms. Our results reveal that the trajectories are displaced by such phase-agnostic stimuli off the oscillator's limit cycle and corralled into a region where oscillation is suppressed, irrespective of the phase at which the stimulus was applied. Our findings suggest the need for broadening theoretical understanding of how complex perturbing waveforms interact with biological oscillators to access their arrhythmic states. In clinical practice, oscillopathies may be treated more effectively with non-traditional stimulus waveforms that obviate the need for phase specificity.

Chance and risk in epilepsy

PURPOSE OF REVIEW

Epilepsy is a dynamical disorder of the brain characterized by sudden, seemingly unpredictable transitions to the ictal state. When and how these transitions occur remain unresolved questions in neurology.

RECENT FINDINGS

Modelling work based on dynamical systems theory proposed that a slow control parameter is necessary to explain the transition between interictal and ictal states. Recently, converging evidence from chronic EEG datasets unravelled the existence of cycles of epileptic brain activity at multiple timescales - circadian, multidien (over multiple days) and circannual - which could reflect cyclical changes in a slow control parameter. This temporal structure of epilepsy has theoretical implications and argues against the conception of seizures as completely random events. The practical significance of cycles in epilepsy is highlighted by their predictive value in computational models for seizure forecasting.

SUMMARY

The canonical randomness of seizures is being reconsidered in light of cycles of brain activity discovered through chronic EEG. This paradigm shift motivates development of next-generation devices to track more closely fluctuations in epileptic brain activity that determine time-varying seizure risk.

Alterations of Neuronal Dynamics as a Mechanism for Cognitive Impairment in Epilepsy

Epilepsy is commonly associated with cognitive and behavioral deficits that dramatically affect the quality of life of patients. In order to identify novel therapeutic strategies aimed at reducing these deficits, it is critical first to understand the mechanisms leading to cognitive impairments in epilepsy. Traditionally, seizures and epileptiform activity in addition to neuronal injury have been considered to be the most significant contributors to cognitive dysfunction. In this review we however highlight the role of a new mechanism: alterations of neuronal dynamics, i.e. the timing at which neurons and networks receive and process neural information. These alterations, caused by the underlying etiologies of epilepsy syndromes, are observed in both animal models and patients in the form of abnormal oscillation patterns in unit firing, local field potentials, and electroencephalogram (EEG). Evidence suggests that such mechanisms significantly contribute to cognitive impairment in epilepsy, independently of seizures and interictal epileptiform activity. Therefore, therapeutic strategies directly targeting neuronal dynamics rather than seizure reduction may significantly benefit the quality of life of patients.

Sharp decrease in the Laplacian matrix rank of phase-space graphs: a potential biomarker in epilepsy

In this paper, phase space reconstruction from stereo-electroencephalography data of ten patients with focal epilepsy forms a series of graphs. Those obtained graphs reflect the transition characteristics of brain dynamical system from pre-seizure to seizure of epilepsy. Interestingly, it is found that the rank of Laplacian matrix of these graphs has a sharp decrease when a seizure is close to happen, which thus might be viewed as a new potential biomarker in epilepsy. In addition, the reliability of this method is numerically verified with a coupled mass neural model. In particular, our simulation suggests that this potential biomarker can play the roles of predictive effect or delayed awareness, depending on the bias current of the Gaussian noise. These results may give new insights into the seizure detection.

FitzHugh-Nagumo oscillators on complex networks mimic epileptic-seizure-related synchronization phenomena

We study patterns of partial synchronization in a network of FitzHugh-Nagumo oscillators with empirical structural connectivity measured in human subjects. We report the spontaneous occurrence of synchronization phenomena that closely resemble the ones seen during epileptic seizures in humans. In order to obtain deeper insights into the interplay between dynamics and network topology, we perform long-term simulations of oscillatory dynamics on different paradigmatic network structures: random networks, regular nonlocally coupled ring networks, ring networks with fractal connectivities, and small-world networks with various rewiring probability. Among these networks, a small-world network with intermediate rewiring probability best mimics the findings achieved with the simulations using the empirical structural connectivity. For the other network topologies, either no spontaneously occurring epileptic-seizure-related synchronization phenomena can be observed in the simulated dynamics, or the overall degree of synchronization remains high throughout the simulation. This indicates that a topology with some balance between regularity and randomness favors the self-initiation and self-termination of episodes of seizure-like strong synchronization.