Phase-resetting curves determine synchronization, phase locking, and clustering in networks of neural oscillators

Dynamic assemblies of parvalbumin interneurons in brain oscillations

Brain oscillations are crucial for perception, memory, and behavior. Parvalbumin-expressing (PV) interneurons are critical for these oscillations, but their population dynamics remain unclear. Using voltage imaging, we simultaneously recorded membrane potentials in up to 26 PV interneurons in vivo during hippocampal ripple oscillations in mice. We found that PV cells generate ripple-frequency rhythms by forming highly dynamic cell assemblies. These assemblies exhibit rapid and significant changes from cycle to cycle, varying greatly in both size and membership. Importantly, this variability is not just random spiking failures of individual neurons. Rather, the activities of other PV cells contain significant information about whether a PV cell spikes or not in a given cycle. This coordination persists without network oscillations, and it exists in subthreshold potentials even when the cells are not spiking. Dynamic assemblies of interneurons may provide a new mechanism to modulate postsynaptic dynamics and impact cognitive functions flexibly and rapidly.

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.

Phase Resetting Curves Determine Stability of Synchrony in One and Two Clusters of Pulse Coupled Oscillators with Delays

Previously, our group used phase response curves under a pulsatile coupling assumption to determine the stability of synchrony within a cluster of neural oscillators and between two clusters of oscillators. The interactions of the within and between cluster terms were considered, demonstrating how an alternating firing pattern between clusters could stabilize within cluster synchrony-even in clusters unable to synchronize themselves in isolation. In addition, criteria were derived for synchrony between two pulse coupled oscillators with synaptic delays. In this study, we update our previous work on one and two clusters of coupled oscillators to include delays and demonstrate the validity of the results using a map of the firing intervals based on the phase resetting curve. We use self-connected neurons to represent clusters and derive conditions under which an oscillator can phase-lock itself with a delayed input. Although this analysis only strictly applies to identical neurons receiving identical synapses from the same number of neurons, the principles are general and can be used to understand how to promote or impede synchrony in physiological networks of neurons. 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.

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.

Patterns of synchronization in 2D networks of inhibitory neurons

Neural firing in many inhibitory networks displays synchronous assembly or clustered firing, in which subsets of neurons fire synchronously, and these subsets may vary with different inputs to, or states of, the network. Most prior analytical and computational modeling of such networks has focused on 1D networks or 2D networks with symmetry (often circular symmetry). Here, we consider a 2D discrete network model on a general torus, where neurons are coupled to two or more nearest neighbors in three directions (horizontal, vertical, and diagonal), and allow different coupling strengths in all directions. Using phase model analysis, we establish conditions for the stability of different patterns of clustered firing behavior in the network. We then apply our results to study how variation of network connectivity and the presence of heterogeneous coupling strengths influence which patterns are stable. We confirm and supplement our results with numerical simulations of biophysical inhibitory neural network models. Our work shows that 2D networks may exhibit clustered firing behavior that cannot be predicted as a simple generalization of a 1D network, and that heterogeneity of coupling can be an important factor in determining which patterns are stable.

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.

Local inhibition in a model of the indirect pathway globus pallidus network slows and deregularizes background firing, but sharpens and synchronizes responses to striatal input

The external segment of globus pallidus (GPe) is a network of oscillatory neurons connected by inhibitory synapses. We studied the intrinsic dynamic and the response to a shared brief inhibitory stimulus in a model GPe network. Individual neurons were simulated using a phase resetting model based on measurements from mouse GPe neurons studied in slices. The neurons showed a broad heterogeneity in their firing rates and in the shapes and sizes of their phase resetting curves. Connectivity in the network was set to match experimental measurements. We generated statistically equivalent neuron heterogeneity in a small-world model, in which 99% of connections were made with near neighbors and 1% at random, and in a model with entirely random connectivity. In both networks, the resting activity was slowed and made more irregular by the local inhibition, but it did not show any periodic pattern. Cross-correlations among neuron pairs were limited to directly connected neurons. When stimulated by a shared inhibitory input, the individual neuron responses separated into two groups: one with a short and stereotyped period of inhibition followed by a transient increase in firing probability, and the other responding with a sustained inhibition. Despite differences in firing rate, the responses of the first group of neurons were of fixed duration and were synchronized across cells.

Universality of stable multi-cluster periodic solutions in a population model of the cell cycle with negative feedback

We study a population model where cells in one part of the cell cycle may affect the progress of cells in another part. If the influence, or feedback, from one part to another is negative, simulations of the model almost always result in multiple temporal clusters formed by groups of cells. We study regions in parameter space where periodic 'k-cyclic' solutions are stable. The regions of stability coincide with sub-triangles on which certain events occur in a fixed order. For boundary sub-triangles with order 'rs1', we prove that the k-cyclic periodic solution is asymptotically stable if the index of the sub-triangle is relatively prime with respect to the number of clusters k and neutrally stable otherwise. For negative linear feedback, we prove that the interior of the parameter set is covered by stable sub-triangles, i.e. a stable k-cyclic solution always exists for some k. We observe numerically that the result also holds for many forms of nonlinear feedback, but may break down in extreme cases.

Emergence of Neuronal Synchronisation in Coupled Areas

One of the most fundamental questions in the field of neuroscience is the emergence of synchronous behaviour in the brain, such as phase, anti-phase, and shift-phase synchronisation. In this work, we investigate how the connectivity between brain areas can influence the phase angle and the neuronal synchronisation. To do this, we consider brain areas connected by means of excitatory and inhibitory synapses, in which the neuron dynamics is given by the adaptive exponential integrate-and-fire model. Our simulations suggest that excitatory and inhibitory connections from one area to another play a crucial role in the emergence of these types of synchronisation. Thus, in the case of unidirectional interaction, we observe that the phase angles of the neurons in the receiver area depend on the excitatory and inhibitory synapses which arrive from the sender area. Moreover, when the neurons in the sender area are synchronised, the phase angle variability of the receiver area can be reduced for some conductance values between the areas. For bidirectional interactions, we find that phase and anti-phase synchronisation can emerge due to excitatory and inhibitory connections. We also verify, for a strong inhibitory-to-excitatory interaction, the existence of silent neuronal activities, namely a large number of excitatory neurons that remain in silence for a long time.

A Hypothesis for Theta Rhythm Frequency Control in CA1 Microcircuits

Computational models of neural circuits with varying levels of biophysical detail have been generated in pursuit of an underlying mechanism explaining the ubiquitous hippocampal theta rhythm. However, within the theta rhythm are at least two types with distinct frequencies associated with different behavioral states, an aspect that must be considered in pursuit of these mechanistic explanations. Here, using our previously developed excitatory-inhibitory network models that generate theta rhythms, we investigate the robustness of theta generation to intrinsic neuronal variability by building a database of heterogeneous excitatory cells and implementing them in our microcircuit model. We specifically investigate the impact of three key "building block" features of the excitatory cell model that underlie our model design: these cells' rheobase, their capacity for post-inhibitory rebound, and their spike-frequency adaptation. We show that theta rhythms at various frequencies can arise dependent upon the combination of these building block features, and we find that the speed of these oscillations are dependent upon the excitatory cells' response to inhibitory drive, as encapsulated by their phase response curves. Taken together, these findings support a hypothesis for theta frequency control that includes two aspects: (i) an internal mechanism that stems from the building block features of excitatory cell dynamics; (ii) an external mechanism that we describe as "inhibition-based tuning" of excitatory cell firing. We propose that these mechanisms control theta rhythm frequencies and underlie their robustness.

Differential contributions of synaptic and intrinsic inhibitory currents to speech segmentation via flexible phase-locking in neural oscillators

Current hypotheses suggest that speech segmentation—the initial division and grouping of the speech stream into candidate phrases, syllables, and phonemes for further linguistic processing—is executed by a hierarchy of oscillators in auditory cortex. Theta (∼3-12 Hz) rhythms play a key role by phase-locking to recurring acoustic features marking syllable boundaries. Reliable synchronization to quasi-rhythmic inputs, whose variable frequency can dip below cortical theta frequencies (down to ∼1 Hz), requires “flexible” theta oscillators whose underlying neuronal mechanisms remain unknown. Using biophysical computational models, we found that the flexibility of phase-locking in neural oscillators depended on the types of hyperpolarizing currents that paced them. Simulated cortical theta oscillators flexibly phase-locked to slow inputs when these inputs caused both (i) spiking and (ii) the subsequent buildup of outward current sufficient to delay further spiking until the next input. The greatest flexibility in phase-locking arose from a synergistic interaction between intrinsic currents that was not replicated by synaptic currents at similar timescales. Flexibility in phase-locking enabled improved entrainment to speech input, optimal at mid-vocalic channels, which in turn supported syllabic-timescale segmentation through identification of vocalic nuclei. Our results suggest that synaptic and intrinsic inhibition contribute to frequency-restricted and -flexible phase-locking in neural oscillators, respectively. Their differential deployment may enable neural oscillators to play diverse roles, from reliable internal clocking to adaptive segmentation of quasi-regular sensory inputs like speech.

Oscillatory activity in auditory cortex is believed to play an important role in auditory and speech processing. One suggested function of these rhythms is to divide the speech stream into candidate phonemes, syllables, words, and phrases, to be matched with learned linguistic templates. This requires brain rhythms to flexibly synchronize with regular acoustic features of the speech stream. How neuronal circuits implement this task remains unknown. In this study, we explored the contribution of inhibitory currents to flexible phase-locking in neuronal theta oscillators, believed to perform initial syllabic segmentation. We found that a combination of specific intrinsic inhibitory currents at multiple timescales, present in a large class of cortical neurons, enabled exceptionally flexible phase-locking, which could be used to precisely segment speech by identifying vowels at mid-syllable. This suggests that the cells exhibiting these currents are a key component in the brain’s auditory and speech processing architecture.

Spatially localized cluster solutions in inhibitory neural networks

Neurons in the inhibitory network of the striatum display cell assembly firing patterns which recent results suggest may consist of spatially compact neural clusters. Previous computational modeling of striatal neural networks has indicated that non-monotonic, distance-dependent coupling may promote spatially localized cluster firing. Here, we identify conditions for the existence and stability of cluster firing solutions in which clusters consist of spatially adjacent neurons in inhibitory neural networks. We consider simple non-monotonic, distance-dependent connectivity schemes in weakly coupled 1-D networks where cells make stronger connections with their kth nearest neighbors on each side and weaker connections with closer neighbors. Using the phase model reduction of the network system, we prove the existence of cluster solutions where neurons that are spatially close together are also synchronized in the same cluster, and find stability conditions for these solutions. Our analysis predicts the long-term behavior for networks of neurons, and we confirm our results by numerical simulations of biophysical neuron network models. Our results demonstrate that an inhibitory network with non-monotonic, distance-dependent connectivity can exhibit cluster solutions where adjacent cells fire together.

Order-indeterminant event-based maps for learning a beat

The process by which humans synchronize to a musical beat is believed to occur through error-correction where an individual's estimates of the period and phase of the beat time are iteratively adjusted to align with an external stimuli. Mathematically, error-correction can be described using a two-dimensional map where convergence to a fixed point corresponds to synchronizing to the beat. In this paper, we show how a neural system, called a beat generator, learns to adapt its oscillatory behavior through error-correction to synchronize to an external periodic signal. We construct a two-dimensional event-based map, which iteratively adjusts an internal parameter of the beat generator to speed up or slow down its oscillatory behavior to bring it into synchrony with the periodic stimulus. The map is novel in that the order of events defining the map are not a priori known. Instead, the type of error-correction adjustment made at each iterate of the map is determined by a sequence of expected events. The map possesses a rich repertoire of dynamics, including periodic solutions and chaotic orbits.

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.

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

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

Delta phase reset predicts conflict-related changes in P3 amplitude and behavior

When multiple competing responses are activated, we respond more slowly than if only one response is activated (response conflict). Conflict-induced slowing is reduced for consecutive high-conflict stimuli, an effect known as conflict adaptation. Verguts and Notebaert's (2009) adaptation by binding theory suggests this is due to Hebbian learning of cognitive control, potentiated by the response of the locus coeruleus norepinephrine (NE) system. Phasic activity of the NE system can potentially be measured non-invasively in humans by recording the P3 component of the event-related potential (ERP), and the P3 is sensitive to conflict adaptation. Bouret and Sara's (2005) network reset theory suggests that phasic NE might functionally reset ongoing large-scale network activity, generating synchronous neural population activity like the P3. To examine the possibility that network reset contributes to conflict effects in the P3, we recorded high-density EEG data while subjects performed a flanker task. As expected, conflict and conflict adaptation modulated P3 amplitudes. Brain-behavior correlation analyses indicated that activity during the rise of the P3 was related to RT and predicted RT differences due to conflict. More importantly, phase of delta oscillations not only predicted reaction time differences between low-conflict and high-conflict conditions, but delta phase reset also predicted the amplitude of the P3. Delta oscillations exhibited dominant peaks both pre and post-stimulus, and delta at stimulus onset predicted the post-stimulus ERP, in particular the N2 and P3. This result bridges human EEG with basic mechanisms suggested by computational neural models and invasive patient recordings, namely that salient cognitive events might reset ongoing oscillations leading to the generation of the phase-locked evoked potential. We conclude that partial phase reset is a cortical mechanism involved in monitoring the environment for unexpected events, and this response contributes to conflict effects in the ERP. These results are in line with theories that phasic NE release might reset ongoing cortical activity, leading to the generation of ERP components like the P3.

Dynamics and sources of response variability and its coordination in visual cortex

The trial-to-trial response variability in sensory cortices and the extent to which this variability can be coordinated among cortical units have strong implications for cortical signal processing. Yet, little is known about the relative contributions and dynamics of defined sources to the cortical response variability and their correlations across cortical units. To fill this knowledge gap, here we obtained and analyzed multisite local field potential (LFP) recordings from visual cortex of turtles in response to repeated naturalistic movie clips and decomposed cortical across-trial LFP response variability into three defined sources, namely, input, network, and local fluctuations. We found that input fluctuations dominate cortical response variability immediately following stimulus onset, whereas network fluctuations dominate the response variability in the steady state during continued visual stimulation. Concurrently, we found that the network fluctuations dominate the correlations of the variability during the ongoing and steady-state epochs, but not immediately following stimulus onset. Furthermore, simulations of various model networks indicated that (i) synaptic time constants, leading to oscillatory activity, and (ii) synaptic clustering and synaptic depression, leading to spatially constrained pockets of coherent activity, are both essential features of cortical circuits to mediate the observed relative contributions and dynamics of input, network, and local fluctuations to the cortical LFP response variability and their correlations across recording sites. In conclusion, these results show how a mélange of multiscale thalamocortical circuit features mediate a complex stimulus-modulated cortical activity that, when naively related to the visual stimulus alone, appears disguised as high and coordinated across-trial response variability.

Itinerant chimeras in an adaptive network of pulse-coupled oscillators

In a network of pulse-coupled oscillators with adaptive coupling, we discover a dynamical regime which we call an "itinerant chimera." Similarly as in classical chimera states, the network splits into two domains, the coherent and the incoherent. The drastic difference is that the composition of the domains is volatile, i.e., the oscillators demonstrate spontaneous switching between the domains. This process can be seen as traveling of the oscillators from one domain to another or as traveling of the chimera core across the network. We explore the basic features of the itinerant chimeras, such as the mean and the variance of the core size, and the oscillators lifetime within the core. We also study the scaling behavior of the system and show that the observed regime is not a finite-size effect but a key feature of the collective dynamics which persists even in large networks.

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.

Memory-induced nonlinear dynamics of excitation in cardiac diseases

Excitable cells, such as cardiac myocytes, exhibit short-term memory, i.e., the state of the cell depends on its history of excitation. Memory can originate from slow recovery of membrane ion channels or from accumulation of intracellular ion concentrations, such as calcium ion or sodium ion concentration accumulation. Here we examine the effects of memory on excitation dynamics in cardiac myocytes under two diseased conditions, early repolarization and reduced repolarization reserve, each with memory from two different sources: slow recovery of a potassium ion channel and slow accumulation of the intracellular calcium ion concentration. We first carry out computer simulations of action potential models described by differential equations to demonstrate complex excitation dynamics, such as chaos. We then develop iterated map models that incorporate memory, which accurately capture the complex excitation dynamics and bifurcations of the action potential models. Finally, we carry out theoretical analyses of the iterated map models to reveal the underlying mechanisms of memory-induced nonlinear dynamics. Our study demonstrates that the memory effect can be unmasked or greatly exacerbated under certain diseased conditions, which promotes complex excitation dynamics, such as chaos. The iterated map models reveal that memory converts a monotonic iterated map function into a nonmonotonic one to promote the bifurcations leading to high periodicity and chaos.

Propagation delays determine neuronal activity and synaptic connectivity patterns emerging in plastic neuronal networks

In plastic neuronal networks, the synaptic strengths are adapted to the neuronal activity. Specifically, spike-timing-dependent plasticity (STDP) is a fundamental mechanism that modifies the synaptic strengths based on the relative timing of pre- and postsynaptic spikes, taking into account the spikes' temporal order. In many studies, propagation delays were neglected to avoid additional dynamic complexity or computational costs. So far, networks equipped with a classic STDP rule typically rule out bidirectional couplings (i.e., either loops or uncoupled states) and are, hence, not able to reproduce fundamental experimental findings. In this review paper, we consider additional features, e.g., extensions of the classic STDP rule or additional aspects like noise, in order to overcome the contradictions between theory and experiment. In addition, we review in detail recent studies showing that a classic STDP rule combined with realistic propagation patterns is able to capture relevant experimental findings. In two coupled oscillatory neurons with propagation delays, bidirectional synapses can be preserved and potentiated. This result also holds for large networks of type-II phase oscillators. In addition, not only the mean of the initial distribution of synaptic weights, but also its standard deviation crucially determines the emergent structural connectivity, i.e., the mean final synaptic weight, the number of two-neuron loops, and the symmetry of the final connectivity pattern. The latter is affected by the firing rates, where more symmetric synaptic configurations emerge at higher firing rates. Finally, we discuss these findings in the context of the computational neuroscience-based development of desynchronizing brain stimulation techniques.

Delay-Induced Multistability and Loop Formation in Neuronal Networks with Spike-Timing-Dependent Plasticity

Spike-timing-dependent plasticity (STDP) adjusts synaptic strengths according to the precise timing of pre- and postsynaptic spike pairs. Theoretical and computational studies have revealed that STDP may contribute to the emergence of a variety of structural and dynamical states in plastic neuronal populations. In this manuscript, we show that by incorporating dendritic and axonal propagation delays in recurrent networks of oscillatory neurons, the asymptotic connectivity displays multistability, where different structures emerge depending on the initial distribution of the synaptic strengths. In particular, we show that the standard deviation of the initial distribution of synaptic weights, besides its mean, determines the main properties of the emergent structural connectivity such as the mean final synaptic weight, the number of two-neuron loops and the symmetry of the final structure. We also show that the firing rates of the neurons affect the evolution of the network, and a more symmetric configuration of the synapses emerges at higher firing rates. We justify the network results based on a two-neuron framework and show how the results translate to large recurrent networks.

Olfaction Warps Visual Time Perception

Our perception of the world builds upon dynamic inputs from multiple senses with different temporal resolutions, and is threaded with the passing of subjective time. How time is extracted from multisensory inputs is scantly known. Utilizing psychophysical testing and electroencephalography, we show in healthy human adults that odors modulate object visibility around critical flicker-fusion frequency (CFF)-the limit at which chromatic flickers become perceived as a stable color-and effectively alter CFF in a congruency-based manner, despite that they afford no clear environmental temporal information. The behavioral gain produced by a congruent relative to an incongruent odor is accompanied by elevated neural oscillatory power around the object's flicker frequency in the right temporal region ~150-300 ms after object onset, and is not mediated by visual awareness. In parallel, odors bias the subjective duration of visual objects without affecting one's temporal sensitivity. These findings point to a neuronal network in the right temporal cortex that executes flexible temporal filtering of upstream visual inputs based on olfactory information. Moreover, they collectively indicate that the very process of sensory integration at the stage of object processing twists time perception, hence casting new insights into the neural timing of multisensory events.

Prediction of memory formation based on absolute electroencephalographic phases in rhinal cortex and hippocampus outperforms prediction based on stimulus-related phase shifts

Absolute (i.e. measured) rhinal and hippocampal phase values are predictive for memory formation. It has been an open question, whether the capability of mediotemporal structures to react to stimulus presentation with phase shifts may be similarly indicative of successful memory formation. We analysed data from 27 epilepsy patients implanted with depth electrodes in the hippocampus and entorhinal cortex, who performed a continuous word recognition task. Electroencephalographic phase information related to the first presentation of repeatedly presented words was used for prediction of subsequent remembering vs. forgetting applying a support vector machine. The capability to predict successful memory formation based on stimulus-related phase shifts was compared to that based on absolute phase values. Average hippocampal phase shifts were larger and rhinal phase shifts were more accumulated for later remembered compared to forgotten trials. Nevertheless, prediction based on absolute phase values clearly outperformed phase shifts and there was no significant increase in prediction accuracies when combining both measures. Our findings indicate that absolute rhinal and hippocampal phases and not stimulus-related phase shifts are most relevant for successful memory formation. Absolute phases possibly affect memory formation via influencing neural membrane potentials and thereby controlling the timing of neural firing.

Phase-locking and bistability in neuronal networks with synaptic depression

We consider a recurrent network of two oscillatory neurons that are coupled with inhibitory synapses. We use the phase response curves of the neurons and the properties of short-term synaptic depression to define Poincaré maps for the activity of the network. The fixed points of these maps correspond to phase-locked modes of the network. Using these maps, we analyze the conditions that allow short-term synaptic depression to lead to the existence of bistable phase-locked, periodic solutions. We show that bistability arises when either the phase response curve of the neuron or the short-term depression profile changes steeply enough. The results apply to any Type I oscillator and we illustrate our findings using the Quadratic Integrate-and-Fire and Morris-Lecar neuron models.

Dichotomous Dynamics in E-I Networks with Strongly and Weakly Intra-connected Inhibitory Neurons

The interconnectivity between excitatory and inhibitory neural networks informs mechanisms by which rhythmic bursts of excitatory activity can be produced in the brain. One such mechanism, Pyramidal Interneuron Network Gamma (PING), relies primarily upon reciprocal connectivity between the excitatory and inhibitory networks, while also including intra-connectivity of inhibitory cells. The causal relationship between excitatory activity and the subsequent burst of inhibitory activity is of paramount importance to the mechanism and has been well studied. However, the role of the intra-connectivity of the inhibitory network, while important for PING, has not been studied in detail, as most analyses of PING simply assume that inhibitory intra-connectivity is strong enough to suppress subsequent firing following the initial inhibitory burst. In this paper we investigate the role that the strength of inhibitory intra-connectivity plays in determining the dynamics of PING-style networks. We show that networks with weak inhibitory intra-connectivity exhibit variations in burst dynamics of both the excitatory and inhibitory cells that are not obtained with strong inhibitory intra-connectivity. Networks with weak inhibitory intra-connectivity exhibit excitatory rhythmic bursts with weak excitatory-to-inhibitory synapses for which classical PING networks would show no rhythmic activity. Additionally, variations in dynamics of these networks as the excitatory-to-inhibitory synaptic weight increases illustrates the important role that consistent pattern formation in the inhibitory cells serves in maintaining organized and periodic excitatory bursts. Finally, motivated by these results and the known diversity of interneurons, we show that a PING-style network with two inhibitory subnetworks, one strongly intra-connected and one weakly intra-connected, exhibits organized and periodic excitatory activity over a larger parameter regime than networks with a homogeneous inhibitory population. Taken together, these results serve to better articulate the role of inhibitory intra-connectivity in generating PING-like rhythms, while also revealing how heterogeneity amongst inhibitory synapses might make such rhythms more robust to a variety of network parameters.

Self-organized emergence of multilayer structure and chimera states in dynamical networks with adaptive couplings

We report the phenomenon of self-organized emergence of hierarchical multilayered structures and chimera states in dynamical networks with adaptive couplings. This process is characterized by a sequential formation of subnetworks (layers) of densely coupled elements, the size of which is ordered in a hierarchical way, and which are weakly coupled between each other. We show that the hierarchical structure causes the decoupling of the subnetworks. Each layer can exhibit either a two-cluster state, a periodic traveling wave, or an incoherent state, and these states can coexist on different scales of subnetwork sizes.

Measurement of phase resetting curves using optogenetic barrage stimuli

BACKGROUND

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

NEW METHOD

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

RESULTS

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

COMPARISON WITH EXISTING METHODS

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

CONCLUSIONS

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

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

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

Memory-Induced Chaos in Cardiac Excitation

Excitable systems display memory, but how memory affects the excitation dynamics of such systems remains to be elucidated. Here we use computer simulation of cardiac action potential models to demonstrate that memory can cause dynamical instabilities that result in complex excitation dynamics and chaos. We develop an iterated map model that correctly describes these dynamics and show that memory converts a monotonic first return map of action potential duration into a nonmonotonic one, resulting in a period-doubling bifurcation route to chaos.

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.

Travelling Wave Pulse Coupled Oscillator (TWPCO) Using a Self-Organizing Scheme for Energy-Efficient Wireless Sensor Networks

Recently, Pulse Coupled Oscillator (PCO)-based travelling waves have attracted substantial attention by researchers in wireless sensor network (WSN) synchronization. Because WSNs are generally artificial occurrences that mimic natural phenomena, the PCO utilizes firefly synchronization of attracting mating partners for modelling the WSN. However, given that sensor nodes are unable to receive messages while transmitting data packets (due to deafness), the PCO model may not be efficient for sensor network modelling. To overcome this limitation, this paper proposed a new scheme called the Travelling Wave Pulse Coupled Oscillator (TWPCO). For this, the study used a self-organizing scheme for energy-efficient WSNs that adopted travelling wave biologically inspired network systems based on phase locking of the PCO model to counteract deafness. From the simulation, it was found that the proposed TWPCO scheme attained a steady state after a number of cycles. It also showed superior performance compared to other mechanisms, with a reduction in the total energy consumption of 25%. The results showed that the performance improved by 13% in terms of data gathering. Based on the results, the proposed scheme avoids the deafness that occurs in the transmit state in WSNs and increases the data collection throughout the transmission states in WSNs.

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.

Intrinsic Cellular Properties and Connectivity Density Determine Variable Clustering Patterns in Randomly Connected Inhibitory Neural Networks

The plethora of inhibitory interneurons in the hippocampus and cortex play a pivotal role in generating rhythmic activity by clustering and synchronizing cell firing. Results of our simulations demonstrate that both the intrinsic cellular properties of neurons and the degree of network connectivity affect the characteristics of clustered dynamics exhibited in randomly connected, heterogeneous inhibitory networks. We quantify intrinsic cellular properties by the neuron's current-frequency relation (IF curve) and Phase Response Curve (PRC), a measure of how perturbations given at various phases of a neurons firing cycle affect subsequent spike timing. We analyze network bursting properties of networks of neurons with Type I or Type II properties in both excitability and PRC profile; Type I PRCs strictly show phase advances and IF curves that exhibit frequencies arbitrarily close to zero at firing threshold while Type II PRCs display both phase advances and delays and IF curves that have a non-zero frequency at threshold. Type II neurons whose properties arise with or without an M-type adaptation current are considered. We analyze network dynamics under different levels of cellular heterogeneity and as intrinsic cellular firing frequency and the time scale of decay of synaptic inhibition are varied. Many of the dynamics exhibited by these networks diverge from the predictions of the interneuron network gamma (ING) mechanism, as well as from results in all-to-all connected networks. Our results show that randomly connected networks of Type I neurons synchronize into a single cluster of active neurons while networks of Type II neurons organize into two mutually exclusive clusters segregated by the cells' intrinsic firing frequencies. Networks of Type II neurons containing the adaptation current behave similarly to networks of either Type I or Type II neurons depending on network parameters; however, the adaptation current creates differences in the cluster dynamics compared to those in networks of Type I or Type II neurons. To understand these results, we compute neuronal PRCs calculated with a perturbation matching the profile of the synaptic current in our networks. Differences in profiles of these PRCs across the different neuron types reveal mechanisms underlying the divergent network dynamics.

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.

Jittering waves in rings of pulse oscillators

Rings of oscillators with delayed pulse coupling are studied analytically, numerically, and experimentally. The basic regimes observed in such rings are rotating waves with constant interspike intervals and phase lags between the neighbors. We show that these rotating waves may destabilize leading to the so-called jittering waves. For these regimes, the interspike intervals are no more equal but form a periodic sequence in time. Analytic criterion for the emergence of jittering waves is derived and confirmed by the numerical and experimental data. The obtained results contribute to the hypothesis that the multijitter instability is universal in systems with pulse coupling.

Phase response curves for models of earthquake fault dynamics

We systematically study effects of external perturbations on models describing earthquake fault dynamics. The latter are based on the framework of the Burridge-Knopoff spring-block system, including the cases of a simple mono-block fault, as well as the paradigmatic complex faults made up of two identical or distinct blocks. The blocks exhibit relaxation oscillations, which are representative for the stick-slip behavior typical for earthquake dynamics. Our analysis is carried out by determining the phase response curves of first and second order. For a mono-block fault, we consider the impact of a single and two successive pulse perturbations, further demonstrating how the profile of phase response curves depends on the fault parameters. For a homogeneous two-block fault, our focus is on the scenario where each of the blocks is influenced by a single pulse, whereas for heterogeneous faults, we analyze how the response of the system depends on whether the stimulus is applied to the block having a shorter or a longer oscillation period.

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.

Resonant Interneurons Can Increase Robustness of Gamma Oscillations

Gamma oscillations are believed to play a critical role in in information processing, encoding, and retrieval. Inhibitory interneuronal network gamma (ING) oscillations may arise from a coupled oscillator mechanism in which individual neurons oscillate or from a population oscillator in which individual neurons fire sparsely and stochastically. All ING mechanisms, including the one proposed herein, rely on alternating waves of inhibition and windows of opportunity for spiking. The coupled oscillator model implemented with Wang-Buzsáki model neurons is not sufficiently robust to heterogeneity in excitatory drive, and therefore intrinsic frequency, to account for in vitro models of ING. Similarly, in a tightly synchronized regime, the stochastic population oscillator model is often characterized by sparse firing, whereas interneurons both in vivo and in vitro do not fire sparsely during gamma, but rather on average every other cycle. We substituted so-called resonator neural models, which exhibit class 2 excitability and postinhibitory rebound (PIR), for the integrators that are typically used. This results in much greater robustness to heterogeneity that actually increases as the average participation in spikes per cycle approximates physiological levels. Moreover, dynamic clamp experiments that show autapse-induced firing in entorhinal cortical interneurons support the idea that PIR can serve as a network gamma mechanism. Furthermore, parvalbumin-positive (PV(+)) cells were much more likely to display both PIR and autapse-induced firing than GAD2(+) cells, supporting the view that PV(+) fast-firing basket cells are more likely to exhibit class 2 excitability than other types of inhibitory interneurons.

SIGNIFICANCE STATEMENT

Gamma oscillations are believed to play a critical role in information processing, encoding, and retrieval. Networks of inhibitory interneurons are thought to be essential for these oscillations. We show that one class of interneurons with an abrupt onset of firing at a threshold frequency may allow more robust synchronization in the presence of noise and heterogeneity. The mechanism for this robustness depends on the intrinsic resonance at this threshold frequency. Moreover, we show experimentally the feasibility of the proposed mechanism and suggest a way to distinguish between this mechanism and another proposed mechanism: that of a stochastic population oscillator independent of the dynamics of individual neurons.

Cooperation and competition of gamma oscillation mechanisms

Oscillations of neuronal activity in different frequency ranges are thought to reflect important aspects of cortical network dynamics. Here we investigate how various mechanisms that contribute to oscillations in neuronal networks may interact. We focus on networks with inhibitory, excitatory, and electrical synapses, where the subnetwork of inhibitory interneurons alone can generate interneuron gamma (ING) oscillations and the interactions between interneurons and pyramidal cells allow for pyramidal-interneuron gamma (PING) oscillations. What type of oscillation will such a network generate? We find that ING and PING oscillations compete: The mechanism generating the higher oscillation frequency "wins"; it determines the frequency of the network oscillation and suppresses the other mechanism. For type I interneurons, the network oscillation frequency is equal to or slightly above the higher of the ING and PING frequencies in corresponding reduced networks that can generate only either of them; if the interneurons belong to the type II class, it is in between. In contrast to ING and PING, oscillations mediated by gap junctions and oscillations mediated by inhibitory synapses may cooperate or compete, depending on the type (I or II) of interneurons and the strengths of the electrical and chemical synapses. We support our computer simulations by a theoretical model that allows a full theoretical analysis of the main results. Our study suggests experimental approaches to deciding to what extent oscillatory activity in networks of interacting excitatory and inhibitory neurons is dominated by ING or PING oscillations and of which class the participating interneurons are.

Clustered Desynchronization from High-Frequency Deep Brain Stimulation

While high-frequency deep brain stimulation is a well established treatment for Parkinson’s disease, its underlying mechanisms remain elusive. Here, we show that two competing hypotheses, desynchronization and entrainment in a population of model neurons, may not be mutually exclusive. We find that in a noisy group of phase oscillators, high frequency perturbations can separate the population into multiple clusters, each with a nearly identical proportion of the overall population. This phenomenon can be understood by studying maps of the underlying deterministic system and is guaranteed to be observed for small noise strengths. When we apply this framework to populations of Type I and Type II neurons, we observe clustered desynchronization at many pulsing frequencies.

While high-frequency deep brain stimulation (DBS) is a decades old treatment for alleviating the motor symptoms Parkinsons disease, the way in which it alleviates these symptoms is not well understood. Making matters more complicated, some experimental results suggest that DBS works by making neurons fire more regularly, while other seemingly contradictory results suggest that DBS works by making neural firing patterns less synchronized. Here we present theoretical and numerical results with the potential to merge these competing hypotheses. For predictable DBS pulsing rates, in the presence of a small amount of noise, a population of neurons will split into distinct clusters, each containing a nearly identical proportion of the overall population. When we observe this clustering phenomenon, on a short time scale, neurons are entrained to high-frequency DBS pulsing, but on a long time scale, they desynchronize predictably.

Emergence and combinatorial accumulation of jittering regimes in spiking oscillators with delayed feedback

Interaction via pulses is common in many natural systems, especially neuronal. In this article we study one of the simplest possible systems with pulse interaction: a phase oscillator with delayed pulsatile feedback. When the oscillator reaches a specific state, it emits a pulse, which returns after propagating through a delay line. The impact of an incoming pulse is described by the oscillator's phase reset curve (PRC). In such a system we discover an unexpected phenomenon: for a sufficiently steep slope of the PRC, a periodic regular spiking solution bifurcates with several multipliers crossing the unit circle at the same parameter value. The number of such critical multipliers increases linearly with the delay and thus may be arbitrary large. This bifurcation is accompanied by the emergence of numerous "jittering" regimes with nonequal interspike intervals (ISIs). Each of these regimes corresponds to a periodic solution of the system with a period roughly proportional to the delay. The number of different "jittering" solutions emerging at the bifurcation point increases exponentially with the delay. We describe the combinatorial mechanism that underlies the emergence of such a variety of solutions. In particular, we show how a periodic solution exhibiting several distinct ISIs can imply the existence of multiple other solutions obtained by rearranging of these ISIs. We show that the theoretical results for phase oscillators accurately predict the behavior of an experimentally implemented electronic oscillator with pulsatile feedback.

Synchrony in stochastically driven neuronal networks with complex topologies

We study the synchronization of a stochastically driven, current-based, integrate-and-fire neuronal model on a preferential-attachment network with scale-free characteristics and high clustering. The synchrony is induced by cascading total firing events where every neuron in the network fires at the same instant of time. We show that in the regime where the system remains in this highly synchronous state, the firing rate of the network is completely independent of the synaptic coupling, and depends solely on the external drive. On the other hand, the ability for the network to maintain synchrony depends on a balance between the fluctuations of the external input and the synaptic coupling strength. In order to accurately predict the probability of repeated cascading total firing events, we go beyond mean-field and treelike approximations and conduct a detailed second-order calculation taking into account local clustering. Our explicit analytical results are shown to give excellent agreement with direct numerical simulations for the particular preferential-attachment network model investigated.

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.

Phase-resetting as a tool of information transmission

Models of information transmission in the brain largely rely on firing rate codes. The abundance of oscillatory activity in the brain suggests that information may be also encoded using the phases of ongoing oscillations. Sensory perception, working memory and spatial navigation have been hypothesized to use phase codes, and cross-frequency coordination and phase synchronization between brain areas have been proposed to gate the flow of information. Phase codes generally require the phase of the oscillations to be reset at specific reference points for consistent coding, and coordination between oscillators requires favorable phase resetting characteristics. Recent evidence supports a role for neural oscillations in providing temporal reference windows that allow for correct parsing of phase-coded information.