Towards a large-scale model of patient-specific epileptic spike-wave discharges

Clinical electroencephalographic (EEG) recordings of the transition into generalised epileptic seizures show a sudden onset of spike-wave dynamics from a low-amplitude irregular background. In addition, non-trivial and variable spatio-temporal dynamics are widely reported in combined EEG/fMRI studies on the scale of the whole cortex. It is unknown whether these characteristics can be accounted for in a large-scale mathematical model with fixed heterogeneous long-range connectivities. Here, we develop a modelling framework with which to investigate such EEG features. We show that a neural field model composed of a few coupled compartments can serve as a low-dimensional prototype for the transition between irregular background dynamics and spike-wave activity. This prototype then serves as a node in a large-scale network with long-range connectivities derived from human diffusion-tensor imaging data. We examine multivariate properties in 42 clinical EEG seizure recordings from 10 patients diagnosed with typical absence epilepsy and 50 simulated seizures from the large-scale model using 10 DTI connectivity sets from humans. The model can reproduce the clinical feature of stereotypy where seizures are more similar within a patient than between patients, essentially creating a patient-specific fingerprint. We propose the approach as a feasible technique for the investigation of patient-specific large-scale epileptic features in space and time.     

A large-scale brain network mechanism for increased seizure propensity in Alzheimer’s disease

People with Alzheimer’s disease (AD) are 6-10 times more likely to develop seizures than the healthy aging population. Leading hypotheses largely consider hyperexcitability of local cortical tissue as primarily responsible for increased seizure prevalence in AD. However, in the general population of people with epilepsy, large-scale brain network organization additionally plays a role in determining seizure likelihood and phenotype. Here, we propose that alterations to large-scale brain network organization seen in AD may contribute to increased seizure likelihood. To test this hypothesis, we combine computational modelling with electrophysiological data using an approach that has proved informative in clinical epilepsy cohorts without AD. EEG was recorded from 21 people with probable AD and 26 healthy controls. At the time of EEG acquisition, all participants were free from seizures. Whole brain functional connectivity derived from source-reconstructed EEG recordings was used to build subject-specific brain network models of seizure transitions. As cortical tissue excitability was increased in the simulations, AD simulations were more likely to transition into seizures than simulations from healthy controls, suggesting an increased group-level probability of developing seizures at a future time for AD participants. We subsequently used the model to assess seizure propensity of different regions across the cortex. We found the most important regions for seizure generation were those typically burdened by amyloid-beta at the early stages of AD, as previously reported by in-vivo and post-mortem staging of amyloid plaques. Analysis of these spatial distributions also give potential insight into mechanisms of increased susceptibility to generalized (as opposed to focal) seizures in AD vs controls. This research suggests avenues for future studies testing patients with seizures, e.g. co-morbid AD/epilepsy patients, and comparisons with PET and MRI scans to relate regional seizure propensity with AD pathologies.     

Epileptic seizures may begin hours in advance of clinical onset: a report of five patients

Mechanisms underlying seizure generation are traditionally thought to act over seconds to minutes before clinical seizure onset. We analyzed continuous 3- to 14-day intracranial EEG recordings from five patients with mesial temporal lobe epilepsy obtained during evaluation for epilepsy surgery. We found localized quantitative EEG changes identifying prolonged bursts of complex epileptiform discharges that became more prevalent 7 hr before seizures and highly localized subclinical seizure-like activity that became more frequent 2 hr prior to seizure onset. Accumulated energy increased in the 50 min before seizure onset, compared to baseline. These observations, from a small number of patients, suggest that epileptic seizures may begin as a cascade of electrophysiological events that evolve over hours and that quantitative measures of preseizure electrical activity could possibly be used to predict seizures far in advance of clinical onset.     

Inflection, canards and excitability threshold in neuronal models

A technique is presented, based on the differential geometry of planar curves, to evaluate the excitability threshold of neuronal models. The aim is to determine regions of the phase plane where solutions to the model equations have zero local curvature, thereby defining a zero-curvature (inflection) set that discerns between sub-threshold and spiking electrical activity. This transition can arise through a Hopf bifurcation, via the so-called canard explosion that happens in an exponentially small parameter variation, and this is typical for a large class of planar neuronal models (FitzHugh–Nagumo, reduced Hodgkin–Huxley), namely, type II neurons (resonators). This transition can also correspond to the crossing of the stable manifold of a saddle equilibrium, in the case of type I neurons (integrators). We compute inflection sets and study how well they approximate the excitability threshold of these neuron models, that is, both in the canard and in the non-canard regime, using tools from invariant manifold theory and singularity theory. With the latter, we investigate the topological changes that inflection sets undergo upon parameter variation. Finally, we show that the concept of inflection set gives a good approximation of the threshold in both the so-called resonator and integrator neuronal cases.     

Multiway analysis of epilepsy tensors

MOTIVATION: The success or failure of an epilepsy surgery depends greatly on the localization of epileptic focus (origin of a seizure). We address the problem of identification of a seizure origin through an analysis of ictal electroencephalogram (EEG), which is proven to be an effective standard in epileptic focus localization. SUMMARY: With a goal of developing an automated and robust way of visual analysis of large amounts of EEG data, we propose a novel approach based on multiway models to study epilepsy seizure structure. Our contributions are 3-fold. First, we construct an Epilepsy Tensor with three modes, i.e. time samples, scales and electrodes, through wavelet analysis of multi-channel ictal EEG. Second, we demonstrate that multiway analysis techniques, in particular parallel factor analysis (PARAFAC), provide promising results in modeling the complex structure of an epilepsy seizure, localizing a seizure origin and extracting artifacts. Third, we introduce an approach for removing artifacts using multilinear subspace analysis and discuss its merits and drawbacks. RESULTS: Ictal EEG analysis of 10 seizures from 7 patients are included in this study. Our results for 8 seizures match with clinical observations in terms of seizure origin and extracted artifacts. On the other hand, for 2 of the seizures, seizure localization is not achieved using an initial trial of PARAFAC modeling. In these cases, first, we apply an artifact removal method and subsequently apply the PARAFAC model on the epilepsy tensor from which potential artifacts have been removed. This method successfully identifies the seizure origin in both cases.     

A method for detecting false bifurcations in dynamical systems: application to neural-field models

In this article, we present a method for tracking changes in curvature of limit cycle solutions that arise due to inflection points. In keeping with previous literature, we term these changes false bifurcations, as they appear to be bifurcations when considering a Poincaré section that is tangent to the solution, but in actual fact the deformation of the solution occurs smoothly as a parameter is varied. These types of solutions arise commonly in electroencephalogram models of absence seizures and correspond to the formation of spikes in these models. Tracking these transitions in parameter space allows regions to be defined corresponding to different types of spike and wave dynamics, that may be of use in clinical neuroscience as a means to classify different subtypes of the more general syndrome.     

Increase trend of correlation and phase synchrony of microwire iEEG before macroseizure onset

Micro/macrowire intracranial EEG (iEEG) signals recorded from implanted micro/macroelectrodes in epileptic patients have received great attention and are considered to include much information of neuron activities in seizure transition compared to scalp EEG from cortical electrodes. Microelectrode is contacted more close to neurons than macroelectrode and it is more sensitive to neuron activity changes than macroelectrode. Microwire iEEG recordings are inevitably advantageous over macrowire iEEG recordings to reveal neuronal mechanisms contributing to the generation of seizures. In this study, we investigate the seizure generation from microwire iEEG recordings and discuss synchronization of microwire iEEGs in four frequency bands: alpha (1−30 Hz), gamma (30−80 Hz), ripple (80–250 Hz), and fast ripple (>250 Hz) via two measures: correlation and phase synchrony. We find that an increase trend of correlation or phase synchrony exists before the macroseizure onset mostly in gamma and ripple bands where the duration of the preictal states varied in different seizures ranging up to a few seconds (minutes). This finding is contrast to the well-known result that a decrease of synchronization in macro domains exists before the macroseizure onset. The finding demonstrates that it is only when the seizure has recruited enough surrounding brain tissue does the signal become strong enough to be observed on the clinical macroelectrode and as a result support the hypothesis of progressive coalescence of microseizure domains. The potential ramifications of such an early detection of microscale seizure activity may open a new window on treatment by making possible disruption of seizure activity before it becomes fully established.     

Widespread EEG Changes Precede Focal Seizures

The process by which the brain transitions into an epileptic seizure is unknown. In this study, we investigated whether the transition to seizure is associated with changes in brain dynamics detectable in the wideband EEG, and whether differences exist across underlying pathologies. Depth electrode ictal EEG recordings from 40 consecutive patients with pharmacoresistant lesional focal epilepsy were low-pass filtered at 500 Hz and sampled at 2,000 Hz. Predefined EEG sections were selected immediately before (immediate preictal), and 30 seconds before the earliest EEG sign suggestive of seizure activity (baseline). Spectral analysis, visual inspection and discrete wavelet transform were used to detect standard (delta, theta, alpha, beta and gamma) and high-frequency bands (ripples and fast ripples). At the group level, each EEG frequency band activity increased significantly from baseline to the immediate preictal section, mostly in a progressive manner and independently of any modification in the state of vigilance. Preictal increases in each frequency band activity were widespread, being observed in the seizure-onset zone and lesional tissue, as well as in remote regions. These changes occurred in all the investigated pathologies (mesial temporal atrophy/sclerosis, local/regional cortical atrophy, and malformations of cortical development), but were more pronounced in mesial temporal atrophy/sclerosis. Our findings indicate that a brain state change with distinctive features, in the form of unidirectional changes across the entire EEG bandwidth, occurs immediately prior to seizure onset. We postulate that these changes might reflect a facilitating state of the brain which enables a susceptible region to generate seizures.     

Persistence of Cortical Sensory Processing during Absence Seizures in Human and an Animal Model: Evidence from EEG and Intracellular Recordings

Absence seizures are caused by brief periods of abnormal synchronized oscillations in the thalamocortical loops, resulting in widespread spike-and-wave discharges (SWDs) in the electroencephalogram (EEG). SWDs are concomitant with a complete or partial impairment of consciousness, notably expressed by an interruption of ongoing behaviour together with a lack of conscious perception of external stimuli. It is largely considered that the paroxysmal synchronizations during the epileptic episode transiently render the thalamocortical system incapable of transmitting primary sensory information to the cortex. Here, we examined in young patients and in the Genetic Absence Epilepsy Rats from Strasbourg (GAERS), a well-established genetic model of absence epilepsy, how sensory inputs are processed in the related cortical areas during SWDs. In epileptic patients, visual event-related potentials (ERPs) were still present in the occipital EEG when the stimuli were delivered during seizures, with a significant increase in amplitude compared to interictal periods and a decrease in latency compared to that measured from non-epileptic subjects. Using simultaneous in vivo EEG and intracellular recordings from the primary somatosensory cortex of GAERS and non-epileptic rats, we found that ERPs and firing responses of related pyramidal neurons to whisker deflection were not significantly modified during SWDs. However, the intracellular subthreshold synaptic responses in somatosensory cortical neurons during seizures had larger amplitude compared to quiescent situations. These convergent findings from human patients and a rodent genetic model show the persistence of cortical responses to sensory stimulations during SWDs, indicating that the brain can still process external stimuli during absence seizures. They also demonstrate that the disruption of conscious perception during absences is not due to an obliteration of information transfer in the thalamocortical system. The possible mechanisms rendering the cortical operation ineffective for conscious perception are discussed, but their definite elucidation will require further investigations.     

Non-linear analysis of the electroencephalogram in Creutzfeldt-Jakob disease

Creutzfeldt-Jakob disease is a rare, neurological, dementing disorder characterised by periodic sharp waves in the electroencephalogram (EEG). Non-linear analysis of these EEG changes may provide insight into the abnormal dynamics of cortical neural networks in this disorder. Babloyantz et al. have suggested that the periodic sharp waves reflect low-dimensional chaotic dynamics in the brain. In the present study this hypothesis was re-examined using newly developed techniques for non-linear time series analysis. We analysed the EEG of a patient with autopsy-proven Creutzfeldt-Jakob disease using the method of non-linear forecasting as introduced by Sugihara and May, and we tested for non-linearity with amplitude-adjusted, phase-randomised surrogate data. Two epochs with generalised periodic sharp waves showed clear evidence for non-linearity. These epochs could be predicted better and further ahead in time than most of the irregular background activity. Testing against cycle-randomised surrogate data and close inspection of the periodograms showed that the non-linearity of the periodic sharp waves may be better explained by quasi-periodicity than by low-dimensional chaos. The EEG further displayed at least one example of a sudden, large qualitative change in the dynamics, highly suggestive of a bifurcation. The presence of quasi-periodicity and bifurcations strongly argues for the use of a non-linear model to describe the EEG in Creutzfeldt-Jakob disease. Received: 28 October 1996 / Accepted in revised form: 8 July 1997     

Application of a dissimilarity index of EEG and its sub-bands on prediction of induced epileptic seizures from rat's EEG signals

ObjectiveEpileptic seizures are defined as manifest of excessive and hyper-synchronous activity of neurons in the cerebral cortex that cause frequent malfunction of the human central nervous system. Therefore, finding precursors and predictors of epileptic seizure is of utmost clinical relevance to reduce the epileptic seizure induced nervous system malfunction consequences. Researchers for this purpose may even guide us to a deep understanding of the seizure generating mechanisms. The goal of this paper is to predict epileptic seizures in epileptic rats.MethodsSeizures were induced in rats using pentylenetetrazole (PTZ) model. EEG signals in interictal, preictal, ictal and postictal periods were then recorded and analyzed to predict epileptic seizures. Epileptic seizures were predicted by calculating an index in consecutive windows of EEG signal and comparing the index with a threshold. In this work, a newly proposed dissimilarity index called Bhattacharyya Based Dissimilarity Index (BBDI), dynamical similarity index and fuzzy similarity index were investigated.ResultsBBDI, dynamical similarity index and fuzzy similarity index were examined on case and control groups and compared to each other. The results show that BBDI outperforms dynamical and fuzzy similarity indices. In order to improve the results, EEG sub-bands were also analyzed. The best result achieved when the proposed dissimilarity index was applied on Delta sub-band that predicts epileptic seizures in all rats with a mean of 299.5 s.ConclusionThe dissimilarity of neural network activity between reference window and present window of EEG signal has a significant increase prior to an epileptic seizure and the proposed dissimilarity index (BBDI) can reveal this variation to predict epileptic seizures. In addition, analyzing EEG sub-bands results in more accurate information about constituent neuronal activities underlying the EEG since certain changes in EEG signal may be amplified when each sub-band is analyzed separately.SignificanceThis paper presents application of a dissimilarity index (BBDI) on EEG signals and its sub-bands to predict PTZ-induced epileptic seizures in rats. Based on the results of this work, BBDI will predict epileptic seizures more accurately and more reliably compared to current indices that increases epileptic patient comfort and improves patient outcomes.     

Complexity of Multi-Channel Electroencephalogram Signal Analysis in Childhood Absence Epilepsy

Absence epilepsy is an important epileptic syndrome in children. Multiscale entropy (MSE), an entropy-based method to measure dynamic complexity at multiple temporal scales, is helpful to disclose the information of brain connectivity. This study investigated the complexity of electroencephalogram (EEG) signals using MSE in children with absence epilepsy. In this research, EEG signals from 19 channels of the entire brain in 21 children aged 5-12 years with absence epilepsy were analyzed. The EEG signals of pre-ictal (before seizure) and ictal states (during seizure) were analyzed by sample entropy (SamEn) and MSE methods. Variations of complexity index (CI), which was calculated from MSE, from the pre-ictal to the ictal states were also analyzed. The entropy values in the pre-ictal state were significantly higher than those in the ictal state. The MSE revealed more differences in analysis compared to the SamEn. The occurrence of absence seizures decreased the CI in all channels. Changes in CI were also significantly greater in the frontal and central parts of the brain, indicating fronto-central cortical involvement of “cortico-thalamo-cortical network” in the occurrence of generalized spike and wave discharges during absence seizures. Moreover, higher sampling frequency was more sensitive in detecting functional changes in the ictal state. There was significantly higher correlation in ictal states in the same patient in different seizures but there were great differences in CI among different patients, indicating that CI changes were consistent in different absence seizures in the same patient but not from patient to patient. This implies that the brain stays in a homogeneous activation state during the absence seizures. In conclusion, MSE analysis is better than SamEn analysis to analyze complexity of EEG, and CI can be used to investigate the functional brain changes during absence seizures.     

Acute unstable depressive syndrome (AUDS) is associated more frequently with epilepsy than major depression

Background  

Depressive disorders are frequent in epilepsy and associated with reduced seizure control. Almost 50% of interictal depressive disorders have to be classified as atypical depressions according to DSM-4 criteria. Research has mainly focused on depressive symptoms in defined populations with epilepsy (e.g., patients admitted to tertiary epilepsy centers). We have chosen the opposite approach. We hypothesized that it is possible to define by clinical means a subgroup of psychiatric patients with higher than expected prevalence of epilepsy and seizures. We hypothesized further that these patients present with an Acute Unstable Depressive Syndrome (AUDS) that does not meet DSM-IV criteria of a Major Depressive Episode (MDE). In a previous publication we have documented that AUDS patients indeed have more often a history of epileptic seizures and abnormal EEG recordings than MDE patients (Vaaler et al. 2009). This study aimed to further classify the differences of depressive symptoms at admittance and follow-up of patients with AUDS and MDE.     

基于二阶C0复杂度的癫痫发作预测

癫痫发作的预测是近年来在临床医学和神经系统科学研究领域中备受关注的问题。如果癫痫发作能够被可靠地预测,则可以提前采取有效的临床预防措施,从而能较大程度地改善癫痫患者的生活质量。文章提出了一种基于二阶C0复杂度的预测算法用于预测癫痫发作。该算法通过分析癫痫患者颅内脑电信号的二阶C0复杂度,利用发作前期复杂度曲线的变化特征预测癫痫发作。作者运用该算法对21组癫痫病人87次发作的临床颅内脑电数据和4组大鼠4次发作的颅内脑电数据进行分析计算,预测准确率分别为94.3%和100%。实验结果表明该算法可以有效地预测癫痫发作,具有潜在的重要临床应用价值。     

Application of recurrence quantification analysis for the automated identification of epileptic EEG signals

Epilepsy is a common neurological disorder that is characterized by the recurrence of seizures. Electroencephalogram (EEG) signals are widely used to diagnose seizures. Because of the non-linear and dynamic nature of the EEG signals, it is difficult to effectively decipher the subtle changes in these signals by visual inspection and by using linear techniques. Therefore, non-linear methods are being researched to analyze the EEG signals. In this work, we use the recorded EEG signals in Recurrence Plots (RP), and extract Recurrence Quantification Analysis (RQA) parameters from the RP in order to classify the EEG signals into normal, ictal, and interictal classes. Recurrence Plot (RP) is a graph that shows all the times at which a state of the dynamical system recurs. Studies have reported significantly different RQA parameters for the three classes. However, more studies are needed to develop classifiers that use these promising features and present good classification accuracy in differentiating the three types of EEG segments. Therefore, in this work, we have used ten RQA parameters to quantify the important features in the EEG signals.These features were fed to seven different classifiers: Support vector machine (SVM), Gaussian Mixture Model (GMM), Fuzzy Sugeno Classifier, K-Nearest Neighbor (KNN), Naive Bayes Classifier (NBC), Decision Tree (DT), and Radial Basis Probabilistic Neural Network (RBPNN). Our results show that the SVM classifier was able to identify the EEG class with an average efficiency of 95.6%, sensitivity and specificity of 98.9% and 97.8%, respectively.     

Characteristic nonlinearities of the 3/s ictal electroencephalogram identified by nonlinear autoregressive analysis

We describe a method for the characterization of electroencephalographic (EEG) signals based on a model which features nonlinear feedback. The characteristic EEG fingerprints obtained through this approach display the time-course of nonlinear interactions, rather than aspects susceptible to standard spectral analysis. Fingerprints of seizure discharges in six patients (five with typical absence seizures, one with complex partial seizures) revealed significant nonlinear interactions. The timing and pattern of these interactions correlated closely with the seizure type. Nonlinear autoregressive (NLAR) analysis is compared with other nonlinear dynamical measures that have been applied to the EEG.     

Identifying an increased risk of epileptic seizures using a multi-feature EEG–ECG classification

Epilepsy, a neurological disorder in which patients suffer from recurring seizures, affects approximately 1% of the world population. In spite of available drug and surgical treatment options, more than 25% of individuals with epilepsy have seizures that are uncontrollable. For these patients with intractable epilepsy, the unpredictability of seizure occurrence underlies an enhanced risk of sudden unexpected death or morbidity. A system that could warn the patient of the impending event or trigger an antiepileptic device would dramatically increase the quality of life for those patients. Here, we proposed a patient-specific algorithm for possible seizure warning using machine learning classification of 34 algorithmic features derived from EEG–ECG recordings. We evaluated our algorithm on unselected and continuous recordings of 12 patients (total of 108 seizures and 3178-h). Good out-of-sample performances were observed around 25% of the patients with an average preictal period around 30 min and independently of the EEG type (scalp or intracranial). Inspection of the most discriminative EEG–ECG features revealed that good classification rates reflected specific physiological precursors, particularly related to certain stages of sleep. From these observations, we conclude that our algorithmic strategy enables a quantitative way to identify “pro-ictal” states with a high risk of seizure generation.     

癫痫病人脑电信号的奇异谱

癫痫是一种常见的神经系统疾患,其唯一客观证据为脑电图的癫痫样发放。在癫痫发作间期,仅有偶发的很难辨别的癫痫样放电,为了正确诊断癫痫病,往往需要医生长时间监测病人的脑电信号,在对脑电信号进行相空间重构,进而对其进行奇异系统分析,发现癫痫病人无论在癫痫发作前、发作中、发作后,其脑电信号的奇异谱曲线不存在噪声平台,明显区别于正常人。是否可以认为脑电信号的奇异谱正代表着大脑的一种基本状态,癫痫患者在未发作时,大脑的基本状态已经处于异常。无论如休,奇异系统分析方法使得可以利用很短的一段脑电数据诊断癫痫。无疑为癫痫病人的临床诊断提供了一条简单、有效的途径。     

Adaptive prediction of epileptic seizures from intracranial recordings

Previous work has demonstrated that some dynamic properties of intracranial EEG signals are indicative of epileptic seizures and hence could be used for prediction in order to realize counter measures. However, most previous studies only investigated predictability via offline analysis of EEG signals as compared to actually predicting seizures in a setting applicable to implantable devices. Here we address this problem, which calls for simple and fast online methods, and based on previous offline analyses we hypothesize that prediction can be further improved when using multiple features to detect the preictal patterns. We propose a simple adaptive online method (an evolving neuro-fuzzy model) to adaptively learn such combined features. The classifier starts out with a simple structure and patient-independent parameters and then grows into a personal seizure predictor as recursive methods tune the model structure and parameters. We apply the adaptive classifier to a publicly available database of intracranial recordings from 21 patients and demonstrate that seizure prediction is improved with our online method as compared to offline non-adaptive techniques. We show that our method is robust with respect to those few model parameters, which are not adapted. Moreover, as we report the performance on data from a publicly available seizure database, our results can serve as a yardstick for future method developments.     

Dynamics on Networks: The Role of Local Dynamics and Global Networks on the Emergence of Hypersynchronous Neural Activity

Graph theory has evolved into a useful tool for studying complex brain networks inferred from a variety of measures of neural activity, including fMRI, DTI, MEG and EEG. In the study of neurological disorders, recent work has discovered differences in the structure of graphs inferred from patient and control cohorts. However, most of these studies pursue a purely observational approach; identifying correlations between properties of graphs and the cohort which they describe, without consideration of the underlying mechanisms. To move beyond this necessitates the development of computational modeling approaches to appropriately interpret network interactions and the alterations in brain dynamics they permit, which in the field of complexity sciences is known as dynamics on networks. In this study we describe the development and application of this framework using modular networks of Kuramoto oscillators. We use this framework to understand functional networks inferred from resting state EEG recordings of a cohort of 35 adults with heterogeneous idiopathic generalized epilepsies and 40 healthy adult controls. Taking emergent synchrony across the global network as a proxy for seizures, our study finds that the critical strength of coupling required to synchronize the global network is significantly decreased for the epilepsy cohort for functional networks inferred from both theta (3–6 Hz) and low-alpha (6–9 Hz) bands. We further identify left frontal regions as a potential driver of seizure activity within these networks. We also explore the ability of our method to identify individuals with epilepsy, observing up to 80 predictive power through use of receiver operating characteristic analysis. Collectively these findings demonstrate that a computer model based analysis of routine clinical EEG provides significant additional information beyond standard clinical interpretation, which should ultimately enable a more appropriate mechanistic stratification of people with epilepsy leading to improved diagnostics and therapeutics.