Differential control of active and silent phases in relaxation models of neuronal rhythms

Firing rate models for gamma oscillations in I-I and E-I networks

Firing rate models for describing the mean-field activities of neuronal ensembles can be used effectively to study network function and dynamics, including synchronization and rhythmicity of excitatory-inhibitory populations. However, traditional Wilson-Cowan-like models, even when extended to include an explicit dynamic synaptic activation variable, are found unable to capture some dynamics such as Interneuronal Network Gamma oscillations (ING). Use of an explicit delay is helpful in simulations at the expense of complicating mathematical analysis. We resolve this issue by introducing a dynamic variable, u, that acts as an effective delay in the negative feedback loop between firing rate (r) and synaptic gating of inhibition (s). In effect, u endows synaptic activation with second order dynamics. With linear stability analysis, numerical branch-tracking and simulations, we show that our r-u-s rate model captures some key qualitative features of spiking network models for ING. We also propose an alternative formulation, a v-u-s model, in which mean membrane potential v satisfies an averaged current-balance equation. Furthermore, we extend the framework to E-I networks. With our six-variable v-u-s model, we demonstrate in firing rate models the transition from Pyramidal-Interneuronal Network Gamma (PING) to ING by increasing the external drive to the inhibitory population without adjusting synaptic weights. Having PING and ING available in a single network, without invoking synaptic blockers, is plausible and natural for explaining the emergence and transition of two different types of gamma oscillations.

Canard-induced complex oscillations in an excitatory network

In Neuroscience, mathematical modelling involving multiple spatial and temporal scales can unveil complex oscillatory activity such as excitable responses to an input current, subthreshold oscillations, spiking or bursting. While the number of slow and fast variables and the geometry of the system determine the type of the complex oscillations, canard structures define boundaries between them. In this study, we use geometric singular perturbation theory to identify and characterise boundaries between different dynamical regimes in multiple-timescale firing rate models of the developing spinal cord. These rate models are either three or four dimensional with state variables chosen within an overall group of two slow and two fast variables. The fast subsystem corresponds to a recurrent excitatory network with fast activity-dependent synaptic depression, and the slow variables represent the cell firing threshold and slow activity-dependent synaptic depression, respectively. We start by demonstrating canard-induced bursting and mixed-mode oscillations in two different three-dimensional rate models. Then, in the full four-dimensional model we show that a canard-mediated slow passage creates dynamics that combine these complex oscillations and give rise to mixed-mode bursting oscillations (MMBOs). We unveil complicated isolas along which MMBOs exist in parameter space. The profile of solutions along each isola undergoes canard-mediated transitions between the sub-threshold regime and the bursting regime; these explosive transitions change the number of oscillations in each regime. Finally, we relate the MMBO dynamics to experimental recordings and discuss their effects on the silent phases of bursting patterns as well as their potential role in creating subthreshold fluctuations that are often interpreted as noise. The mathematical framework used in this paper is relevant for modelling multiple timescale dynamics in excitable systems.

Divisive gain modulation enables flexible and rapid entrainment in a neocortical microcircuit model

Neocortical circuits exhibit a rich dynamic repertoire, and their ability to achieve entrainment (adjustment of their frequency to match the input frequency) is thought to support many cognitive functions and indicate functional flexibility. Although previous studies have explored the influence of various circuit properties on this phenomenon, the role of divisive gain modulation (or divisive inhibition) is unknown. This gain control mechanism is thought to be delivered mainly by the soma-targeting interneurons in neocortical microcircuits. In this study, we use a neural mass model of the neocortical microcircuit (extended Wilson-Cowan model) featuring both soma-targeting and dendrite-targeting interneuronal subpopulations to investigate the role of divisive gain modulation in entrainment. Our results demonstrate that the presence of divisive inhibition in the microcircuit, as delivered by the soma-targeting interneurons, enables its entrainment to a wider range of input frequencies. Divisive inhibition also promotes a faster entrainment, with the microcircuit needing less time to converge to the fully entrained state. We suggest that divisive inhibition, working alongside subtractive inhibition, allows for more adaptive oscillatory responses in neocortical circuits and, thus, supports healthy brain functioning.NEW & NOTEWORTHY We introduce a computational neocortical microcircuit model that features two inhibitory neural populations, with one providing subtractive and the other divisive inhibition to the excitatory population. We demonstrate that divisive inhibition widens the range of input frequencies to which the microcircuit can become entrained and diminishes the time needed to reach full entrainment. We suggest that divisive inhibition enables more adaptive oscillatory activity, with important implications for both normal and pathological brain function.

Mechanisms underlying sharpening of visual response dynamics with familiarity

Experience-dependent modifications of synaptic connections are thought to change patterns of network activities and stimulus tuning with learning. However, only a few studies explored how synaptic plasticity shapes the response dynamics of cortical circuits. Here, we investigated the mechanism underlying sharpening of both stimulus selectivity and response dynamics with familiarity observed in monkey inferotemporal cortex. Broadening the distribution of activities and stronger oscillations in the response dynamics after learning provide evidence for synaptic plasticity in recurrent connections modifying the strength of positive feedback. Its interplay with slow negative feedback via firing rate adaptation is critical in sharpening response dynamics. Analysis of changes in temporal patterns also enables us to disentangle recurrent and feedforward synaptic plasticity and provides a measure for the strengths of recurrent synaptic plasticity. Overall, this work highlights the importance of analyzing changes in dynamics as well as network patterns to further reveal the mechanisms of visual learning.

Encoding information into autonomously bursting neural network with pairs of time-delayed pulses

Biological neural networks with many plastic synaptic connections can store external input information in the map of synaptic weights as a form of unsupervised learning. However, the same neural network often produces dramatic reverberating events in which many neurons fire almost simultaneously – a phenomenon coined as ‘population burst.’ The autonomous bursting activity is a consequence of the delicate balance between recurrent excitation and self-inhibition; as such, any periodic sequences of burst-generating stimuli delivered even at a low frequency (~1 Hz) can easily suppress the entire network connectivity. Here we demonstrate that ‘Δt paired-pulse stimulation’, can be a novel way for encoding spatially-distributed high-frequency (~10 Hz) information into such a system without causing a complete suppression. The encoded memory can be probed simply by delivering multiple probing pulses and then estimating the precision of the arrival times of the subsequent evoked recurrent bursts.

Propensity for Bistability of Bursting and Silence in the Leech Heart Interneuron

The coexistence of neuronal activity regimes has been reported under normal and pathological conditions. Such multistability could enhance the flexibility of the nervous system and has many implications for motor control, memory, and decision making. Multistability is commonly promoted by neuromodulation targeting specific membrane ionic currents. Here, we investigated how modulation of different ionic currents could affect the neuronal propensity for bistability. We considered a leech heart interneuron model. It exhibits bistability of bursting and silence in a narrow range of the leak current parameters, conductance (g_leak) and reversal potential (E_leak). We assessed the propensity for bistability of the model by using bifurcation diagrams. On the diagram (g_leak, E_leak), we mapped bursting and silent regimes. For the canonical value of E_leak we determined the range of g_leak which supported the bistability. We use this range as an index of propensity for bistability. We investigated how this index was affected by alterations of ionic currents. We systematically changed their conductances, one at a time, and built corresponding bifurcation diagrams in parameter planes of the maximal conductance of a given current and the leak conductance. We found that conductance of only one current substantially affected the index of propensity; the increase of the maximal conductance of the hyperpolarization-activated cationic current increased the propensity index. The second conductance with the strongest effect was the conductance of the low-threshold fast Ca²⁺ current; its reduction increased the propensity index although the effect was about two times smaller in magnitude. Analyzing the model with both changes applied simultaneously, we found that the diagram (g_leak, E_leak) showed a progressively expanded area of bistability of bursting and silence.

Phasic Firing and Coincidence Detection by Subthreshold Negative Feedback: Divisive or Subtractive or, Better, Both

Phasic neurons typically fire only for a fast-rising input, say at the onset of a step current, but not for steady or slow inputs, a property associated with type III excitability. Phasic neurons can show extraordinary temporal precision for phase locking and coincidence detection. Exemplars are found in the auditory brain stem where precise timing is used in sound localization. Phasicness at the cellular level arises from a dynamic, voltage-gated, negative feedback that can be recruited subthreshold, preventing the neuron from reaching spike threshold if the voltage does not rise fast enough. We consider two mechanisms for phasicness: a low threshold potassium current (subtractive mechanism) and a sodium current with subthreshold inactivation (divisive mechanism). We develop and analyze three reduced models with either divisive or subtractive mechanisms or both to gain insight into the dynamical mechanisms for the potentially high temporal precision of type III-excitable neurons. We compare their firing properties and performance for a range of stimuli. The models have characteristic non-monotonic input-output relations, firing rate vs. input intensity, for either stochastic current injection or Poisson-timed excitatory synaptic conductance trains. We assess performance according to precision of phase-locking and coincidence detection by the models' responses to repetitive packets of unitary excitatory synaptic inputs with more or less temporal coherence. We find that each mechanism contributes features but best performance is attained if both are present. The subtractive mechanism confers extraordinary precision for phase locking and coincidence detection but only within a restricted parameter range when the divisive mechanism of sodium inactivation is inoperative. The divisive mechanism guarantees robustness of phasic properties, without compromising excitability, although with somewhat less precision. Finally, we demonstrate that brief transient inhibition if properly timed can enhance the reliability of firing.

Gain control through divisive inhibition prevents abrupt transition to chaos in a neural mass model

Experimental results suggest that there are two distinct mechanisms of inhibition in cortical neuronal networks: subtractive and divisive inhibition. They modulate the input-output function of their target neurons either by increasing the input that is needed to reach maximum output or by reducing the gain and the value of maximum output itself, respectively. However, the role of these mechanisms on the dynamics of the network is poorly understood. We introduce a novel population model and numerically investigate the influence of divisive inhibition on network dynamics. Specifically, we focus on the transitions from a state of regular oscillations to a state of chaotic dynamics via period-doubling bifurcations. The model with divisive inhibition exhibits a universal transition rate to chaos (Feigenbaum behavior). In contrast, in an equivalent model without divisive inhibition, transition rates to chaos are not bounded by the universal constant (non-Feigenbaum behavior). This non-Feigenbaum behavior, when only subtractive inhibition is present, is linked to the interaction of bifurcation curves in the parameter space. Indeed, searching the parameter space showed that such interactions are impossible when divisive inhibition is included. Therefore, divisive inhibition prevents non-Feigenbaum behavior and, consequently, any abrupt transition to chaos. The results suggest that the divisive inhibition in neuronal networks could play a crucial role in keeping the states of order and chaos well separated and in preventing the onset of pathological neural dynamics.

Determining the contributions of divisive and subtractive feedback in the Hodgkin-Huxley model

The Hodgkin-Huxley (HH) model is the basis for numerous neural models. There are two negative feedback processes in the HH model that regulate rhythmic spiking. The first is an outward current with an activation variable n that has an opposite influence to the excitatory inward current and therefore provides subtractive negative feedback. The other is the inactivation of an inward current with an inactivation variable h that reduces the amount of positive feedback and therefore provides divisive feedback. Rhythmic spiking can be obtained with either negative feedback process, so we ask what is gained by having two feedback processes. We also ask how the different negative feedback processes contribute to spiking. We show that having two negative feedback processes makes the HH model more robust to changes in applied currents and conductance densities than models that possess only one negative feedback variable. We also show that the contributions made by the subtractive and divisive feedback variables are not static, but depend on time scales and conductance values. In particular, they contribute differently to the dynamics in Type I versus Type II neurons.

Scale-free bursting in human cortex following hypoxia at birth

The human brain is fragile in the face of oxygen deprivation. Even a brief interruption of metabolic supply at birth challenges an otherwise healthy neonatal cortex, leading to a cascade of homeostatic responses. During recovery from hypoxia, cortical activity exhibits a period of highly irregular electrical fluctuations known as burst suppression. Here we show that these bursts have fractal properties, with power-law scaling of burst sizes across a remarkable 5 orders of magnitude and a scale-free relationship between burst sizes and durations. Although burst waveforms vary greatly, their average shape converges to a simple form that is asymmetric at long time scales. Using a simple computational model, we argue that this asymmetry reflects activity-dependent changes in the excitatory-inhibitory balance of cortical neurons. Bursts become more symmetric following the resumption of normal activity, with a corresponding reorganization of burst scaling relationships. These findings place burst suppression in the broad class of scale-free physical processes termed crackling noise and suggest that the resumption of healthy activity reflects a fundamental reorganization in the relationship between neuronal activity and its underlying metabolic constraints.

Variety of alternative stable phase-locking in networks of electrically coupled relaxation oscillators

We studied the dynamics of a large-scale model network comprised of oscillating electrically coupled neurons. Cells are modeled as relaxation oscillators with short duty cycle, so they can be considered either as models of pacemaker cells, spiking cells with fast regenerative and slow recovery variables or firing rate models of excitatory cells with synaptic depression or cellular adaptation.

It was already shown that electrically coupled relaxation oscillators exhibit not only synchrony but also anti-phase behavior if electrical coupling is weak. We show that a much wider spectrum of spatiotemporal patterns of activity can emerge in a network of electrically coupled cells as a result of switching from synchrony, produced by short external signals of different spatial profiles. The variety of patterns increases with decreasing rate of neuronal firing (or duty cycle) and with decreasing strength of electrical coupling. We study also the effect of network topology - from all-to-all – to pure ring connectivity, where only the closest neighbors are coupled.

We show that the ring topology promotes anti-phase behavior as compared to all-to-all coupling. It also gives rise to a hierarchical organization of activity: during each of the main phases of a given pattern cells fire in a particular sequence determined by the local connectivity. We have analyzed the behavior of the network using geometric phase plane methods and we give heuristic explanations of our findings.

Our results show that complex spatiotemporal activity patterns can emerge due to the action of stochastic or sensory stimuli in neural networks without chemical synapses, where each cell is equally coupled to others via gap junctions. This suggests that in developing nervous systems where only electrical coupling is present such a mechanism can lead to the establishment of proto-networks generating premature multiphase oscillations whereas the subsequent emergence of chemical synapses would later stabilize generated patterns.

From neural plate to cortical arousal-a neuronal network theory of sleep derived from in vitro "model" systems for primordial patterns of spontaneous bioelectric activity in the vertebrate central nervous system

In the early 1960s intrinsically generated widespread neuronal discharges were discovered to be the basis for the earliest motor behavior throughout the animal kingdom. The pattern generating system is in fact programmed into the developing nervous system, in a regionally specific manner, already at the early neural plate stage. Such rhythmically modulated phasic bursts were next discovered to be a general feature of developing neural networks and, largely on the basis of experimental interventions in cultured neural tissues, to contribute significantly to their morpho-physiological maturation. In particular, the level of spontaneous synchronized bursting is homeostatically regulated, and has the effect of constraining the development of excessive network excitability. After birth or hatching, this "slow-wave" activity pattern becomes sporadically suppressed in favor of sensory oriented "waking" behaviors better adapted to dealing with environmental contingencies. It nevertheless reappears periodically as "sleep" at several species-specific points in the diurnal/nocturnal cycle. Although this "default" behavior pattern evolves with development, its essential features are preserved throughout the life cycle, and are based upon a few simple mechanisms which can be both experimentally demonstrated and simulated by computer modeling. In contrast, a late onto- and phylogenetic aspect of sleep, viz., the intermittent "paradoxical" activation of the forebrain so as to mimic waking activity, is much less well understood as regards its contribution to brain development. Some recent findings dealing with this question by means of cholinergically induced "aroused" firing patterns in developing neocortical cell cultures, followed by quantitative electrophysiological assays of immediate and longterm sequelae, will be discussed in connection with their putative implications for sleep ontogeny.

Call it sleep -- what animals without backbones can tell us about the phylogeny of intrinsically generated neuromotor rhythms during early development

A comprehensive overview is presented of the literature dealing with the development of sleep-like motility and neuronal activity patterns in non-vertebrate animals. it has been established that spontaneous, periodically modulated, neurogenic bursts of movement appear to be a universal feature of prenatal behavior. New empirical data are presented showing that such' seismic sleep' or 'rapid-body-movement' bursts in cuttlefish persist for some time after birth. Extensive ontogenetic research in both vertebrates and non-vertebrates is thus essential before current hypotheses about the phylogeny of motorically active sleep-like states can be taken seriously.

No phylogeny without ontogeny: a comparative and developmental search for the sources of sleep-like neural and behavioral rhythms

A comprehensive review is presented of reported aspects and putative mechanisms of sleep-like motility rhythms throughout the animal kingdom. It is proposed that 'rapid eye movement (REM) sleep' be regarded as a special case of a distinct but much broader category of behavior, 'rapid body movement (RBM) sleep', defined by intrinsically-generated and apparently non-purposive movements. Such a classification completes a 2 × 2 matrix defined by the axes sleep versus waking and active versus quiet. Although 'paradoxical' arousal of forebrain electrical activity is restricted to warm-blooded vertebrates, we urge that juvenile or even infantile stages of development be investigated in cold-blooded animals, in view of the many reports of REM-like spontaneous motility (RBMs) in a wide range of species during sleep. The neurophysiological bases for motorically active sleep at the brainstem level and for slow-wave sleep in the forebrain appear to be remarkably similar, and to be subserved in both cases by a primitive diffuse mode of neuronal organization. Thus, the spontaneous synchronous burst discharges which are characteristics of the sleeping brain can be readily simulated even by highly unstructured neural network models. Neuromotor discharges during active sleep appear to reflect a hierarchy of simple relaxation oscillation mechanisms, spanning a wide range of spike-dependent relaxation times, whereas the periodic alternation of active and quiet sleep states more likely results from the entrainment of intrinsic cellular rhythms and/or from activity-dependent homeostatic changes in network excitability.

Modulating the precision of recurrent bursts in cultured neural networks

Synchronized bursts are a very common feature in biological neural networks, and they play an important role in various brain functions and neurological diseases. This Letter investigates "recurrent synchronized bursts" induced by a single pulse stimulation in cultured networks of rat cortical neurons. We look at how the precision in their arrival times can be modified by a noble time-delayed stimulation protocol, which we term as "Δt training." The emergence of recurrent bursts and the change of the precision in their arrival times can be explained by the stochastic resonance of a damped, subthreshold, neural oscillation.

Explicit maps to predict activation order in multiphase rhythms of a coupled cell network

We present a novel extension of fast-slow analysis of clustered solutions to coupled networks of three cells, allowing for heterogeneity in the cells' intrinsic dynamics. In the model on which we focus, each cell is described by a pair of first-order differential equations, which are based on recent reduced neuronal network models for respiratory rhythmogenesis. Within each pair of equations, one dependent variable evolves on a fast time scale and one on a slow scale. The cells are coupled with inhibitory synapses that turn on and off on the fast time scale. In this context, we analyze solutions in which cells take turns activating, allowing any activation order, including multiple activations of two of the cells between successive activations of the third. Our analysis proceeds via the derivation of a set of explicit maps between the pairs of slow variables corresponding to the non-active cells on each cycle. We show how these maps can be used to determine the order in which cells will activate for a given initial condition and how evaluation of these maps on a few key curves in their domains can be used to constrain the possible activation orders that will be observed in network solutions. Moreover, under a small set of additional simplifying assumptions, we collapse the collection of maps into a single 2D map that can be computed explicitly. From this unified map, we analytically obtain boundary curves between all regions of initial conditions producing different activation patterns.

Stochastic synchronization of neuronal populations with intrinsic and extrinsic noise

We extend the theory of noise-induced phase synchronization to the case of a neural master equation describing the stochastic dynamics of an ensemble of uncoupled neuronal population oscillators with intrinsic and extrinsic noise. The master equation formulation of stochastic neurodynamics represents the state of each population by the number of currently active neurons, and the state transitions are chosen so that deterministic Wilson-Cowan rate equations are recovered in the mean-field limit. We apply phase reduction and averaging methods to a corresponding Langevin approximation of the master equation in order to determine how intrinsic noise disrupts synchronization of the population oscillators driven by a common extrinsic noise source. We illustrate our analysis by considering one of the simplest networks known to generate limit cycle oscillations at the population level, namely, a pair of mutually coupled excitatory (E) and inhibitory (I) subpopulations. We show how the combination of intrinsic independent noise and extrinsic common noise can lead to clustering of the population oscillators due to the multiplicative nature of both noise sources under the Langevin approximation. Finally, we show how a similar analysis can be carried out for another simple population model that exhibits limit cycle oscillations in the deterministic limit, namely, a recurrent excitatory network with synaptic depression; inclusion of synaptic depression into the neural master equation now generates a stochastic hybrid system.

Quantifying the relative contributions of divisive and subtractive feedback to rhythm generation

Biological systems are characterized by a high number of interacting components. Determining the role of each component is difficult, addressed here in the context of biological oscillations. Rhythmic behavior can result from the interplay of positive feedback that promotes bistability between high and low activity, and slow negative feedback that switches the system between the high and low activity states. Many biological oscillators include two types of negative feedback processes: divisive (decreases the gain of the positive feedback loop) and subtractive (increases the input threshold) that both contribute to slowly move the system between the high- and low-activity states. Can we determine the relative contribution of each type of negative feedback process to the rhythmic activity? Does one dominate? Do they control the active and silent phase equally? To answer these questions we use a neural network model with excitatory coupling, regulated by synaptic depression (divisive) and cellular adaptation (subtractive feedback). We first attempt to apply standard experimental methodologies: either passive observation to correlate the variations of a variable of interest to system behavior, or deletion of a component to establish whether a component is critical for the system. We find that these two strategies can lead to contradictory conclusions, and at best their interpretive power is limited. We instead develop a computational measure of the contribution of a process, by evaluating the sensitivity of the active (high activity) and silent (low activity) phase durations to the time constant of the process. The measure shows that both processes control the active phase, in proportion to their speed and relative weight. However, only the subtractive process plays a major role in setting the duration of the silent phase. This computational method can be used to analyze the role of negative feedback processes in a wide range of biological rhythms.

As modern experimental techniques uncover new components in biological systems and describe their mutual interactions, the problem of determining the contribution of each component becomes critical. The many feedback loops created by these interactions can lead to oscillatory behavior. Examples of oscillations in biology include the cell cycle, circadian rhythms, the electrical activity of excitable cells, and predator-prey systems. While we understand how negative feedback loops can cause oscillations, when multiple feedback loops are present it becomes difficult to identify the dominant mechanism(s), if any. We address the problem of establishing the relative contribution of a feedback process using a biological oscillator model for which oscillations are controlled by two types of slow negative feedback. To determine which is the dominant process, we first use standard experimental methodologies: either passive observation to correlate a variable's behavior to system activity, or deletion of a component to establish whether that component is critical for the system. We find that these methods have limited applicability to the determination of the dominant process. We then develop a new quantitative measure of the contribution of each process to the oscillations. This computational method can be extended to a wide variety of oscillatory systems.

Neurophysiological and computational principles of cortical rhythms in cognition

Synchronous rhythms represent a core mechanism for sculpting temporal coordination of neural activity in the brain-wide network. This review focuses on oscillations in the cerebral cortex that occur during cognition, in alert behaving conditions. Over the last two decades, experimental and modeling work has made great strides in elucidating the detailed cellular and circuit basis of these rhythms, particularly gamma and theta rhythms. The underlying physiological mechanisms are diverse (ranging from resonance and pacemaker properties of single cells to multiple scenarios for population synchronization and wave propagation), but also exhibit unifying principles. A major conceptual advance was the realization that synaptic inhibition plays a fundamental role in rhythmogenesis, either in an interneuronal network or in a reciprocal excitatory-inhibitory loop. Computational functions of synchronous oscillations in cognition are still a matter of debate among systems neuroscientists, in part because the notion of regular oscillation seems to contradict the common observation that spiking discharges of individual neurons in the cortex are highly stochastic and far from being clocklike. However, recent findings have led to a framework that goes beyond the conventional theory of coupled oscillators and reconciles the apparent dichotomy between irregular single neuron activity and field potential oscillations. From this perspective, a plethora of studies will be reviewed on the involvement of long-distance neuronal coherence in cognitive functions such as multisensory integration, working memory, and selective attention. Finally, implications of abnormal neural synchronization are discussed as they relate to mental disorders like schizophrenia and autism.

Mechanism for the universal pattern of activity in developing neuronal networks

Spontaneous episodic activity is a fundamental mode of operation of developing networks. Surprisingly, the duration of an episode of activity correlates with the length of the silent interval that precedes it, but not with the interval that follows. Here we use a modeling approach to explain this characteristic, but thus far unexplained, feature of developing networks. Because the correlation pattern is observed in networks with different structures and components, a satisfactory model needs to generate the right pattern of activity regardless of the details of network architecture or individual cell properties. We thus developed simple models incorporating excitatory coupling between heterogeneous neurons and activity-dependent synaptic depression. These models robustly generated episodic activity with the correct correlation pattern. The correlation pattern resulted from episodes being triggered at random levels of recovery from depression while they terminated around the same level of depression. To explain this fundamental difference between episode onset and termination, we used a mean field model, where only average activity and average level of recovery from synaptic depression are considered. In this model, episode onset is highly sensitive to inputs. Thus noise resulting from random coincidences in the spike times of individual neurons led to the high variability at episode onset and to the observed correlation pattern. This work further shows that networks with widely different architectures, different cell types, and different functions all operate according to the same general mechanism early in their development.

Properties and mechanisms of spontaneous activity in the embryonic chick hindbrain

Spontaneous activity regulates many aspects of central nervous system development. We demonstrate that in the embryonic chick hindbrain, spontaneous activity is expressed between embryonic days (E) 6-9. Over this period the frequency of activity decreases significantly, although the events maintain a consistent rhythm on the timescale of minutes. At E6, the activity is pharmacologically dependent on serotonin, nACh, GABA(A), and glycine input, but not on muscarinic, glutamatergic, or GABA(B) receptor activation. It also depends on gap junctions, t-type calcium channels and TTX-sensitive ion channels. In intact spinal cord-hindbrain preparations, E6 spontaneous events originate in the spinal cord and propagate into lateral hindbrain tissue; midline activity follows the appearance of lateral activity. However, the spinal cord is not required for hindbrain activity. There are two invariant points of origin of activity along the midline, both within the caudal group of serotonin-expressing cell bodies; one point is caudal to the nV exit point while the other is caudal to the nVII exit point. Additional caudal midline points of origin are seen in a minority of cases. Using immunohistochemistry, we show robust differentiation of the serotonergic raphe near the midline at E6, and extensive fiber tracts expressing GAD65/67 and the nAChR in lateral areas; this suggests that the medial activity is dependent on serotonergic neuron activation, while lateral activity requires other transmitters. Although there are differences between species, this activity is highly conserved between mouse and chick, suggesting that developmental event(s) within the hindbrain are dependent on expression of this spontaneous activity.

Noise-induced transitions in slow wave neuronal dynamics

Many neuronal systems exhibit slow random alternations and sudden switches in activity states. Models with noisy relaxation dynamics (oscillatory, excitable or bistable) account for these temporal, slow wave, patterns and the fluctuations within states. The noise-induced transitions in a relaxation dynamics are analogous to escape by a particle in a slowly changing double-well potential. In this formalism, we obtain semi-analytically the first and second order statistical properties: the distributions of the slow process at the transitions and the temporal correlations of successive switching events. We find that the temporal correlations can be used to help distinguish among biophysical mechanisms for the slow negative feedback, such as divisive or subtractive. We develop our results in the context of models for cellular pacemaker neurons; they also apply to mean-field models for spontaneously active networks with slow wave dynamics.

Control of oscillation periods and phase durations in half-center central pattern generators: a comparative mechanistic analysis

Central pattern generators (CPGs) consisting of interacting groups of neurons drive a variety of repetitive, rhythmic behaviors in invertebrates and vertebrates, such as arise in locomotion, respiration, mastication, scratching, and so on. These CPGs are able to generate rhythmic activity in the absence of afferent feedback or rhythmic inputs. However, functionally relevant CPGs must adaptively respond to changing demands, manifested as changes in oscillation period or in relative phase durations in response to variations in non-patterned inputs or drives. Although many half-center CPG models, composed of symmetric units linked by reciprocal inhibition yet varying in their intrinsic cellular properties, have been proposed, the precise oscillatory mechanisms operating in most biological CPGs remain unknown. Using numerical simulations and phase-plane analysis, we comparatively investigated how the intrinsic cellular features incorporated in different CPG models, such as subthreshold activation based on a slowly inactivating persistent sodium current, adaptation based on slowly activating calcium-dependent potassium current, or post-inhibitory rebound excitation, can contribute to the control of oscillation period and phase durations in response to changes in excitatory external drive to one or both half-centers. Our analysis shows that both the sensitivity of oscillation period to alterations of excitatory drive and the degree to which the duration of each phase can be separately controlled depend strongly on the intrinsic cellular mechanisms involved in rhythm generation and phase transitions. In particular, the CPG formed from units incorporating a slowly inactivating persistent sodium current shows the greatest range of oscillation periods and the greatest degree of independence in phase duration control by asymmetric inputs. These results are explained based on geometric analysis of the phase plane structures corresponding to the dynamics for each CPG type, which in particular helps pinpoint the roles of escape and release from synaptic inhibition in the effects we find.

Correlation analysis a tool for comparing relaxation-type models to experimental data

We describe a new technique for comparing mathematical models to the biological systems that are described. This technique is appropriate for systems that produce relaxation oscillations or bursting oscillations, and takes advantage of noise that is inherent to all biological systems. Both types of oscillations are composed of active phases of activity followed by silent phases, repeating periodically. The presence of noise adds variability to the durations of the different phases. The central idea of the technique is that the active phase duration may be correlated with either/both the previous or next silent phase duration, and the resulting correlation pattern provides information about the dynamic structure of the system. Correlation patterns can easily be determined by making scatter plots and applying correlation analysis to the cluster of data points. This could be done both with experimental data and with model simulation data. If the model correlation pattern is in general agreement with the experimental data, then this adds support for the validity of the model. Otherwise, the model must be corrected. While this tool is only one test of many required to validate a mathematical model, it is easy to implement and is noninvasive.

Multi-stability and pattern-selection in oscillatory networks with fast inhibition and electrical synapses

A model or hybrid network consisting of oscillatory cells interconnected by inhibitory and electrical synapses may express different stable activity patterns without any change of network topology or parameters, and switching between the patterns can be induced by specific transient signals. However, little is known of properties of such signals. In the present study, we employ numerical simulations of neural networks of different size composed of relaxation oscillators, to investigate switching between in-phase (IP) and anti-phase (AP) activity patterns. We show that the time windows of susceptibility to switching between the patterns are similar in 2-, 4- and 6-cell fully-connected networks. Moreover, in a network (N = 4, 6) expressing a given AP pattern, a stimulus with a given profile consisting of depolarizing and hyperpolarizing signals sent to different subpopulations of cells can evoke switching to another AP pattern. Interestingly, the resulting pattern encodes the profile of the switching stimulus. These results can be extended to different network architectures. Indeed, relaxation oscillators are not only models of cellular pacemakers, bursting or spiking, but are also analogous to firing-rate models of neural activity. We show that rules of switching similar to those found for relaxation oscillators apply to oscillating circuits of excitatory cells interconnected by electrical synapses and cross-inhibition. Our results suggest that incoming information, arriving in a proper time window, may be stored in an oscillatory network in the form of a specific spatio-temporal activity pattern which is expressed until new pertinent information arrives.

Fluctuation-driven rhythmogenesis in an excitatory neuronal network with slow adaptation

We study an excitatory all-to-all coupled network of N spiking neurons with synaptically filtered background noise and slow activity-dependent hyperpolarization currents. Such a system exhibits noise-induced burst oscillations over a range of values of the noise strength (variance) and level of cell excitability. Since both of these quantities depend on the rate of background synaptic inputs, we show how noise can provide a mechanism for increasing the robustness of rhythmic bursting and the range of burst frequencies. By exploiting a separation of time scales we also show how the system dynamics can be reduced to low-dimensional mean field equations in the limit N --> infinity. Analysis of the bifurcation structure of the mean field equations provides insights into the dynamical mechanisms for initiating and terminating the bursts.

Dynamical characteristics common to neuronal competition models

Models implementing neuronal competition by reciprocally inhibitory populations are widely used to characterize bistable phenomena such as binocular rivalry. We find common dynamical behavior in several models of this general type, which differ in their architecture in the form of their gain functions, and in how they implement the slow process that underlies alternating dominance. We focus on examining the effect of the input strength on the rate (and existence) of oscillations. In spite of their differences, all considered models possess similar qualitative features, some of which we report here for the first time. Experimentally, dominance durations have been reported to decrease monotonically with increasing stimulus strength (such as Levelt's "Proposition IV"). The models predict this behavior; however, they also predict that at a lower range of input strength dominance durations increase with increasing stimulus strength. The nonmonotonic dependency of duration on stimulus strength is common to both deterministic and stochastic models. We conclude that additional experimental tests of Levelt's Proposition IV are needed to reconcile models and perception.