Journal of Mathematical Neuroscience最新文献

In 1972-1973 Wilson and Cowan introduced a mathematical model of the population dynamics of synaptically coupled excitatory and inhibitory neurons in the neocortex. The model dealt only with the mean numbers of activated and quiescent excitatory and inhibitory neurons, and said nothing about fluctuations and correlations of such activity. However, in 1997 Ohira and Cowan, and then in 2007-2009 Buice and Cowan introduced Markov models of such activity that included fluctuation and correlation effects. Here we show how both models can be used to provide a quantitative account of the population dynamics of neocortical activity.We first describe how the Markov models account for many recent measurements of the resting or spontaneous activity of the neocortex. In particular we show that the power spectrum of large-scale neocortical activity has a Brownian motion baseline, and that the statistical structure of the random bursts of spiking activity found near the resting state indicates that such a state can be represented as a percolation process on a random graph, called directed percolation.Other data indicate that resting cortex exhibits pair correlations between neighboring populations of cells, the amplitudes of which decay slowly with distance, whereas stimulated cortex exhibits pair correlations which decay rapidly with distance. Here we show how the Markov model can account for the behavior of the pair correlations.Finally we show how the 1972-1973 Wilson-Cowan equations can account for recent data which indicates that there are at least two distinct modes of cortical responses to stimuli. In mode 1 a low intensity stimulus triggers a wave that propagates at a velocity of about 0.3 m/s, with an amplitude that decays exponentially. In mode 2 a high intensity stimulus triggers a larger response that remains local and does not propagate to neighboring regions.

The tools of weakly coupled phase oscillator theory have had a profound impact on the neuroscience community, providing insight into a variety of network behaviours ranging from central pattern generation to synchronisation, as well as predicting novel network states such as chimeras. However, there are many instances where this theory is expected to break down, say in the presence of strong coupling, or must be carefully interpreted, as in the presence of stochastic forcing. There are also surprises in the dynamical complexity of the attractors that can robustly appear-for example, heteroclinic network attractors. In this review we present a set of mathematical tools that are suitable for addressing the dynamics of oscillatory neural networks, broadening from a standard phase oscillator perspective to provide a practical framework for further successful applications of mathematics to understanding network dynamics in neuroscience.

Sensory input to the lamprey central pattern generator (CPG) for locomotion is known to have a significant role in modulating lamprey swimming. Lamprey CPGs are known to have the ability to entrain to a bending stimulus, that is, in the presence of a rhythmic signal, the CPG will change its frequency to match the stimulus frequency. Bending experiments in which the lamprey spinal cord has been removed and mechanically bent back and forth at a single point have been used to determine the range of frequencies that can entrain the CPG rhythm. First, we model the lamprey locomotor CPG as a chain of neural oscillators with three classes of neurons and sinusoidal forcing representing edge cell input. We derive a phase model using the connections described in the neural model. This results in a simpler model yet maintains some properties of the neural model. For both the neural model and the derived phase model, entrainment ranges are computed for forcing at different points along the chain while varying both intersegmental coupling strength and the coupling strength between the forcer and chain. Entrainment ranges for chains with nonuniform intersegmental coupling asymmetry are larger when forcing is applied to the middle of the chain than when it is applied to either end, a result that is qualitatively similar to the experimental results. In the limit of weak coupling in the chain, the entrainment results of the neural model approach the entrainment results for the derived phase model. Both biological experiments and the robustness of non-monotonic entrainment ranges as a function of the forcing position across different classes of CPG models with nonuniform asymmetric coupling suggest that a specific property of the intersegmental coupling of the CPG is key to entrainment.

We investigate the dynamics of unidirectional semi-infinite chains of type-I oscillators that are periodically forced at their root node, as an archetype of wave generation in neural networks. In previous studies, numerical simulations based on uniform forcing have revealed that trajectories approach a traveling wave in the far-downstream, large time limit. While this phenomenon seems typical, it is hardly anticipated because the system does not exhibit any of the crucial properties employed in available proofs of existence of traveling waves in lattice dynamical systems. Here, we give a full mathematical proof of generation under uniform forcing in a simple piecewise affine setting for which the dynamics can be solved explicitly. In particular, our analysis proves existence, global stability, and robustness with respect to perturbations of the forcing, of families of waves with arbitrary period/wave number in some range, for every value of the parameters in the system.

A fundamental question concerning the way the visual world is represented in our brain is how a cortical cell responds when its classical receptive field contains a plurality of stimuli. Two opposing models have been proposed. In the response-averaging model, the neuron responds with a weighted average of all individual stimuli. By contrast, in the probability-mixing model, the cell responds to a plurality of stimuli as if only one of the stimuli were present. Here we apply the probability-mixing and the response-averaging model to leaky integrate-and-fire neurons, to describe neuronal behavior based on observed spike trains. We first estimate the parameters of either model using numerical methods, and then test which model is most likely to have generated the observed data. Results show that the parameters can be successfully estimated and the two models are distinguishable using model selection.

The original neural field model of Wilson and Cowan is often interpreted as the averaged behaviour of a network of switch like neural elements with a distribution of switch thresholds, giving rise to the classic sigmoidal population firing-rate function so prevalent in large scale neuronal modelling. In this paper we explore the effects of such threshold noise without recourse to averaging and show that spatial correlations can have a strong effect on the behaviour of waves and patterns in continuum models. Moreover, for a prescribed spatial covariance function we explore the differences in behaviour that can emerge when the underlying stationary distribution is changed from Gaussian to non-Gaussian. For travelling front solutions, in a system with exponentially decaying spatial interactions, we make use of an interface approach to calculate the instantaneous wave speed analytically as a series expansion in the noise strength. From this we find that, for weak noise, the spatially averaged speed depends only on the choice of covariance function and not on the shape of the stationary distribution. For a system with a Mexican-hat spatial connectivity we further find that noise can induce localised bump solutions, and using an interface stability argument show that there can be multiple stable solution branches.

We study a population of spiking neurons which are subject to independent noise processes and a strong common time-dependent input. We show that the response of output spikes to independent noise shapes information transmission of such populations even when information transmission properties of single neurons are left unchanged. In particular, we consider two Poisson models in which independent noise either (i) adds and deletes spikes (AD model) or (ii) shifts spike times (STS model). We show that in both models suprathreshold stochastic resonance (SSR) can be observed, where the information transmitted by a neural population is increased with addition of independent noise. In the AD model, the presence of the SSR effect is robust and independent of the population size or the noise spectral statistics. In the STS model, the information transmission properties of the population are determined by the spectral statistics of the noise, leading to a strongly increased effect of SSR in some regimes, or an absence of SSR in others. Furthermore, we observe a high-pass filtering of information in the STS model that is absent in the AD model. We quantify information transmission by means of the lower bound on the mutual information rate and the spectral coherence function. To this end, we derive the signal-output cross-spectrum, the output power spectrum, and the cross-spectrum of two spike trains for both models analytically.

We investigate the propagation of probabilistic uncertainty through the action potential mechanism in nerve cells. Using the Hodgkin-Huxley (H-H) model and Stochastic Collocation on Sparse Grids, we obtain an accurate probabilistic interpretation of the deterministic dynamics of the transmembrane potential and gating variables. Using Sobol indices, out of the 11 uncertain parameters in the H-H model, we unravel two main uncertainty sources, which account for more than 90 % of the fluctuations in neuronal responses, and have a direct biophysical interpretation. We discuss how this interesting feature of the H-H model allows one to reduce greatly the probabilistic degrees of freedom in uncertainty quantification analyses, saving CPU time in numerical simulations and opening possibilities for probabilistic generalisation of other deterministic models of great importance in physiology and mathematical neuroscience.

Simple-spike synchrony between Purkinje cells projecting to a common neuron in the deep cerebellar nucleus is emerging as an important factor in the encoding of output information from cerebellar cortex. A phenomenon known as stochastic synchronization happens when uncoupled oscillators synchronize due to correlated inputs. Stochastic synchronization is a viable mechanism through which simple-spike synchrony could be generated, but it has received scarce attention, perhaps because the presence of feedforward inhibition in the input to Purkinje cells makes insights difficult. This paper presents a method to account for feedforward inhibition so the usual mathematical approaches to stochastic synchronization can be applied. The method consists in finding a single Phase Response Curve, called the equivalent PRC, that accounts for the effects of both excitatory inputs and delayed feedforward inhibition from molecular layer interneurons. The results suggest that a theory of stochastic synchronization for the case of feedforward inhibition may not be necessary, since this case can be approximately reduced to the case of inputs characterized by a single PRC. Moreover, feedforward inhibition could in many situations increase the level of synchrony experienced by Purkinje cells.