An investigation of EEG dynamics in an animal model of temporal lobe epilepsy using the maximum Lyapunov exponent

Directed Connectivity Analysis of the Neuro-Cardio- and Respiratory Systems Reveals Novel Biomarkers of Susceptibility to SUDEP

Sudden unexpected death in epilepsy (SUDEP) is the leading cause of epilepsy-related mortality and its pathophysiological mechanisms remain unknown. Goal: We set to record and analyze for the first time concurrent electroencephalographic (EEG), electrocardiographic (ECG), and unrestrained whole-body plethysmographic (Pleth) signals from control (WT - wild type) and SUDEP-prone mice (KO- knockout Kcna1 animal model). Methods: Employing multivariate autoregressive models (MVAR) we measured all tri-organ effective directional interactions by the generalized partial directed coherence (GPDC) in the frequency domain over time (hours). Results: When compared to the control (WT) animals, the SUDEP-prone (KO) animals exhibited (p < 0.001) reduced afferent and efferent interactions between the heart and the brain over the full frequency spectrum (0-200Hz), enhanced efferent interactions from the brain to the lungs and from the heart to the lungs at high (>90 Hz) frequencies (especially during periods with seizure activity), and decreased feedback from the lungs to the brain at low (<40 Hz) frequencies. Conclusions: These results show that impairment in the afferent and efferent pathways in the holistic neuro-cardio-respiratory network could lead to SUDEP, and effective connectivity measures and their dynamics could serve as novel biomarkers of susceptibility to SUDEP and seizures respectively.

EEG Analysis in Structural Focal Epilepsy Using the Methods of Nonlinear Dynamics (Lyapunov Exponents, Lempel-Ziv Complexity, and Multiscale Entropy)

This paper analyzes a case with the patient having focal structural epilepsy by processing electroencephalogram (EEG) fragments containing the "sharp wave" pattern of brain activity. EEG signals were recorded using 21 channels. Based on the fact that EEG signals are time series, an approach has been developed for their analysis using nonlinear dynamics tools: calculating the Lyapunov exponent's spectrum, multiscale entropy, and Lempel-Ziv complexity. The calculation of the first Lyapunov exponent is carried out by three methods: Wolf, Rosenstein, and Sano-Sawada, to obtain reliable results. The seven Lyapunov exponent spectra are calculated by the Sano-Sawada method. For the observed patient, studies showed that with medical treatment, his condition did not improve, and as a result, it was recommended to switch from conservative treatment to surgical. The obtained results of the patient's EEG study using the indicated nonlinear dynamics methods are in good agreement with the medical report and MRI data. The approach developed for the analysis of EEG signals by nonlinear dynamics methods can be applied for early detection of structural changes.

Evaluation of dynamic scaling of growing interfaces in EEG fluctuations of seizures in animal model of temporal lobe epilepsy

Epileptic seizures, as a dynamic phenomenon in brain behavior, obey a scale-free behavior, frequently analyzed by electrical activity recording. This recording can be seen as a surface that roughens with time. Dynamic scaling studies roughening processes or growing interfaces. In this theory, a set of exponents -obtained from scale invariance properties- characterize rough interfaces growth. The aim of the present study was to investigate scaling behavior in EEG time series fluctuations of a chemical animal model of temporal lobe epilepsy, with dynamic scaling to detect changes on seizure onset. We analyzed local variables in different sampling intervals and estimated rough, scaling and dynamic exponents. Results exhibited long-range correlations in interictal activity. Results of renormalization and data collapsing confirmed that each epoch of EEG fluctuations for interictal, preictal and postictal collapse in a curve in different scales, each segment independently; remarkably, we found non-scaling behavior in seizures epochs. Data for the different sampling intervals for ictal period do not collapse in one curve, which implies that ictal activity does not exhibit the same scaling behavior than the other epochs. Statistical significant differences of growth exponent were found between interictal and ictal segment, while for scaling exponent, significant differences were found between interictal and postictal segment. These results confirm the potential of scaling exponents as characteristic parameters to detect changes on seizure onset, which suggests their use as inputs for analysis methods for seizure detection in long-term recordings, while changes in growth exponent are potentially useful for prediction purposes.

Dynamic trajectory of multiple single-unit activity during working memory task in rats

Working memory plays an important role in complex cognitive tasks. A popular theoretical view is that transient properties of neuronal dynamics underlie cognitive processing. The question raised here as to how the transient dynamics evolve in working memory. To address this issue, we investigated the multiple single-unit activity dynamics in rat medial prefrontal cortex (mPFC) during a Y-maze working memory task. The approach worked by reconstructing state space from delays of the original single-unit firing rate variables, which were further analyzed using kernel principal component analysis (KPCA). Then the neural trajectories were obtained to visualize the multiple single-unit activity. Furthermore, the maximal Lyapunov exponent (MLE) was calculated to quantitatively evaluate the neural trajectories during the working memory task. The results showed that the neuronal activity produced stable and reproducible neural trajectories in the correct trials while showed irregular trajectories in the incorrect trials, which may establish a link between the neurocognitive process and behavioral performance in working memory. The MLEs significantly increased during working memory in the correctly performed trials, indicating an increased divergence of the neural trajectories. In the incorrect trials, the MLEs were nearly zero and remained unchanged during the task. Taken together, the trial-specific neural trajectory provides an effective way to track the instantaneous state of the neuronal population during the working memory task and offers valuable insights into working memory function. The MLE describes the changes of neural dynamics in working memory and may reflect different neuronal population states in working memory.

Influence of neuropathology on convection-enhanced delivery in the rat hippocampus

Local drug delivery techniques, such as convention-enhanced delivery (CED), are promising novel strategies for delivering therapeutic agents otherwise limited by systemic toxicity and blood-brain-barrier restrictions. CED uses positive pressure to deliver infusate homogeneously into interstitial space, but its distribution is dependent upon appropriate tissue targeting and underlying neuroarchitecture. To investigate effects of local tissue pathology and associated edema on infusate distribution, CED was applied to the hippocampi of rats that underwent electrically-induced, self-sustaining status epilepticus (SE), a prolonged seizure. Infusion occurred 24 hours post-SE, using a macromolecular tracer, the magnetic resonance (MR) contrast agent gadolinium chelated with diethylene triamine penta-acetic acid and covalently attached to albumin (Gd-albumin). High-resolution T1- and T2-relaxation-weighted MR images were acquired at 11.1 Tesla in vivo prior to infusion to generate baseline contrast enhancement images and visualize morphological changes, respectively. T1-weighted imaging was repeated post-infusion to visualize final contrast-agent distribution profiles. Histological analysis was performed following imaging to characterize injury. Infusions of Gd-albumin into injured hippocampi resulted in larger distribution volumes that correlated with increased injury severity, as measured by hyperintense regions seen in T2-weighted images and corresponding histological assessments of neuronal degeneration, myelin degradation, astrocytosis, and microglial activation. Edematous regions included the CA3 hippocampal subfield, ventral subiculum, piriform and entorhinal cortex, amygdalar nuclei, middle and laterodorsal/lateroposterior thalamic nuclei. This study demonstrates MR-visualized injury processes are reflective of cellular alterations that influence local distribution volume, and provides a quantitative basis for the planning of local therapeutic delivery strategies in pathological brain regions.

Complexity and multifractality of neuronal noise in mouse and human hippocampal epileptiform dynamics

Fractal methods offer an invaluable means of investigating turbulent nonlinearity in non-stationary biomedical recordings from the brain. Here, we investigate properties of complexity (i.e. the correlation dimension, maximum Lyapunov exponent, 1/f(γ) noise and approximate entropy) and multifractality in background neuronal noise-like activity underlying epileptiform transitions recorded at the intracellular and local network scales from two in vitro models: the whole-intact mouse hippocampus and lesional human hippocampal slices. Our results show evidence for reduced dynamical complexity and multifractal signal features following transition to the ictal epileptiform state. These findings suggest that pathological breakdown in multifractal complexity coincides with loss of signal variability or heterogeneity, consistent with an unhealthy ictal state that is far from the equilibrium of turbulent yet healthy fractal dynamics in the brain. Thus, it appears that background noise-like activity successfully captures complex and multifractal signal features that may, at least in part, be used to classify and identify brain state transitions in the healthy and epileptic brain, offering potential promise for therapeutic neuromodulatory strategies for afflicted patients suffering from epilepsy and other related neurological disorders.

Kernel duration and modulation gain in a coupled oscillator model and their implications on the progression of seizures

The coupled oscillator model has previously been used for the simulation of neuronal activities in in vitro rat hippocampal slice seizure data and the evaluation of seizure suppression algorithms. Each model unit can be described as either an oscillator which can generate action potential spike trains without inputs, or a threshold-based unit. With the change of only one parameter, each unit can either be an oscillator or a threshold-based spiking unit. This would eliminate the need of a new set of equations for each type of unit. Previous analysis has suggested that long kernel duration and imbalance of inhibitory feedback can cause the system to intermittently transition into and out of ictal activities. The state transitions of seizure-like events were investigated here; specifically, how the system excitability may change when the system underwent transitions in the preictal and postictal processes. Analysis showed that the area of the excitation kernel is positively correlated with the mean firing rate of ictal activity. The kernel duration is also correlated to the amount of ictal activity. The transition into ictal involved the escape from the saddle point foci in the state space trajectory identified using Newton's method.

RESPONSIVE NEUROMODULATORS BASED ON ARTIFICIAL NEURAL NETWORKS USED TO CONTROL SEIZURE-LIKE EVENTS IN A COMPUTATIONAL MODEL OF EPILEPSY

Deep brain stimulation (DBS) has been noted for its potential to suppress epileptic seizures. To date, DBS has achieved mixed results as a therapeutic approach to seizure control. Using a computational model, we demonstrate that high-complexity, biologically-inspired responsive neuromodulation is superior to periodic forms of neuromodulation (responsive and non-responsive) such as those implemented in DBS, as well as neuromodulation using random and random repetitive-interval stimulation. We configured radial basis function (RBF) networks to generate outputs modeling interictal time series recorded from rodent hippocampal slices that were perfused with low Mg2+/high K+solution. We then compared the performance of RBF-based interictal modulation, periodic biphasic-pulse modulation, random modulation and random repetitive modulation on a cognitive rhythm generator (CRG) model of spontaneous seizure-like events (SLEs), testing efficacy of SLE control. A statistically significant improvement in SLE mitigation for the RBF interictal modulation case versus the periodic and random cases was observed, suggesting that the use of biologically-inspired neuromodulators may achieve better results for the purpose of electrical control of seizures in a clinical setting.

Seizure prediction and its applications

Epilepsy is characterized by intermittent, paroxysmal, hypersynchronous electrical activity that may remain localized and/or spread and severely disrupt the brain's normal multitask and multiprocessing function. Epileptic seizures are the hallmarks of such activity. The ability to issue warnings in real time of impending seizures may lead to novel diagnostic tools and treatments for epilepsy. Applications may range from a warning to the patient to avert seizure-associated injuries, to automatic timely administration of an appropriate stimulus. Seizure prediction could become an integral part of the treatment of epilepsy through neuromodulation, especially in the new generation of closed-loop seizure control systems.

Complexity in neuronal noise depends on network interconnectivity

"Noise," or noise-like activity (NLA), defines background electrical membrane potential fluctuations at the cellular level of the nervous system, comprising an important aspect of brain dynamics. Using whole-cell voltage recordings from fast-spiking stratum oriens interneurons and stratum pyramidale neurons located in the CA3 region of the intact mouse hippocampus, we applied complexity measures from dynamical systems theory (i.e., 1/f(γ) noise and correlation dimension) and found evidence for complexity in neuronal NLA, ranging from high- to low-complexity dynamics. Importantly, these high- and low-complexity signal features were largely dependent on gap junction and chemical synaptic transmission. Progressive neuronal isolation from the surrounding local network via gap junction blockade (abolishing gap junction-dependent spikelets) and then chemical synaptic blockade (abolishing excitatory and inhibitory post-synaptic potentials), or the reverse order of these treatments, resulted in emergence of high-complexity NLA dynamics. Restoring local network interconnectivity via blockade washout resulted in resolution to low-complexity behavior. These results suggest that the observed increase in background NLA complexity is the result of reduced network interconnectivity, thereby highlighting the potential importance of the NLA signal to the study of network state transitions arising in normal and abnormal brain dynamics (such as in epilepsy, for example).

Early MR diffusion and relaxation changes in the parahippocampal gyrus precede the onset of spontaneous seizures in an animal model of chronic limbic epilepsy

Structural changes in limbic regions are often observed in individuals with temporal lobe epilepsy (TLE) and in animal models. However, the brain structural changes during the evolution into epilepsy remain largely unknown. Therefore, the purpose of this study was to define the temporal changes in limbic structures after experimental status epilepticus (SE) during the latency period of epileptogenesis in vivo, with quantitative diffusion tensor imaging (DTI) and T2 relaxometry in an animal model of chronic TLE. A pair of fifty micron electrodes was implanted into the ventral hippocampus in twelve male adult rats. Self-sustaining SE was induced with electrical stimulation in eleven rats. Three rats served as age-matched controls. In vivo diffusion tensor and T2 magnetic resonance imaging (MRI) was performed at 11.1 Tesla, pre- and post-implantation of electrodes and 3, 5, 7, 10, 20, 40 and 60 days post-SE to assess structural changes. Spontaneous seizures were identified with continuous time-locked video-monitoring. Following imaging in vivo, fixed, excised brains were MR imaged at 17.6 Tesla. Subsequently, histological analysis was correlated with MRI results. Following SE, 8/11 injured rats developed spontaneous seizures. Unique to these 8 rats, early T2, diffusivity and anisotropy changes were observed in vivo within the parahippocampal gyrus (contralateral) and fimbria (bilateral). In excised brains, bilateral increase in anisotropy was observed in the dentate gyrus, corresponding to mossy fiber sprouting as determined by Timm staining. Using T2 relaxometry and DTI, specific transient and long-term structural changes were observed only in rats that developed spontaneous limbic seizures.