Effects of conduction delays on the existence and stability of one to one phase locking between two pulse-coupled oscillators

A mean field theory for pulse-coupled neural oscillators based on the spike time response curve

A mean field method for pulse-coupled oscillators with delays used a self-connected oscillator to represent a synchronous cluster of N - 1 oscillators and a single oscillator assumed to be perturbed from the cluster. A periodic train of biexponential conductance input was divided into a tonic and a phasic component representing the mean field input. A single cycle of the phasic conductance from the cluster was applied to the single oscillator embedded in the tonic component at different phases to measure the change in the cycle length in which the perturbation was initiated, that is, the first-order phase response curve (PRC), and the second-order PRC in the following cycle. A homogeneous network of 100 biophysically calibrated inhibitory interneurons with either shunting or hyperpolarizing inhibition tested the predictive power of the method. A self-consistency criterion predicted the oscillation frequency of the network from the PRCs as a function of the synaptic delay. The major determinant of the stability of synchrony was the sign of the slope of the first-order PRC of the single oscillator in response to an input from the self-connected cluster at a phase corresponding to the delay value. For most short delays, first-order PRCs correctly predicted the frequency and stability of simulated network activity. However, considering the second-order PRC improved the frequency prediction and resolved an incorrect prediction of stability of global synchrony at delays close to the free running period of single neurons in which a discontinuity in the PRC precluded existence of 1:1 self-locking.NEW & NOTEWORTHY A mean field theory for synchrony in neural networks in which neurons are generally above threshold in the mean-driven regime is developed to extend and complement mean field theory previously developed by others for neurons that are generally below threshold in the fluctuation-driven regime. This work extends phase response curve theory as applied to high-frequency oscillations in networks with synaptic inputs that are not short with respect to the network period.

Existence and stability criteria for global synchrony and for synchrony in two alternating clusters of pulse-coupled oscillators updated to include conduction delays

Phase Response Curves (PRCs) have been useful in determining and analyzing various phase-locking modes in networks of oscillators under pulse-coupling assumptions, as reviewed in Mathematical Biosciences, 226:77-96, 2010. Here, we update that review to include progress since 2010 on pulse coupled oscillators with conduction delays. We then present original results that extend the derivation of the criteria for stability of global synchrony in networks of pulse-coupled oscillators to include conduction delays. We also incorporate conduction delays to extend previous studies that showed how an alternating firing pattern between two synchronized clusters could enforce within-cluster synchrony, even for clusters unable to synchronize themselves in isolation. To obtain these results, we used self-connected neurons to represent clusters. These results greatly extend the applicability of the stability analyses to networks of pulse-coupled oscillators since conduction delays are ubiquitous and strongly impact the stability of synchrony. Although these analyses only strictly apply to identical oscillators with identical connections to other oscillators, the principles are general and suggest how to promote or impede synchrony in physiological networks of neurons, for example. Heterogeneity can be interpreted as a form of frozen noise, and approximate synchrony can be sustained despite heterogeneity. The pulse-coupled oscillator model can not only be used to describe biological neuronal networks but also cardiac pacemakers, lasers, fireflies, artificial neural networks, social self-organization, and wireless sensor networks.

Existence and Stability Criteria for Global Synchrony and for Synchrony in two Alternating Clusters of Pulse-Coupled Oscillators Updated to Include Conduction Delays

Phase Response Curves (PRCs) have been useful in determining and analyzing various phase-locking modes in networks of oscillators under pulse-coupling assumptions, as reviewed in Mathematical Biosciences, 226:77-96, 2010. Here, we update that review to include progress since 2010 on pulse coupled oscillators with conduction delays. We then present original results that extend the derivation of the criteria for stability of global synchrony in networks of pulse-coupled oscillators to include conduction delays. We also incorporate conduction delays to extend previous studies that showed how an alternating firing pattern between two synchronized clusters could enforce within cluster synchrony, even for clusters unable to synchronize themselves in isolation. To obtain these results, we used self-connected neurons to represent clusters. These results greatly extend the applicability of the stability analyses to networks of pulse-coupled oscillators since conduction delays are ubiquitous and strongly impact the stability of synchrony. Although these analyses only strictly apply to identical oscillators with identical connections to other oscillators, the principles are general and suggest how to promote or impede synchrony in physiological networks of neurons, for example. Heterogeneity can be interpreted as a form of frozen noise, and approximate synchrony can be sustained despite heterogeneity. The pulse-coupled oscillator model can not only be used to describe biological neuronal networks but also cardiac pacemakers, lasers, fireflies, artificial neural networks, social self-organization, and wireless sensor networks.

AMS Subject Classification

37N25, 39A06, 39A30, 92B25, 92C20.

Synchronization of delayed coupled neurons with multiple synaptic connections

Synchronization is a key feature of the brain dynamics and is necessary for information transmission across brain regions and in higher brain functions like cognition, learning and memory. Experimental findings demonstrated that in cortical microcircuits there are multiple synapses between pairs of connected neurons. Synchronization of neurons in the presence of multiple synaptic connections may be relevant for optimal learning and memory, however, its effect on the dynamics of the neurons is not adequately studied. Here, we address the question that how changes in the strength of the synaptic connections and transmission delays between neurons impact synchronization in a two-neuron system with multiple synapses. To this end, we analytically and computationally investigated synchronization dynamics by considering both phase oscillator model and conductance-based Hodgkin-Huxley (HH) model. Our results show that symmetry/asymmetry of feedforward and feedback connections crucially determines stability of the phase locking of the system based on the strength of connections and delays. In both models, the two-neuron system with multiple synapses achieves in-phase synchrony in the presence of small and large delays, whereas an anti-phase synchronization state is favored for median delays. Our findings can expand the understanding of the functional role of multisynaptic contacts in neuronal synchronization and may shed light on the dynamical consequences of pathological multisynaptic connectivity in a number of brain disorders.

Synaptic dependence of dynamic regimes when coupling neural populations

In this article we focus on the study of the collective dynamics of neural networks. The analysis of two recent models of coupled "next-generation" neural mass models allows us to observe different global mean dynamics of large neural populations. These models describe the mean dynamics of all-to-all coupled networks of quadratic integrate-and-fire spiking neurons. In addition, one of these models considers the influence of the synaptic adaptation mechanism on the macroscopic dynamics. We show how both models are related through a parameter and we study the evolution of the dynamics when switching from one model to the other by varying that parameter. Interestingly, we have detected three main dynamical regimes in the coupled models: Rössler-type (funnel type), bursting-type, and spiking-like (oscillator-type) dynamics. This result opens the question of which regime is the most suitable for realistic simulations of large neural networks and shows the possibility of the emergence of chaotic collective dynamics when synaptic adaptation is very weak.

Delay-dependent transitions of phase synchronization and coupling symmetry between neurons shaped by spike-timing-dependent plasticity

Synchronization plays a key role in learning and memory by facilitating the communication between neurons promoted by synaptic plasticity. Spike-timing-dependent plasticity (STDP) is a form of synaptic plasticity that modifies the strength of synaptic connections between neurons based on the coincidence of pre- and postsynaptic spikes. In this way, STDP simultaneously shapes the neuronal activity and synaptic connectivity in a feedback loop. However, transmission delays due to the physical distance between neurons affect neuronal synchronization and the symmetry of synaptic coupling. To address the question that how transmission delays and STDP can jointly determine the emergent pairwise activity-connectivity patterns, we studied phase synchronization properties and coupling symmetry between two bidirectionally coupled neurons using both phase oscillator and conductance-based neuron models. We show that depending on the range of transmission delays, the activity of the two-neuron motif can achieve an in-phase/anti-phase synchronized state and its connectivity can attain a symmetric/asymmetric coupling regime. The coevolutionary dynamics of the neuronal system and the synaptic weights due to STDP stabilizes the motif in either one of these states by transitions between in-phase/anti-phase synchronization states and symmetric/asymmetric coupling regimes at particular transmission delays. These transitions crucially depend on the phase response curve (PRC) of the neurons, but they are relatively robust to the heterogeneity of transmission delays and potentiation-depression imbalance of the STDP profile.

Interneuronal network model of theta-nested fast oscillations predicts differential effects of heterogeneity, gap junctions and short term depression for hyperpolarizing versus shunting inhibition

Theta and gamma oscillations in the hippocampus have been hypothesized to play a role in the encoding and retrieval of memories. Recently, it was shown that an intrinsic fast gamma mechanism in medial entorhinal cortex can be recruited by optogenetic stimulation at theta frequencies, which can persist with fast excitatory synaptic transmission blocked, suggesting a contribution of interneuronal network gamma (ING). We calibrated the passive and active properties of a 100-neuron model network to capture the range of passive properties and frequency/current relationships of experimentally recorded PV+ neurons in the medial entorhinal cortex (mEC). The strength and probabilities of chemical and electrical synapses were also calibrated using paired recordings, as were the kinetics and short-term depression (STD) of the chemical synapses. Gap junctions that contribute a noticeable fraction of the input resistance were required for synchrony with hyperpolarizing inhibition; these networks exhibited theta-nested high frequency oscillations similar to the putative ING observed experimentally in the optogenetically-driven PV-ChR2 mice. With STD included in the model, the network desynchronized at frequencies above ~200 Hz, so for sufficiently strong drive, fast oscillations were only observed before the peak of the theta. Because hyperpolarizing synapses provide a synchronizing drive that contributes to robustness in the presence of heterogeneity, synchronization decreases as the hyperpolarizing inhibition becomes weaker. In contrast, networks with shunting inhibition required non-physiological levels of gap junctions to synchronize using conduction delays within the measured range.

Noise-induced switching in an oscillator with pulse delayed feedback: A discrete stochastic modeling approach

We study the dynamics of an oscillatory system with pulse delayed feedback and noise of two types: (i) phase noise acting on the oscillator and (ii) stochastic fluctuations of the feedback delay. Using an event-based approach, we reduce the system dynamics to a stochastic discrete map. For weak noise, we find that the oscillator fluctuates around a deterministic state, and we derive an autoregressive model describing the system dynamics. For stronger noise, the oscillator demonstrates noise-induced switching between various deterministic states; our theory provides a good estimate of the switching statistics in the linear limit. We show that the robustness of the system toward this switching is strikingly different depending on the type of noise. We compare the analytical results for linear coupling to numerical simulations of nonlinear coupling and find that the linear model also provides a qualitative explanation for the differences in robustness to both types of noise. Moreover, phase noise drives the system toward higher frequencies, while stochastic delays do not, and we relate this effect to our theoretical results.

Temporal patterns of synchrony in a pyramidal-interneuron gamma (PING) network

Synchronization in neural systems plays an important role in many brain functions. Synchronization in the gamma frequency band (30-100 Hz) is involved in a variety of cognitive phenomena; abnormalities of the gamma synchronization are found in schizophrenia and autism spectrum disorder. Frequently, the strength of synchronization is not high, and synchronization is intermittent even on short time scales (few cycles of oscillations). That is, the network exhibits intervals of synchronization followed by intervals of desynchronization. Neural circuit dynamics may show different distributions of desynchronization durations even if the synchronization strength is fixed. We use a conductance-based neural network exhibiting pyramidal-interneuron gamma rhythm to study the temporal patterning of synchronized neural oscillations. We found that changes in the synaptic strength (as well as changes in the membrane kinetics) can alter the temporal patterning of synchrony. Moreover, we found that the changes in the temporal pattern of synchrony may be independent of the changes in the average synchrony strength. Even though the temporal patterning may vary, there is a tendency for dynamics with short (although potentially numerous) desynchronizations, similar to what was observed in experimental studies of neural synchronization in the brain. Recent studies suggested that the short desynchronizations dynamics may facilitate the formation and the breakup of transient neural assemblies. Thus, the results of this study suggest that changes of synaptic strength may alter the temporal patterning of the gamma synchronization as to make the neural networks more efficient in the formation of neural assemblies and the facilitation of cognitive phenomena.

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.

Phase response theory explains cluster formation in sparsely but strongly connected inhibitory neural networks and effects of jitter due to sparse connectivity

We show how to predict whether a neural network will exhibit global synchrony (a one-cluster state) or a two-cluster state based on the assumption of pulsatile coupling and critically dependent upon the phase response curve (PRC) generated by the appropriate perturbation from a partner cluster. Our results hold for a monotonically increasing (meaning longer delays as the phase increases) PRC, which likely characterizes inhibitory fast-spiking basket and cortical low-threshold-spiking interneurons in response to strong inhibition. Conduction delays stabilize synchrony for this PRC shape, whereas they destroy two-cluster states, the former by avoiding a destabilizing discontinuity and the latter by approaching it. With conduction delays, stronger coupling strength can promote a one-cluster state, so the weak coupling limit is not applicable here. We show how jitter can destabilize global synchrony but not a two-cluster state. Local stability of global synchrony in an all-to-all network does not guarantee that global synchrony can be observed in an appropriately scaled sparsely connected network; the basin of attraction can be inferred from the PRC and must be sufficiently large. Two-cluster synchrony is not obviously different from one-cluster synchrony in the presence of noise and may be the actual substrate for oscillations observed in the local field potential (LFP) and the electroencephalogram (EEG) in situations where global synchrony is not possible. Transitions between cluster states may change the frequency of the rhythms observed in the LFP or EEG. Transitions between cluster states within an inhibitory subnetwork may allow more effective recruitment of pyramidal neurons into the network rhythm. NEW & NOTEWORTHY We show that jitter induced by sparse connectivity can destabilize global synchrony but not a two-cluster state with two smaller clusters firing alternately. On the other hand, conduction delays stabilize synchrony and destroy two-cluster states. These results hold if each cluster exhibits a phase response curve similar to one that characterizes fast-spiking basket and cortical low-threshold-spiking cells for strong inhibition. Either a two-cluster or a one-cluster state might provide the oscillatory substrate for neural computations.

Globally attracting synchrony in a network of oscillators with all-to-all inhibitory pulse coupling

The synchronization tendencies of networks of oscillators have been studied intensely. We assume a network of all-to-all pulse-coupled oscillators in which the effect of a pulse is independent of the number of oscillators that simultaneously emit a pulse and the normalized delay (the phase resetting) is a monotonically increasing function of oscillator phase with the slope everywhere less than 1 and a value greater than 2φ-1, where φ is the normalized phase. Order switching cannot occur; the only possible solutions are globally attracting synchrony and cluster solutions with a fixed firing order. For small conduction delays, we prove the former stable and all other possible attractors nonexistent due to the destabilizing discontinuity of the phase resetting at a phase of 0.

Stimulus-dependent synchronization in delayed-coupled neuronal networks

Time delay is a general feature of all interactions. Although the effects of delayed interaction are often neglected when the intrinsic dynamics is much slower than the coupling delay, they can be crucial otherwise. We show that delayed coupled neuronal networks support transitions between synchronous and asynchronous states when the level of input to the network changes. The level of input determines the oscillation period of neurons and hence whether time-delayed connections are synchronizing or desynchronizing. We find that synchronizing connections lead to synchronous dynamics, whereas desynchronizing connections lead to out-of-phase oscillations in network motifs and to frustrated states with asynchronous dynamics in large networks. Since the impact of a neuronal network to downstream neurons increases when spikes are synchronous, networks with delayed connections can serve as gatekeeper layers mediating the firing transfer to other regions. This mechanism can regulate the opening and closing of communicating channels between cortical layers on demand.

Zero-lag synchronization despite inhomogeneities in a relay system

A novel proposal for the zero-lag synchronization of the delayed coupled neurons, is to connect them indirectly via a third relay neuron. In this study, we develop a Poincaré map to investigate the robustness of the synchrony in such a relay system against inhomogeneity in the neurons and synaptic parameters. We show that when the inhomogeneity does not violate the symmetry of the system, synchrony is maintained and in some cases inhomogeneity enhances synchrony. On the other hand if the inhomogeneity breaks the symmetry of the system, zero lag synchrony can not be preserved. In this case we give analytical results for the phase lag of the spiking of the neurons in the stable state.

Cross-frequency synchronization of oscillators with time-delayed coupling

We carry out theoretical and experimental studies of cross-frequency synchronization of two pulse oscillators with time-delayed coupling. In the theoretical part of the paper we utilize the concept of phase resetting curves and analyze the system dynamics in the case of weak coupling. We construct a Poincaré map and obtain the synchronization zones in the parameter space for m:n synchronization. To challenge the theoretical results we designed an electronic circuit implementing the coupled oscillators and studied its dynamics experimentally. We show that the developed theory predicts dynamical properties of the realistic system, including location of the synchronization zones and bifurcations inside them.

Effects of synaptic plasticity on phase and period locking in a network of two oscillatory neurons

We study the effects of synaptic plasticity on the determination of firing period and relative phases in a network of two oscillatory neurons coupled with reciprocal inhibition. We combine the phase response curves of the neurons with the short-term synaptic plasticity properties of the synapses to define Poincaré maps for the activity of an oscillatory network. Fixed points of these maps correspond to the phase-locked modes of the network. These maps allow us to analyze the dependence of the resulting network activity on the properties of network components. Using a combination of analysis and simulations, we show how various parameters of the model affect the existence and stability of phase-locked solutions. We find conditions on the synaptic plasticity profiles and the phase response curves of the neurons for the network to be able to maintain a constant firing period, while varying the phase of locking between the neurons or vice versa. A generalization to cobwebbing for two-dimensional maps is also discussed.

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.

Direct connections assist neurons to detect correlation in small amplitude noises

We address a question on the effect of common stochastic inputs on the correlation of the spike trains of two neurons when they are coupled through direct connections. We show that the change in the correlation of small amplitude stochastic inputs can be better detected when the neurons are connected by direct excitatory couplings. Depending on whether intrinsic firing rate of the neurons is identical or slightly different, symmetric or asymmetric connections can increase the sensitivity of the system to the input correlation by changing the mean slope of the correlation transfer function over a given range of input correlation. In either case, there is also an optimum value for synaptic strength which maximizes the sensitivity of the system to the changes in input correlation.

Synchronization of delayed coupled neurons in presence of inhomogeneity

In principle, two directly coupled limit cycle oscillators can overcome mismatch in intrinsic rates and match their frequencies, but zero phase lag synchronization is just achievable in the limit of zero mismatch, i.e., with identical oscillators. Delay in communication, on the other hand, can exert phase shift in the activity of the coupled oscillators. In this study, we address the question of how phase locked, and in particular zero phase lag synchronization, can be achieved for a heterogeneous system of two delayed coupled neurons. We have analytically studied the possibility of inphase synchronization and near inphase synchronization when the neurons are not identical or the connections are not exactly symmetric. We have shown that while any single source of inhomogeneity can violate isochronous synchrony, multiple sources of inhomogeneity can compensate for each other and maintain synchrony. Numeric studies on biologically plausible models also support the analytic results.

Controlling the oscillation phase through precisely timed closed-loop optogenetic stimulation: a computational study

Dynamic oscillatory coherence is believed to play a central role in flexible communication between brain circuits. To test this communication-through-coherence hypothesis, experimental protocols that allow a reliable control of phase-relations between neuronal populations are needed. In this modeling study, we explore the potential of closed-loop optogenetic stimulation for the control of functional interactions mediated by oscillatory coherence. The theory of non-linear oscillators predicts that the efficacy of local stimulation will depend not only on the stimulation intensity but also on its timing relative to the ongoing oscillation in the target area. Induced phase-shifts are expected to be stronger when the stimulation is applied within specific narrow phase intervals. Conversely, stimulations with the same or even stronger intensity are less effective when timed randomly. Stimulation should thus be properly phased with respect to ongoing oscillations (in order to optimally perturb them) and the timing of the stimulation onset must be determined by a real-time phase analysis of simultaneously recorded local field potentials (LFPs). Here, we introduce an electrophysiologically calibrated model of Channelrhodopsin 2 (ChR2)-induced photocurrents, based on fits holding over two decades of light intensity. Through simulations of a neural population which undergoes coherent gamma oscillations—either spontaneously or as an effect of continuous optogenetic driving—we show that precisely-timed photostimulation pulses can be used to shift the phase of oscillation, even at transduction rates smaller than 25%. We consider then a canonic circuit with two inter-connected neural populations oscillating with gamma frequency in a phase-locked manner. We demonstrate that photostimulation pulses applied locally to a single population can induce, if precisely phased, a lasting reorganization of the phase-locking pattern and hence modify functional interactions between the two populations.

Synchronization by elastic neuronal latencies

Psychological and physiological considerations entail that formation and functionality of neuronal cell assemblies depend upon synchronized repeated activation such as zero-lag synchronization. Several mechanisms for the emergence of this phenomenon have been suggested, including the global network quantity, the greatest common divisor of neuronal circuit delay loops. However, they require strict biological prerequisites such as precisely matched delays and connectivity, and synchronization is represented as a stationary mode of activity instead of a transient phenomenon. Here we show that the unavoidable increase in neuronal response latency to ongoing stimulation serves as a nonuniform gradual stretching of neuronal circuit delay loops. This apparent nuisance is revealed to be an essential mechanism in various types of neuronal time controllers, where synchronization emerges as a transient phenomenon and without predefined precisely matched synaptic delays. These findings are described in an experimental procedure where conditioned stimulations were enforced on a circuit of neurons embedded within a large-scale network of cortical cells in vitro, and are corroborated and extended by simulations of circuits composed of Hodgkin-Huxley neurons with time-dependent latencies. These findings announce a cortical time scale for time controllers based on tens of microseconds stretching of neuronal circuit delay loops per spike. They call for a reexamination of the role of the temporal periodic mode in brain functionality using advanced in vitro and in vivo experiments.

When Long-Range Zero-Lag Synchronization is Feasible in Cortical Networks

Many studies have reported long-range synchronization of neuronal activity between brain areas, in particular in the beta and gamma bands with frequencies in the range of 14–30 and 40–80 Hz, respectively. Several studies have reported synchrony with zero phase lag, which is remarkable considering the synaptic and conduction delays inherent in the connections between distant brain areas. This result has led to many speculations about the possible functional role of zero-lag synchrony, such as for neuronal communication, attention, memory, and feature binding. However, recent studies using recordings of single-unit activity and local field potentials report that neuronal synchronization may occur with non-zero phase lags. This raises the questions whether zero-lag synchrony can occur in the brain and, if so, under which conditions. We used analytical methods and computer simulations to investigate which connectivity between neuronal populations allows or prohibits zero-lag synchrony. We did so for a model where two oscillators interact via a relay oscillator. Analytical results and computer simulations were obtained for both type I Mirollo–Strogatz neurons and type II Hodgkin–Huxley neurons. We have investigated the dynamics of the model for various types of synaptic coupling and importantly considered the potential impact of Spike-Timing Dependent Plasticity (STDP) and its learning window. We confirm previous results that zero-lag synchrony can be achieved in this configuration. This is much easier to achieve with Hodgkin–Huxley neurons, which have a biphasic phase response curve, than for type I neurons. STDP facilitates zero-lag synchrony as it adjusts the synaptic strengths such that zero-lag synchrony is feasible for a much larger range of parameters than without STDP.

Impact of adaptation currents on synchronization of coupled exponential integrate-and-fire neurons

The ability of spiking neurons to synchronize their activity in a network depends on the response behavior of these neurons as quantified by the phase response curve (PRC) and on coupling properties. The PRC characterizes the effects of transient inputs on spike timing and can be measured experimentally. Here we use the adaptive exponential integrate-and-fire (aEIF) neuron model to determine how subthreshold and spike-triggered slow adaptation currents shape the PRC. Based on that, we predict how synchrony and phase locked states of coupled neurons change in presence of synaptic delays and unequal coupling strengths. We find that increased subthreshold adaptation currents cause a transition of the PRC from only phase advances to phase advances and delays in response to excitatory perturbations. Increased spike-triggered adaptation currents on the other hand predominantly skew the PRC to the right. Both adaptation induced changes of the PRC are modulated by spike frequency, being more prominent at lower frequencies. Applying phase reduction theory, we show that subthreshold adaptation stabilizes synchrony for pairs of coupled excitatory neurons, while spike-triggered adaptation causes locking with a small phase difference, as long as synaptic heterogeneities are negligible. For inhibitory pairs synchrony is stable and robust against conduction delays, and adaptation can mediate bistability of in-phase and anti-phase locking. We further demonstrate that stable synchrony and bistable in/anti-phase locking of pairs carry over to synchronization and clustering of larger networks. The effects of adaptation in aEIF neurons on PRCs and network dynamics qualitatively reflect those of biophysical adaptation currents in detailed Hodgkin-Huxley-based neurons, which underscores the utility of the aEIF model for investigating the dynamical behavior of networks. Our results suggest neuronal spike frequency adaptation as a mechanism synchronizing low frequency oscillations in local excitatory networks, but indicate that inhibition rather than excitation generates coherent rhythms at higher frequencies.

Synchronization of neuronal spiking in the brain is related to cognitive functions, such as perception, attention, and memory. It is therefore important to determine which properties of neurons influence their collective behavior in a network and to understand how. A prominent feature of many cortical neurons is spike frequency adaptation, which is caused by slow transmembrane currents. We investigated how these adaptation currents affect the synchronization tendency of coupled model neurons. Using the efficient adaptive exponential integrate-and-fire (aEIF) model and a biophysically detailed neuron model for validation, we found that increased adaptation currents promote synchronization of coupled excitatory neurons at lower spike frequencies, as long as the conduction delays between the neurons are negligible. Inhibitory neurons on the other hand synchronize in presence of conduction delays, with or without adaptation currents. Our results emphasize the utility of the aEIF model for computational studies of neuronal network dynamics. We conclude that adaptation currents provide a mechanism to generate low frequency oscillations in local populations of excitatory neurons, while faster rhythms seem to be caused by inhibition rather than excitation.

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.