Synchronous and asynchronous bursting states: role of intrinsic neural dynamics

Macroscopic Gamma Oscillation With Bursting Neuron Model Under Stochastic Fluctuation

Gamma oscillations are thought to play a role in information processing in the brain. Bursting neurons, which exhibit periodic clusters of spiking activity, are a type of neuron that are thought to contribute largely to gamma oscillations. However, little is known about how the properties of bursting neurons affect the emergence of gamma oscillation, its waveforms, and its synchronized characteristics, especially when subjected to stochastic fluctuations. In this study, we proposed a bursting neuron model that can analyze the bursting ratio and the phase response function. Then we theoretically analyzed the neuronal population dynamics composed of bursting excitatory neurons, mixed with inhibitory neurons. The bifurcation analysis of the equivalent Fokker-Planck equation exhibits three types of gamma oscillations of unimodal firing, bimodal firing in the inhibitory population, and bimodal firing in the excitatory population under different interaction strengths. The analyses of the macroscopic phase response function by the adjoint method of the Fokker-Planck equation revealed that the inhibitory doublet facilitates synchronization of the high-frequency oscillations. When we keep the strength of interactions constant, decreasing the bursting ratio of the individual neurons increases the relative high-gamma component of the populational phase-coupling functions. This also improves the ability of the neuronal population model to synchronize with faster oscillatory input. The analytical frameworks in this study provide insight into nontrivial dynamics of the population of bursting neurons, which further suggest that bursting neurons have an important role in rhythmic activities.

Bayesian Estimation of Phase Dynamics Based on Partially Sampled Spikes Generated by Realistic Model Neurons

A dynamic system showing stable rhythmic activity can be represented by the dynamics of phase oscillators. This would provide a useful mathematical framework through which one can understand the system's dynamic properties. A recent study proposed a Bayesian approach capable of extracting the underlying phase dynamics directly from time-series data of a system showing rhythmic activity. Here we extended this method to spike data that otherwise provide only limited phase information. To determine how this method performs with spike data, we applied it to simulated spike data generated by a realistic neuronal network model. We then compared the estimated dynamics obtained based on the spike data with the dynamics theoretically derived from the model. The method successfully extracted the modeled phase dynamics, particularly the interaction function, when the amount of available data was sufficiently large. Furthermore, the method was able to infer synaptic connections based on the estimated interaction function. Thus, the method was found to be applicable to spike data and practical for understanding the dynamic properties of rhythmic neural systems.

Neural dynamics and information representation in microcircuits of motor cortex

The brain has to analyze and respond to external events that can change rapidly from time to time, suggesting that information processing by the brain may be essentially dynamic rather than static. The dynamical features of neural computation are of significant importance in motor cortex that governs the process of movement generation and learning. In this paper, we discuss these features based primarily on our recent findings on neural dynamics and information coding in the microcircuit of rat motor cortex. In fact, cortical neurons show a variety of dynamical behavior from rhythmic activity in various frequency bands to highly irregular spike firing. Of particular interest are the similarity and dissimilarity of the neuronal response properties in different layers of motor cortex. By conducting electrophysiological recordings in slice preparation, we report the phase response curves (PRCs) of neurons in different cortical layers to demonstrate their layer-dependent synchronization properties. We then study how motor cortex recruits task-related neurons in different layers for voluntary arm movements by simultaneous juxtacellular and multiunit recordings from behaving rats. The results suggest an interesting difference in the spectrum of functional activity between the superficial and deep layers. Furthermore, the task-related activities recorded from various layers exhibited power law distributions of inter-spike intervals (ISIs), in contrast to a general belief that ISIs obey Poisson or Gamma distributions in cortical neurons. We present a theoretical argument that this power law of in vivo neurons may represent the maximization of the entropy of firing rate with limited energy consumption of spike generation. Though further studies are required to fully clarify the functional implications of this coding principle, it may shed new light on information representations by neurons and circuits in motor cortex.

Higher-order spike triggered analysis of neural oscillators

For the purpose of elucidating the neural coding process based on the neural excitability mechanism, researchers have recently investigated the relationship between neural dynamics and the spike triggered stimulus ensemble (STE). Ermentrout et al. analytically derived the relational equation between the phase response curve (PRC) and the spike triggered average (STA). The STA is the first cumulant of the STE. However, in order to understand the neural function as the encoder more explicitly, it is necessary to elucidate the relationship between the PRC and higher-order cumulants of the STE. In this paper, we give a general formulation to relate the PRC and the nth moment of the STE. By using this formulation, we derive a relational equation between the PRC and the spike triggered covariance (STC), which is the covariance of the STE. We show the effectiveness of the relational equation through numerical simulations and use the equation to identify the feature space of the rat hippocampal CA1 pyramidal neurons from their PRCs. Our result suggests that the hippocampal CA1 pyramidal neurons oscillating in the theta frequency range are commonly sensitive to inputs composed of theta and gamma frequency components.

The role of voltage dependence of the NMDA receptor in cellular and network oscillation

Unraveling the mechanisms underlying oscillatory behavior is critical for understanding normal and pathological brain processes. Here we used electrophysiology in mouse neocortical slices and principles of nonlinear dynamics to demonstrate how an increase in the N-methyl-d-aspartic acid receptor (NMDAR) conductance can create a nonlinear whole-cell current-voltage (I-V) relationship which leads to changes in cellular stability. We discovered two behaviorally and morphologically distinct pyramidal cell populations. Under control conditions, both cell types responded to depolarizing current injection with regular spiking patterns. However, upon NMDAR activation, an intrinsic oscillatory (IO) cell type (n = 44) showed a nonlinear whole-cell I-V relationship, intrinsic voltage-dependent oscillations plus amplification of alternating input current, and these properties persisted after disabling action potential generation with tetrodotoxin (TTX). The other non-oscillatory (NO) neuronal population (n = 24) demonstrated none of these behaviors. Simultaneous intra- and extracellular recordings demonstrated the NMDAR's capacity to promote low-frequency seizure-like network oscillations via its effects on intrinsic neuronal properties. The two pyramidal cell types demonstrated different relationships with network oscillation--the IO cells were leaders that were activated early in the population activity cycle while the activation of the NO cell type was distributed across network bursts. The properties of IO neurons disappeared in a low-magnesium environment where the voltage dependence of the receptor is abolished; concurrently, the cellular contribution to network oscillation switched to synchronous firing. Thus, depending upon the efficacy of NMDAR in altering the linearity of the whole-cell I-V relationship, the two cell populations played different roles in sustaining network oscillation.

Dampening of hyperexcitability in CA1 pyramidal neurons by polyunsaturated fatty acids acting on voltage-gated ion channels

A ketogenic diet is an alternative treatment of epilepsy in infants. The diet, rich in fat and low in carbohydrates, elevates the level of polyunsaturated fatty acids (PUFAs) in plasma. These substances have therefore been suggested to contribute to the anticonvulsive effect of the diet. PUFAs modulate the properties of a range of ion channels, including K and Na channels, and it has been hypothesized that these changes may be part of a mechanistic explanation of the ketogenic diet. Using computational modelling, we here study how experimentally observed PUFA-induced changes of ion channel activity affect neuronal excitability in CA1, in particular responses to synaptic input of high synchronicity. The PUFA effects were studied in two pathological models of cellular hyperexcitability associated with epileptogenesis. We found that experimentally derived PUFA modulation of the A-type K (K_A) channel, but not the delayed-rectifier K channel, restored healthy excitability by selectively reducing the response to inputs of high synchronicity. We also found that PUFA modulation of the transient Na channel was effective in this respect if the channel's steady-state inactivation was selectively affected. Furthermore, PUFA-induced hyperpolarization of the resting membrane potential was an effective approach to prevent hyperexcitability. When the combined effect of PUFA on the K_A channel, the Na channel, and the resting membrane potential, was simulated, a lower concentration of PUFA was needed to restore healthy excitability. We therefore propose that one explanation of the beneficial effect of PUFAs lies in its simultaneous action on a range of ion-channel targets. Furthermore, this work suggests that a pharmacological cocktail acting on the voltage dependence of the Na-channel inactivation, the voltage dependences of K_A channels, and the resting potential can be an effective treatment of epilepsy.

Natural movies evoke spike trains with low spike time variability in cat primary visual cortex

Neuronal responses in primary visual cortex have been found to be highly variable. This has led to the widespread notion that neuronal responses have to be averaged over large numbers of neurons to obtain suitably invariant responses that can be used to reliably encode or represent external stimuli. However, it is possible that the high variability of neuronal responses may result from the use of simple, artificial stimuli and that the visual cortex may respond differently to dynamic, naturalistic images. To investigate this question, we recorded the responses of primary visual cortical neurons in the anesthetized cat under stimulation with time-varying natural movies. We found that cortical neurons on the whole exhibited a high degree of spike count variability, but a surprisingly low degree of spike time variability. The spike count variability was further reduced when all but the first spike in a burst were removed. We also found that responses exhibiting low spike time variability exhibited low spike count variability, suggesting that rate coding and temporal coding might be more compatible than previously thought. In addition, we found the spike time variability to be significantly lower when stimulated by natural movies as compared with stimulation using drifting gratings. Our results indicate that response variability in primary visual cortex is stimulus dependent and significantly lower than previous measurements have indicated.

Influences of membrane properties on phase response curve and synchronization stability in a model globus pallidus neuron

The activity patterns of the globus pallidus (GPe) and subthalamic nucleus (STN) are closely associated with motor function and dysfunction in the basal ganglia. In the pathological state caused by dopamine depletion, the STN-GPe network exhibits rhythmic synchronous activity accompanied by rebound bursts in the STN. Therefore, the mechanism of activity transition is a key to understand basal ganglia functions. As synchronization in GPe neurons could induce pathological STN rebound bursts, it is important to study how synchrony is generated in the GPe. To clarify this issue, we applied the phase-reduction technique to a conductance-based GPe neuronal model in order to derive the phase response curve (PRC) and interaction function between coupled GPe neurons. Using the PRC and interaction function, we studied how the steady-state activity of the GPe network depends on intrinsic membrane properties, varying ionic conductances on the membrane. We noted that a change in persistent sodium current, fast delayed rectifier Kv3 potassium current, M-type potassium current and small conductance calcium-dependent potassium current influenced the PRC shape and the steady state. The effect of those currents on the PRC shape could be attributed to extension of the firing period and reduction of the phase response immediately after an action potential. In particular, the slow potassium current arising from the M-type potassium and the SK current was responsible for the reduction of the phase response. These results suggest that the membrane property modulation controls synchronization/asynchronization in the GPe and the pathological pattern of STN-GPe activity.

A putatively novel form of spontaneous coordination in neural activity

We simultaneously recorded local field potentials from three sites along the olfactory-entorhinal axis in rats lightly anesthetized with isoflurane, as part of another experiment. While analyzing the initial data from that experiment with spectrograms, we discovered a potentially novel form of correlated neural activity, with near-simultaneous occurrence across the three widely separated brain sites. After validating their existence further, we named these events Synchronous Frequency Bursts (SFBs). Here we report our initial investigations into their properties and their potential functional significance. In Experiment 1, we found that SFBs have highly regular properties, consisting of brief (approximately 250 ms), high amplitude bursts of LFP energy spanning frequency ranges from the delta band (1-4 Hz) to at least the low gamma band (30-50 Hz). SFBs occurred almost simultaneously across recording sites, usually with onsets <25 ms apart, and there was no clear pattern of temporal leading or lagging among the sites. While the SFBs had fairly typical, exponentially decaying power spectral density plots, their coherence structure was unusual, with high peaks in several narrow frequency ranges and little coherence in other bands. In Experiment 2, we found that SFBs occurred far more often under light anesthesia than deeper anesthetic states, and were especially prevalent as the animals regained consciousness. Finally, in Experiment 3 we showed that SFBs occur simultaneously at a significant rate across brain sites from putatively different functional subsystems--olfactory versus motor pathways. We suggest that SFBs do not carry information per se, but rather, play a role in coordinating activity in different frequency bands, potentially brain-wide, as animals progress from sleep or anesthesia toward full consciousness.

Complex evolution of spike patterns during burst propagation through feed-forward networks

Stable signal transmission is crucial for information processing by the brain. Synfire-chains, defined as feed-forward networks of spiking neurons, are a well-studied class of circuit structure that can propagate a packet of single spikes while maintaining a fixed packet profile. Here, we studied the stable propagation of spike bursts, rather than single spike activities, in a feed-forward network of a general class of excitable bursting neurons. In contrast to single spikes, bursts can propagate stably without converging to any fixed profiles. Spike timings of bursts continue to change cyclically or irregularly during propagation depending on intrinsic properties of the neurons and the coupling strength of the network. To find the conditions under which bursts lose fixed profiles, we propose an analysis based on timing shifts of burst spikes similar to the phase response analysis of limit-cycle oscillators.

Interplay between a phase response curve and spike-timing-dependent plasticity leading to wireless clustering

A phase response curve (PRC) characterizes the signal transduction between oscillators such as neurons on a fixed network in a minimal manner, while spike-timing-dependent plasiticity (STDP) characterizes the way of rewiring networks in an activity-dependent manner. This paper demonstrates that these two key properties both related to the interaction times of oscillators work synergetically to carve functionally useful circuits. STDP working on neurons that prefer asynchrony converts the initial asynchronous firing to clustered firing with synchrony within a cluster. They get synchronized within a cluster despite their preference to asynchrony because STDP selectively disrupts intracluster connections, which we call wireless clustering. Our PRC analysis reveals a triad mechanism: the network structure affects how the PRC is read out to determine the synchrony tendency, the synchrony tendency affects how the STDP works, and STDP affects the network structure, closing the loop.