Epileptic seizures in a heterogeneous excitatory network with short-term plasticity

The neuronal and synaptic dynamics underlying post-inhibitory rebound burst related to major depressive disorder in the lateral habenula neuron model

A burst behavior observed in the lateral habenula (LHb) neuron related to major depressive disorder has attracted much attention. The burst is induced from silence by the excitatory N-methyl-D-aspartate (NMDA) synapse or by the inhibitory stimulation, i.e., a post-inhibitory rebound (PIR) burst, which has not been explained clearly. In the present paper, the neuronal and synaptic dynamics for the PIR burst are acquired in a theoretical neuron model. At first, dynamic cooperations between the fast rise of inhibitory γ-aminobutyric acid (GABA) synapse, slow rise of NMDA synapse, and T-type calcium current to evoke the PIR burst are obtained. Similar to the inhibitory pulse stimulation, fast rising GABA current can reduce the membrane potential to a level low enough to de-inactivate the low threshold T-type calcium current to evoke a PIR spike, which can enhance the slow rising NMDA current activated at a time before or after the PIR spike. The NMDA current following the PIR spike exhibits slow decay to induce multiple spikes to form the PIR burst. Such results present a theoretical explanation and a candidate for the PIR burst in real LHb neurons. Then, the dynamical mechanism for the PIR spike mediated by the T-type calcium channel is obtained. At large conductance of T-type calcium channel, the resting state corresponds to a stable focus near Hopf bifurcation and exhibits an "uncommon" threshold curve with membrane potential much lower than the resting membrane potential. Inhibitory modulation induces membrane potential decreased to run across the threshold curve to evoke the PIR spike. At small conductance of the T-type calcium channel, a stable node appears and manifests a common threshold curve with higher membrane potential, resulting in non-PIR phenomenon. The results present the dynamic cooperations between neuronal dynamics and fast/slow dynamics of different synapses for the PIR burst observed in the LHb neuron, which is helpful for the modulations to major depressive disorder.

The co-activation patterns of multiple brain regions in Juvenile Myoclonic Epilepsy

Juvenile myoclonic epilepsy (JME) as an idiopathic generalized epilepsy has been studied by many advanced neuroimaging techniques to elucidate its neuroanatomical basis and pathophysiological mechanisms. In this paper, we used co-activation patterns (CAPs) to explore the differences of dynamic brain activity changes in resting state between JME patients and healthy controls. 27 cases JME patients and 27 cases healthy of fMRI data were collected. The structural image data of the subjects were analyzed by voxel-based morphological analysis, and the regions with gray matter volume atrophy and high voxel were selected as the regions of interest. Further, the mean disease duration was used as boundary to divide the patients' data into the below-average time and the above-average time groups, which were defined as patient disease duration groups. And these data were used to construct CAPs and to compare changes in brain dynamics. It was found that the number of patterns occurrences and the possibility of switching between patterns were smaller than those in the healthy control, which indicated patients with damage to brain regions. For the patient time control group, the number of patterns occurrences and the possibility of switching between patterns were similar, while there was linear regression between the three values and disease duration. Collectively, this study provides important evidence for revealing the key brain regions of JME by studying the transformation between CAPs. Future studies could investigate the effects of receiving treatment on patient dynamic brain activity.

The seizure classification of focal epilepsy based on the network motif analysis

Due to the complexity of focal epilepsy and its risk for transiting to the generalized epilepsy, the development of reliable classification methods to accurately predict and classify focal and generalized seizures is critical for the clinical management of patients with epilepsy. In order to holistically understand the seizure propagation behavior of focal epilepsy, we propose a three-node motif reduced network by respectively simplifying the focal region, surrounding healthy region and their critical regions as the single node. Because three-node motif can richly characterize information evolutions, the motif analysis method could comprehensively investigate the seizure behavior of focal epilepsy. Firstly, we define a new seizure propagation marker value to capture the seizure onsets and intensity. Based on the three-node motif analysis, it is shown that the focal seizure and spreading can be categorized as inhibitory seizure, focal seizure, focal-critical seizure and generalized seizures, respectively. The four types of seizures correspond to specific modal types respectively, reflecting the strong correlation between seizure behavior and information flow evolution. In addition, it is found that the intensity difference of outflow and inflow information from the critical node (connection heterogeneity) and the excitability of the critical node significantly affected the distribution and transition of the four seizure types. In particular, the method of local linear stability analysis also verifies the effectiveness of four types of seizures classification. In sum, this paper computationally confirms the complex dynamic behavior of focal seizures, and the study of criticality is helpful to propose novel seizure control strategies.

Nonlinear mechanism for the enhanced bursting activities induced by fast inhibitory autapse and reduced activities by fast excitatory autapse

The paradoxical phenomena that excitatory modulation does not enhance but reduces or inhibitory modulation not suppresses but promotes neural firing activities have attracted increasing attention. In the present study, paradoxical phenomena induced by both fast excitatory and inhibitory autapses in a "Fold/Big Homoclinic" bursting are simulated, and the corresponding nonlinear and biophysical mechanisms are presented. Firstly, the enhanced conductance of excitatory autapse induces the number of spikes per burst and firing rate reduced, while the enhanced inhibitory autapse cause both indicators increased. Secondly, with fast-slow variable dissection, the burst of bursting is identified to locate between a fold bifurcation and a big saddle-homoclinic orbit bifurcation of the fast subsystem. Enhanced excitatory or inhibitory autapses cannot induce changes of both bifurcation points, i.e., burst width. However, width of slow variable between two successive spikes within a burst becomes wider for the excitatory autapse and narrower for the inhibitory autapse, resulting in the less and more spikes per burst, respectively. Last, the autaptic current of fast autapse mainly plays a role during the peak of action potential, differing from the slow autaptic current with exponential decay, which can play roles following the peak of action potential. The fast excitatory autaptic current enhances the amplitude of the action potential and reduces the repolarization of the action potential to lengthen the interspike interval (ISI) of the spiking of the fast subsystem, resulting in the wide width of slow variable between successive spikes. The fast inhibitory autaptic current reduces the amplitude of action potential and ISI of spiking, resulting in narrow width of slow variable. The novel example of the paradoxical responses for both fast modulations and nonlinear mechanism extend the contents of neurodynamics, which presents potential functions of the fast autapse.

Correlation Analysis of Synchronization Type and Degree in Respiratory Neural Network

Pre-Bötzinger complex (PBC) is a necessary condition for the generation of respiratory rhythm. Due to the existence of synaptic gaps, delay plays a key role in the synchronous operation of coupled neurons. In this study, the relationship between synchronization and correlation degree is established for the first time by using ISI bifurcation and correlation coefficient, and the relationship between synchronization and correlation degree is discussed under the conditions of no delay, symmetric delay, and asymmetric delay. The results show that the phase synchronization of two coupling PBCs is closely related to the weak correlation, that is, the weak phase synchronization may occur under the condition of incomplete synchronization. Moreover, the time delay and coupling strength are controlled in the modified PBC network model, which not only reveals the law of PBC firing transition but also reveals the complex synchronization behavior in the coupled chaotic neurons. Especially, when the two coupled neurons are nonidentical, the complete synchronization will disappear. These results fully reveal the dynamic behavior of the PBC neural system, which is helpful to explore the signal transmission and coding of PBC neurons and provide theoretical value for further understanding respiratory rhythm.

Deep-layer motif method for estimating information flow between EEG signals

Accurate identification for the information flow between epileptic seizure signals is the key to construct the directional epileptic brain network which can be used to localize epileptic focus. In this paper, our concern is on how to improve the direction identification of information flow and also investigate how it can be cut off or weakened. In view of this, we propose the deep-layer motif method. Based on the directional index (DI) estimation using permutation conditional mutual information, the effectiveness of the proposed deep-layer motif method is numerically assessed with the coupled mass neural model. Furthermore, we investigate the robustness of this method in considering the interference of autaptic coupling, time delay and short-term plasticity. Results show that compared to the simple 1-layer motif method, the 2nd- and 3rd-layer motif methods have the dominant enhancement effects for the direction identification. In particular, deep-layer motif method possesses good anti-jamming performance and good robustness in calculating DI. In addition, we investigate the effect of deep brain stimulation (DBS) on the information flow. It is found that this deep-layer motif method is still superior to the single-layer motif method in direction identification and is robust to weak DBS. However, the high-frequency strong DBS can effectively decrease the DI suggesting the weakened information flow. These results may give new insights into the seizure detection and control.

Complex dynamics of hair bundle of auditory nervous system (I): spontaneous oscillations and two cases of steady states

The hair bundles of inner hair cells in the auditory nervous exhibit spontaneous oscillations, which is the prerequisite for an important auditory function to enhance the sensitivity of inner ear to weak sounds, otoacoustic emission. In the present paper, the dynamics of spontaneous oscillations and relationships to steady state are acquired in a two-dimensional model with fast variable X (displacement of hair bundles) and slow variable X a . The spontaneous oscillations are derived from negative stiffness modulated by two biological factors (S and D) and are identified to appear in multiple two-dimensional parameter planes. In (S, D) plane, comprehensive bifurcations including 4 types of codimension-2 bifurcation and 5 types of codimension-1 bifurcation related to the spontaneous oscillations are acquired. The spontaneous oscillations are surrounded by supercritical and subcritical Hopf bifurcation curves, and outside of the curves are two cases of steady state. Case-1 and Case-2 steady states exhibit Z-shaped (coexistence of X) and N-shaped (coexistence of X a ) X-nullclines, respectively. In (S, D) plane, left and right to the spontaneous oscillations are two subcases of Case-1, which exhibit the stable equilibrium point locating on the upper and lower branches of X-nullcline, respectively, resembling that of the neuron. Lower to the spontaneous oscillations are 3 subcases of Case-2 from left to right, which manifest stable equilibrium point locating on left, middle, and right branches of X-nullcline, respectively, differing from that of the neuron. The phase plane for steady state is divided into four parts by nullclines, which manifest different vector fields. The phase trajectory of transient behavior beginning from a phase point in the four regions to the stable equilibrium point exhibits different dynamics determined by the vector fields, which is the basis to identify dynamical mechanism of complex forced oscillations induced by external signal. The results present comprehensive viewpoint and deep understanding for dynamics of the spontaneous oscillations and steady states of hair bundles, which can be used to well explain the experimental observations and to modulate functions of spontaneous oscillations.

Extracting the transition network of epileptic seizure onset

In presurgical monitoring, focal seizure onset is visually assessed from intracranial electroencephalogram (EEG), typically based on the selection of channels that show the strongest changes in amplitude and frequency. As epileptic seizure dynamics is increasingly considered to reflect changes in potentially distributed neural networks, it becomes important to also assess the interrelationships between channels. We propose a workflow to quantitatively extract the nodes and edges contributing to the seizure onset using an across-seizure scoring. We propose a quantification of the consistency of EEG channel contributions to seizure onset within a patient. The workflow is exemplified using recordings from patients with different degrees of seizure-onset consistency.