Early detection of epileptic seizures based on parameter identification of neural mass model

Computational Evidence for a Competitive Thalamocortical Model of Spikes and Spindle Activity in Rolandic Epilepsy

Rolandic epilepsy (RE) is the most common idiopathic focal childhood epilepsy syndrome, characterized by sleep-activated epileptiform spikes and seizures and cognitive deficits in school age children. Recent evidence suggests that this disease may be caused by disruptions to the Rolandic thalamocortical circuit, resulting in both an abundance of epileptiform spikes and a paucity of sleep spindles in the Rolandic cortex during non-rapid eye movement sleep (NREM); electrographic features linked to seizures and cognitive symptoms, respectively. The neuronal mechanisms that support the competitive shared thalamocortical circuitry between pathological epileptiform spikes and physiological sleep spindles are not well-understood. In this study we introduce a computational thalamocortical model for the sleep-activated epileptiform spikes observed in RE. The cellular and neuronal circuits of this model incorporate recent experimental observations in RE, and replicate the electrophysiological features of RE. Using this model, we demonstrate that: (1) epileptiform spikes can be triggered and promoted by either a reduced NMDA current or h-type current; and (2) changes in inhibitory transmission in the thalamic reticular nucleus mediates an antagonistic dynamic between epileptiform spikes and spindles. This work provides the first computational model that both recapitulates electrophysiological features and provides a mechanistic explanation for the thalamocortical switch between the pathological and physiological electrophysiological rhythms observed during NREM sleep in this common epileptic encephalopathy.

A New Neural Mass Model Driven Method and Its Application in Early Epileptic Seizure Detection

OBJECTIVE

Despite numerous neural computational models proposed to explain physiological and pathological mechanisms of brain activity, a large gap remains between theory and application of the models. Building on the successful application of data-driven methods in epileptic seizure detection, we aim to build a bridge between data and models in this paper.

METHODS

We first propose a novel model-driven seizure detection method based on dynamic features in epileptic EEGs, where the rationale for dynamic features in epileptic EEGs can be clarified in theory by characterizing the variation of parameters of the model. Then we apply the proposed D&F-model-driven method to the problem of early epileptic seizure detection, where the evolution of model parameters selected and optimized by the proposed method is measured and used to detect the starting point of the seizure.

RESULTS

Numerical results on two open EEG databases demonstrate that our proposed method does a good job of early epileptic seizure detection. The average detection sensitivity, false positive rate and early detection period attain 100%, 0.1/h, and 7.1 s respectively.

CONCLUSION

This paper provides a strategy to characterize EEG signals using a NMM-related method and the model parameters optimized by real EEG may then serve as features in their own right for early seizure detection.

SIGNIFICANCE

An useful attempt to early detect epileptic seizures by combining the neural mass model with data analysis.

Seizure tracking of epileptic EEGs using a model-driven approach

OBJECTIVE

As a chronic neurological disorder, epilepsy is characterized by recurrent and unprovoked epileptic seizures that can disrupt the normal neuro-biologic, cognitive, psychological conditions of patients. Therefore, it is worthwhile to give a detailed account of how the epileptic EEG evolves during a period of seizure so that an effective control can be guided for epileptic patients in clinics.

APPROACH

Considering the successful application of the neural mass model (NMM) in exploring the insights into brain activities for epilepsy, in this paper, we aim to construct a model-driven approach to track the development of seizures using epileptic EEGs. We first propose a new time-delay Wendling model with sub-populations (TD-W-SP model) with respect to three aspects of improvements. Then we introduce a model-driven seizure tracking approach, where a model training method is designed based on extracted features from epileptic EEGs and a tracking index is defined as a function of the trained model parameters.

MAIN RESULTS

Numerical results on eight patients on CHB-MIT database demonstrate that our proposed method performs well in simulating epileptic-like EEGs as well as tracking the evolution of three stages (that is, from pre-ictal to ictal and from ictal to post-ictal) during a period of epileptic seizure.

SIGNIFICANCE

A useful attempt to track epileptic seizures by combining the NMM with the data analysis.

Automated epileptic seizure detection based on break of excitation/inhibition balance

Physiological models are attractive for seizure detection, as their parameters are related to physiological meanings. We propose an algorithm to early detect epileptic seizures based on automatic estimation of average synaptic gains (excitatory Ae, slow and fast inhibitory B and G) by combining clinical data with a neural mass model. Three indices (Ae/B, Ae/G and Ae/(B + G)), all related to excitation/inhibition balance, were calculated and used as cues to detect seizures. A simple thresholding method was employed. We evaluated the algorithm against the manual scoring of a human expert on intracranial EEG samples from 23 patients suffering from different types of epilepsy. Best performance was achieved using Ae/(B + G) as a cue, i.e. excitation/(slow + fast) inhibition, on temporal lobe epilepsy (TLE) patients. A leave-one-out cross-validation showed that the algorithm achieved 92.98% sensitivity for TLE patients. The median false positive rate was 0.16 per hour, and median detection delay was 14.5 s. Of interest, the threshold values determined by a leave-one-out cross-validation did nearly not vary among TLE patients, suggesting a general excitation/inhibition balance baseline in TLE patients. The same approach could be used with other types of epilepsy by adapting the neural mass model to these types.

Seizure evolution can be characterized as path through synaptic gain space of a neural mass model

Physiologically based models could facilitate better understanding of mechanisms underlying epileptic seizures. In this paper, we attempt to reveal the dynamic evolution of intracranial EEG activity during epileptic seizures based on synaptic gain identification procedure of a neural mass model. The distribution of average excitatory, slow and fast inhibitory synaptic gain in the parameter space and their temporal evolution, i.e., the path through the model parameter space, were analyzed in thirty seizures from ten temporal lobe epileptic patients. Results showed that the synaptic gain values located roughly on a plane before seizure onset, dispersed during seizure and returned to the plane when seizure terminated. Cluster analysis was performed on seizure paths and demonstrated consistency in synaptic gain evolution across different seizures from the individual patient. Furthermore, two patient groups were identified, each one corresponding to a specific synaptic gain evolution in the parameter space during a seizure. Results were validated by a bootstrapping approach based on comparison with random paths. The differences in the path revealed variations in EEG dynamics for patients despite showing identical seizure onset pattern. Our approach may have the potential to classify the epileptic patients into subgroups based on different mechanisms revealed by subtle changes in synaptic gains and further enable more robust decisions regarding treatment strategy.

Dynamics underlying interictal to ictal transition in temporal lobe epilepsy: insights from a neural mass model

We propose an approach that combines a neural mass model and clinical intracranial electroencephalographic (iEEG) recordings to explore the potential pathophysiological mechanisms (at the neuronal population level) of ictogenesis. Thirty iEEG recordings from 10 temporal lobe epilepsy (TLE) patients around seizure onset were investigated. Physiologically meaningful parameters [average excitatory (A_e ), slow (B), and fast (G) inhibitory synaptic gain] were identified during interictal to ictal transition. Four ratios (A_e /G, A_e /B, A_e /(B + G), and B/G) were derived from these parameters, and their evolution over time was analyzed. The excitation/inhibition ratio increased around seizure onset and decreased before seizure offset, indicating the impairment and re-emergence of excitation/inhibition balance around seizure onset and before seizure offset, respectively. Moreover, the slow inhibition may have an earlier effect on excitation/inhibition imbalance. We confirm the decrease in excitation/inhibition ratio upon seizure termination in human temporal lobe epilepsy, as revealed by optogenetic approaches both in vivo in animal models and in vitro. The increase in excitation/inhibition ratio around seizure occurrence could be an indicator to detect seizures.

Particle swarm optimization algorithm based parameters estimation and control of epileptiform spikes in a neural mass model

This paper proposes an epilepsy detection and closed-loop control strategy based on Particle Swarm Optimization (PSO) algorithm. The proposed strategy can effectively suppress the epileptic spikes in neural mass models, where the epileptiform spikes are recognized as the biomarkers of transitions from the normal (interictal) activity to the seizure (ictal) activity. In addition, the PSO algorithm shows capabilities of accurate estimation for the time evolution of key model parameters and practical detection for all the epileptic spikes. The estimation effects of unmeasurable parameters are improved significantly compared with unscented Kalman filter. When the estimated excitatory-inhibitory ratio exceeds a threshold value, the epileptiform spikes can be inhibited immediately by adopting the proportion-integration controller. Besides, numerical simulations are carried out to illustrate the effectiveness of the proposed method as well as the potential value for the model-based early seizure detection and closed-loop control treatment design.

Tracking slow modulations in synaptic gain using dynamic causal modelling: validation in epilepsy

In this work we propose a proof of principle that dynamic causal modelling can identify plausible mechanisms at the synaptic level underlying brain state changes over a timescale of seconds. As a benchmark example for validation we used intracranial electroencephalographic signals in a human subject. These data were used to infer the (effective connectivity) architecture of synaptic connections among neural populations assumed to generate seizure activity. Dynamic causal modelling allowed us to quantify empirical changes in spectral activity in terms of a trajectory in parameter space — identifying key synaptic parameters or connections that cause observed signals. Using recordings from three seizures in one patient, we considered a network of two sources (within and just outside the putative ictal zone). Bayesian model selection was used to identify the intrinsic (within-source) and extrinsic (between-source) connectivity. Having established the underlying architecture, we were able to track the evolution of key connectivity parameters (e.g., inhibitory connections to superficial pyramidal cells) and test specific hypotheses about the synaptic mechanisms involved in ictogenesis. Our key finding was that intrinsic synaptic changes were sufficient to explain seizure onset, where these changes showed dissociable time courses over several seconds. Crucially, these changes spoke to an increase in the sensitivity of principal cells to intrinsic inhibitory afferents and a transient loss of excitatory–inhibitory balance.

•

We propose a framework to characterise slow dynamical changes in the brain.

•

Dynamical causal modelling finds the most likely connectivity among two brain areas.

•

The synaptic weights defining these connections are tracked in time.

•

We analyse brain activity of an epileptic subject, at the focus and just outside it.

•

We point to modulations of synaptic connections as responsible of the seizure.

On the modelling of seizure dynamics

This scientific commentary refers to ‘On the nature of seizure dynamics’, by V. Jirsa et al. (doi:10.1093/brain/awu133).

Functional brain connectivity from EEG in epilepsy: seizure prediction and epileptogenic focus localization

Today, neuroimaging techniques are frequently used to investigate the integration of functionally specialized brain regions in a network. Functional connectivity, which quantifies the statistical dependencies among the dynamics of simultaneously recorded signals, allows to infer the dynamical interactions of segregated brain regions. In this review we discuss how the functional connectivity patterns obtained from intracranial and scalp electroencephalographic (EEG) recordings reveal information about the dynamics of the epileptic brain and can be used to predict upcoming seizures and to localize the seizure onset zone. The added value of extracting information that is not visibly identifiable in the EEG data using functional connectivity analysis is stressed. Despite the fact that many studies have showed promising results, we must conclude that functional connectivity analysis has not made its way into clinical practice yet.