A new approach for determining phase response curves reveals that Purkinje cells can act as perfect integrators

How neuronal morphology impacts the synchronisation state of neuronal networks

The biophysical properties of neurons not only affect how information is processed within cells, they can also impact the dynamical states of the network. Specifically, the cellular dynamics of action-potential generation have shown relevance for setting the (de)synchronisation state of the network. The dynamics of tonically spiking neurons typically fall into one of three qualitatively distinct types that arise from distinct mathematical bifurcations of voltage dynamics at the onset of spiking. Accordingly, changes in ion channel composition or even external factors, like temperature, have been demonstrated to switch network behaviour via changes in the spike onset bifurcation and hence its associated dynamical type. A thus far less addressed modulator of neuronal dynamics is cellular morphology. Based on simplified and anatomically realistic mathematical neuron models, we show here that the extent of dendritic arborisation has an influence on the neuronal dynamical spiking type and therefore on the (de)synchronisation state of the network. Specifically, larger dendritic trees prime neuronal dynamics for in-phase-synchronised or splayed-out activity in weakly coupled networks, in contrast to cells with otherwise identical properties yet smaller dendrites. Our biophysical insights hold for generic multicompartmental classes of spiking neuron models (from ball-and-stick-type to anatomically reconstructed models) and establish a connection between neuronal morphology and the susceptibility of neural tissue to synchronisation in health and disease.

Basket to Purkinje Cell Inhibitory Ephaptic Coupling Is Abolished in Episodic Ataxia Type 1

Dominantly inherited missense mutations of the KCNA1 gene, which encodes the KV1.1 potassium channel subunit, cause Episodic Ataxia type 1 (EA1). Although the cerebellar incoordination is thought to arise from abnormal Purkinje cell output, the underlying functional deficit remains unclear. Here we examine synaptic and non-synaptic inhibition of Purkinje cells by cerebellar basket cells in an adult mouse model of EA1. The synaptic function of basket cell terminals was unaffected, despite their intense enrichment for KV1.1-containing channels. In turn, the phase response curve quantifying the influence of basket cell input on Purkine cell output was maintained. However, ultra-fast non-synaptic ephaptic coupling, which occurs in the cerebellar 'pinceau' formation surrounding the axon initial segment of Purkinje cells, was profoundly reduced in EA1 mice in comparison with their wild type littermates. The altered temporal profile of basket cell inhibition of Purkinje cells underlines the importance of Kv1.1 channels for this form of signalling, and may contribute to the clinical phenotype of EA1.

Firing rate-dependent phase responses of Purkinje cells support transient oscillations

Both spike rate and timing can transmit information in the brain. Phase response curves (PRCs) quantify how a neuron transforms input to output by spike timing. PRCs exhibit strong firing-rate adaptation, but its mechanism and relevance for network output are poorly understood. Using our Purkinje cell (PC) model, we demonstrate that the rate adaptation is caused by rate-dependent subthreshold membrane potentials efficiently regulating the activation of Na+ channels. Then, we use a realistic PC network model to examine how rate-dependent responses synchronize spikes in the scenario of reciprocal inhibition-caused high-frequency oscillations. The changes in PRC cause oscillations and spike correlations only at high firing rates. The causal role of the PRC is confirmed using a simpler coupled oscillator network model. This mechanism enables transient oscillations between fast-spiking neurons that thereby form PC assemblies. Our work demonstrates that rate adaptation of PRCs can spatio-temporally organize the PC input to cerebellar nuclei.

The effects of chloride dynamics on substantia nigra pars reticulata responses to pallidal and striatal inputs

As a rodent basal ganglia (BG) output nucleus, the substantia nigra pars reticulata (SNr) is well positioned to impact behavior. SNr neurons receive GABAergic inputs from the striatum (direct pathway) and globus pallidus (GPe, indirect pathway). Dominant theories of action selection rely on these pathways' inhibitory actions. Yet, experimental results on SNr responses to these inputs are limited and include excitatory effects. Our study combines experimental and computational work to characterize, explain, and make predictions about these pathways. We observe diverse SNr responses to stimulation of SNr-projecting striatal and GPe neurons, including biphasic and excitatory effects, which our modeling shows can be explained by intracellular chloride processing. Our work predicts that ongoing GPe activity could tune the SNr operating mode, including its responses in decision-making scenarios, and GPe output may modulate synchrony and low-frequency oscillations of SNr neurons, which we confirm using optogenetic stimulation of GPe terminals within the SNr.

How to correctly quantify neuronal phase-response curves from noisy recordings

At the level of individual neurons, various coding properties can be inferred from the input-output relationship of a cell. For small inputs, this relation is captured by the phase-response curve (PRC), which measures the effect of a small perturbation on the timing of the subsequent spike. Experimentally, however, an accurate experimental estimation of PRCs is challenging. Despite elaborate measurement efforts, experimental PRC estimates often cannot be related to those from modeling studies. In particular, experimental PRCs rarely resemble the characteristic theoretical PRC expected close to spike initiation, which is indicative of the underlying spike-onset bifurcation. Here, we show for conductance-based model neurons that the correspondence between theoretical and measured phase-response curve is lost when the stimuli used for the estimation are too large. In this case, the derived phase-response curve is distorted beyond recognition and takes on a generic shape that reflects the measurement protocol and masks the spike-onset bifurcation. We discuss how to identify appropriate stimulus strengths for perturbation and noise-stimulation methods, which permit to estimate PRCs that reliably reflect the spike-onset bifurcation - a task that is particularly difficult if a lower bound for the stimulus amplitude is dictated by prominent intrinsic neuronal noise.

Macroscopic phase resetting-curves determine oscillatory coherence and signal transfer in inter-coupled neural circuits

Macroscopic oscillations of different brain regions show multiple phase relationships that are persistent across time and have been implicated in routing information. While multiple cellular mechanisms influence the network oscillatory dynamics and structure the macroscopic firing motifs, one of the key questions is to identify the biophysical neuronal and synaptic properties that permit such motifs to arise. A second important issue is how the different neural activity coherence states determine the communication between the neural circuits. Here we analyse the emergence of phase-locking within bidirectionally delayed-coupled spiking circuits in which global gamma band oscillations arise from synaptic coupling among largely excitable neurons. We consider both the interneuronal (ING) and the pyramidal-interneuronal (PING) population gamma rhythms and the inter coupling targeting the pyramidal or the inhibitory neurons. Using a mean-field approach together with an exact reduction method, we reduce each spiking network to a low dimensional nonlinear system and derive the macroscopic phase resetting-curves (mPRCs) that determine how the phase of the global oscillation responds to incoming perturbations. This is made possible by the use of the quadratic integrate-and-fire model together with a Lorentzian distribution of the bias current. Depending on the type of gamma (PING vs. ING), we show that incoming excitatory inputs can either speed up the macroscopic oscillation (phase advance; type I PRC) or induce both a phase advance and a delay (type II PRC). From there we determine the structure of macroscopic coherence states (phase-locking) of two weakly synaptically-coupled networks. To do so we derive a phase equation for the coupled system which links the synaptic mechanisms to the coherence states of the system. We show that a synaptic transmission delay is a necessary condition for symmetry breaking, i.e. a non-symmetric phase lag between the macroscopic oscillations. This potentially provides an explanation to the experimentally observed variety of gamma phase-locking modes. Our analysis further shows that symmetry-broken coherence states can lead to a preferred direction of signal transfer between the oscillatory networks where this directionality also depends on the timing of the signal. Hence we suggest a causal theory for oscillatory modulation of functional connectivity between cortical circuits.

Large scale brain oscillations emerge from synaptic interactions within neuronal circuits. Over the past years, such macroscopic rhythms have been suggested to play a crucial role in routing the flow of information across cortical regions, resulting in a functional connectome. The underlying mechanism is cortical oscillations that bind together following a well-known motif called phase-locking. While there is significant experimental support for multiple phase-locking modes in the brain, it is still unclear what is the underlying mechanism that permits macroscopic rhythms to phase lock. In the present paper we take up with this issue, and to show that, one can study the emergent macroscopic phase-locking within the mathematical framework of weakly coupled oscillators. We find that under synaptic delays, fully symmetrically coupled networks can display symmetry-broken states of activity, where one network starts to lead in phase the second (also sometimes known as stuttering states). When we analyse how incoming transient signals affect the coupled system, we find that in the symmetry-broken state, the effect depends strongly on which network is targeted (the leader or the follower) as well as the timing of the input. Hence we show how the dynamics of the emergent phase-locked activity imposes a functional directionality on how signals are processed. We thus offer clarification on the synaptic and circuit properties responsible for the emergence of multiple phase-locking patterns and provide support for its functional implication in information transfer.

Temporal adjustment of short calls according to a partner during vocal turn-taking in Japanese macaques

Turn-taking is a common feature in human speech, and is also seen in the communication of other primate species. However, evidence of turn-taking in vocal exchanges within a short time frame is still scarce in nonhuman primates. This study investigated whether dynamic adjustment during turn-taking in short calls exists in Japanese macaques Macaca fuscata. We observed exchanges of short calls such as grunts, girneys, and short, low coos during social interactions in a free-ranging group of Japanese macaques. We found that the median gap between the turns of two callers was 250 ms. Call intervals varied among individuals, suggesting that call intervals were not fixed among individuals. Solo call intervals were shorter than call intervals interrupted by responses from partners (i.e., exchanges) and longer than those between the partner's reply and the reply to that call, indicating that the monkeys did not just repeat calls at certain intervals irrespective of the social situation. The differences in call intervals during exchanged and solo call sequences were explained by the response interval of the partner, suggesting an adjustment of call timing according to the tempo of the partner's call utterance. These findings suggest that monkeys display dynamic temporal adjustment in a short time window, which is comparable with turn-taking in human speech.

Physiology of the cerebellum

The cerebellum is a central brain structure deeply integrated into major loops with the cerebral cortex, brainstem, and spinal cord. The cerebellum shows a complex regional organization consisting of modules with sagittal orientation. The cerebellum takes part in motor control and its lesions cause a movement incoordination syndrome called ataxia. Recent observations also imply involvement of the cerebellum in cognition and executive control, with an impact on pathologies like dyslexia and autism. The cerebellum operates as a forward controller learning to predict the precise timing of correlated events. The physiologic mechanisms of cerebellar functioning are still the object of intense research. The signals entering the cerebellum through the mossy fibers are processed in the granular layer and transmitted to Purkinje cells, while a collateral pathway activates the deep cerebellar nuclei (DCN). Purkinje cells in turn inhibit DCN, so that the cerebellar cortex operates as a side loop controlling the DCN. Learning is now known to occur through synaptic plasticity at multiple synapses in the granular layer, molecular layer, and DCN, extending the original concept of the Motor Learning Theory that predicted a single form of plasticity at the synapse between parallel fibers and Purkinje cells under the supervision of climbing fibers deriving from the inferior olive. Coordination derives from the precise regulation of timing and gain in the different cerebellar modules. The investigation of cerebellar dynamics using advanced physiologic recordings and computational models is now providing new clues on how the cerebellar network performs its internal computations.

Synaptic Activation of a Detailed Purkinje Cell Model Predicts Voltage-Dependent Control of Burst-Pause Responses in Active Dendrites

The dendritic processing in cerebellar Purkinje cells (PCs), which integrate synaptic inputs coming from hundreds of thousands granule cells and molecular layer interneurons, is still unclear. Here we have tested a leading hypothesis maintaining that the significant PC output code is represented by burst-pause responses (BPRs), by simulating PC responses in a biophysically detailed model that allowed to systematically explore a broad range of input patterns. BPRs were generated by input bursts and were more prominent in Zebrin positive than Zebrin negative (Z+ and Z-) PCs. Different combinations of parallel fiber and molecular layer interneuron synapses explained type I, II and III responses observed in vivo. BPRs were generated intrinsically by Ca-dependent K channel activation in the somato-dendritic compartment and the pause was reinforced by molecular layer interneuron inhibition. BPRs faithfully reported the duration and intensity of synaptic inputs, such that synaptic conductance tuned the number of spikes and release probability tuned their regularity in the millisecond range. Interestingly, the burst and pause of BPRs depended on the stimulated dendritic zone reflecting the different input conductance and local engagement of voltage-dependent channels. Multiple local inputs combined their actions generating complex spatio-temporal patterns of dendritic activity and BPRs. Thus, local control of intrinsic dendritic mechanisms by synaptic inputs emerges as a fundamental PC property in activity regimens characterized by bursting inputs from granular and molecular layer neurons.

Effect of Phase Response Curve Shape and Synaptic Driving Force on Synchronization of Coupled Neuronal Oscillators

The role of the phase response curve (PRC) shape on the synchrony of synaptically coupled oscillating neurons is examined. If the PRC is independent of the phase, because of the synaptic form of the coupling, synchrony is found to be stable for both excitatory and inhibitory coupling at all rates, whereas the antisynchrony becomes stable at low rates. A faster synaptic rise helps extend the stability of antisynchrony to higher rates. If the PRC is not constant but has a profile like that of a leaky integrate-and-fire model, then, in contrast to the earlier reports that did not include the voltage effects, mutual excitation could lead to stable synchrony provided the synaptic reversal potential is below the voltage level the neuron would have reached in the absence of the interaction and threshold reset. This level is controlled by the applied current and the leakage parameters. Such synchrony is contingent on significant phase response (that would result, for example, by a sharp PRC jump) occurring during the synaptic rising phase. The rising phase, however, does not contribute significantly if it occurs before the voltage spike reaches its peak. Then a stable near-synchronous state can still exist between type 1 PRC neurons if the PRC shows a left skewness in its shape. These results are examined comprehensively using perfect integrate-and-fire, leaky integrate-and-fire, and skewed PRC shapes under the assumption of the weakly coupled oscillator theory applied to synaptically coupled neuron models.

Qualitative changes in phase-response curve and synchronization at the saddle-node-loop bifurcation

Prominent changes in neuronal dynamics have previously been attributed to a specific switch in onset bifurcation, the Bogdanov-Takens (BT) point. This study unveils another, relevant and so far underestimated transition point: the saddle-node-loop bifurcation, which can be reached by several parameters, including capacitance, leak conductance, and temperature. This bifurcation turns out to induce even more drastic changes in synchronization than the BT transition. This result arises from a direct effect of the saddle-node-loop bifurcation on the limit cycle and hence spike dynamics. In contrast, the BT bifurcation exerts its immediate influence upon the subthreshold dynamics and hence only indirectly relates to spiking. We specifically demonstrate that the saddle-node-loop bifurcation (i) ubiquitously occurs in planar neuron models with a saddle node on invariant cycle onset bifurcation, and (ii) results in a symmetry breaking of the system's phase-response curve. The latter entails an increase in synchronization range in pulse-coupled oscillators, such as neurons. The derived bifurcation structure is of interest in any system for which a relaxation limit is admissible, such as Josephson junctions and chemical oscillators.

Phase-Dependent Modulation of Oscillatory Phase and Synchrony by Long-Lasting Depolarizing Inputs in Central Neurons

Oscillatory neural activities have been implicated in various types of information processing in the CNS. The procerebral (PC) lobe of the land mollusk Limax valentianus shows an ongoing oscillatory local field potential (LFP). Olfactory input increases both the frequency and spatial synchrony of the LFP oscillation by a nitric oxide (NO)-mediated mechanism, but how NO modulates the activity in a specific manner has been unclear. In the present study, we used electrical stimulation and NO uncaging to systematically analyze the response of the LFP oscillation and found phase-dependent effects on phase shifting and synchrony. The neurons that presumably release NO in the PC lobe preferentially fired at phases in which NO has a synchronizing effect, suggesting that the timing of NO release is regulated to induce a stereotyped response to natural sensory stimuli. The phase-response curve (PRC) describes the timing dependence of responses of an oscillatory system to external input. PRCs are usually constructed by recording the temporal shifts of the neural activity in response to brief electrical pulses. However, NO evokes a long-lasting depolarization persisting for several cycles of oscillation. The phase-response relationship obtained by NO stimulation was approximately the integral of the PRC. A similar relationship was also shown for regular firing of mouse cerebellar Purkinje cells receiving step depolarization, suggesting the generality of the results to oscillatory neural systems with highly distinct properties. These results indicate novel dynamic effects of long-lasting inputs on network oscillation and synchrony, which are based on simple and ubiquitous mechanisms.

Inverse Stochastic Resonance in Cerebellar Purkinje Cells

Purkinje neurons play an important role in cerebellar computation since their axons are the only projection from the cerebellar cortex to deeper cerebellar structures. They have complex internal dynamics, which allow them to fire spontaneously, display bistability, and also to be involved in network phenomena such as high frequency oscillations and travelling waves. Purkinje cells exhibit type II excitability, which can be revealed by a discontinuity in their f-I curves. We show that this excitability mechanism allows Purkinje cells to be efficiently inhibited by noise of a particular variance, a phenomenon known as inverse stochastic resonance (ISR). While ISR has been described in theoretical models of single neurons, here we provide the first experimental evidence for this effect. We find that an adaptive exponential integrate-and-fire model fitted to the basic Purkinje cell characteristics using a modified dynamic IV method displays ISR and bistability between the resting state and a repetitive activity limit cycle. ISR allows the Purkinje cell to operate in different functional regimes: the all-or-none toggle or the linear filter mode, depending on the variance of the synaptic input. We propose that synaptic noise allows Purkinje cells to quickly switch between these functional regimes. Using mutual information analysis, we demonstrate that ISR can lead to a locally optimal information transfer between the input and output spike train of the Purkinje cell. These results provide the first experimental evidence for ISR and suggest a functional role for ISR in cerebellar information processing.

How neurons generate output spikes in response to various combinations of inputs is a central issue in contemporary neuroscience. Due to their large dendritic tree and complex intrinsic properties, cerebellar Purkinje cells are an important model system to study this input-output transformation. Here we examine how noise can change the parameters of this transformation. In experiments we found that spike generation in Purkinje cells can be efficiently inhibited by noise of a particular amplitude. This effect is called inverse stochastic resonance (ISR) and has previously been described only in theoretical models of neurons. We explain the mechanism underlying ISR using a simple model matching the properties of experimentally characterized Purkinje cells. We found that ISR is present in Purkinje cells when the mean input current is near threshold for spike generation. ISR can be explained by the co-existence of resting and spiking solutions of the simple model. Changes of the input noise variance change the lifetime of these resting and spiking states, suggesting a mechanism for a tunable filter with long time constants implemented by a Purkinje cell population in the cerebellum. Finally, ISR leads to locally optimal information transfer from the input to the output of a Purkinje cell.

Time-invariant feed-forward inhibition of Purkinje cells in the cerebellar cortex in vivo

KEY POINTS

We performed extracellular recording of pairs of interneuron-Purkinje cells in vivo. A single interneuron produces a substantial, short-lasting, inhibition of Purkinje cells. Feed-forward inhibition is associated with characteristic asymmetric cross-correlograms. In vivo, Purkinje cell spikes only depend on the most recent synaptic activity.

ABSTRACT

Cerebellar molecular layer interneurons are considered to control the firing rate and spike timing of Purkinje cells. However, interactions between these cell types are largely unexplored in vivo. Using tetrodes, we performed simultaneous extracellular recordings of neighbouring Purkinje cells and molecular layer interneurons, presumably basket cells, in adult rats in vivo. The high levels of afferent synaptic activity encountered in vivo yield irregular spiking and reveal discharge patterns characteristic of feed-forward inhibition, thus suggesting an overlap of the afferent excitatory inputs between Purkinje cells and basket cells. Under conditions of intense background synaptic inputs, interneuron spikes exert a short-lasting inhibitory effect, delaying the following Purkinje cell spike by an amount remarkably independent of the Purkinje cell firing cycle. This effect can be explained by the short memory time of the Purkinje cell potential as a result of the intense incoming synaptic activity. Finally, we found little evidence for any involvement of the interneurons that we recorded with the cerebellar high-frequency oscillations promoting Purkinje cell synchrony. The rapid interactions between interneurons and Purkinje cells might be of particular importance in fine motor control because the inhibitory action of interneurons on Purkinje cells leads to deep cerebellar nuclear disinhibition and hence increased cerebellar output.

Enhanced Sensitivity to Rapid Input Fluctuations by Nonlinear Threshold Dynamics in Neocortical Pyramidal Neurons

The way in which single neurons transform input into output spike trains has fundamental consequences for network coding. Theories and modeling studies based on standard Integrate-and-Fire models implicitly assume that, in response to increasingly strong inputs, neurons modify their coding strategy by progressively reducing their selective sensitivity to rapid input fluctuations. Combining mathematical modeling with in vitro experiments, we demonstrate that, in L5 pyramidal neurons, the firing threshold dynamics adaptively adjust the effective timescale of somatic integration in order to preserve sensitivity to rapid signals over a broad range of input statistics. For that, a new Generalized Integrate-and-Fire model featuring nonlinear firing threshold dynamics and conductance-based adaptation is introduced that outperforms state-of-the-art neuron models in predicting the spiking activity of neurons responding to a variety of in vivo-like fluctuating currents. Our model allows for efficient parameter extraction and can be analytically mapped to a Generalized Linear Model in which both the input filter—describing somatic integration—and the spike-history filter—accounting for spike-frequency adaptation—dynamically adapt to the input statistics, as experimentally observed. Overall, our results provide new insights on the computational role of different biophysical processes known to underlie adaptive coding in single neurons and support previous theoretical findings indicating that the nonlinear dynamics of the firing threshold due to Na⁺-channel inactivation regulate the sensitivity to rapid input fluctuations.

Over the last decades, a variety of simplified spiking models have been shown to achieve a surprisingly high performance in predicting the neuronal responses to in vitro somatic current injections. Because of the complex adaptive behavior featured by cortical neurons, this success is however restricted to limited stimulus ranges: model parameters optimized for a specific input regime are often inappropriate to describe the response to input currents with different statistical properties. In the present study, a new spiking neuron model is introduced that captures single-neuron computation over a wide range of input statistics and explains different aspects of the neuronal dynamics within a single framework. Our results indicate that complex forms of single neuron adaptation are mediated by the nonlinear dynamics of the firing threshold and that the input-output transformation performed by cortical pyramidal neurons can be intuitively understood in terms of an enhanced Generalized Linear Model in which both the input filter and the spike-history filter adapt to the input statistics.

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.

Dynamics of on-off neural firing patterns and stochastic effects near a sub-critical Hopf bifurcation

On-off firing patterns, in which repetition of clusters of spikes are interspersed with epochs of subthreshold oscillations or quiescent states, have been observed in various nervous systems, but the dynamics of this event remain unclear. Here, we report that on-off firing patterns observed in three experimental models (rat sciatic nerve subject to chronic constrictive injury, rat CA1 pyramidal neuron, and rabbit blood pressure baroreceptor) appeared as an alternation between quiescent state and burst containing multiple period-1 spikes over time. Burst and quiescent state had various durations. The interspike interval (ISI) series of on-off firing pattern was suggested as stochastic using nonlinear prediction and autocorrelation function. The resting state was changed to a period-1 firing pattern via on-off firing pattern as the potassium concentration, static pressure, or depolarization current was changed. During the changing process, the burst duration of on-off firing pattern increased and the duration of the quiescent state decreased. Bistability of a limit cycle corresponding to period-1 firing and a focus corresponding to resting state was simulated near a sub-critical Hopf bifurcation point in the deterministic Morris-Lecar (ML) model. In the stochastic ML model, noise-induced transitions between the coexisting regimes formed an on-off firing pattern, which closely matched that observed in the experiment. In addition, noise-induced exponential change in the escape rate from the focus, and noise-induced coherence resonance were identified. The distinctions between the on-off firing pattern and stochastic firing patterns generated near three other types of bifurcations of equilibrium points, as well as other viewpoints on the dynamics of on-off firing pattern, are discussed. The results not only identify the on-off firing pattern as noise-induced stochastic firing pattern near a sub-critical Hopf bifurcation point, but also offer practical indicators to discriminate bifurcation types and neural excitability types.

On the firing rate dependency of the phase response curve of rat Purkinje neurons in vitro

Synchronous spiking during cerebellar tasks has been observed across Purkinje cells: however, little is known about the intrinsic cellular mechanisms responsible for its initiation, cessation and stability. The Phase Response Curve (PRC), a simple input-output characterization of single cells, can provide insights into individual and collective properties of neurons and networks, by quantifying the impact of an infinitesimal depolarizing current pulse on the time of occurrence of subsequent action potentials, while a neuron is firing tonically. Recently, the PRC theory applied to cerebellar Purkinje cells revealed that these behave as phase-independent integrators at low firing rates, and switch to a phase-dependent mode at high rates. Given the implications for computation and information processing in the cerebellum and the possible role of synchrony in the communication with its post-synaptic targets, we further explored the firing rate dependency of the PRC in Purkinje cells. We isolated key factors for the experimental estimation of the PRC and developed a closed-loop approach to reliably compute the PRC across diverse firing rates in the same cell. Our results show unambiguously that the PRC of individual Purkinje cells is firing rate dependent and that it smoothly transitions from phase independent integrator to a phase dependent mode. Using computational models we show that neither channel noise nor a realistic cell morphology are responsible for the rate dependent shift in the phase response curve.

A spiking network model of cerebellar Purkinje cells and molecular layer interneurons exhibiting irregular firing

While the anatomy of the cerebellar microcircuit is well-studied, how it implements cerebellar function is not understood. A number of models have been proposed to describe this mechanism but few emphasize the role of the vast network Purkinje cells (PKJs) form with the molecular layer interneurons (MLIs)—the stellate and basket cells. We propose a model of the MLI-PKJ network composed of simple spiking neurons incorporating the major anatomical and physiological features. In computer simulations, the model reproduces the irregular firing patterns observed in PKJs and MLIs in vitro and a shift toward faster, more regular firing patterns when inhibitory synaptic currents are blocked. In the model, the time between PKJ spikes is shown to be proportional to the amount of feedforward inhibition from an MLI on average. The two key elements of the model are: (1) spontaneously active PKJs and MLIs due to an endogenous depolarizing current, and (2) adherence to known anatomical connectivity along a parasagittal strip of cerebellar cortex. We propose this model to extend previous spiking network models of the cerebellum and for further computational investigation into the role of irregular firing and MLIs in cerebellar learning and function.

Predicting the responses of repetitively firing neurons to current noise

We used phase resetting methods to predict firing patterns of rat subthalamic nucleus (STN) neurons when their rhythmic firing was densely perturbed by noise. We applied sequences of contiguous brief (0.5–2 ms) current pulses with amplitudes drawn from a Gaussian distribution (10–100 pA standard deviation) to autonomously firing STN neurons in slices. Current noise sequences increased the variability of spike times with little or no effect on the average firing rate. We measured the infinitesimal phase resetting curve (PRC) for each neuron using a noise-based method. A phase model consisting of only a firing rate and PRC was very accurate at predicting spike timing, accounting for more than 80% of spike time variance and reliably reproducing the spike-to-spike pattern of irregular firing. An approximation for the evolution of phase was used to predict the effect of firing rate and noise parameters on spike timing variability. It quantitatively predicted changes in variability of interspike intervals with variation in noise amplitude, pulse duration and firing rate over the normal range of STN spontaneous rates. When constant current was used to drive the cells to higher rates, the PRC was altered in size and shape and accurate predictions of the effects of noise relied on incorporating these changes into the prediction. Application of rate-neutral changes in conductance showed that changes in PRC shape arise from conductance changes known to accompany rate increases in STN neurons, rather than the rate increases themselves. Our results show that firing patterns of densely perturbed oscillators cannot readily be distinguished from those of neurons randomly excited to fire from the rest state. The spike timing of repetitively firing neurons may be quantitatively predicted from the input and their PRCs, even when they are so densely perturbed that they no longer fire rhythmically.

Most neurons receive thousands of synaptic inputs per second. Each of these may be individually weak but collectively they shape the temporal pattern of firing by the postsynaptic neuron. If the postsynaptic neuron fires repetitively, its synaptic inputs need not directly trigger action potentials, but may instead control the timing of action potentials that would occur anyway. The phase resetting curve encapsulates the influence of an input on the timing of the next action potential, depending on its time of arrival. We measured the phase resetting curves of neurons in the subthalamic nucleus and used them to accurately predict the timing of action potentials in a phase model subjected to complex input patterns. A simple approximation to the phase model accurately predicted the changes in firing pattern evoked by dense patterns of noise pulses varying in amplitude and pulse duration, and by changes in firing rate. We also showed that the phase resetting curve changes systematically with changes in total neuron conductance, and doing so predicts corresponding changes in firing pattern. Our results indicate that the phase model may accurately represent the temporal integration of complex patterns of input to repetitively firing neurons.

In vivo cardiac phase response curve elucidates human respiratory heart rate variability

Recovering interaction of endogenous rhythms from observations is challenging, especially if a mathematical model explaining the behaviour of the system is unknown. The decisive information for successful reconstruction of the dynamics is the sensitivity of an oscillator to external influences, which is quantified by its phase response curve. Here we present a technique that allows the extraction of the phase response curve from a non-invasive observation of a system consisting of two interacting oscillators--in this case heartbeat and respiration--in its natural environment and under free-running conditions. We use this method to obtain the phase-coupling functions describing cardiorespiratory interactions and the phase response curve of 17 healthy humans. We show for the first time the phase at which the cardiac beat is susceptible to respiratory drive and extract the respiratory-related component of heart rate variability. This non-invasive method for the determination of phase response curves of coupled oscillators may find application in many scientific disciplines.

Coupled oscillator dynamics of vocal turn-taking in monkeys

Cooperation is central to human communication. The foundation of cooperative verbal communication is taking turns to speak. Given the universality of turn-taking, it is natural to ask how it evolved. We used marmoset monkeys to explore whether another primate species exhibits cooperative vocal communication by taking turns. Marmosets share with humans a cooperative breeding strategy and volubility. Cooperative care behaviors are thought to scaffold prosocial cognitive processes. Moreover, marmosets and other callitrichid primates are very vocal and readily exchange vocalizations with conspecifics. By measuring the natural statistics of marmoset vocal exchanges, we observed that they take turns in extended sequences and show that this vocal turn-taking has as its foundation dynamics characteristic of coupled oscillators--one that is similar to the dynamics proposed for human conversational turn-taking. As marmoset monkeys are on a different branch of the evolutionary tree that led to humans, our data demonstrate convergent evolution of vocal cooperation. Perhaps more importantly, our data offer a plausible alternative scenario to "gestural origin" hypotheses for how human cooperative vocal communication could have evolved.

Multipulse phase resetting curves

In this paper, we introduce and study systematically, in terms of phase response curves, the effect of dual-pulse excitation on the dynamics of an autonomous oscillator. Specifically, we test the deviations from linear summation of phase advances resulting from two small perturbations. We analytically derive a correction term, which generally appears for oscillators whose intrinsic dimensionality is >1. The nonlinear correction term is found to be proportional to the square of the perturbation. We demonstrate this effect in the Stuart-Landau model and in various higher dimensional neuronal models. This deviation from the superposition principle needs to be taken into account in studies of networks of pulse-coupled oscillators. Further, this deviation could be used in the verification of oscillator models via a dual-pulse excitation.

Spectral reconstruction of phase response curves reveals the synchronization properties of mouse globus pallidus neurons

The propensity of a neuron to synchronize is captured by its infinitesimal phase response curve (iPRC). Determining whether an iPRC is biphasic, meaning that small depolarizing perturbations can actually delay the next spike, if delivered at appropriate phases, is a daunting experimental task because negative lobes in the iPRC (unlike positive ones) tend to be small and may be occluded by the normal discharge variability of a neuron. To circumvent this problem, iPRCs are commonly derived from numerical models of neurons. Here, we propose a novel and natural method to estimate the iPRC by direct estimation of its spectral modes. First, we show analytically that the spectral modes of the iPRC of an arbitrary oscillator are readily measured by applying weak harmonic perturbations. Next, applying this methodology to biophysical neuronal models, we show that a low-dimensional spectral reconstruction is sufficient to capture the structure of the iPRC. This structure was preserved reasonably well even with added physiological scale jitter in the neuronal models. To validate the methodology empirically, we applied it first to a low-noise electronic oscillator with a known design and then to cortical pyramidal neurons, recorded in whole cell configuration, that are known to possess a monophasic iPRC. Finally, using the methodology in conjunction with perforated-patch recordings from pallidal neurons, we show, in contrast to recent modeling studies, that these neurons have biphasic somatic iPRCs. Biphasic iPRCs would cause lateral somatically targeted pallidal inhibition to desynchronize pallidal neurons, providing a plausible explanation for their lack of synchrony in vivo.

An integrator circuit in cerebellar cortex

The brain builds dynamic models of the body and the outside world to predict the consequences of actions and stimuli. A well-known example is the oculomotor integrator, which anticipates the position-dependent elasticity forces acting on the eye ball by mathematically integrating over time oculomotor velocity commands. Many models of neural integration have been proposed, based on feedback excitation, lateral inhibition or intrinsic neuronal nonlinearities. We report here that a computational model of the cerebellar cortex, a structure thought to implement dynamic models, reveals a hitherto unrecognized integrator circuit. In this model, comprising Purkinje cells, molecular layer interneurons and parallel fibres, Purkinje cells were able to generate responses lasting more than 10 s, to which both neuronal and network mechanisms contributed. Activation of the somatic fast sodium current by subthreshold voltage fluctuations was able to maintain pulse-evoked graded persistent activity, whereas lateral inhibition among Purkinje cells via recurrent axon collaterals further prolonged the responses to step and sine wave stimulation. The responses of Purkinje cells decayed with a time-constant whose value depended on their baseline spike rate, with integration vanishing at low (< 1 per s) and high rates (> 30 per s). The model predicts that the apparently fast circuit of the cerebellar cortex may control the timing of slow processes without having to rely on sensory feedback. Thus, the cerebellar cortex may contain an adaptive temporal integrator, with the sensitivity of integration to the baseline spike rate offering a potential mechanism of plasticity of the response time-constant.

Hippocampal CA1 pyramidal neurons exhibit type 1 phase-response curves and type 1 excitability

Phase-resetting properties of neurons determine their functionality as integrators (type 1) vs. resonators (type 2), as well as their synchronization tendencies. We introduce a novel, bias-correction method to estimate the infinitesimal phase-resetting curve (iPRC) and confirm type 1 excitability in hippocampal pyramidal CA1 neurons in vitro by two independent methods. First, PRCs evoked using depolarizing pulses consisted only of advances, consistent with type 1. Second, the frequency/current (f/I) plots showed no minimum frequency, again consistent with type 1. Type 1 excitability was also confirmed by the absence of a resonant peak in the interspike interval histograms derived from the f/I data. The PRC bias correction assumed that the distribution of noisy phase resetting is truncated, because an input cannot advance a spike to a point in time before the input (the causal limit) and successfully removed the statistical bias for delays in the null PRC in response to zero-magnitude input by computing the phase resetting as the mean of the untruncated distribution. The PRC for depolarization peaked at late phases and decreased to zero by the end of the cycle, whereas delays observed in response to hyperpolarization increased monotonically. The bias correction did not affect this difference in shape, which was due instead to the causal limit obscuring the iPRC for depolarization but not hyperpolarization. Our results suggest that weak periodic hyperpolarizing drive can theoretically entrain CA1 pyramidal neurons at any phase but that strong excitation will preferentially phase-lock them with zero time lag.

Phase response curves of subthalamic neurons measured with synaptic input and current injection

Infinitesimal phase response curves (iPRCs) provide a simple description of the response of repetitively firing neurons and may be used to predict responses to any pattern of synaptic input. Their simplicity makes them useful for understanding the dynamics of neurons when certain conditions are met. For example, the sizes of evoked phase shifts should scale linearly with stimulus strength, and the form of the iPRC should remain relatively constant as firing rate varies. We measured the PRCs of rat subthalamic neurons in brain slices using corticosubthalamic excitatory postsynaptic potentials (EPSPs; mediated by both AMPA- and NMDA-type receptors) and injected current pulses and used them to calculate the iPRC. These were relatively insensitive to both the size of the stimulus and the cell's firing rate, suggesting that the iPRC can predict the response of subthalamic nucleus cells to extrinsic inputs. However, the iPRC calculated using EPSPs differed from that obtained using current pulses. EPSPs (normalized for charge) were much more effective at altering the phase of subthalamic neurons than current pulses. The difference was not attributable to the extended time course of NMDA receptor-mediated currents, being unaffected by blockade of NMDA receptors. The iPRC provides a good description of subthalamic neurons' response to input, but iPRCs are best estimated using synaptic inputs rather than somatic current injection.

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.

Identifying type I excitability using dynamics of stochastic neural firing patterns

The stochastic firing patterns are simulated near saddle-node bifurcation on an invariant cycle corresponding to type I excitability in stochastic Morris-Lecar model. In absence of external periodic signal, the stochastic firing manifests continuous distribution in ISI histogram (ISIH), whose amplitude at first increases sharply and then decreases exponentially. In presence of the external periodic signal, stochastic firing patterns appear as two cases of integer multiple firing with multiple discrete peaks in ISIH. One manifests perfect exponential decay in all peaks and the other imperfect exponential decay except a lower first peak. These stochastic firing patterns simulated with or without external periodic signal can be demonstrated in the experiments on rat hippocampal CA1 pyramidal neurons. The exponential decay laws in the multiple peaks are also acquired using probability analysis method. The perfect decay law is determined by the independent characteristic within the firing while the imperfect decay law is from the inhibitory effect. In addition, the stochastic firing patterns corresponding to type I excitability are compared to those of type II excitability. The results not only reveal the dynamics of stochastic firing patterns with or without external signal corresponding to type I excitability, but also provide practical indicators to availably identify type I excitability.

Phase response theory extended to nonoscillatory network components

New tools for analysis of oscillatory networks using phase response theory (PRT) under the assumption of pulsatile coupling have been developed steadily since the 1980s, but none have yet allowed for analysis of mixed systems containing nonoscillatory elements. This caveat has excluded the application of PRT to most real systems, which are often mixed. We show that a recently developed tool, the functional phase resetting curve (fPRC), provides a serendipitous benefit: it allows incorporation of nonoscillatory elements into systems of oscillators where PRT can be applied. We validate this method in a model system of neural oscillators and a biological system, the pyloric network of crustacean decapods.

Neural circuits of the cerebellum: hypothesis for function

The rapid growth of cerebellar research is going to clarify several aspects of cellular and circuit physiology. However, the concepts about cerebellar mechanisms of function are still largely related to clinical observations and to models elaborated before the last discoveries appeared. In this paper, the major issues are revisited, suggesting that previous concepts can now be refined and modified. The cerebellum is fundamentally involved in timing and in controlling the ordered and precise execution of motor sequences. The fast reaction of the cerebellum to the inputs is sustained by specific cellular mechanisms ensuring precision on the millisecond scale. These include burst-burst reconversion in the granular layer and instantaneous frequency modulation on the 100-Hz band in Purkinje and deep cerebellar nuclei cells. Precisely timed signals can be used for perceptron operations in Purkinje cells and to establish appropriate correlations with climbing fiber signals inducing learning at parallel fiber synapses. In the granular layer, plasticity turns out to be instrumental to timing, providing a conceptual solution to the discrepancy between cerebellar learning and timing. The granular layer sub-circuit can be tuned by long-term synaptic plasticity and synaptic inhibition to delay the incoming signals over a 100-ms range. For longer sequences, large circuit sections can be entrained into coherent activity in 100-ms cycles. These dynamic aspects, which have not been accounted for by original theories, could in fact represent the essence of cerebellar functioning. It is suggested that the cerebellum can, in this way, operate the realignment of temporally incongruent signals, allowing their binding and pattern recognition in Purkinje cells. The demonstration of these principles, their behavioral relevance and their relationship with internal model theories represent a challenge for future cerebellar research.

Short conduction delays cause inhibition rather than excitation to favor synchrony in hybrid neuronal networks of the entorhinal cortex

How stable synchrony in neuronal networks is sustained in the presence of conduction delays is an open question. The Dynamic Clamp was used to measure phase resetting curves (PRCs) for entorhinal cortical cells, and then to construct networks of two such neurons. PRCs were in general Type I (all advances or all delays) or weakly type II with a small region at early phases with the opposite type of resetting. We used previously developed theoretical methods based on PRCs under the assumption of pulsatile coupling to predict the delays that synchronize these hybrid circuits. For excitatory coupling, synchrony was predicted and observed only with no delay and for delays greater than half a network period that cause each neuron to receive an input late in its firing cycle and almost immediately fire an action potential. Synchronization for these long delays was surprisingly tight and robust to the noise and heterogeneity inherent in a biological system. In contrast to excitatory coupling, inhibitory coupling led to antiphase for no delay, very short delays and delays close to a network period, but to near-synchrony for a wide range of relatively short delays. PRC-based methods show that conduction delays can stabilize synchrony in several ways, including neutralizing a discontinuity introduced by strong inhibition, favoring synchrony in the case of noisy bistability, and avoiding an initial destabilizing region of a weakly type II PRC. PRCs can identify optimal conduction delays favoring synchronization at a given frequency, and also predict robustness to noise and heterogeneity.

Individual oscillators, such as pendulum-based clocks and fireflies, can spontaneously organize into a coherent, synchronized entity with a common frequency. Neurons can oscillate under some circumstances, and can synchronize their firing both within and across brain regions. Synchronized assemblies of neurons are thought to underlie cognitive functions such as recognition, recall, perception and attention. Pathological synchrony can lead to epilepsy, tremor and other dynamical diseases, and synchronization is altered in most mental disorders. Biological neurons synchronize despite conduction delays, heterogeneous circuit composition, and noise. In biological experiments, we built simple networks in which two living neurons could interact via a computer in real time. The computer precisely controlled the nature of the connectivity and the length of the communication delays. We characterized the synchronization tendencies of individual, isolated oscillators by measuring how much a single input delivered by the computer transiently shortened or lengthened the cycle period of the oscillation. We then used this information to correctly predict the strong dependence of the coordination pattern of the firing of the component neurons on the length of the communication delays. Upon this foundation, we can begin to build a theory of the basic principles of synchronization in more complex brain circuits.

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.

Inclusion of noise in iterated firing time maps based on the phase response curve

The infinitesimal phase response curve (PRC) of a neural oscillator to a weak input is a powerful predictor of network dynamics; however, many networks have strong coupling and require direct measurement of the PRC for strong inputs under the assumption of pulsatile coupling. We incorporate measured noise levels in firing time maps constructed from PRCs to predict phase-locked modes of activity, phase difference, and locking strength in 78 heterogeneous hybrid networks of 2 neurons constructed using the dynamic clamp. We show that noise may either destroy or stabilize a phase-locked mode of activity.