Class-II neurons display a higher degree of stochastic synchronization than class-I neurons

The Effects of Background Noise on a Biophysical Model of Olfactory Bulb Mitral Cells

The spiking activity of mitral cells (MC) in the olfactory bulb is a key attribute in olfactory sensory information processing to downstream cortical areas. A more detailed understanding of the modulation of MC spike statistics could shed light on mechanistic studies of olfactory bulb circuits and olfactory coding. We study the spike response of a recently developed single-compartment biophysical MC model containing seven known ionic currents and calcium dynamics subject to constant current input with background white noise. We observe rich spiking dynamics even with constant current input, including multimodal peaks in the interspike interval distribution (ISI). Although weak-to-moderate background noise for a fixed current input does not change the firing rate much, the spike dynamics can change dramatically, exhibiting non-monotonic spike variability not commonly observed in standard neuron models. We explain these dynamics with a phenomenological model of the ISI probability density function. Our study clarifies some of the complexities of MC spiking dynamics.

On the interpretation of Dirac δ pulses in differential equations for phase oscillators

In this note, we discuss the usage of the Dirac δ function in models of phase oscillators with pulsatile inputs. Many authors use a product of the delta function and the phase response curve in the right-hand side of an ordinary differential equation to describe the discontinuous phase dynamics in such systems. We point out that this notation has to be treated with care as it is ambiguous. We argue that the presumably most canonical interpretation does not lead to the intended behavior in many cases.

Shunting Inhibition Improves Synchronization in Heterogeneous Inhibitory Interneuronal Networks with Type 1 Excitability Whereas Hyperpolarizing Inhibition Is Better for Type 2 Excitability

All-to-all homogeneous networks of inhibitory neurons synchronize completely under the right conditions; however, many modeling studies have shown that biological levels of heterogeneity disrupt synchrony. Our fundamental scientific question is “how can neurons maintain partial synchrony in the presence of heterogeneity and noise?” A particular subset of strongly interconnected interneurons, the PV+ fast-spiking (FS) basket neurons, are strongly implicated in γ oscillations and in phase locking of nested γ oscillations to theta. Their excitability type apparently varies between brain regions: in CA1 and the dentate gyrus they have type 1 excitability, meaning that they can fire arbitrarily slowly, whereas in the striatum and cortex they have type 2 excitability, meaning that there is a frequency thresh old below which they cannot sustain repetitive firing. We constrained the models to study the effect of excitability type (more precisely bifurcation type) in isolation from all other factors. We use sparsely connected, heterogeneous, noisy networks with synaptic delays to show that synchronization properties, namely the resistance to suppression and the strength of theta phase to γ amplitude coupling, are strongly dependent on the pairing of excitability type with the type of inhibition. Shunting inhibition performs better for type 1 and hyperpolarizing inhibition for type 2. γ Oscillations and their nesting within theta oscillations are thought to subserve cognitive functions like memory encoding and recall; therefore, it is important to understand the contribution of intrinsic properties to these rhythms.

Synchronization of stochastic mean field networks of Hodgkin-Huxley neurons with noisy channels

In this work we are interested in a mathematical model of the collective behavior of a fully connected network of finitely many neurons, when their number and when time go to infinity. We assume that every neuron follows a stochastic version of the Hodgkin-Huxley model, and that pairs of neurons interact through both electrical and chemical synapses, the global connectivity being of mean field type. When the leak conductance is strictly positive, we prove that if the initial voltages are uniformly bounded and the electrical interaction between neurons is strong enough, then, uniformly in the number of neurons, the whole system synchronizes exponentially fast as time goes to infinity, up to some error controlled by (and vanishing with) the channels noise level. Moreover, we prove that if the random initial condition is exchangeable, on every bounded time interval the propagation of chaos property for this system holds (regardless of the interaction intensities). Combining these results, we deduce that the nonlinear McKean-Vlasov equation describing an infinite network of such neurons concentrates, as time goes to infinity, around the dynamics of a single Hodgkin-Huxley neuron with chemical neurotransmitter channels. Our results are illustrated and complemented with numerical simulations.

A theoretical framework for analyzing coupled neuronal networks: Application to the olfactory system

Determining how synaptic coupling within and between regions is modulated during sensory processing is an important topic in neuroscience. Electrophysiological recordings provide detailed information about neural spiking but have traditionally been confined to a particular region or layer of cortex. Here we develop new theoretical methods to study interactions between and within two brain regions, based on experimental measurements of spiking activity simultaneously recorded from the two regions. By systematically comparing experimentally-obtained spiking statistics to (efficiently computed) model spike rate statistics, we identify regions in model parameter space that are consistent with the experimental data. We apply our new technique to dual micro-electrode array in vivo recordings from two distinct regions: olfactory bulb (OB) and anterior piriform cortex (PC). Our analysis predicts that: i) inhibition within the afferent region (OB) has to be weaker than the inhibition within PC, ii) excitation from PC to OB is generally stronger than excitation from OB to PC, iii) excitation from PC to OB and inhibition within PC have to both be relatively strong compared to presynaptic inputs from OB. These predictions are validated in a spiking neural network model of the OB–PC pathway that satisfies the many constraints from our experimental data. We find when the derived relationships are violated, the spiking statistics no longer satisfy the constraints from the data. In principle this modeling framework can be adapted to other systems and be used to investigate relationships between other neural attributes besides network connection strengths. Thus, this work can serve as a guide to further investigations into the relationships of various neural attributes within and across different regions during sensory processing.

Sensory processing is known to span multiple regions of the nervous system. However, electrophysiological recordings during sensory processing have traditionally been limited to a single region or brain layer. With recent advances in experimental techniques, recorded spiking activity from multiple regions simultaneously is feasible. However, other important quantities— such as inter-region connection strengths—cannot yet be measured. Here, we develop new theoretical tools to leverage data obtained by recording from two different brain regions simultaneously. We address the following questions: what are the crucial neural network attributes that enable sensory processing across different regions, and how are these attributes related to one another? With a novel theoretical framework to efficiently calculate spiking statistics, we can characterize a high dimensional parameter space that satisfies data constraints. We apply our results to the olfactory system to make specific predictions about effective network connectivity. Our framework relies on incorporating relatively easy-to-measure quantities to predict hard-to-measure interactions across multiple brain regions. Because this work is adaptable to other systems, we anticipate it will be a valuable tool for analysis of other larger scale brain recordings.

Measurement of phase resetting curves using optogenetic barrage stimuli

BACKGROUND

The phase resetting curve (PRC) is a primary measure of a rhythmically firing neuron's responses to synaptic input, quantifying the change in phase of the firing oscillation as a function of the input phase. PRCs provide information about whether neurons will synchronize due to synaptic coupling or shared input. However, PRC estimation has been limited to in vitro preparations where stable intracellular recordings can be obtained and background activity is minimal, and new methods are required for in vivo applications.

NEW METHOD

We estimated PRCs using dense optogenetic stimuli and extracellular spike recording. Autonomously firing neurons in substantia nigra pars reticulata (SNr) of Thy1-channelrhodopsin 2 (ChR2) transgenic mice were stimulated with random barrages of light pulses, and PRCs were determined using multiple linear regression.

RESULTS

The PRCs obtained were type-I, showing only phase advances in response to depolarizing input, and generally sloped upward from early to late phases. Secondary PRCs, indicating the effect on the subsequent ISI, showed phase delays primarily for stimuli arriving at late phases. Phase models constructed from the optogenetic PRCs accounted for a large fraction of the variance in ISI length and provided a good approximation of the spike-triggered average stimulus.

COMPARISON WITH EXISTING METHODS

Compared to methods based on intracellular current injection, the new method sacrifices some temporal resolution. However, it should be much more widely applicable in vivo, because only extracellular recording and optogenetic stimulation are required.

CONCLUSIONS

These results demonstrate PRC estimation using methods suitable for in vivo applications.

Stimulus Phase Locking of Cortical Oscillations for Rhythmic Tone Sequences in Rats

Humans can rapidly detect regular patterns (i.e., within few cycles) without any special attention to the acoustic environment. This suggests that human sensory systems are equipped with a powerful mechanism for automatically predicting forthcoming stimuli to detect regularity. It has recently been hypothesized that the neural basis of sensory predictions exists for not only what happens (predictive coding) but also when a particular stimulus occurs (predictive timing). Here, we hypothesize that the phases of neural oscillations are critical in predictive timing, and these oscillations are modulated in a band-specific manner when acoustic patterns become predictable, i.e., regular. A high-density microelectrode array (10 × 10 within 4 × 4 mm²) was used to characterize spatial patterns of band-specific oscillations when a random-tone sequence was switched to a regular-tone sequence. Increasing the regularity of the tone sequence enhanced phase locking in a band-specific manner, notwithstanding the type of the regular sound pattern. Gamma-band phase locking increased immediately after the transition from random to regular sequences, while beta-band phase locking gradually evolved with time after the transition. The amplitude of the tone-evoked response, in contrast, increased with frequency separation with respect to the prior tone, suggesting that the evoked-response amplitude encodes sequence information on a local scale, i.e., the local order of tones. The phase locking modulation spread widely over the auditory cortex, while the amplitude modulation was confined around the activation foci. Thus, our data suggest that oscillatory phase plays a more important role than amplitude in the neuronal detection of tone sequence regularity, which is closely related to predictive timing. Furthermore, band-specific contributions may support recent theories that gamma oscillations encode bottom-up prediction errors, whereas beta oscillations are involved in top-down prediction.

Coding Properties of Three Intrinsically Distinct Retinal Ganglion Cells under Periodic Stimuli: A Computational Study

As the sole output neurons in the retina, ganglion cells play significant roles in transforming visual information into spike trains, and then transmitting them to the higher visual centers. However, coding strategies that retinal ganglion cells (RGCs) adopt to accomplish these processes are not completely clear yet. To clarify these issues, we investigate the coding properties of three types of RGCs (repetitive spiking, tonic firing, and phasic firing) by two different measures (spike-rate and spike-latency). Model results show that for periodic stimuli, repetitive spiking RGC and tonic RGC exhibit similar spike-rate patterns. Their spike- rates decrease gradually with increased stimulus frequency, moreover, variation of stimulus amplitude would change the two RGCs' spike-rate patterns. For phasic RGC, it activates strongly at medium levels of frequency when the stimulus amplitude is low. While if high stimulus amplitude is applied, phasic RGC switches to respond strongly at low frequencies. These results suggest that stimulus amplitude is a prominent factor in regulating RGCs in encoding periodic signals. Similar conclusions can be drawn when analyzes spike-latency patterns of the three RGCs. More importantly, the above phenomena can be accurately reproduced by Hodgkin's three classes of neurons, indicating that RGCs can perform the typical three classes of firing dynamics, depending on the distinctions of ion channel densities. Consequently, model results from the three RGCs may be not specific, but can also applicable to neurons in other brain regions which exhibit part(s) or all of the Hodgkin's three excitabilities.

Interplay between excitability type and distributions of neuronal connectivity determines neuronal network synchronization

While the interplay between neuronal excitability properties and global properties of network topology is known to affect network propensity for synchronization, it is not clear how detailed characteristics of these properties affect spatiotemporal pattern formation. Here we study mixed networks, composed of neurons having type I and/or type II phase response curves, with varying distributions of local and random connections and show that not only average network properties, but also the connectivity distribution statistics, significantly affect network synchrony. Namely, we study networks with fixed networkwide properties, but vary the number of random connections that nodes project. We show that varying node excitability (type I vs type II) influences network synchrony most dramatically for systems with long-tailed distributions of the number of random connections per node. This indicates that a cluster of even a few highly rewired cells with a high propensity for synchronization can alter the degree of synchrony in the network as a whole. We show this effect generally on a network of coupled Kuramoto oscillators and investigate the impact of this effect more thoroughly in pulse-coupled networks of biophysical neurons.

Optimization of noise-induced synchronization of oscillator networks

We investigate common-noise-induced synchronization between two identical networks of coupled phase oscillators exhibiting fully locked collective oscillations. Using the collective phase description method for fully locked oscillators, we demonstrate that two noninteracting networks of coupled phase oscillators can exhibit in-phase synchronization between the networks when driven by weak common noise. We derive the Lyapunov exponent characterizing the relaxation time for synchronization and develop a method of obtaining the optimal input pattern of common noise to achieve fast synchronization. We illustrate the theory using three representative networks with heterogeneous, global, and local coupling. The theoretical results are validated by direct numerical simulations.

Criterion for noise-induced synchronization: Application to colloidal alignment

Colloidal bodies of irregular shape rotate as they descend under gravity in solution. This rotational response provides a means of bringing a dispersion of identical bodies into a synchronized rotation with the same orientation using programed forcing. We use the notion of statistical entropy to derive bounds on the rate of synchronization. These bounds apply generally to dynamical systems with stable periodic motion with a phase ϕ(t), when subjected to an impulsive perturbation. The impulse causes a change of phase expressible as a phase map ψ(ϕ). We derive an upper limit on the average change of entropy 〈ΔH〉 in terms of this phase map; when this limit is negative, alignment must occur. For systems that have achieved a low entropy, the 〈ΔH〉 approaches this upper limit.

Stochastic Synchronization in Purkinje Cells with Feedforward Inhibition Could Be Studied with Equivalent Phase-Response Curves

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.

Synchronization properties of heterogeneous neuronal networks with mixed excitability type

We study the synchronization of neuronal networks with dynamical heterogeneity, showing that network structures with the same propensity for synchronization (as quantified by master stability function analysis) may develop dramatically different synchronization properties when heterogeneity is introduced with respect to neuronal excitability type. Specifically, we investigate networks composed of neurons with different types of phase response curves (PRCs), which characterize how oscillating neurons respond to excitatory perturbations. Neurons exhibiting type 1 PRC respond exclusively with phase advances, while neurons exhibiting type 2 PRC respond with either phase delays or phase advances, depending on when the perturbation occurs. We find that Watts-Strogatz small world networks transition to synchronization gradually as the proportion of type 2 neurons increases, whereas scale-free networks may transition gradually or rapidly, depending upon local correlations between node degree and excitability type. Random placement of type 2 neurons results in gradual transition to synchronization, whereas placement of type 2 neurons as hubs leads to a much more rapid transition, showing that type 2 hub cells easily "hijack" neuronal networks to synchronization. These results underscore the fact that the degree of synchronization observed in neuronal networks is determined by a complex interplay between network structure and the dynamical properties of individual neurons, indicating that efforts to recover structural connectivity from dynamical correlations must in general take both factors into account.

Effect of heterogeneity and noise on cross frequency phase-phase and phase-amplitude coupling

Cross-frequency coupling is hypothesized to play a functional role in neural computation. We apply phase resetting theory to two types of cross-frequency coupling that can occur when a slower oscillator periodically forces one or more oscillators: phase-phase coupling, in which the two oscillations are phase-locked, and phase-amplitude coupling, in which the amplitude of the driven oscillation is modulated. Our first result is that the shape of the phase resetting curve predicts the tightness of locking to a pulsatile forcing periodic input at any ratio of forced to intrinsic period; the tightness of the locking decreases as the ratio increases. Theoretical expressions were obtained for the probability density of the phases for a population of heterogeneous oscillators or a noisy single oscillator. Results were confirmed using two types of simulated networks and experiments on hippocampal CA1 neurons. Theoretical expressions were also obtained and confirmed for the probability density of N spike times within a single cycle of low frequency forcing. The second result is a suggested mechanism for phase-amplitude coupling in which progressive desynchronization leads to decreasing amplitude during a low frequency forcing cycle. Network simulations confirmed the theoretical viability of this mechanism, and that it generalizes to more diffuse input.

Noise-induced synchronization of oscillatory convection and its optimization

We investigate common-noise-induced phase synchronization between uncoupled identical Hele-Shaw cells exhibiting oscillatory convection. Using the phase description method for oscillatory convection, we demonstrate that the uncoupled systems of oscillatory Hele-Shaw convection can exhibit in-phase synchronization when driven by weak common noise. We derive the Lyapunov exponent determining the relaxation time for the synchronization, and develop a method for obtaining the optimal spatial pattern of the common noise to achieve synchronization. The theoretical results are confirmed by direct numerical simulations.

Stimulus phase locking of cortical oscillation for auditory stream segregation in rats

The phase of cortical oscillations contains rich information and is valuable for encoding sound stimuli. Here we hypothesized that oscillatory phase modulation, instead of amplitude modulation, is a neural correlate of auditory streaming. Our behavioral evaluation provided compelling evidences for the first time that rats are able to organize auditory stream. Local field potentials (LFPs) were investigated in the cortical layer IV or deeper in the primary auditory cortex of anesthetized rats. In response to ABA- sequences with different inter-tone intervals and frequency differences, neurometric functions were characterized with phase locking as well as the band-specific amplitude evoked by test tones. Our results demonstrated that under large frequency differences and short inter-tone intervals, the neurometric function based on stimulus phase locking in higher frequency bands, particularly the gamma band, could better describe van Noorden's perceptual boundary than the LFP amplitude. Furthermore, the gamma-band neurometric function showed a build-up-like effect within around 3 seconds from sequence onset. These findings suggest that phase locking and amplitude have different roles in neural computation, and support our hypothesis that temporal modulation of cortical oscillations should be considered to be neurophysiological mechanisms of auditory streaming, in addition to forward suppression, tonotopic separation, and multi-second adaptation.

Amplitude and phase-locking adaptation of neural oscillation in the rat auditory cortex in response to tone sequence

Sensory adaptation allows stimulus sensitivity to be dynamically modulated according to stimulus statistics and plays pivotal roles in efficient neural computation. Here, it is hypothesized that in the auditory cortex, phase locking of local field potentials (LFPs) to test tones exhibits an adaptation property, i.e., phase-locking adaptation, which is distinct from the amplitude adaptation of oscillatory components. Series of alternating tone sequences were applied in which the inter-tone interval (ITI) and frequency difference (ΔF) between successive tones were varied. Then, adaptation was characterized by the temporal evolution of the band-specific amplitude and phase locking evoked by the test tones. Differences as well as similarities were revealed between amplitude and phase-locking adaptations. First, both amplitude and phase-locking adaptations were enhanced by short ITIs and small ΔFs. Second, the amplitude adaptation was more effective in a higher frequency band, while the phase-locking adaptation was more effective in a lower frequency band. Third, as with the adaptation of multiunit activities (MUAs), the amplitude adaptation occurred mainly within a second, while the phase-locking showed multi-second adaptation specifically in the gamma band for short ITI and small ΔF conditions. Fourth, the amplitude adaptation and phase-locking adaptation were co-modulated in a within-second time scale, while this co-modulation was not observed in a multi-second time scale. These findings suggest that the amplitude and phase-locking adaptations have different mechanisms and functions. The phase-locking adaptation is likely to play more crucial roles in encoding a temporal structure of stimulus than the amplitude adaptation.

Impact of neuronal heterogeneity on correlated colored noise-induced synchronization

Synchronization plays an important role in neural signal processing and transmission. Many hypotheses have been proposed to explain the origin of neural synchronization. In recent years, correlated noise-induced synchronization has received support from many theoretical and experimental studies. However, many of these prior studies have assumed that neurons have identical biophysical properties and that their inputs are well modeled by white noise. In this context, we use colored noise to induce synchronization between oscillators with heterogeneity in both phase-response curves and frequencies. In the low noise limit, we derive novel analytical theory showing that the time constant of colored noise influences correlated noise-induced synchronization and that oscillator heterogeneity can limit synchronization. Surprisingly, however, heterogeneous oscillators may synchronize better than homogeneous oscillators given low input correlations. We also find resonance of oscillator synchronization to colored noise inputs when firing frequencies diverge. Collectively, these results prove robust for both relatively high noise regimes and when applied to biophysically realistic spiking neuron models, and further match experimental recordings from acute brain slices.

The Hodgkin-Huxley heritage: from channels to circuits

The Hodgkin-Huxley studies of the action potential, published 60 years ago, are a central pillar of modern neuroscience research, ranging from molecular investigations of the structural basis of ion channel function to the computational implications at circuit level. In this Symposium Review, we aim to demonstrate the ongoing impact of Hodgkin's and Huxley's ideas. The Hodgkin-Huxley model established a framework in which to describe the structural and functional properties of ion channels, including the mechanisms of ion permeation, selectivity, and gating. At a cellular level, the model is used to understand the conditions that control both the rate and timing of action potentials, essential for neural encoding of information. Finally, the Hodgkin-Huxley formalism is central to computational neuroscience to understand both neuronal integration and circuit level information processing, and how these mechanisms might have evolved to minimize energy cost.

Stimulus features, resetting curves, and the dependence on adaptation

We derive a formula that relates the spike-triggered covariance (STC) to the phase resetting curve (PRC) of a neural oscillator. We use this to show how changes in the shape of the PRC alter the sensitivity of the neuron to different stimulus features, which are the eigenvectors of the STC. We compute the PRC and STC for some biophysical models. We compare the STCs and their spectral properties for a two-parameter family of PRCs. Surprisingly, the skew of the PRC has a larger effect on the spectrum and shape of the STC than does the bimodality of the PRC (which plays a large role in synchronization properties). Finally, we relate the STC directly to the spike-triggered average and apply this theory to an olfactory bulb mitral cell recording.

The shape of phase-resetting curves in oscillators with a saddle node on an invariant circle bifurcation

We introduce a simple two-dimensional model that extends the Poincaré oscillator so that the attracting limit cycle undergoes a saddle node bifurcation on an invariant circle (SNIC) for certain parameter values. Arbitrarily close to this bifurcation, the phase-resetting curve (PRC) continuously depends on parameters, where its shape can be not only primarily positive or primarily negative but also nearly sinusoidal. This example system shows that one must be careful inferring anything about the bifurcation structure of the oscillator from the shape of its PRC.

Intrinsic heterogeneity in oscillatory dynamics limits correlation-induced neural synchronization

Synchronous neural oscillations are found throughout the brain and are thought to contribute to neural coding and the propagation of activity. Several proposed mechanisms of synchronization have gained support through combined theoretical and experimental investigation, including mechanisms based on coupling and correlated input. Here, we ask how correlation-induced synchrony is affected by physiological heterogeneity across neurons. To address this question, we examined cell-to-cell differences in phase-response curves (PRCs), which characterize the response of periodically firing neurons to weak perturbations. Using acute slice electrophysiology, we measured PRCs across a single class of principal neurons capable of sensory-evoked oscillations in vivo: the olfactory bulb mitral cells (MCs). Periodically firing MCs displayed a broad range of PRCs, each of which was well fit by a simple three-parameter model. MCs also displayed differences in firing rate-current relationships and in preferred firing rate ranges. Both the observed PRC heterogeneity and moderate firing rate differences (∼10 Hz) separately reduced the maximum correlation-induced synchrony between MCs by up to 25-30%. Simulations further demonstrated that these components of heterogeneity alone were sufficient to account for the difference in synchronization among heterogeneous vs. homogeneous populations in vitro. Within this simulation framework, independent modulation of specific PRC features additionally revealed which aspects of PRC heterogeneity most strongly impact correlation-induced synchronization. Finally, we demonstrated good agreement of novel mathematical theory with our experimental and simulation results, providing a theoretical basis for the influence of heterogeneity on correlation-induced neural synchronization.

Optimal phase response functions for fast pulse-coupled synchronization in wireless sensor networks

Synchronization is crucial to wireless sensor networks. Recently a pulse-coupled synchronization strategy that emulates biological pulse-coupled agents has been used to achieve this goal. We propose to optimize the phase response function such that synchronization rate is maximized. Since the synchronization rate is increased independently of transmission power, energy consumption is reduced, hence extending the life of battery-powered sensor networks. Comparison with existing phase response functions confirms the effectiveness of the method.

Efficient estimation of phase-response curves via compressive sensing

The phase-response curve (PRC), relating the phase shift of an oscillator to external perturbation, is an important tool to study neurons and their population behavior. It can be experimentally estimated by measuring the phase changes caused by probe stimuli. These stimuli, usually short pulses or continuous noise, have a much wider frequency spectrum than that of neuronal dynamics. This makes the experimental data high dimensional while the number of data samples tends to be small. Current PRC estimation methods have not been optimized for efficiently discovering the relevant degrees of freedom from such data. We propose a systematic and efficient approach based on a recently developed signal processing theory called compressive sensing (CS). CS is a framework for recovering sparsely constructed signals from undersampled data and is suitable for extracting information about the PRC from finite but high-dimensional experimental measurements. We illustrate how the CS algorithm can be translated into an estimation scheme and demonstrate that our CS method can produce good estimates of the PRCs with simulated and experimental data, especially when the data size is so small that simple approaches such as naive averaging fail. The tradeoffs between degrees of freedom vs. goodness-of-fit were systematically analyzed, which help us to understand better what part of the data has the most predictive power. Our results illustrate that finite sizes of neuroscientific data in general compounded by large dimensionality can hamper studies of the neural code and suggest that CS is a good tool for overcoming this challenge.

A-current and type I/type II transition determine collective spiking from common input

The mechanisms and impact of correlated, or synchronous, firing among pairs and groups of neurons are under intense investigation throughout the nervous system. A ubiquitous circuit feature that can give rise to such correlations consists of overlapping, or common, inputs to pairs and populations of cells, leading to common spike train responses. Here, we use computational tools to study how the transfer of common input currents into common spike outputs is modulated by the physiology of the recipient cells. We focus on a key conductance, g(A), for the A-type potassium current, which drives neurons between "type II" excitability (low g(A)), and "type I" excitability (high g(A)). Regardless of g(A), cells transform common input fluctuations into a tendency to spike nearly simultaneously. However, this process is more pronounced at low g(A) values. Thus, for a given level of common input, type II neurons produce spikes that are relatively more correlated over short time scales. Over long time scales, the trend reverses, with type II neurons producing relatively less correlated spike trains. This is because these cells' increased tendency for simultaneous spiking is balanced by an anticorrelation of spikes at larger time lags. These findings extend and interpret prior findings for phase oscillators to conductance-based neuron models that cover both oscillatory (superthreshold) and subthreshold firing regimes. We demonstrate a novel implication for neural signal processing: downstream cells with long time constants are selectively driven by type I cell populations upstream and those with short time constants by type II cell populations. Our results are established via high-throughput numerical simulations and explained via the cells' filtering properties and nonlinear dynamics.

Robustness, variability, phase dependence, and longevity of individual synaptic input effects on spike timing during fluctuating synaptic backgrounds: a modeling study of globus pallidus neuron phase response properties

A neuron's phase response curve (PRC) shows how inputs arriving at different times during the spike cycle differentially affect the timing of subsequent spikes. Using a full morphological model of a globus pallidus (GP) neuron, we previously demonstrated that dendritic conductances shape the PRC in a spike frequency-dependent manner, suggesting different functional roles of perisomatic and distal dendritic synapses in the control of patterned network activity. In the present study we extend this analysis to examine the impact of physiologically realistic high conductance states on somatic and dendritic PRCs and the time course of spike train perturbations. First, we found that average somatic and dendritic PRCs preserved their shapes and spike frequency dependence when the model was driven by spatially-distributed, stochastic conductance inputs rather than tonic somatic current. However, responses to inputs during specific synaptic backgrounds often deviated substantially from the average PRC. Therefore, we analyzed the interactions of PRC stimuli with transient fluctuations in the synaptic background on a trial-by-trial basis. We found that the variability in responses to PRC stimuli and the incidence of stimulus-evoked added or skipped spikes were stimulus-phase-dependent and reflected the profile of the average PRC, suggesting commonality in the underlying mechanisms. Clear differences in the relation between the phase of input and variability of spike response between dendritic and somatic inputs indicate that these regions generally represent distinct dynamical subsystems of synaptic integration with respect to influencing the stability of spike time attractors generated by the overall synaptic conductance.

Cellular and circuit mechanisms maintain low spike co-variability and enhance population coding in somatosensory cortex

The responses of cortical neurons are highly variable across repeated presentations of a stimulus. Understanding this variability is critical for theories of both sensory and motor processing, since response variance affects the accuracy of neural codes. Despite this influence, the cellular and circuit mechanisms that shape the trial-to-trial variability of population responses remain poorly understood. We used a combination of experimental and computational techniques to uncover the mechanisms underlying response variability of populations of pyramidal (E) cells in layer 2/3 of rat whisker barrel cortex. Spike trains recorded from pairs of E-cells during either spontaneous activity or whisker deflected responses show similarly low levels of spiking co-variability, despite large differences in network activation between the two states. We developed network models that show how spike threshold non-linearities dilute E-cell spiking co-variability during spontaneous activity and low velocity whisker deflections. In contrast, during high velocity whisker deflections, cancelation mechanisms mediated by feedforward inhibition maintain low E-cell pairwise co-variability. Thus, the combination of these two mechanisms ensure low E-cell population variability over a wide range of whisker deflection velocities. Finally, we show how this active decorrelation of population variability leads to a drastic increase in the population information about whisker velocity. The prevalence of spiking non-linearities and feedforward inhibition in the nervous system suggests that the mechanisms for low network variability presented in our study may generalize throughout the brain.

Temperature-modulated synchronization transition in coupled neuronal oscillators

We study two firing properties to characterize the activities of a neuron: frequency-current (f-I) curves and phase response curves (PRCs), with variation in the intrinsic temperature scaling parameter (μ) controlling the opening and closing of ionic channels. We show a peak of the firing frequency for small μ in a class I neuron with the I value immediately after the saddle-node bifurcation, which is entirely different from previous experimental reports as well as model studies. The PRC takes a type II form on a logarithmic f-I curve when μ is small. Then, we analyze the synchronization phenomena in a two-neuron network using the phase-reduction method. We find common μ-dependent transition and bifurcation of synchronizations, regardless of the values of I. Such results give us helpful insight into synchronizations tuned with a sinusoidal-wave temperature modulation on neurons.

Balanced synaptic input shapes the correlation between neural spike trains

Stimulus properties, attention, and behavioral context influence correlations between the spike times produced by a pair of neurons. However, the biophysical mechanisms that modulate these correlations are poorly understood. With a combined theoretical and experimental approach, we show that the rate of balanced excitatory and inhibitory synaptic input modulates the magnitude and timescale of pairwise spike train correlation. High rate synaptic inputs promote spike time synchrony rather than long timescale spike rate correlations, while low rate synaptic inputs produce opposite results. This correlation shaping is due to a combination of enhanced high frequency input transfer and reduced firing rate gain in the high input rate state compared to the low state. Our study extends neural modulation from single neuron responses to population activity, a necessary step in understanding how the dynamics and processing of neural activity change across distinct brain states.

Correlation transfer from basal ganglia to thalamus in Parkinson's disease

Spike trains from neurons in the basal ganglia of parkinsonian primates show increased pairwise correlations, oscillatory activity, and burst rate compared to those from neurons recorded during normal brain activity. However, it is not known how these changes affect the behavior of downstream thalamic neurons. To understand how patterns of basal ganglia population activity may affect thalamic spike statistics, we study pairs of model thalamocortical (TC) relay neurons receiving correlated inhibitory input from the internal segment of the globus pallidus (GPi), a primary output nucleus of the basal ganglia. We observe that the strength of correlations of TC neuron spike trains increases with the GPi correlation level, and bursty firing patterns such as those seen in the parkinsonian GPi allow for stronger transfer of correlations than do firing patterns found under normal conditions. We also show that the T-current in the TC neurons does not significantly affect correlation transfer, despite its pronounced effects on spiking. Oscillatory firing patterns in GPi are shown to affect the timescale at which correlations are best transferred through the system. To explain this last result, we analytically compute the spike count correlation coefficient for oscillatory cases in a reduced point process model. Our analysis indicates that the dependence of the timescale of correlation transfer is robust to different levels of input spike and rate correlations and arises due to differences in instantaneous spike correlations, even when the long timescale rhythmic modulations of neurons are identical. Overall, these results show that parkinsonian firing patterns in GPi do affect the transfer of correlations to the thalamus.

Correlation transfer in stochastically driven neural oscillators over long and short time scales

In the absence of synaptic coupling, two or more neural oscillators may become synchronized by virtue of the statistical correlations in their noisy input streams. Recent work has shown that the degree of correlation transfer from input currents to output spikes depends not only on intrinsic oscillator dynamics, but also on the length of the observation window over which the correlation is calculated. In this paper we use stochastic phase reduction and regular perturbations to derive the correlation of the total phase elapsed over long time scales, a quantity that provides a convenient proxy for the spike count correlation. Over short time scales, we derive the spike count correlation directly using straightforward probabilistic reasoning applied to the density of the phase difference. Our approximations show that output correlation scales with the autocorrelation of the phase resetting curve over long time scales. We also find a concise expression for the influence of the shape of the phase resetting curve on the initial slope of the output correlation over short time scales. These analytic results together with numerical simulations provide new intuitions for the recent counterintuitive finding that type I oscillators transfer correlations more faithfully than do type II over long time scales, while the reverse holds true for the better understood case of short time scales.

Measurement of infinitesimal phase response curves from noisy real neurons

We sought to measure infinitesimal phase response curves (iPRCs) from rat hippocampal CA1 pyramidal neurons. It is difficult to measure iPRCs from noisy neurons because of the dilemma that either the linearity or the signal-to-noise ratio of responses to external perturbations must be sacrificed. To overcome this difficulty, we used an iPRC measurement model formulated as the Langevin phase equation (LPE) to extract iPRCs in the Bayesian scheme. We then simultaneously verified the effectiveness of the measurement model and the reliability of the estimated iPRCs by demonstrating that LPEs with the estimated iPRCs could predict the stochastic behaviors of the same neurons, whose iPRCs had been measured, when they were perturbed by periodic stimulus currents. Our results suggest that the LPE is an effective model for real oscillating neurons and that many theoretical frameworks based on it may be applicable to real nerve systems.

Optimal phase response curves for stochastic synchronization of limit-cycle oscillators by common Poisson noise

We consider optimization of phase response curves for stochastic synchronization of noninteracting limit-cycle oscillators by common Poisson impulsive signals. The optimal functional shape for sufficiently weak signals is sinusoidal, but can differ for stronger signals. By solving the Euler-Lagrange equation associated with the minimization of the Lyapunov exponent characterizing synchronization efficiency, the optimal phase response curve is obtained. We show that the optimal shape mutates from a sinusoid to a sawtooth as the constraint on its squared amplitude is varied.

Robust entrainment of circadian oscillators requires specific phase response curves

The circadian clocks keeping time in many living organisms rely on self-sustained biochemical oscillations entrained by external cues, such as light, to the 24-h cycle induced by Earth's rotation. However, environmental cues are unreliable due to the variability of habitats, weather conditions, or cue-sensing mechanisms among individuals. A tempting hypothesis is that circadian clocks have evolved so as to be robust to fluctuations in the signal that entrains them. To support this hypothesis, we analyze the synchronization behavior of weakly and periodically forced oscillators in terms of their phase response curve (PRC), which measures phase changes induced by a perturbation applied at different times of the cycle. We establish a general relationship between the robustness of key entrainment properties, such as stability and oscillator phase, on the one hand, and the shape of the PRC as characterized by a specific curvature or the existence of a dead zone, on the other hand. The criteria obtained are applied to computational models of circadian clocks and account for the disparate robustness properties of various forcing schemes. Finally, the analysis of PRCs measured experimentally in several organisms strongly suggests a case of convergent evolution toward an optimal strategy for maintaining a clock that is accurate and robust to environmental fluctuations.

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.

Three synaptic components contributing to robust network synchronization

Robust synchronized activity is widely observed in real neural systems. However, a mechanism for robust synchronization that can be understood analytically, and has a clear physical meaning, remains elusive. This paper considers such a mechanism by formalizing three synaptic components contributing to robust synchronization in networks with heterogeneous external drive currents and conductance-based synapses. The first component arises from the assumption that the aggregate post-synaptic potential of a neuron decays more if it fires later within a spike volley. The second component results because neurons with smaller drives have reached a lower membrane potential at the time when the volley of inputs arrives than that reached by neurons with larger drives. The third component arises from the assumption that neurons firing later in the previous volley have had less time to integrate their drives than neurons firing earlier have had, again causing a lower membrane potential at the time when the volley of inputs arrives. Because of the voltage-dependent properties of synaptic currents, either of the last two components will cause smaller inhibitions for the later-firing neurons if the synapses are inhibitory. This smaller inhibition causes the later-firing neurons to fire earlier in the next cycle, thereby forcing them toward synchrony. With these three synaptic components, we discuss the relationship between the robustness of the synchrony and the parameters, search for the optimal parameter set for the robust network synchronization of a certain frequency band, and demonstrate the key role of the voltage-dependent properties of synaptic currents in robust or stable synchronization.

Coupled noisy spiking neurons as velocity-controlled oscillators in a model of grid cell spatial firing

One of the two primary classes of models of grid cell spatial firing uses interference between oscillators at dynamically modulated frequencies. Generally, these models are presented in terms of idealized oscillators (modeled as sinusoids), which differ from biological oscillators in multiple important ways. Here we show that two more realistic, noisy neural models (Izhikevich's simple model and a biophysical model of an entorhinal cortex stellate cell) can be successfully used as oscillators in a model of this type. When additive noise is included in the models such that uncoupled or sparsely coupled cells show realistic interspike interval variance, both synaptic and gap-junction coupling can synchronize networks of cells to produce comparatively less variable network-level oscillations. We show that the frequency of these oscillatory networks can be controlled sufficiently well to produce stable grid cell spatial firing on the order of at least 2-5 min, despite the high noise level. Our results suggest that the basic principles of oscillatory interference models work with more realistic models of noisy neurons. Nevertheless, a number of simplifications were still made and future work should examine increasingly realistic models.

Phase-response curves and synchronized neural networks

We review the principal assumptions underlying the application of phase-response curves (PRCs) to synchronization in neuronal networks. The PRC measures how much a given synaptic input perturbs spike timing in a neural oscillator. Among other applications, PRCs make explicit predictions about whether a given network of interconnected neurons will synchronize, as is often observed in cortical structures. Regarding the assumptions of the PRC theory, we conclude: (i) The assumption of noise-tolerant cellular oscillations at or near the network frequency holds in some but not all cases. (ii) Reduced models for PRC-based analysis can be formally related to more realistic models. (iii) Spike-rate adaptation limits PRC-based analysis but does not invalidate it. (iv) The dependence of PRCs on synaptic location emphasizes the importance of improving methods of synaptic stimulation. (v) New methods can distinguish between oscillations that derive from mutual connections and those arising from common drive. (vi) It is helpful to assume linear summation of effects of synaptic inputs; experiments with trains of inputs call this assumption into question. (vii) Relatively subtle changes in network structure can invalidate PRC-based predictions. (viii) Heterogeneity in the preferred frequencies of component neurons does not invalidate PRC analysis, but can annihilate synchronous activity.

Synchronization of uncoupled oscillators by common gamma impulses: From phase locking to noise-induced synchronization

Nonlinear oscillators can mutually synchronize when they are driven by common external impulses. Two important scenarios are (i) synchronization resulting from phase locking of each oscillator to regular periodic impulses and (ii) noise-induced synchronization caused by the Poisson random impulses, but their difference has not been fully quantified. Here, we analyze a pair of uncoupled oscillators subject to common random impulses with gamma-distributed intervals, which can be smoothly interpolated between the regular periodic and the random Poisson impulses. Their dynamics are characterized by phase distributions, frequency detuning, Lyapunov exponents, and information-theoretic measures, which clearly reveal the differences between the two synchronization scenarios.

A comparison of methods to determine neuronal phase-response curves

The phase-response curve (PRC) is an important tool to determine the excitability type of single neurons which reveals consequences for their synchronizing properties. We review five methods to compute the PRC from both model data and experimental data and compare the numerically obtained results from each method. The main difference between the methods lies in the reliability which is influenced by the fluctuations in the spiking data and the number of spikes available for analysis. We discuss the significance of our results and provide guidelines to choose the best method based on the available data.

Phase response curve analysis of a full morphological globus pallidus neuron model reveals distinct perisomatic and dendritic modes of synaptic integration

Synchronization of globus pallidus (GP) neurons and cortically entrained oscillations between GP and other basal ganglia nuclei are key features of the pathophysiology of Parkinson's disease. Phase response curves (PRCs), which tabulate the effects of phasic inputs within a neuron's spike cycle on output spike timing, are efficient tools for predicting the emergence of synchronization in neuronal networks and entrainment to periodic input. In this study we apply physiologically realistic synaptic conductance inputs to a full morphological GP neuron model to determine the phase response properties of the soma and different regions of the dendritic tree. We find that perisomatic excitatory inputs delivered throughout the interspike interval advance the phase of the spontaneous spike cycle yielding a type I PRC. In contrast, we demonstrate that distal dendritic excitatory inputs can either delay or advance the next spike depending on whether they occur early or late in the spike cycle. We find this latter pattern of responses, summarized by a biphasic (type II) PRC, was a consequence of dendritic activation of the small conductance calcium-activated potassium current, SK. We also evaluate the spike-frequency dependence of somatic and dendritic PRC shapes, and we demonstrate the robustness of our results to variations of conductance densities, distributions, and kinetic parameters. We conclude that the distal dendrite of GP neurons embodies a distinct dynamical subsystem that could promote synchronization of pallidal networks to excitatory inputs. These results highlight the need to consider different effects of perisomatic and dendritic inputs in the control of network behavior.

Amplification of asynchronous inhibition-mediated synchronization by feedback in recurrent networks

Synchronization of 30–80 Hz oscillatory activity of the principle neurons in the olfactory bulb (mitral cells) is believed to be important for odor discrimination. Previous theoretical studies of these fast rhythms in other brain areas have proposed that principle neuron synchrony can be mediated by short-latency, rapidly decaying inhibition. This phasic inhibition provides a narrow time window for the principle neurons to fire, thus promoting synchrony. However, in the olfactory bulb, the inhibitory granule cells produce long lasting, small amplitude, asynchronous and aperiodic inhibitory input and thus the narrow time window that is required to synchronize spiking does not exist. Instead, it has been suggested that correlated output of the granule cells could serve to synchronize uncoupled mitral cells through a mechanism called “stochastic synchronization”, wherein the synchronization arises through correlation of inputs to two neural oscillators. Almost all work on synchrony due to correlations presumes that the correlation is imposed and fixed. Building on theory and experiments that we and others have developed, we show that increased synchrony in the mitral cells could produce an increase in granule cell activity for those granule cells that share a synchronous group of mitral cells. Common granule cell input increases the input correlation to the mitral cells and hence their synchrony by providing a positive feedback loop in correlation. Thus we demonstrate the emergence and temporal evolution of input correlation in recurrent networks with feedback. We explore several theoretical models of this idea, ranging from spiking models to an analytically tractable model.

Neurons in many parts of the brain fire spikes rhythmically and synchronously in many behaviorally and functionally relevant contexts. There are many mechanisms for producing oscillatory synchronization between populations of biological oscillators. One way to produce synchrony is that the population of oscillators receives common correlated input. In this paper, we study a population of oscillating neurons (mitral cells) that are not directly coupled to each other but receive broadband correlated input from a second population of neurons (granule cells). The granule cell population, in turn, receives inputs from the mitral cells; hence, the mitral and granule cells are reciprocally connected. Correlated input to the oscillating mitral cells produces tighter synchrony in the activity of the mitral cell population. We hypothesize that this increased mitral cell synchrony will evoke greater activity in specific groups of granule cells and that these specific granule cells, in turn, become the source of the correlated input to the mitral cells. That is, the synchronous input from the mitral cells increases the fraction of correlated feedback. Thus, we close the correlation loop. We show through analysis and simulations that this feedback mechanism can lead to the spontaneous appearance of highly synchronous activity within the mitral cells. We show that there is good experimental support for this mechanism in the circuitry of the olfactory bulb. We speculate that such mechanisms could also arise in other parts of the brain.

Coupling regularizes individual units in noisy populations

The regularity of a noisy system can modulate in various ways. It is well known that coupling in a population can lower the variability of the entire network; the collective activity is more regular. Here, we show that diffusive (reciprocal) coupling of two simple Ornstein-Uhlenbeck (O-U) processes can regularize the individual, even when it is coupled to a noisier process. In cellular networks, the regularity of individual cells is important when a select few play a significant role. The regularizing effect of coupling surprisingly applies also to general nonlinear noisy oscillators. However, unlike with the O-U process, coupling-induced regularity is robust to different kinds of coupling. With two coupled noisy oscillators, we derive an asymptotic formula assuming weak noise and coupling for the variance of the period (i.e., spike times) that accurately captures this effect. Moreover, we find that reciprocal coupling can regularize the individual period of higher dimensional oscillators such as the Morris-Lecar and Brusselator models, even when coupled to noisier oscillators. Coupling can have a counterintuitive and beneficial effect on noisy systems. These results have implications for the role of connectivity with noisy oscillators and the modulation of variability of individual oscillators.

Time scales of spike-train correlation for neural oscillators with common drive

We examine the effect of the phase-resetting curve on the transfer of correlated input signals into correlated output spikes in a class of neural models receiving noisy superthreshold stimulation. We use linear-response theory to approximate the spike correlation coefficient in terms of moments of the associated exit time problem and contrast the results for type I vs type II models and across the different time scales over which spike correlations can be assessed. We find that, on long time scales, type I oscillators transfer correlations much more efficiently than type II oscillators. On short time scales this trend reverses, with the relative efficiency switching at a time scale that depends on the mean and standard deviation of input currents. This switch occurs over time scales that could be exploited by downstream circuits.

Evaluation of the oscillatory interference model of grid cell firing through analysis and measured period variance of some biological oscillators

Models of the hexagonally arrayed spatial activity pattern of grid cell firing in the literature generally fall into two main categories: continuous attractor models or oscillatory interference models. Burak and Fiete (2009, PLoS Comput Biol) recently examined noise in two continuous attractor models, but did not consider oscillatory interference models in detail. Here we analyze an oscillatory interference model to examine the effects of noise on its stability and spatial firing properties. We show analytically that the square of the drift in encoded position due to noise is proportional to time and inversely proportional to the number of oscillators. We also show there is a relatively fixed breakdown point, independent of many parameters of the model, past which noise overwhelms the spatial signal. Based on this result, we show that a pair of oscillators are expected to maintain a stable grid for approximately t = 5µ³/(4πσ)² seconds where µ is the mean period of an oscillator in seconds and σ² its variance in seconds². We apply this criterion to recordings of individual persistent spiking neurons in postsubiculum (dorsal presubiculum) and layers III and V of entorhinal cortex, to subthreshold membrane potential oscillation recordings in layer II stellate cells of medial entorhinal cortex and to values from the literature regarding medial septum theta bursting cells. All oscillators examined have expected stability times far below those seen in experimental recordings of grid cells, suggesting the examined biological oscillators are unfit as a substrate for current implementations of oscillatory interference models. However, oscillatory interference models can tolerate small amounts of noise, suggesting the utility of circuit level effects which might reduce oscillator variability. Further implications for grid cell models are discussed.

For many animals, including rats, accurate spatial memory over relatively large areas is important in order to find food and shelter. Just as unique points in time can be efficiently represented by combinations of repeating elements like hours, days, and months, points in space can be represented as combinations of elements that repeat at different spatial scales. Just such a code has been identified in the brains of rats and it shows an intriguing triangular spacing of encoded locations. Two different explanations have been developed as to what general mechanism in the brain might be able to generate this unusual code. However, to date there is not conclusive experimental evidence indicating whether either of the two explanations is correct. Here we show in detail that one of the explanations, called oscillatory interference, has specific requirements regarding the amount of variability in the system that implements it. We then report data experimentally examining candidate systems to evaluate their levels of noise. The large amount of noise that we find presents a challenge to the currently suggested biological implementations of oscillatory interference, but it does not provide support for the alternative explanation.

Type-II phase resetting curve is optimal for stochastic synchrony

The phase-resetting curve (PRC) describes the response of a neural oscillator to small perturbations in membrane potential. Its usefulness for predicting the dynamics of weakly coupled deterministic networks has been well characterized. However, the inputs to real neurons may often be more accurately described as barrages of synaptic noise. Effective connectivity between cells may thus arise in the form of correlations between the noisy input streams. We use constrained optimization and perturbation methods to prove that the PRC shape determines susceptibility to synchrony among otherwise uncoupled noise-driven neural oscillators. PRCs can be placed into two general categories: type-I PRCs are non-negative, while type-II PRCs have a large negative region. Here we show that oscillators with type-II PRCs receiving common noisy input synchronize more readily than those with type-I PRCs.

Synchronization dynamics of two coupled neural oscillators receiving shared and unshared noisy stimuli

The response of neurons to external stimuli greatly depends on the intrinsic dynamics of the network. Here, the intrinsic dynamics are modeled as coupling and the external input is modeled as shared and unshared noise. We assume the neurons are repetitively firing action potentials (i.e., neural oscillators), are weakly and identically coupled, and the external noise is weak. Shared noise can induce bistability between the synchronous and anti-phase states even though the anti-phase state is the only stable state in the absence of noise. We study the Fokker-Planck equation of the system and perform an asymptotic reduction rho(0). The rho(0) solution is more computationally efficient than both the Monte Carlo simulations and the 2D Fokker-Planck solver, and agrees remarkably well with the full system with weak noise and weak coupling. With moderate noise and coupling, rho(0) is still qualitatively correct despite the small noise and coupling assumption in the asymptotic reduction. Our phase model accurately predicts the behavior of a realistic synaptically coupled Morris-Lecar system.

Reliability, synchrony and noise

The brain is noisy. Neurons receive tens of thousands of highly fluctuating inputs and generate spike trains that appear highly irregular. Much of this activity is spontaneous - uncoupled to overt stimuli or motor outputs - leading to questions about the functional impact of this noise. Although noise is most often thought of as disrupting patterned activity and interfering with the encoding of stimuli, recent theoretical and experimental work has shown that noise can play a constructive role - leading to increased reliability or regularity of neuronal firing in single neurons and across populations. These results raise fundamental questions about how noise can influence neural function and computation.