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 → ∞. Analysis of the bifurcation structure of the mean field equations provides insights into the dynamical mechanisms for initiating and terminating the bursts

Membrane Properties and the Balance between Excitation and Inhibition Control Gamma-Frequency Oscillations Arising from Feedback Inhibition

Computational studies as well as in vivo and in vitro results have shown that many cortical neurons fire in a highly irregular manner and at low average firing rates. These patterns seem to persist even when highly rhythmic signals are recorded by local field potential electrodes or other methods that quantify the summed behavior of a local population. Models of the 30–80 Hz gamma rhythm in which network oscillations arise through ‘stochastic synchrony’ capture the variability observed in the spike output of single cells while preserving network-level organization. We extend upon these results by constructing model networks constrained by experimental measurements and using them to probe the effect of biophysical parameters on network-level activity. We find in simulations that gamma-frequency oscillations are enabled by a high level of incoherent synaptic conductance input, similar to the barrage of noisy synaptic input that cortical neurons have been shown to receive in vivo. This incoherent synaptic input increases the emergent network frequency by shortening the time scale of the membrane in excitatory neurons and by reducing the temporal separation between excitation and inhibition due to decreased spike latency in inhibitory neurons. These mechanisms are demonstrated in simulations and in vitro current-clamp and dynamic-clamp experiments. Simulation results further indicate that the membrane potential noise amplitude has a large impact on network frequency and that the balance between excitatory and inhibitory currents controls network stability and sensitivity to external inputs

The Dynamic Brain in Action: Cortical Oscillations and Coordination Dynamics

Cortical oscillations are electrical activities with rhythmic and/or repetitive nature generated spontaneously and in response to stimuli. Study of cortical oscillations has become an area of converging interests since the last two decades and has deepened our understanding of its physiological basis across different behavioral states. Experimental and modeling work has taught us that there is a wide diversity of cellular and circuit mechanisms underlying the generation of cortical rhythms. A wildly diverse set of functions has pertained to synchronous oscillations but their significance in cognition should be better appraised in the more general framework of correlation between spike times of neurons. Oscillations are the core mechanism in adjusting neuronal interactions and shaping temporal coordination of neural activity. In the first part of this thesis, we review essential feature of cortical oscillations in membrane potentials and local field potentials recorded from turtle ex vivo preparation. Then we develop a simple computational model that reproduces the observed features. This modeling investigation suggests a plausible underlying mechanism for rhythmogenesis through cellular and circuit properties. The second part of the thesis is about temporal coordination dynamics quantified by signal and noise correlations. Here, again, we present a computational model to show how temporal coordination and synchronous oscillations can be sewn together. More importantly, what biophysical ingrediants are necessary for a network to reproduce the observed coordination dynamics

One-Dimensional Population Density Approaches to Recurrently Coupled Networks of Neurons with Noise

Mean-field systems have been previously derived for networks of coupled, two-dimensional, integrate-and-fire neurons such as the Izhikevich, adapting exponential (AdEx) and quartic integrate and fire (QIF), among others. Unfortunately, the mean-field systems have a degree of frequency error and the networks analyzed often do not include noise when there is adaptation. Here, we derive a one-dimensional partial differential equation (PDE) approximation for the marginal voltage density under a first order moment closure for coupled networks of integrate-and-fire neurons with white noise inputs. The PDE has substantially less frequency error than the mean-field system, and provides a great deal more information, at the cost of analytical tractability. The convergence properties of the mean-field system in the low noise limit are elucidated. A novel method for the analysis of the stability of the asynchronous tonic firing solution is also presented and implemented. Unlike previous attempts at stability analysis with these network types, information about the marginal densities of the adaptation variables is used. This method can in principle be applied to other systems with nonlinear partial differential equations.Comment: 26 Pages, 6 Figure

Emergent bursting and synchrony in computer simulations of neuronal cultures

Experimental studies of neuronal cultures have revealed a wide variety of spiking network activity ranging from sparse, asynchronous firing to distinct, network-wide synchronous bursting. However, the functional mechanisms driving these observed firing patterns are not well understood. In this work, we develop an in silico network of cortical neurons based on known features of similar in vitro networks. The activity from these simulations is found to closely mimic experimental data. Furthermore, the strength or degree of network bursting is found to depend on a few parameters: the density of the culture, the type of synaptic connections, and the ratio of excitatory to inhibitory connections. Network bursting gradually becomes more prominent as either the density, the fraction of long range connections, or the fraction of excitatory neurons is increased. Interestingly, biologically prevalent values of parameters result in networks that are at the transition between strong bursting and sparse firing. Using principal components analysis, we show that a large fraction of the variance in firing rates is captured by the first component for bursting networks. These results have implications for understanding how information is encoded at the population level as well as for why certain network parameters are ubiquitous in cortical tissue

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

An astrocyte-dependent mechanism for neuronal rhythmogenesis

Communication between neurons rests on their capacity to change their firing pattern to encode different messages. For several vital functions, such as respiration and mastication, neurons need to generate a rhythmic firing pattern. Here we show in the rat trigeminal sensori-motor circuit for mastication that this ability depends on regulation of the extracellular Ca2+ concentration ([Ca2+]e) by astrocytes. In this circuit, astrocytes respond to sensory stimuli that induce neuronal rhythmic activity, and their blockade with a Ca2+ chelator prevents neurons from generating a rhythmic bursting pattern. This ability is restored by adding S100b, an astrocytic Ca2+-binding protein, to the extracellular space, while application of an anti-S100b antibody prevents generation of rhythmic activity. These results indicate that astrocytes regulate a fundamental neuronal property: the capacity to change firing pattern. These findings may have broad implications for many other neural networks whose functions depend on the generation of rhythmic activity

Noise-induced behaviors in neural mean field dynamics

The collective behavior of cortical neurons is strongly affected by the presence of noise at the level of individual cells. In order to study these phenomena in large-scale assemblies of neurons, we consider networks of firing-rate neurons with linear intrinsic dynamics and nonlinear coupling, belonging to a few types of cell populations and receiving noisy currents. Asymptotic equations as the number of neurons tends to infinity (mean field equations) are rigorously derived based on a probabilistic approach. These equations are implicit on the probability distribution of the solutions which generally makes their direct analysis difficult. However, in our case, the solutions are Gaussian, and their moments satisfy a closed system of nonlinear ordinary differential equations (ODEs), which are much easier to study than the original stochastic network equations, and the statistics of the empirical process uniformly converge towards the solutions of these ODEs. Based on this description, we analytically and numerically study the influence of noise on the collective behaviors, and compare these asymptotic regimes to simulations of the network. We observe that the mean field equations provide an accurate description of the solutions of the network equations for network sizes as small as a few hundreds of neurons. In particular, we observe that the level of noise in the system qualitatively modifies its collective behavior, producing for instance synchronized oscillations of the whole network, desynchronization of oscillating regimes, and stabilization or destabilization of stationary solutions. These results shed a new light on the role of noise in shaping collective dynamics of neurons, and gives us clues for understanding similar phenomena observed in biological networks

Mean-field equations for stochastic firing-rate neural fields with delays: Derivation and noise-induced transitions

In this manuscript we analyze the collective behavior of mean-field limits of large-scale, spatially extended stochastic neuronal networks with delays. Rigorously, the asymptotic regime of such systems is characterized by a very intricate stochastic delayed integro-differential McKean-Vlasov equation that remain impenetrable, leaving the stochastic collective dynamics of such networks poorly understood. In order to study these macroscopic dynamics, we analyze networks of firing-rate neurons, i.e. with linear intrinsic dynamics and sigmoidal interactions. In that case, we prove that the solution of the mean-field equation is Gaussian, hence characterized by its two first moments, and that these two quantities satisfy a set of coupled delayed integro-differential equations. These equations are similar to usual neural field equations, and incorporate noise levels as a parameter, allowing analysis of noise-induced transitions. We identify through bifurcation analysis several qualitative transitions due to noise in the mean-field limit. In particular, stabilization of spatially homogeneous solutions, synchronized oscillations, bumps, chaotic dynamics, wave or bump splitting are exhibited and arise from static or dynamic Turing-Hopf bifurcations. These surprising phenomena allow further exploring the role of noise in the nervous system.Comment: Updated to the latest version published, and clarified the dependence in space of Brownian motion