Effect of sharp jumps at the edges of phase response curves on synchronization of electrically coupled neuronal oscillators

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.

Variety of alternative stable phase-locking in networks of electrically coupled relaxation oscillators

We studied the dynamics of a large-scale model network comprised of oscillating electrically coupled neurons. Cells are modeled as relaxation oscillators with short duty cycle, so they can be considered either as models of pacemaker cells, spiking cells with fast regenerative and slow recovery variables or firing rate models of excitatory cells with synaptic depression or cellular adaptation.

It was already shown that electrically coupled relaxation oscillators exhibit not only synchrony but also anti-phase behavior if electrical coupling is weak. We show that a much wider spectrum of spatiotemporal patterns of activity can emerge in a network of electrically coupled cells as a result of switching from synchrony, produced by short external signals of different spatial profiles. The variety of patterns increases with decreasing rate of neuronal firing (or duty cycle) and with decreasing strength of electrical coupling. We study also the effect of network topology - from all-to-all – to pure ring connectivity, where only the closest neighbors are coupled.

We show that the ring topology promotes anti-phase behavior as compared to all-to-all coupling. It also gives rise to a hierarchical organization of activity: during each of the main phases of a given pattern cells fire in a particular sequence determined by the local connectivity. We have analyzed the behavior of the network using geometric phase plane methods and we give heuristic explanations of our findings.

Our results show that complex spatiotemporal activity patterns can emerge due to the action of stochastic or sensory stimuli in neural networks without chemical synapses, where each cell is equally coupled to others via gap junctions. This suggests that in developing nervous systems where only electrical coupling is present such a mechanism can lead to the establishment of proto-networks generating premature multiphase oscillations whereas the subsequent emergence of chemical synapses would later stabilize generated patterns.

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.