Phase-response curves give the responses of neurons to transient inputs

Revealing the mechanism and function underlying pairwise temporal coupling in collective motion

Coordinated motion in animal groups has predominantly been studied with a focus on spatial interactions, such as how individuals position and orient themselves relative to one another. Temporal aspects have, by contrast, received much less attention. Here, by studying pairwise interactions in juvenile zebrafish (Danio rerio)-including using immersive volumetric virtual reality (VR) with which we can directly test models of social interactions in situ-we reveal that there exists a rhythmic out-of-phase (i.e., an alternating) temporal coordination dynamic. We find that reciprocal (bi-directional) feedback is both necessary and sufficient to explain this emergent coupling. Beyond a mechanistic understanding, we find, both from VR experiments and analysis of freely swimming pairs, that temporal coordination considerably improves spatial responsiveness, such as to changes in the direction of motion of a partner. Our findings highlight the synergistic role of spatial and temporal coupling in facilitating effective communication between individuals on the move.

Conditioning and pseudoconditioning differently change intrinsic excitability of inhibitory interneurons in the neocortex

Many studies indicate a broad role of various classes of GABAergic interneurons in the processes related to learning. However, little is known about how the learning process affects intrinsic excitability of specific classes of interneurons in the neocortex. To determine this, we employed a simple model of conditional learning in mice where vibrissae stimulation was used as a conditioned stimulus and a tail shock as an unconditioned one. In vitro whole-cell patch-clamp recordings showed an increase in intrinsic excitability of low-threshold spiking somatostatin-expressing interneurons (SST-INs) in layer 4 (L4) of the somatosensory (barrel) cortex after the conditioning paradigm. In contrast, pseudoconditioning reduced intrinsic excitability of SST-LTS, parvalbumin-expressing interneurons (PV-INs), and vasoactive intestinal polypeptide-expressing interneurons (VIP-INs) with accommodating pattern in L4 of the barrel cortex. In general, increased intrinsic excitability was accompanied by narrowing of action potentials (APs), whereas decreased intrinsic excitability coincided with AP broadening. Altogether, these results show that both conditioning and pseudoconditioning lead to plastic changes in intrinsic excitability of GABAergic interneurons in a cell-specific manner. In this way, changes in intrinsic excitability can be perceived as a common mechanism of learning-induced plasticity in the GABAergic system.

A framework for macroscopic phase-resetting curves for generalised spiking neural networks

Brain rhythms emerge from synchronization among interconnected spiking neurons. Key properties of such rhythms can be gleaned from the phase-resetting curve (PRC). Inferring the PRC and developing a systematic phase reduction theory for large-scale brain rhythms remains an outstanding challenge. Here we present a theoretical framework and methodology to compute the PRC of generic spiking networks with emergent collective oscillations. We adopt a renewal approach where neurons are described by the time since their last action potential, a description that can reproduce the dynamical feature of many cell types. For a sufficiently large number of neurons, the network dynamics are well captured by a continuity equation known as the refractory density equation. We develop an adjoint method for this equation giving a semi-analytical expression of the infinitesimal PRC. We confirm the validity of our framework for specific examples of neural networks. Our theoretical framework can link key biological properties at the individual neuron scale and the macroscopic oscillatory network properties. Beyond spiking networks, the approach is applicable to a broad class of systems that can be described by renewal processes.

The formation of oscillatory neuronal assemblies at the network level has been hypothesized to be fundamental to many cognitive and motor functions. One prominent tool to understand the dynamics of oscillatory activity response to stimuli, and hence the neural code for which it is a substrate, is a nonlinear measure called Phase-Resetting Curve (PRC). At the network scale, the PRC defines the measure of how a given synaptic input perturbs the timing of next upcoming volley of spike assemblies: either advancing or delaying this timing. As a further application, one can use PRCs to make unambiguous predictions about whether communicating networks of neurons will phase-lock as it is often observed across the cortical areas and what would be this stable phase-configuration: synchronous, asynchronous or with asymmetric phase-shifts. The latter configuration also implies a preferential flow of information form the leading network to the follower, thereby giving causal signatures of directed functional connectivity. Because of the key position of the PRC in studying synchrony, information flow and entrainment to external forcing, it is crucial to move toward a theory that allows to compute the PRCs of network-wide oscillations not only for a restricted class of models, as has been done in the past, but to network descriptions that are generalized and can reflect flexibly single cell properties. In this manuscript, we tackle this issue by showing how the PRC for network oscillations can be computed using the adjoint systems of partial differential equations that define the dynamics of the neural activity density.

Temperature elevations can induce switches to homoclinic action potentials that alter neural encoding and synchronization

Almost seventy years after the discovery of the mechanisms of action potential generation, some aspects of their computational consequences are still not fully understood. Based on mathematical modeling, we here explore a type of action potential dynamics – arising from a saddle-node homoclinic orbit bifurcation - that so far has received little attention. We show that this type of dynamics is to be expected by specific changes in common physiological parameters, like an elevation of temperature. Moreover, we demonstrate that it favours synchronization patterns in networks – a feature that becomes particularly prominent when system parameters change such that homoclinic spiking is induced. Supported by in-vitro hallmarks for homoclinic spikes in the rodent brain, we hypothesize that the prevalence of homoclinic spikes in the brain may be underestimated and provide a missing link between the impact of biophysical parameters on abrupt transitions between asynchronous and synchronous states of electrical activity in the brain.

The intrinsic dynamics of neurons, in particular the generation action potentials, can impact neural network states and processes of encoding information. The authors demonstrate how the elevation of temperature induces a type of action potential dynamics that favors synchronization patterns in neural networks.

Developing control-theoretic objectives for large-scale brain dynamics and cognitive enhancement

The development of technologies for brain stimulation provides a means for scientists and clinicians to directly actuate the brain and nervous system. Brain stimulation has shown intriguing potential in terms of modifying particular symptom clusters in patients and behavioral characteristics of subjects. The stage is thus set for optimization of these techniques and the pursuit of more nuanced stimulation objectives, including the modification of complex cognitive functions such as memory and attention. Control theory and engineering will play a key role in the development of these methods, guiding computational and algorithmic strategies for stimulation. In particular, realizing this goal will require new development of frameworks that allow for controlling not only brain activity, but also latent dynamics that underlie neural computation and information processing. In the current opinion, we review recent progress in brain stimulation and outline challenges and potential research pathways associated with exogenous control of cognitive function.

Phase response analyses support a relaxation oscillator model of locomotor rhythm generation in Caenorhabditis elegans

Neural circuits coordinate with muscles and sensory feedback to generate motor behaviors appropriate to an animal's environment. In C. elegans, the mechanisms by which the motor circuit generates undulations and modulates them based on the environment are largely unclear. We quantitatively analyzed C. elegans locomotion during free movement and during transient optogenetic muscle inhibition. Undulatory movements were highly asymmetrical with respect to the duration of bending and unbending during each cycle. Phase response curves induced by brief optogenetic inhibition of head muscles showed gradual increases and rapid decreases as a function of phase at which the perturbation was applied. A relaxation oscillator model based on proprioceptive thresholds that switch the active muscle moment was developed and is shown to quantitatively agree with data from free movement, phase responses, and previous results for gait adaptation to mechanical loadings. Our results suggest a neuromuscular mechanism underlying C. elegans motor pattern generation within a compact circuit.

Acetylcholine-gated current translates wake neuronal firing rate information into a spike timing-based code in Non-REM sleep, stabilizing neural network dynamics during memory consolidation

Sleep is critical for memory consolidation, although the exact mechanisms mediating this process are unknown. Combining reduced network models and analysis of in vivo recordings, we tested the hypothesis that neuromodulatory changes in acetylcholine (ACh) levels during non-rapid eye movement (NREM) sleep mediate stabilization of network-wide firing patterns, with temporal order of neurons’ firing dependent on their mean firing rate during wake. In both reduced models and in vivo recordings from mouse hippocampus, we find that the relative order of firing among neurons during NREM sleep reflects their relative firing rates during prior wake. Our modeling results show that this remapping of wake-associated, firing frequency-based representations is based on NREM-associated changes in neuronal excitability mediated by ACh-gated potassium current. We also show that learning-dependent reordering of sequential firing during NREM sleep, together with spike timing-dependent plasticity (STDP), reconfigures neuronal firing rates across the network. This rescaling of firing rates has been reported in multiple brain circuits across periods of sleep. Our model and experimental data both suggest that this effect is amplified in neural circuits following learning. Together our data suggest that sleep may bias neural networks from firing rate-based towards phase-based information encoding to consolidate memories.

We show that neuromodulatory changes during non-rapid eye movement (NREM) sleep generate stable spike timing relationships between neurons, the ordering of which reflects the neurons’ relative firing rates during wake. Learning-dependent ordering of firing in the hippocampus during NREM, acting in tandem with spike timing-dependent plasticity, reconfigures neuronal firing rates across the hippocampal network. This “rescaling” of neuronal firing rates has recently been reported in multiple brain circuits across periods of sleep. Together, our results suggest that the brain is remapping frequency-biased representations of information formed during wake into timing biased-representations during NREM sleep.

Interspike interval correlations in neuron models with adaptation and correlated noise

The generation of neural action potentials (spikes) is random but nevertheless may result in a rich statistical structure of the spike sequence. In particular, contrary to the popular renewal assumption of theoreticians, the intervals between adjacent spikes are often correlated. Experimentally, different patterns of interspike-interval correlations have been observed and computational studies have identified spike-frequency adaptation and correlated noise as the two main mechanisms that can lead to such correlations. Analytical studies have focused on the single cases of either correlated (colored) noise or adaptation currents in combination with uncorrelated (white) noise. For low-pass filtered noise or adaptation, the serial correlation coefficient can be approximated as a single geometric sequence of the lag between the intervals, providing an explanation for some of the experimentally observed patterns. Here we address the problem of interval correlations for a widely used class of models, multidimensional integrate-and-fire neurons subject to a combination of colored and white noise sources and a spike-triggered adaptation current. Assuming weak noise, we derive a simple formula for the serial correlation coefficient, a sum of two geometric sequences, which accounts for a large class of correlation patterns. The theory is confirmed by means of numerical simulations in a number of special cases including the leaky, quadratic, and generalized integrate-and-fire models with colored noise and spike-frequency adaptation. Furthermore we study the case in which the adaptation current and the colored noise share the same time scale, corresponding to a slow stochastic population of adaptation channels; we demonstrate that our theory can account for a nonmonotonic dependence of the correlation coefficient on the channel’s time scale. Another application of the theory is a neuron driven by network-noise-like fluctuations (green noise). We also discuss the range of validity of our weak-noise theory and show that by changing the relative strength of white and colored noise sources, we can change the sign of the correlation coefficient. Finally, we apply our theory to a conductance-based model which demonstrates its broad applicability.

The elementary processing units in the central nervous system are neurons that transmit information by short electrical pulses, so called action potentials or spikes. The generation of the action potential is a random process that can be shaped by correlated fluctuations (colored noise) and by adaptation. A consequence of these two ubiquitous features is that the successive time intervals between spikes, the interspike intervals, are not independent but correlated. As these correlations can significantly improve information transmission and weak-signal detection, it is an important task to develop analytical approaches to these statistics for well-established computational models. Here we present a theory of interval correlations for a widely used class of integrate-and-fire models endowed with an adaptation mechanism and subject to correlated fluctuations. We demonstrate which patterns of interval correlations can be expected from the interplay of colored noise, adaptation and intrinsic nonlinear dynamics.

On the Role of Arkypallidal and Prototypical Neurons for Phase Transitions in the External Pallidum

The external pallidum (globus pallidus pars externa [GPe]) plays a central role for basal ganglia functions and dynamics and, consequently, has been included in most computational studies of the basal ganglia. These studies considered the GPe as a homogeneous neural population. However, experimental studies have shown that the GPe contains at least two distinct cell types (prototypical and arkypallidal cells). In this work, we provide in silico insight into how pallidal heterogeneity modulates dynamic regimes inside the GPe and how they affect the GPe response to oscillatory input. We derive a mean-field model of the GPe system from a microscopic spiking neural network of recurrently coupled prototypical and arkypallidal neurons. Using bifurcation analysis, we examine the influence of dopamine-dependent changes of intrapallidal connectivity on the GPe dynamics. We find that increased self-inhibition of prototypical cells can induce oscillations, whereas increased inhibition of prototypical cells by arkypallidal cells leads to the emergence of a bistable regime. Furthermore, we show that oscillatory input to the GPe, arriving from striatum, leads to characteristic patterns of cross-frequency coupling observed at the GPe. Based on these findings, we propose two different hypotheses of how dopamine depletion at the GPe may lead to phase-amplitude coupling between the parkinsonian beta rhythm and a GPe-intrinsic γ rhythm. Finally, we show that these findings generalize to realistic spiking neural networks of sparsely coupled Type I excitable GPe neurons.SIGNIFICANCE STATEMENT Our work provides (1) insight into the theoretical implications of a dichotomous globus pallidus pars externa (GPe) organization, and (2) an exact mean-field model that allows for future investigations of the relationship between GPe spiking activity and local field potential fluctuations. We identify the major phase transitions that the GPe can undergo when subject to static or periodic input and link these phase transitions to the emergence of synchronized oscillations and cross-frequency coupling in the basal ganglia. Because of the close links between our model and experimental findings on the structure and dynamics of prototypical and arkypallidal cells, our results can be used to guide both experimental and computational studies on the role of the GPe for basal ganglia dynamics in health and disease.

Efficient and robust coding in heterogeneous recurrent networks

Cortical networks show a large heterogeneity of neuronal properties. However, traditional coding models have focused on homogeneous populations of excitatory and inhibitory neurons. Here, we analytically derive a class of recurrent networks of spiking neurons that close to optimally track a continuously varying input online, based on two assumptions: 1) every spike is decoded linearly and 2) the network aims to reduce the mean-squared error between the input and the estimate. From this we derive a class of predictive coding networks, that unifies encoding and decoding and in which we can investigate the difference between homogeneous networks and heterogeneous networks, in which each neurons represents different features and has different spike-generating properties. We find that in this framework, 'type 1' and 'type 2' neurons arise naturally and networks consisting of a heterogeneous population of different neuron types are both more efficient and more robust against correlated noise. We make two experimental predictions: 1) we predict that integrators show strong correlations with other integrators and resonators are correlated with resonators, whereas the correlations are much weaker between neurons with different coding properties and 2) that 'type 2' neurons are more coherent with the overall network activity than 'type 1' neurons.

Evidence of two modes of spiking evoked in human firing motoneurones by Ia afferent electrical stimulation

Neurone firing behaviour is a result of complex interaction between synaptic inputs and cellular intrinsic properties. Intriguing firing behaviour, delayed spiking, was shown in some neurones, in particular, in cat neocortical neurones and rat pyramidal hippocampal neurones. In contrast, the similar spiking mode was not reported for animal spinal motoneurones. In the present study, an attempt was made to look for possible evidence of delayed spiking in human motoneurones firing within the low-frequency, sub-primary range, characteristic for voluntary muscle contractions and postural tasks. Forty-seven firing motor units (MUs) were analyzed in ten experiments on three muscles (the flexor carpi ulnaris, the tibialis anterior, and the abductor pollicis brevis) in four healthy humans. Single MUs were activated by gentle voluntary muscle contractions. MU peri-stimulus time histograms, durations of inter-spike intervals, and motoneurone excitability changes within a target interspike interval were analyzed. It was found that during testing the firing motoneurone excitability by small, transient excitatory Ia afferent volley, depending firstly on volley timing within a target interspike interval and excitatory volley strength, the same motoneurone displayed either the direct short-latency response (the H-reflex) or the delayed response (with prolonged and variable latency). Thus, the findings, for the first time, provide evidence for a possibility of two modes of spiking in firing motoneurones. Methods of the estimation of delayed responses and their possible functional significance are discussed. It is emphasized that, for understanding of this issue, the integration of data from studies on experimental animals and humans is desirable.

Interactions of multiple rhythms in a biophysical network of neurons

Neural oscillations, including rhythms in the beta1 band (12-20 Hz), are important in various cognitive functions. Often neural networks receive rhythmic input at frequencies different from their natural frequency, but very little is known about how such input affects the network's behavior. We use a simplified, yet biophysical, model of a beta1 rhythm that occurs in the parietal cortex, in order to study its response to oscillatory inputs. We demonstrate that a cell has the ability to respond at the same time to two periodic stimuli of unrelated frequencies, firing in phase with one, but with a mean firing rate equal to that of the other. We show that this is a very general phenomenon, independent of the model used. We next show numerically that the behavior of a different cell, which is modeled as a high-dimensional dynamical system, can be described in a surprisingly simple way, owing to a reset that occurs in the state space when the cell fires. The interaction of the two cells leads to novel combinations of properties for neural dynamics, such as mode-locking to an input without phase-locking to it.

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.

Global phase-amplitude description of oscillatory dynamics via the parameterization method

In this paper, we use the parameterization method to provide a complete description of the dynamics of an n-dimensional oscillator beyond the classical phase reduction. The parameterization method allows us, via efficient algorithms, to obtain a parameterization of the attracting invariant manifold of the limit cycle in terms of the phase-amplitude variables. The method has several advantages. It provides analytically a Fourier-Taylor expansion of the parameterization up to any order, as well as a simplification of the dynamics that allows for a numerical globalization of the manifolds. Thus, one can obtain the local and global isochrons and isostables, including the slow attracting manifold, up to high accuracy, which offer a geometrical portrait of the oscillatory dynamics. Furthermore, it provides straightforwardly the infinitesimal phase and amplitude response functions, that is, the extended infinitesimal phase and amplitude response curves, which monitor the phase and amplitude shifts beyond the asymptotic state. Thus, the methodology presented yields an accurate description of the phase dynamics for perturbations not restricted to the limit cycle but to its attracting invariant manifold. Finally, we explore some strategies to reduce the dimension of the dynamics, including the reduction of the dynamics to the slow stable submanifold. We illustrate our methods by applying them to different three-dimensional single neuron and neural population models in neuroscience.

Indirect pathway control of firing rate and pattern in the substantia nigra pars reticulata

Unitary pallido-nigral synaptic currents were measured using optogenetic stimulation, which activated up to three unitary synaptic inputs to each substantia nigra pars reticulata (SNr) cell. Episodic barrages of synaptic conductances were generated based on in vivo firing patterns of globus pallidus pars externa (GPe) cells and applied to SNr cells using conductance clamp. Barrage inputs were compared to continuous step conductances with the same mean. Barrage inputs and steps both slowed SNr neuron firing and produced disinhibition responses seen in peristimulus histograms. Barrages were less effective than steps at producing inhibition and disinhibition responses. Barrages, but not steps, produced irregular firing during the inhibitory response. Phase models of SNr neurons were constructed from their phase-resetting curves. The phase models reproduced the inhibition and disinhibition responses to the same inputs applied to the neurons. The disinhibition response did not require rebound currents but arose from reset of the cells' oscillation. The differences in firing rate and irregularity in response to barrage and step inhibition resulted from the high sensitivity of SNr neurons to inhibition at late phases in their intrinsic oscillation. During step inhibition, cells continued rhythmic firing at a reduced rate. During barrages, brief bouts of intense inhibition stalled the cells' phase evolution late in their cycle, close to firing, and even very brief respites from inhibition rapidly released single action potentials. The SNr cell firing pattern reflected the fine structure of the synaptic barrage from GPe, as well as its onset and offset.NEW & NOTEWORTHY The pallido-nigral pathway connects the striatum to spontaneously active basal ganglia output neurons in the substantia nigra. Each substantia nigra neuron receives powerful inhibitory synaptic connections from a small group of globus pallidus cells and may fire during pauses in pallidal activity. Despite lacking any hyperpolarization-activated rebound currents, they fire quickly to even brief pauses in the pallido-nigral inhibition. The mechanism of their rapid disinhibitory response is explained by features of their phase-resetting curves.

Synchronization of Electrically Coupled Resonate-and-Fire Neurons

Electrical coupling between neurons is broadly present across brain areas and is typically assumed to synchronize network activity. However, intrinsic properties of the coupled cells can complicate this simple picture. Many cell types with electrical coupling show a diversity of post-spike subthreshold fluctuations, often linked to subthreshold resonance, which are transmitted through electrical synapses in addition to action potentials. Using the theory of weakly coupled oscillators, we explore the effect of both subthreshold and spike-mediated coupling on synchrony in small networks of electrically coupled resonate-and-fire neurons, a hybrid neuron model with damped subthreshold oscillations and a range of post-spike voltage dynamics. We calculate the phase response curve using an extension of the adjoint method that accounts for the discontinuous post-spike reset rule. We find that both spikes and subthreshold fluctuations can jointly promote synchronization. The subthreshold contribution is strongest when the voltage exhibits a significant post-spike elevation in voltage, or plateau potential. Additionally, we show that the geometry of trajectories approaching the spiking threshold causes a "reset-induced shear" effect that can oppose synchrony in the presence of network asymmetry, despite having no effect on the phase-locking of symmetrically coupled pairs.

Frequency-dependent entrainment of striatal fast-spiking interneurons

Striatal fast-spiking interneurons (FSIs) fire in variable-length runs of action potentials at 20-200 spikes/s separated by pauses. In vivo, or with fluctuating applied current, both runs and pauses become briefer and more variable. During runs, spikes are entrained specifically to gamma-frequency components of the input fluctuations. We stimulated parvalbumin-expressing striatal FSIs in mouse brain slices with broadband noise currents added to direct current steps and measured spike entrainment across all frequencies. As the constant current level was increased, FSIs produced longer runs and showed sharper frequency tuning, with best entrainment at the stimulus frequency matching their intrarun firing rate. We separated the contributions of previous spikes from that of the fluctuating stimulus, revealing a strong contribution of previous action potentials to gamma-frequency entrainment. In contrast, after subtraction of the effect inherited from the previous spike, the remaining stimulus contribution to spike generation was less sharply tuned, showing a larger contribution of lower frequencies. The frequency specificity of entrainment within a run was reproduced with a phase resetting model based on experimentally measured phase resetting curves of the same FSIs. In the model, broadly tuned phase entrainment for the first spike in a run evolved into sharply tuned gamma entrainment over the next few spikes. The data and modeling results indicate that for FSIs firing in brief runs and pauses firing within runs is entrained by gamma-frequency components of the input, whereas the onset timing of runs may be sensitive to a wider range of stimulus frequency components.NEW & NOTEWORTHY Specific types of neurons entrain their spikes to particular oscillation frequencies in their synaptic input. This entrainment is commonly understood in terms of the subthreshold voltage response, but how this translates to spiking is not clear. We show that in striatal fast-spiking interneurons, entrainment to gamma-frequency input depends on rhythmic spike runs and is explained by the phase resetting curve, whereas run initiation can be triggered by a broad range of input frequencies.

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.

Population dynamics and entrainment of basal ganglia pacemakers are shaped by their dendritic arbors

The theory of phase oscillators is an essential tool for understanding population dynamics of pacemaking neurons. GABAergic pacemakers in the substantia nigra pars reticulata (SNr), a main basal ganglia (BG) output nucleus, receive inputs from the direct and indirect pathways at distal and proximal regions of their dendritic arbors, respectively. We combine theory, optogenetic stimulation and electrophysiological experiments in acute brain slices to ask how dendritic properties impact the propensity of the various inputs, arriving at different locations along the dendrite, to recruit or entrain SNr pacemakers. By combining cable theory with sinusoidally-modulated optogenetic activation of either proximal somatodendritic regions or the entire somatodendritic arbor of SNr neurons, we construct an analytical model that accurately fits the empirically measured somatic current response to inputs arising from illuminating the soma and various portions of the dendritic field. We show that the extent of the dendritic tree that is illuminated generates measurable and systematic differences in the pacemaker’s phase response curve (PRC), causing a shift in its peak. Finally, we show that the divergent PRCs correctly predict differences in two major features of the collective dynamics of SNr neurons: the fidelity of population responses to sudden step-like changes in inputs; and the phase latency at which SNr neurons are entrained by rhythmic stimulation, which can occur in the BG under both physiological and pathophysiological conditions. Our novel method generates measurable and physiologically meaningful spatial effects, and provides the first empirical demonstration of how the collective responses of SNr pacemakers are determined by the transmission properties of their dendrites. SNr dendrites may serve to delay distal striatal inputs so that they impinge on the spike initiation zone simultaneously with pallidal and subthalamic inputs in order to guarantee a fair competition between the influence of the monosynaptic direct- and polysynaptic indirect pathways.

The substantia nigra pars reticulata (SNr) is a main output nucleus of the basal ganglia (BG), where inputs from the competing direct and indirect pathways converge onto the same neurons. Interestingly, these inputs are differentially distributed with direct and indirect pathway projections arriving at distal and proximal regions of the dendritic arbor, respectively. We employ a novel method combining theory with electrophysiological experiments and optogenetics to study the distinct effects of inputs arriving at different locations along the dendrite. Our approach represents a useful compromise between complexity and reduction in modelling. Our work addresses the question of high fidelity encoding of inputs by networks of neurons in the new context of pacemaking neurons, which are driven to fire by their intrinsic dynamics rather than by a network state. We provide the first empirical demonstration that dendritic delays can introduce latencies in the responses of a population of neurons that are commensurate with synaptic delays, suggesting a new role for SNr dendrites with implications for BG function.

Predicting responses to inhibitory synaptic input in substantia nigra pars reticulata neurons

The changes in firing probability produced by a synaptic input are usually visualized using the poststimulus time histogram (PSTH). It would be useful if postsynaptic firing patterns could be predicted from patterns of afferent synaptic activation, but attempts to predict the PSTH from synaptic potential waveforms using reasoning based on voltage trajectory and spike threshold have not been successful, especially for inhibitory inputs. We measured PSTHs for substantia nigra pars reticulata (SNr) neurons inhibited by optogenetic stimulation of striato-nigral inputs or by matching artificial inhibitory conductances applied by dynamic clamp. The PSTH was predicted by a model based on each SNr cell's phase-resetting curve (PRC). Optogenetic activation of striato-nigral input or artificial synaptic inhibition produced a PSTH consisting of an initial depression of firing followed by oscillatory increases and decreases repeating at the SNr cell's baseline firing rate. The phase resetting model produced PSTHs closely resembling the cell data, including the primary pause in firing and the oscillation. Key features of the PSTH, including the onset rate and duration of the initial inhibitory phase, and the subsequent increase in firing probability could be explained from the characteristic shape of the SNr cell's PRC. The rate of damping of the late oscillation was explained by the influence of asynchronous phase perturbations producing firing rate jitter and wander. Our results demonstrate the utility of phase-resetting models as a general method for predicting firing in spontaneously active neurons and their value in interpretation of the striato-nigral PSTH. NEW & NOTEWORTHY The coupling of patterned presynaptic input to sequences of postsynaptic firing is a Gordian knot, complicated by the multidimensionality of neuronal state and the diversity of potential initial states. Even so, it is fundamental for even the simplest understanding of network dynamics. We show that a simple phase-resetting model constructed from experimental measurements can explain and predict the sequence of spike rate changes following synaptic inhibition of an oscillating basal ganglia output neuron.

Predicting the response of striatal spiny neurons to sinusoidal input

Spike-timing effects of small-amplitude sinusoidal currents were measured in mouse striatal spiny neurons firing repetitively. Spike-timing reliability varied with the stimulus frequency. For frequencies near the cell's firing rate, the cells altered firing rate to match the stimulus and became phase locked to it. The stimulus phase of firing during lock depended on the stimulus frequency relative to the cell's unperturbed firing rate. Interspike intervals during sinusoidal stimulation were predicted using an iterative map constructed from the cells' phase-resetting curve. Variability of interspike intervals was reduced by stimulation at all frequencies higher than about half the cell's unperturbed rate, and interspike intervals were accurately predicted by the map. Long sequences of spike times were predicted by iterating on the map. The accuracy of that prediction varied with frequency. Spike time predictability was highest near and during phase lock. The map predicted the phase of firing on the input and its dependence on stimulus frequency. Prediction errors, when they occurred, were of two kinds: unpredicted variation in interspike interval from intrinsic cell noise and accumulation of prediction errors from previous interspike intervals. Each type of prediction error arose from a different mechanism, and their impact was also predicted from the phase model. When two oscillatory input currents were presented simultaneously, striatal neurons responded selectively to only one of them, the one closest in frequency to the cell's unperturbed firing rate. Their spike times encoded the frequency and phase of that single oscillatory input.NEW & NOTEWORTHY During repetitive firing, the timing of action potentials is determined by the interaction between the input and voltage-sensitive currents throughout the interspike interval. This interaction is encapsulated in the neuron's phase-resetting curve. The phase-resetting curve predicted spike timing to small sinusoidal currents over a wide range of stimulus frequencies. Firing patterns were most sensitive to oscillatory components near the cell's own firing rate, even in the presence of noise and other inputs.

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.

Dendritic and Axonal Propagation Delays Determine Emergent Structures of Neuronal Networks with Plastic Synapses

Spike-timing-dependent plasticity (STDP) modifies synaptic strengths based on the relative timing of pre- and postsynaptic spikes. The temporal order of spikes turned out to be crucial. We here take into account how propagation delays, composed of dendritic and axonal delay times, may affect the temporal order of spikes. In a minimal setting, characterized by neglecting dendritic and axonal propagation delays, STDP eliminates bidirectional connections between two coupled neurons and turns them into unidirectional connections. In this paper, however, we show that depending on the dendritic and axonal propagation delays, the temporal order of spikes at the synapses can be different from those in the cell bodies and, consequently, qualitatively different connectivity patterns emerge. In particular, we show that for a system of two coupled oscillatory neurons, bidirectional synapses can be preserved and potentiated. Intriguingly, this finding also translates to large networks of type-II phase oscillators and, hence, crucially impacts on the overall hierarchical connectivity patterns of oscillatory neuronal networks.

The Ca2+-activated chloride channel anoctamin-2 mediates spike-frequency adaptation and regulates sensory transmission in thalamocortical neurons

Neuronal firing patterns, which are crucial for determining the nature of encoded information, have been widely studied; however, the molecular identity and cellular mechanisms of spike-frequency adaptation are still not fully understood. Here we show that spike-frequency adaptation in thalamocortical (TC) neurons is mediated by the Ca²⁺-activated Cl⁻ channel (CACC) anoctamin-2 (ANO2). Knockdown of ANO2 in TC neurons results in significantly reduced spike-frequency adaptation along with increased tonic spiking. Moreover, thalamus-specific knockdown of ANO2 increases visceral pain responses. These results indicate that ANO2 contributes to reductions in spike generation in highly activated TC neurons and thereby restricts persistent information transmission.

Spike-frequency adaptation in thalamocortical (TC) neurons is important for sensory transmission though the underlying mechanisms are not fully understood. Here, the authors identify a role for the calcium-activated chloride channel, ANO2, in mediating TC spiking adaptations and visceral pain response.

Dopamine Neurons Change the Type of Excitability in Response to Stimuli

The dynamics of neuronal excitability determine the neuron's response to stimuli, its synchronization and resonance properties and, ultimately, the computations it performs in the brain. We investigated the dynamical mechanisms underlying the excitability type of dopamine (DA) neurons, using a conductance-based biophysical model, and its regulation by intrinsic and synaptic currents. Calibrating the model to reproduce low frequency tonic firing results in N-methyl-D-aspartate (NMDA) excitation balanced by γ-Aminobutyric acid (GABA)-mediated inhibition and leads to type I excitable behavior characterized by a continuous decrease in firing frequency in response to hyperpolarizing currents. Furthermore, we analyzed how excitability type of the DA neuron model is influenced by changes in the intrinsic current composition. A subthreshold sodium current is necessary for a continuous frequency decrease during application of a negative current, and the low-frequency "balanced" state during simultaneous activation of NMDA and GABA receptors. Blocking this current switches the neuron to type II characterized by the abrupt onset of repetitive firing. Enhancing the anomalous rectifier Ih current also switches the excitability to type II. Key characteristics of synaptic conductances that may be observed in vivo also change the type of excitability: a depolarized γ-Aminobutyric acid receptor (GABAR) reversal potential or co-activation of α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptors (AMPARs) leads to an abrupt frequency drop to zero, which is typical for type II excitability. Coactivation of N-methyl-D-aspartate receptors (NMDARs) together with AMPARs and GABARs shifts the type I/II boundary toward more hyperpolarized GABAR reversal potentials. To better understand how altering each of the aforementioned currents leads to changes in excitability profile of DA neuron, we provide a thorough dynamical analysis. Collectively, these results imply that type I excitability in dopamine neurons might be important for low firing rates and fine-tuning basal dopamine levels, while switching excitability to type II during NMDAR and AMPAR activation may facilitate a transient increase in dopamine concentration, as type II neurons are more amenable to synchronization by mutual excitation.

Mathematical Frameworks for Oscillatory Network Dynamics in Neuroscience

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

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.

Neurons as oscillators

Regularly spiking neurons can be described as oscillators. In this article we review some of the insights gained from this conceptualization and their relevance for systems neuroscience. First, we explain how a regularly spiking neuron can be viewed as an oscillator and how the phase-response curve (PRC) describes the response of the neuron's spike times to small perturbations. We then discuss the meaning of the PRC for a single neuron's spiking behavior and review the PRCs measured from a variety of neurons in a range of spiking regimes. Next, we show how the PRC can be related to a number of common measures used to quantify neuronal firing, such as the spike-triggered average and the peristimulus histogram. We further show that the response of a neuron to correlated inputs depends on the shape of the PRC. We then explain how the PRC of single neurons can be used to predict neural network behavior. Given the PRC, conduction delays, and the waveform and time course of the synaptic potentials, it is possible to predict neural population behavior such as synchronization. The PRC also allows us to quantify the robustness of the synchronization to heterogeneity and noise. We finally ask how to combine the measured PRCs and the predictions based on PRC to further the understanding of systems neuroscience. As an example, we discuss how the change of the PRC by the neuromodulator acetylcholine could lead to a destabilization of cortical network dynamics. Although all of these studies are grounded in mathematical abstractions that do not strictly hold in biology, they provide good estimates for the emergence of the brain's network activity from the properties of individual neurons. The study of neurons as oscillators can provide testable hypotheses and mechanistic explanations for systems neuroscience.

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.

Neural Cross-Frequency Coupling: Connecting Architectures, Mechanisms, and Functions

Neural oscillations are ubiquitously observed in the mammalian brain, but it has proven difficult to tie oscillatory patterns to specific cognitive operations. Notably, the coupling between neural oscillations at different timescales has recently received much attention, both from experimentalists and theoreticians. We review the mechanisms underlying various forms of this cross-frequency coupling. We show that different types of neural oscillators and cross-frequency interactions yield distinct signatures in neural dynamics. Finally, we associate these mechanisms with several putative functions of cross-frequency coupling, including neural representations of multiple environmental items, communication over distant areas, internal clocking of neural processes, and modulation of neural processing based on temporal predictions.

Estimation of the phase response curve from Parkinsonian tremor

Phase response curves (PRCs), characterizing the response of an oscillator to weak external perturbation, have been estimated from a broad range of biological oscillators, including single neurons in vivo. PRC estimates, in turn, provide an intuitive insight into how oscillatory systems become entrained and how they can be desynchronized. Here, we explore the application of PRC theory to the case of Parkinsonian tremor. Initial attempts to establish a causal effect of subthreshold transcranial magnetic stimulation applied to primary motor cortex on the filtered tremor phase were unsuccessful. We explored the possible explanations of this and demonstrate that assumptions made when estimating the PRC in a traditional setting, such as a single neuron, are not arbitrary when applied to the case of tremor PRC estimation. We go on to extract the PRC of Parkinsonian tremor using an iterative method that requires varying the definition of the tremor cycle and estimating the PRC at multiple peristimulus time samples. Justification for this method is supported by estimates of PRC from simulated single neuron data. We provide an approach to estimating confidence limits for tremor PRC and discuss the interpretational caveats introduced by tremor harmonics and the intrinsic variability of the tremor's period.

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.

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.

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.

Subthreshold membrane currents confer distinct tuning properties that enable neurons to encode the integral or derivative of their input

Neurons rely on action potentials, or spikes, to encode information. But spikes can encode different stimulus features in different neurons. We show here through simulations and experiments how neurons encode the integral or derivative of their input based on the distinct tuning properties conferred upon them by subthreshold currents. Slow-activating subthreshold inward (depolarizing) current mediates positive feedback control of subthreshold voltage, sustaining depolarization and allowing the neuron to spike on the basis of its integrated stimulus waveform. Slow-activating subthreshold outward (hyperpolarizing) current mediates negative feedback control of subthreshold voltage, truncating depolarization and forcing the neuron to spike on the basis of its differentiated stimulus waveform. Depending on its direction, slow-activating subthreshold current cooperates or competes with fast-activating inward current during spike initiation. This explanation predicts that sensitivity to the rate of change of stimulus intensity differs qualitatively between integrators and differentiators. This was confirmed experimentally in spinal sensory neurons that naturally behave as specialized integrators or differentiators. Predicted sensitivity to different stimulus features was confirmed by covariance analysis. Integration and differentiation, which are themselves inverse operations, are thus shown to be implemented by the slow feedback mediated by oppositely directed subthreshold currents expressed in different neurons.

Examination of rhythmicity of extracellularly recorded neurons in the entorhinal cortex

A number of studies have examined the theta-rhythmic modulation of neuronal firing in the hippocampal circuit. For extracellular recordings, this is often done by examining spectral properties of the spike-time autocorrelogram, most significantly, for validating the presence or absence of theta modulation across species. These techniques can show significant rhythmicity for high firing rate, highly rhythmic neurons; however, they are substantially biased by several factors including the peak firing rate of the neuron, the amount of time spent in the neuron's receptive field, and other temporal properties of the rhythmicity such as cycle-skipping. These limitations make it difficult to examine rhythmic modulation in neurons with low firing rates or when an animal has short dwell times within the firing field and difficult to compare rhythmicity under disparate experimental conditions when these factors frequently differ. Here, we describe in detail the challenges that researchers face when using these techniques and apply our findings to recent recordings from bat entorhinal grid cells, suggesting that they may have lacked enough data to examine theta rhythmicity robustly. We describe a more sensitive and statistically rigorous method using maximum likelihood estimation (MLE) of a parametric model of the lags within the autocorrelation window, which helps to alleviate some of the problems of traditional methods and was also unable to detect rhythmicity in bat grid cells. Using large batteries of simulated data, we explored the boundaries for which the MLE technique and the theta index can detect rhythmicity. The MLE technique is less sensitive to many features of the autocorrelogram and provides a framework for statistical testing to detect rhythmicity as well as changes in rhythmicity in individual sessions providing a substantial improvement over previous methods.

Modulation of a Single Neuron Has State-Dependent Actions on Circuit Dynamics,

We use computational models to examine the circumstances under which modulation of a single neuron can alter the dynamics of an entire circuit. We show that under some circumstances, the circuit's behavior is robust to the neuromodulation of a single hub neuron, while under other circumstances the same neuromodulatory action of a single neuron can produce a large variety of circuit outcomes.

Cover Figure

Adaptation and shunting inhibition leads to pyramidal/interneuron gamma with sparse firing of pyramidal cells

Gamma oscillations are a prominent phenomenon related to a number of brain functions. Data show that individual pyramidal neurons can fire at rate below gamma with the population showing clear gamma oscillations and synchrony. In one kind of idealized model of such weak gamma, pyramidal neurons fire in clusters. Here we provide a theory for clustered gamma PING rhythms with strong inhibition and weaker excitation. Our simulations of biophysical models show that the adaptation of pyramidal neurons coupled with their low firing rate leads to cluster formation. A partially analytic study of a canonical model shows that the phase response curves with a near zero flat region, caused by the presence of the slow adaptive current, are the key to the formation of clusters. Furthermore we examine shunting inhibition and show that clusters become robust and generic.

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.

Solving Winfree's puzzle: the isochrons in the FitzHugh-Nagumo model

We consider the FitzHugh-Nagumo model, an example of a system with two time scales for which Winfree was unable to determine the overall structure of the isochrons. An isochron is the set of all points in the basin of an attracting periodic orbit that converge to this periodic orbit with the same asymptotic phase. We compute the isochrons as one-dimensional parametrised curves with a method based on the continuation of suitable two-point boundary value problems. This allows us to present in detail the geometry of how the basin of attraction is foliated by isochrons. They exhibit extreme sensitivity and feature sharp turns, which is why Winfree had difficulties finding them. We observe that the sharp turns and sensitivity of the isochrons are associated with the slow-fast nature of the FitzHugh-Nagumo system; more specifically, it occurs near its repelling (unstable) slow manifold.

Dynamics of the exponential integrate-and-fire model with slow currents and adaptation

In order to properly capture spike-frequency adaptation with a simplified point-neuron model, we study approximations of Hodgkin-Huxley (HH) models including slow currents by exponential integrate-and-fire (EIF) models that incorporate the same types of currents. We optimize the parameters of the EIF models under the external drive consisting of AMPA-type conductance pulses using the current-voltage curves and the van Rossum metric to best capture the subthreshold membrane potential, firing rate, and jump size of the slow current at the neuron’s spike times. Our numerical simulations demonstrate that, in addition to these quantities, the approximate EIF-type models faithfully reproduce bifurcation properties of the HH neurons with slow currents, which include spike-frequency adaptation, phase-response curves, critical exponents at the transition between a finite and infinite number of spikes with increasing constant external drive, and bifurcation diagrams of interspike intervals in time-periodically forced models. Dynamics of networks of HH neurons with slow currents can also be approximated by corresponding EIF-type networks, with the approximation being at least statistically accurate over a broad range of Poisson rates of the external drive. For the form of external drive resembling realistic, AMPA-like synaptic conductance response to incoming action potentials, the EIF model affords great savings of computation time as compared with the corresponding HH-type model. Our work shows that the EIF model with additional slow currents is well suited for use in large-scale, point-neuron models in which spike-frequency adaptation is important.

Effect of phase response curve skew on synchronization with and without conduction delays

A central problem in cortical processing including sensory binding and attentional gating is how neurons can synchronize their responses with zero or near-zero time lag. For a spontaneously firing neuron, an input from another neuron can delay or advance the next spike by different amounts depending upon the timing of the input relative to the previous spike. This information constitutes the phase response curve (PRC). We present a simple graphical method for determining the effect of PRC shape on synchronization tendencies and illustrate it using type 1 PRCs, which consist entirely of advances (delays) in response to excitation (inhibition). We obtained the following generic solutions for type 1 PRCs, which include the pulse-coupled leaky integrate and fire model. For pairs with mutual excitation, exact synchrony can be stable for strong coupling because of the stabilizing effect of the causal limit region of the PRC in which an input triggers a spike immediately upon arrival. However, synchrony is unstable for short delays, because delayed inputs arrive during a refractory period and cannot trigger an immediate spike. Right skew destabilizes antiphase and enables modes with time lags that grow as the conduction delay is increased. Therefore, right skew favors near synchrony at short conduction delays and a gradual transition between synchrony and antiphase for pairs coupled by mutual excitation. For pairs with mutual inhibition, zero time lag synchrony is stable for conduction delays ranging from zero to a substantial fraction of the period for pairs. However, for right skew there is a preferred antiphase mode at short delays. In contrast to mutual excitation, left skew destabilizes antiphase for mutual inhibition so that synchrony dominates at short delays as well. These pairwise synchronization tendencies constrain the synchronization properties of neurons embedded in larger networks.

Status epilepticus as a system disturbance: is status epilepticus due to synchronization or desynchronization?

The traditional view of seizure activity is one in which there is extreme hypersynchronization. Although what is meant by hypersynchronization is rarely explicitly and fully defined, it can be understood to imply large numbers of neurons firing together essentially simultaneously. In this discussion we explore the possibility that seizures-both self-terminating and sustained in status-are not purely synchronous in time or in space. We investigate the alternative possibility that much seizure activity represents spatiotemporal desynchronization. Furthermore, we discuss the possibility that, in contrast to canonical views of epileptic activity, a high degree of synchronization is a prerequisite for termination of the seizure rather than a marker of early and ongoing seizure activity. These ideas will be discussed with reference to results from our collaborative group based on microelectrode recordings in patients with epilepsy as well as to the many studies done by others in both patients and animal models. Finally, we will explore implications for these hypotheses in the treatment of patients with epilepsy and in status epilepticus.

Local linear approximation of the Jacobian matrix better captures phase resetting of neural limit cycle oscillators

One effect of any external perturbations, such as presynaptic inputs, received by limit cycle oscillators when they are part of larger neural networks is a transient change in their firing rate, or phase resetting. A brief external perturbation moves the figurative point outside the limit cycle, a geometric perturbation that we mapped into a transient change in the firing rate, or a temporal phase resetting. In order to gain a better qualitative understanding of the link between the geometry of the limit cycle and the phase resetting curve (PRC), we used a moving reference frame with one axis tangent and the others normal to the limit cycle. We found that the stability coefficients associated with the unperturbed limit cycle provided good quantitative predictions of both the tangent and the normal geometric displacements induced by external perturbations. A geometric-to-temporal mapping allowed us to correctly predict the PRC while preserving the intuitive nature of this geometric approach.

Interaction function of oscillating coupled neurons

Large scale simulations of electrically coupled neuronal oscillators often employ the phase coupled oscillator paradigm to understand and predict network behavior. We study the nature of the interaction between such coupled oscillators using weakly coupled oscillator theory. By employing piecewise linear approximations for phase response curves and voltage time courses and parametrizing their shapes, we compute the interaction function for all such possible shapes and express it in terms of discrete Fourier modes. We find that reasonably good approximation is achieved with four Fourier modes that comprise of both sine and cosine terms.