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

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.

Physiological features of parvalbumin-expressing GABAergic interneurons contributing to high-frequency oscillations in the cerebral cortex

Parvalbumin-expressing (PV+) inhibitory interneurons drive gamma oscillations (30-80 Hz), which underlie higher cognitive functions. In this review, we discuss two groups/aspects of fundamental properties of PV+ interneurons. In the first group (dubbed Before Axon), we list properties representing optimal synaptic integration in PV+ interneurons designed to support fast oscillations. For example: [i] Information can neither enter nor leave the neocortex without the engagement of fast PV+ -mediated inhibition; [ii] Voltage responses in PV+ interneuron dendrites integrate linearly to reduce impact of the fluctuations in the afferent drive; and [iii] Reversed somatodendritic Rm gradient accelerates the time courses of synaptic potentials arriving at the soma. In the second group (dubbed After Axon), we list morphological and biophysical properties responsible for (a) short synaptic delays, and (b) efficient postsynaptic outcomes. For example: [i] Fast-spiking ability that allows PV+ interneurons to outpace other cortical neurons (pyramidal neurons). [ii] Myelinated axon (which is only found in the PV+ subclass of interneurons) to secure fast-spiking at the initial axon segment; and [iii] Inhibitory autapses - autoinhibition, which assures brief biphasic voltage transients and supports postinhibitory rebounds. Recent advent of scientific tools, such as viral strategies to target PV cells and the ability to monitor PV cells via in vivo imaging during behavior, will aid in defining the role of PV cells in the CNS. Given the link between PV+ interneurons and cognition, in the future, it would be useful to carry out physiological recordings in the PV+ cell type selectively and characterize if and how psychiatric and neurological diseases affect initiation and propagation of electrical signals in this cortical sub-circuit. Voltage imaging may allow fast recordings of electrical signals from many PV+ interneurons simultaneously.

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.

Automatic Adaptation of Model Neurons and Connections to Build Hybrid Circuits with Living Networks

Hybrid circuits built by creating mono- or bi-directional interactions among living cells and model neurons and synapses are an effective way to study neuron, synaptic and neural network dynamics. However, hybrid circuit technology has been largely underused in the context of neuroscience studies mainly because of the inherent difficulty in implementing and tuning this type of interactions. In this paper, we present a set of algorithms for the automatic adaptation of model neurons and connections in the creation of hybrid circuits with living neural networks. The algorithms perform model time and amplitude scaling, real-time drift adaptation, goal-driven synaptic and model tuning/calibration and also automatic parameter mapping. These algorithms have been implemented in RTHybrid, an open-source library that works with hard real-time constraints. We provide validation examples by building hybrid circuits in a central pattern generator. The results of the validation experiments show that the proposed dynamic adaptation facilitates closed-loop communication among living and artificial model neurons and connections, and contributes to characterize system dynamics, achieve control, automate experimental protocols and extend the lifespan of the preparations.

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.

RTHybrid: A Standardized and Open-Source Real-Time Software Model Library for Experimental Neuroscience

Closed-loop technologies provide novel ways of online observation, control and bidirectional interaction with the nervous system, which help to study complex non-linear and partially observable neural dynamics. These protocols are often difficult to implement due to the temporal precision required when interacting with biological components, which in many cases can only be achieved using real-time technology. In this paper we introduce RTHybrid (www.github.com/GNB-UAM/RTHybrid), a free and open-source software that includes a neuron and synapse model library to build hybrid circuits with living neurons in a wide variety of experimental contexts. In an effort to encourage the standardization of real-time software technology in neuroscience research, we compared different open-source real-time operating system patches, RTAI, Xenomai 3 and Preempt-RT, according to their performance and usability. RTHybrid has been developed to run over Linux operating systems supporting both Xenomai 3 and Preempt-RT real-time patches, and thus allowing an easy implementation in any laboratory. We report a set of validation tests and latency benchmarks for the construction of hybrid circuits using this library. With this work we want to promote the dissemination of standardized, user-friendly and open-source software tools developed for open- and closed-loop experimental neuroscience.

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.

Anticipation via canards in excitable systems

Neurons can anticipate incoming signals by exploiting a physiological mechanism that is not well understood. This article offers a novel explanation on how a receiver neuron can predict the sender's dynamics in a unidirectionally-coupled configuration, in which both sender and receiver follow the evolution of a multi-scale excitable system. We present a novel theoretical viewpoint based on a mathematical object, called canard, to explain anticipation in excitable systems. We provide a numerical approach, which allows to determine the transient effects of canards. To demonstrate the general validity of canard-mediated anticipation in the context of excitable systems, we illustrate our framework in two examples, a multi-scale radio-wave circuit (the van der Pol model) that inspired a caricature neuronal model (the FitzHugh-Nagumo model) and a biophysical neuronal model (a 2-dimensional reduction of the Hodgkin-Huxley model), where canards act as messengers to the senders' prediction. We also propose an experimental paradigm that would enable experimental neuroscientists to validate our predictions. We conclude with an outlook to possible fascinating research avenues to further unfold the mechanisms underpinning anticipation. We envisage that our approach can be employed by a wider class of excitable systems with appropriate theoretical extensions.

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.

Modeling of synchronization behavior of bursting neurons at nonlinearly coupled dynamical networks

Synchronization behaviors of bursting neurons coupled through electrical and dynamic chemical synapses are investigated. The Izhikevich model is used with random and small world network of bursting neurons. Various currents which consist of diffusive electrical and time-delayed dynamic chemical synapses are used in the simulations to investigate the influences of synaptic currents and couplings on synchronization behavior of bursting neurons. The effects of parameters, such as time delay, inhibitory synaptic strengths, and decay time on synchronization behavior are investigated. It is observed that in random networks with no delay, bursting synchrony is established with the electrical synapse alone, single spiking synchrony is observed with hybrid coupling. In small world network with no delay, periodic bursting behavior with multiple spikes is observed when only chemical and only electrical synapse exist. Single-spike and multiple-spike bursting are established with hybrid couplings. A decrease in the synchronization measure is observed with zero time delay, as the decay time is increased in random network. For synaptic delays which are above active phase period, synchronization measure increases with an increase in synaptic strength and time delay in small world network. However, in random network, it increases with only an increase in synaptic strength.

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.

Unitary synaptic connections among substantia nigra pars reticulata neurons

Neurons in substantia nigra pars reticulata (SNr) are synaptically coupled by local axon collaterals, providing a potential mechanism for local signal processing. Because SNr neurons fire spontaneously, these synapses are constantly active. To investigate their properties, we recorded spontaneous inhibitory postsynaptic currents (sIPSCs) from SNr neurons in brain slices, in which afferents from upstream nuclei are severed, and the cells fire rhythmically. The sIPSC trains contained a mixture of periodic and aperiodic events. Autocorrelation analysis of sIPSC trains showed that a majority of cells had one to four active unitary inputs. The properties of the unitary IPSCs (uIPSCs) were analyzed for cells with one unitary input, using a model of periodic presynaptic firing and stochastic synaptic transmission. The inferred presynaptic firing rates and coefficient of variation of interspike intervals (ISIs) corresponded well with direct measurements of spiking in SNr neurons. Methods were developed to estimate the success probability, amplitude distributions, and kinetics of the uIPSCs, while removing the contribution from aperiodic sIPSCs. The sIPSC amplitudes were not increased upon release from halorhodopsin silencing, suggesting that most synapses were not depressed at the spontaneous firing rate. Gramicidin perforated-patch recordings indicated that the average reversal potential of spontaneous inhibitory postsynaptic potentials was -64 mV. Because of the change in driving force across the ISI, the unitary inputs are predicted to have a larger postsynaptic impact when they arrive late in the ISI. Simulations of network activity suggest that this very sparse inhibitory coupling may act to desynchronize the activity of SNr neurons while having only a small effect on firing rate.

Slow noise in the period of a biological oscillator underlies gradual trends and abrupt transitions in phasic relationships in hybrid neural networks

In order to study the ability of coupled neural oscillators to synchronize in the presence of intrinsic as opposed to synaptic noise, we constructed hybrid circuits consisting of one biological and one computational model neuron with reciprocal synaptic inhibition using the dynamic clamp. Uncoupled, both neurons fired periodic trains of action potentials. Most coupled circuits exhibited qualitative changes between one-to-one phase-locking with fairly constant phasic relationships and phase slipping with a constant progression in the phasic relationships across cycles. The phase resetting curve (PRC) and intrinsic periods were measured for both neurons, and used to construct a map of the firing intervals for both the coupled and externally forced (PRC measurement) conditions. For the coupled network, a stable fixed point of the map predicted phase locking, and its absence produced phase slipping. Repetitive application of the map was used to calibrate different noise models to simultaneously fit the noise level in the measurement of the PRC and the dynamics of the hybrid circuit experiments. Only a noise model that added history-dependent variability to the intrinsic period could fit both data sets with the same parameter values, as well as capture bifurcations in the fixed points of the map that cause switching between slipping and locking. We conclude that the biological neurons in our study have slowly-fluctuating stochastic dynamics that confer history dependence on the period. Theoretical results to date on the behavior of ensembles of noisy biological oscillators may require re-evaluation to account for transitions induced by slow noise dynamics.

The past, present, and future of real-time control in cellular electrophysiology

For over 60 years, real-time control has been an important technique in the study of excitable cells. Two such control-based technologies are reviewed here. First, voltage-clamp methods revolutionized the study of excitable cells. In this family of techniques, membrane potential is controlled, allowing one to parameterize a powerful class of models that describe the voltage-current relationship of cell membranes simply, flexibly, and accurately. Second, dynamic-clamp methods allow the addition of new, "virtual" membrane mechanisms to living cells. Dynamic clamp allows researchers unprecedented ways of testing computationally based hypotheses in biological preparations. The review ends with predictions of how control-based technologies will be improved and adapted for new uses in the near future.

Intercellular delay regulates the collective period of repressively coupled gene regulatory oscillator networks

Most biological rhythms are generated by a population of cellular oscillators coupled through intercellular signaling. Recent experimental evidence shows that the collective period may differ significantly from the autonomous period in the presence of intercellular delays. The phenomenon has been investigated using delay-coupled phase oscillators, but the proposed phase model contains no direct biological mechanism, which may weaken the model's reliability in unraveling biophysical principles. Based on a published gene regulatory oscillator model, we analyze the collective period of delay-coupled biological oscillators using the multivariable harmonic balance technique. We prove that, in contradiction to the common intuition that the collective period increases linearly with the coupling delay, the collective period turns out to be a periodic function of the intercellular delay. More surprisingly, the collective period may even decrease with the intercellular delay when the delay resides in certain regions. The collective period is given in a closed-form in terms of biochemical reaction constants and thus provides biological insights as well as guidance in synthetic-biological-oscillator design. Simulation results are given based on a segmentation clock model to confirm the theoretical predictions.

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.

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

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

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.