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

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

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.     

Analysis of State-Dependent Transitions in Frequency and Long-Distance Coordination in a Model Oscillatory Cortical Circuit

Changes in behavioral state are typically accompanied by changes in the frequency and spatial coordination of rhythmic activity in the neocortex. In this article, we analyze the effects of neuromodulation on ionic conductances in an oscillating cortical circuit model. The model consists of synaptically-coupled excitatory and inhibitory neurons and supports rhythmic activity in the alpha, beta, and gamma ranges. We find that the effects of neuromodulation on ionic conductances are, by themselves, sufficient to induce transitions between synchronous gamma and beta rhythms and asynchronous alpha rhythms. Moreover, these changes are consistent with changes in behavioral state, with the rhythm transitioning from the slower alpha to the faster gamma and beta as arousal increases. We also observe that it is the same set of underlying intrinsic and network mechanisms that appear to be simultaneously responsible for both the observed transitions between the rhythm types and between their synchronization properties. Spike time response curves (STRCs) are used to study the relationship between the transitions in rhythm and the underlying biophysics.     

Effective models of periodically driven networks

Circadian rhythms are governed by a highly coupled, complex network of genes. Due to feedback within the network, any modification of the system's state requires coherent changes in several nodes. A model of the underlying network is necessary to compute these modifications. We use an effective modeling approach for this task. Rather than inferred biochemical interactions, our method utilizes microarray data from a group of mutants for its construction. With simulated data, we develop an effective model for a circadian network in a peripheral tissue, subject to driving by the suprachiasmatic nucleus, the mammalian pacemaker. The effective network can predict time-dependent gene expression levels in other mutants.     

Dynamics of recurrent neural networks with delayed unreliable synapses: metastable clustering

The influence of unreliable synapses on the dynamic properties of a neural network is investigated for a homogeneous integrate-and-fire network with delayed inhibitory synapses. Numerical and analytical calculations show that the network relaxes to a state with dynamic clusters of identical size which permanently exchange neurons. We present analytical results for the number of clusters and their distribution of firing times which are determined by the synaptic properties. The number of possible configurations increases exponentially with network size. In addition to states with a maximal number of clusters, metastable ones with a smaller number of clusters survive for an exponentially large time scale. An externally excited cluster survives for some time, too, thus clusters may encode information.     

Synchronized oscillations in the visual cortex — a synergetic model

We present an oscillator network model for the synchronization of oscillatory neuronal activity underlying visual processing. The single neuron is modeled by means of a limit cycle oscillator with an eigenfrequency corresponding to visual stimulation. The eigenfrequency may be time dependent. The mutual coupling strengths are unsymmetrical and activity dependent, and they scatter within the network. Synchronized clusters (groups) of neurons emerge in the network due to the visual stimulation. The different clusters correspond to different visual stimuli. There is no limitation of the number of stimuli. Distinct clusters do not perturb each other, although the coupling strength between all model neurons is of the same order of magnitude. Our analysis is not restricted to weak coupling strength. The scatter of the couplings causes shifts of the cluster frequencies. The model's behavior is compared with the experimental findings. The coupling mechanism is extended in order to model the influence of bicucullin upon the neural network. We additionally investigate repulsive couplings, which lead to constant phase differences between clusters of the same frequency. Finally, we consider the problem of selective attention from the viewpoint of our model.     

Synchronized oscillations in the visual cortex — a synergetic model

 We present an oscillator network model for the synchronization of oscillatory neuronal activity underlying visual processing. The single neuron is modeled by means of a limit cycle oscillator with an eigenfrequency corresponding to visual stimulation. The eigenfrequency may be time dependent. The mutual coupling strengths are unsymmetrical and activity dependent, and they scatter within the network. Synchronized clusters (groups) of neurons emerge in the network due to the visual stimulation. The different clusters correspond to different visual stimuli. There is no limitation of the number of stimuli. Distinct clusters do not perturb each other, although the coupling strength between all model neurons is of the same order of magnitude. Our analysis is not restricted to weak coupling strength. The scatter of the couplings causes shifts of the cluster frequencies. The model’s behavior is compared with the experimental findings. The coupling mechanism is extended in order to model the influence of bicucullin upon the neural network. We additionally investigate repulsive couplings, which lead to constant phase differences between clusters of the same frequency. Finally, we consider the problem of selective attention from the viewpoint of our model. Received: 15 February 1995/Accepted in revised form: 18 July 1995     

Minimal size of cell assemblies coordinated by gamma oscillations

In networks of excitatory and inhibitory neurons with mutual synaptic coupling, specific drive to sub-ensembles of cells often leads to gamma-frequency (25-100 Hz) oscillations. When the number of driven cells is too small, however, the synaptic interactions may not be strong or homogeneous enough to support the mechanism underlying the rhythm. Using a combination of computational simulation and mathematical analysis, we study the breakdown of gamma rhythms as the driven ensembles become too small, or the synaptic interactions become too weak and heterogeneous. Heterogeneities in drives or synaptic strengths play an important role in the breakdown of the rhythms; nonetheless, we find that the analysis of homogeneous networks yields insight into the breakdown of rhythms in heterogeneous networks. In particular, if parameter values are such that in a homogeneous network, it takes several gamma cycles to converge to synchrony, then in a similar, but realistically heterogeneous network, synchrony breaks down altogether. This leads to the surprising conclusion that in a network with realistic heterogeneity, gamma rhythms based on the interaction of excitatory and inhibitory cell populations must arise either rapidly, or not at all. For given synaptic strengths and heterogeneities, there is a (soft) lower bound on the possible number of cells in an ensemble oscillating at gamma frequency, based simply on the requirement that synaptic interactions between the two cell populations be strong enough. This observation suggests explanations for recent experimental results concerning the modulation of gamma oscillations in macaque primary visual cortex by varying spatial stimulus size or attention level, and for our own experimental results, reported here, concerning the optogenetic modulation of gamma oscillations in kainate-activated hippocampal slices. We make specific predictions about the behavior of pyramidal cells and fast-spiking interneurons in these experiments.     

Mathematical framework for place coding in the auditory system

In the auditory system, tonotopy is postulated to be the substrate for a place code, where sound frequency is encoded by the location of the neurons that fire during the stimulus. Though conceptually simple, the computations that allow for the representation of intensity and complex sounds are poorly understood. Here, a mathematical framework is developed in order to define clearly the conditions that support a place code. To accommodate both frequency and intensity information, the neural network is described as a space with elements that represent individual neurons and clusters of neurons. A mapping is then constructed from acoustic space to neural space so that frequency and intensity are encoded, respectively, by the location and size of the clusters. Algebraic operations -addition and multiplication- are derived to elucidate the rules for representing, assembling, and modulating multi-frequency sound in networks. The resulting outcomes of these operations are consistent with network simulations as well as with electrophysiological and psychophysical data. The analyses show how both frequency and intensity can be encoded with a purely place code, without the need for rate or temporal coding schemes. The algebraic operations are used to describe loudness summation and suggest a mechanism for the critical band. The mathematical approach complements experimental and computational approaches and provides a foundation for interpreting data and constructing models.     

The generation of phase differences and frequency changes in a network model of inferior olive subthreshold oscillations

It is commonly accepted that the Inferior Olive (IO) provides a timing signal to the cerebellum. Stable subthreshold oscillations in the IO can facilitate accurate timing by phase-locking spikes to the peaks of the oscillation. Several theoretical models accounting for the synchronized subthreshold oscillations have been proposed, however, two experimental observations remain an enigma. The first is the observation of frequent alterations in the frequency of the oscillations. The second is the observation of constant phase differences between simultaneously recorded neurons. In order to account for these two observations we constructed a canonical network model based on anatomical and physiological data from the IO. The constructed network is characterized by clustering of neurons with similar conductance densities, and by electrical coupling between neurons. Neurons inside a cluster are densely connected with weak strengths, while neurons belonging to different clusters are sparsely connected with stronger connections. We found that this type of network can robustly display stable subthreshold oscillations. The overall frequency of the network changes with the strength of the inter-cluster connections, and phase differences occur between neurons of different clusters. Moreover, the phase differences provide a mechanistic explanation for the experimentally observed propagating waves of activity in the IO. We conclude that the architecture of the network of electrically coupled neurons in combination with modulation of the inter-cluster coupling strengths can account for the experimentally observed frequency changes and the phase differences.     

Synchronization properties of networks of electrically coupled neurons in the presence of noise and heterogeneities

We investigate how synchrony can be generated or induced in networks of electrically coupled integrate-and-fire neurons subject to noisy and heterogeneous inputs. Using analytical tools, we find that in a network under constant external inputs, synchrony can appear via a Hopf bifurcation from the asynchronous state to an oscillatory state. In a homogeneous net work, in the oscillatory state all neurons fire in synchrony, while in a heterogeneous network synchrony is looser, many neurons skipping cycles of the oscillation. If the transmission of action potentials via the electrical synapses is effectively excitatory, the Hopf bifurcation is supercritical, while effectively inhibitory transmission due to pronounced hyperpolarization leads to a subcritical bifurcation. In the latter case, the network exhibits bistability between an asynchronous state and an oscillatory state where all the neurons fire in synchrony. Finally we show that for time-varying external inputs, electrical coupling enhances the synchronization in an asynchronous network via a resonance at the firing-rate frequency.

Dynamics of Sparsely Connected Networks of Excitatory and Inhibitory Spiking Neurons

The dynamics of networks of sparsely connected excitatory and inhibitory integrate-and-fire neurons are studied analytically. The analysis reveals a rich repertoire of states, including synchronous states in which neurons fire regularly; asynchronous states with stationary global activity and very irregular individual cell activity; and states in which the global activity oscillates but individual cells fire irregularly, typically at rates lower than the global oscillation frequency. The network can switch between these states, provided the external frequency, or the balance between excitation and inhibition, is varied. Two types of network oscillations are observed. In the fast oscillatory state, the network frequency is almost fully controlled by the synaptic time scale. In the slow oscillatory state, the network frequency depends mostly on the membrane time constant. Finite size effects in the asynchronous state are also discussed.     

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

How stable synchrony in neuronal networks is sustained in the presence of conduction delays is an open question. The Dynamic Clamp was used to measure phase resetting curves (PRCs) for entorhinal cortical cells, and then to construct networks of two such neurons. PRCs were in general Type I (all advances or all delays) or weakly type II with a small region at early phases with the opposite type of resetting. We used previously developed theoretical methods based on PRCs under the assumption of pulsatile coupling to predict the delays that synchronize these hybrid circuits. For excitatory coupling, synchrony was predicted and observed only with no delay and for delays greater than half a network period that cause each neuron to receive an input late in its firing cycle and almost immediately fire an action potential. Synchronization for these long delays was surprisingly tight and robust to the noise and heterogeneity inherent in a biological system. In contrast to excitatory coupling, inhibitory coupling led to antiphase for no delay, very short delays and delays close to a network period, but to near-synchrony for a wide range of relatively short delays. PRC-based methods show that conduction delays can stabilize synchrony in several ways, including neutralizing a discontinuity introduced by strong inhibition, favoring synchrony in the case of noisy bistability, and avoiding an initial destabilizing region of a weakly type II PRC. PRCs can identify optimal conduction delays favoring synchronization at a given frequency, and also predict robustness to noise and heterogeneity.     

Synchronous Clusters in a Noisy Inhibitory Neural Network

We study the stability and information encoding capacity of synchronized states in a neuronal network model that represents part of thalamic circuitry. Our model neurons have a Hodgkin-Huxley-type low-threshold calcium channel, display postinhibitory rebound, and are connected via GABAergic inhibitory synapses.We find that there is a threshold in synaptic strength, _c, below which there are no stable spiking network states. Above threshold the stable spiking state is a cluster state, where different groups of neurons fire consecutively, and each neuron fires with the same cluster each time. Weak noise destabilizes this state, but stronger noise drives the system into a different, self-organized, stochastically synchronized state. Neuronal firing is still organized in clusters, but individual neurons can hop from cluster to cluster. Noise can actually induce and sustain such a state below the threshold of synaptic strength. We do find a qualitative difference in the firing patterns between small (10 neurons) and large (1000 neurons) networks.We determine the information content of the spike trains in terms of two separate contributions: the spike-time jitter around cluster firing times, and the hopping from cluster to cluster. We quantify the information loss due to temporally correlated interspike intervals. Recent experiments on the locust olfactory system and striatal neurons suggest that the nervous system may actually use these two channels to encode separate and unique information.     

Cortical network modeling: analytical methods for firing rates and some properties of networks of LIF neurons.

The circuitry of cortical networks involves interacting populations of excitatory (E) and inhibitory (I) neurons whose relationships are now known to a large extent. Inputs to E- and I-cells may have their origins in remote or local cortical areas. We consider a rudimentary model involving E- and I-cells. One of our goals is to test an analytic approach to finding firing rates in neural networks without using a diffusion approximation and to this end we consider in detail networks of excitatory neurons with leaky integrate-and-fire (LIF) dynamics. A simple measure of synchronization, denoted by S(q), where q is between 0 and 100 is introduced. Fully connected E-networks have a large tendency to become dominated by synchronously firing groups of cells, except when inputs are relatively weak. We observed random or asynchronous firing in such networks with diverse sets of parameter values. When such firing patterns were found, the analytical approach was often able to accurately predict average neuronal firing rates. We also considered several properties of E-E networks, distinguishing several kinds of firing pattern. Included were those with silences before or after periods of intense activity or with periodic synchronization. We investigated the occurrence of synchronized firing with respect to changes in the internal excitatory postsynaptic potential (EPSP) magnitude in a network of 100 neurons with fixed values of the remaining parameters. When the internal EPSP size was less than a certain value, synchronization was absent. The amount of synchronization then increased slowly as the EPSP amplitude increased until at a particular EPSP size the amount of synchronization abruptly increased, with S(5) attaining the maximum value of 100%. We also found network frequency transfer characteristics for various network sizes and found a linear dependence of firing frequency over wide ranges of the external afferent frequency, with non-linear effects at lower input frequencies. The theory may also be applied to sparsely connected networks, whose firing behaviour was found to change abruptly as the probability of a connection passed through a critical value. The analytical method was also found to be useful for a feed-forward excitatory network and a network of excitatory and inhibitory neurons.     

Cell Type-Specific Separation of Subicular Principal Neurons during Network Activities

The hippocampal output structure, the subiculum, expresses two major memory relevant network rhythms, sharp wave ripple and gamma frequency oscillations. To this date, it remains unclear how the two distinct types of subicular principal cells, intrinsically bursting and regular spiking neurons, participate in these two network rhythms. Using concomitant local field potential and intracellular recordings in an in vitro mouse model that allows the investigation of both network rhythms, we found a cell type-specific segregation of principal neurons into participating intrinsically bursting and non-participating regular spiking cells. However, if regular spiking cells were kept at a more depolarized level, they did participate in a specific manner, suggesting a potential bimodal working model dependent on the level of excitation. Furthermore, intrinsically bursting and regular spiking cells exhibited divergent intrinsic membrane and synaptic properties in the active network. Thus, our results suggest a cell-type-specific segregation of principal cells into two separate groups during network activities, supporting the idea of two parallel streams of information processing within the subiculum.     

Network Adaptation Improves Temporal Representation of Naturalistic Stimuli in Drosophila Eye: II Mechanisms

Retinal networks must adapt constantly to best present the ever changing visual world to the brain. Here we test the hypothesis that adaptation is a result of different mechanisms at several synaptic connections within the network. In a companion paper (Part I), we showed that adaptation in the photoreceptors (R1–R6) and large monopolar cells (LMC) of the Drosophila eye improves sensitivity to under-represented signals in seconds by enhancing both the amplitude and frequency distribution of LMCs'' voltage responses to repeated naturalistic contrast series. In this paper, we show that such adaptation needs both the light-mediated conductance and feedback-mediated synaptic conductance. A faulty feedforward pathway in histamine receptor mutant flies speeds up the LMC output, mimicking extreme light adaptation. A faulty feedback pathway from L2 LMCs to photoreceptors slows down the LMC output, mimicking dark adaptation. These results underline the importance of network adaptation for efficient coding, and as a mechanism for selectively regulating the size and speed of signals in neurons. We suggest that concert action of many different mechanisms and neural connections are responsible for adaptation to visual stimuli. Further, our results demonstrate the need for detailed circuit reconstructions like that of the Drosophila lamina, to understand how networks process information.     

Modeling the Seasonal Adaptation of Circadian Clocks by Changes in the Network Structure of the Suprachiasmatic Nucleus

The dynamics of circadian rhythms needs to be adapted to day length changes between summer and winter. It has been observed experimentally, however, that the dynamics of individual neurons of the suprachiasmatic nucleus (SCN) does not change as the seasons change. Rather, the seasonal adaptation of the circadian clock is hypothesized to be a consequence of changes in the intercellular dynamics, which leads to a phase distribution of electrical activity of SCN neurons that is narrower in winter and broader during summer. Yet to understand this complex intercellular dynamics, a more thorough understanding of the impact of the network structure formed by the SCN neurons is needed. To that effect, we propose a mathematical model for the dynamics of the SCN neuronal architecture in which the structure of the network plays a pivotal role. Using our model we show that the fraction of long-range cell-to-cell connections and the seasonal changes in the daily rhythms may be tightly related. In particular, simulations of the proposed mathematical model indicate that the fraction of long-range connections between the cells adjusts the phase distribution and consequently the length of the behavioral activity as follows: dense long-range connections during winter lead to a narrow activity phase, while rare long-range connections during summer lead to a broad activity phase. Our model is also able to account for the experimental observations indicating a larger light-induced phase-shift of the circadian clock during winter, which we show to be a consequence of higher synchronization between neurons. Our model thus provides evidence that the variations in the seasonal dynamics of circadian clocks can in part also be understood and regulated by the plasticity of the SCN network structure.