Sparsely synchronized neuronal oscillations

Frequency-domain order parameters for the burst and spike synchronization transitions of bursting neurons

We are interested in characterization of synchronization transitions of bursting neurons in the frequency domain. Instantaneous population firing rate (IPFR) [Formula: see text], which is directly obtained from the raster plot of neural spikes, is often used as a realistic collective quantity describing population activities in both the computational and the experimental neuroscience. For the case of spiking neurons, a realistic time-domain order parameter, based on [Formula: see text], was introduced in our recent work to characterize the spike synchronization transition. Unlike the case of spiking neurons, the IPFR [Formula: see text] of bursting neurons exhibits population behaviors with both the slow bursting and the fast spiking timescales. For our aim, we decompose the IPFR [Formula: see text] into the instantaneous population bursting rate [Formula: see text] (describing the bursting behavior) and the instantaneous population spike rate [Formula: see text] (describing the spiking behavior) via frequency filtering, and extend the realistic order parameter to the case of bursting neurons. Thus, we develop the frequency-domain bursting and spiking order parameters which are just the bursting and spiking "coherence factors" [Formula: see text] and [Formula: see text] of the bursting and spiking peaks in the power spectral densities of [Formula: see text] and [Formula: see text] (i.e., "signal to noise" ratio of the spectral peak height and its relative width). Through calculation of [Formula: see text] and [Formula: see text], we obtain the bursting and spiking thresholds beyond which the burst and spike synchronizations break up, respectively. Consequently, it is shown in explicit examples that the frequency-domain bursting and spiking order parameters may be usefully used for characterization of the bursting and the spiking transitions, respectively.

Bursting reverberation as a multiscale neuronal network process driven by synaptic depression-facilitation

Neuronal networks can generate complex patterns of activity that depend on membrane properties of individual neurons as well as on functional synapses. To decipher the impact of synaptic properties and connectivity on neuronal network behavior, we investigate the responses of neuronal ensembles from small (5-30 cells in a restricted sphere) and large (acute hippocampal slice) networks to single electrical stimulation: in both cases, a single stimulus generated a synchronous long-lasting bursting activity. While an initial spike triggered a reverberating network activity that lasted 2-5 seconds for small networks, we found here that it lasted only up to 300 milliseconds in slices. To explain this phenomena present at different scales, we generalize the depression-facilitation model and extracted the network time constants. The model predicts that the reverberation time has a bell shaped relation with the synaptic density, revealing that the bursting time cannot exceed a maximum value. Furthermore, before reaching its maximum, the reverberation time increases sub-linearly with the synaptic density of the network. We conclude that synaptic dynamics and connectivity shape the mean burst duration, a property present at various scales of the networks. Thus bursting reverberation is a property of sufficiently connected neural networks, and can be generated by collective depression and facilitation of underlying functional synapses.

Noise-induced burst and spike synchronizations in an inhibitory small-world network of subthreshold bursting neurons

We are interested in noise-induced firings of subthreshold neurons which may be used for encoding environmental stimuli. Noise-induced population synchronization was previously studied only for the case of global coupling, unlike the case of subthreshold spiking neurons. Hence, we investigate the effect of complex network architecture on noise-induced synchronization in an inhibitory population of subthreshold bursting Hindmarsh-Rose neurons. For modeling complex synaptic connectivity, we consider the Watts-Strogatz small-world network which interpolates between regular lattice and random network via rewiring, and investigate the effect of small-world connectivity on emergence of noise-induced population synchronization. Thus, noise-induced burst synchronization (synchrony on the slow bursting time scale) and spike synchronization (synchrony on the fast spike time scale) are found to appear in a synchronized region of the [Formula: see text]-[Formula: see text] plane ([Formula: see text]: synaptic inhibition strength and [Formula: see text]: noise intensity). As the rewiring probability [Formula: see text] is decreased from 1 (random network) to 0 (regular lattice), the region of spike synchronization shrinks rapidly in the [Formula: see text]-[Formula: see text] plane, while the region of the burst synchronization decreases slowly. We separate the slow bursting and the fast spiking time scales via frequency filtering, and characterize the noise-induced burst and spike synchronizations by employing realistic order parameters and statistical-mechanical measures introduced in our recent work. Thus, the bursting and spiking thresholds for the burst and spike synchronization transitions are determined in terms of the bursting and spiking order parameters, respectively. Furthermore, we also measure the degrees of burst and spike synchronizations in terms of the statistical-mechanical bursting and spiking measures, respectively.

Oscillations emerging from noise-driven steady state in networks with electrical synapses and subthreshold resonance

Oscillations play a critical role in cognitive phenomena and have been observed in many brain regions. Experimental evidence indicates that classes of neurons exhibit properties that could promote oscillations, such as subthreshold resonance and electrical gap junctions. Typically, these two properties are studied separately but it is not clear which is the dominant determinant of global network rhythms. Our aim is to provide an analytical understanding of how these two effects destabilize the fluctuation-driven state, in which neurons fire irregularly, and lead to an emergence of global synchronous oscillations. Here we show how the oscillation frequency is shaped by single neuron resonance, electrical and chemical synapses.The presence of both gap junctions and subthreshold resonance are necessary for the emergence of oscillations. Our results are in agreement with several experimental observations such as network responses to oscillatory inputs and offer a much-needed conceptual link connecting a collection of disparate effects observed in networks.

Self-consistent determination of the spike-train power spectrum in a neural network with sparse connectivity

A major source of random variability in cortical networks is the quasi-random arrival of presynaptic action potentials from many other cells. In network studies as well as in the study of the response properties of single cells embedded in a network, synaptic background input is often approximated by Poissonian spike trains. However, the output statistics of the cells is in most cases far from being Poisson. This is inconsistent with the assumption of similar spike-train statistics for pre- and postsynaptic cells in a recurrent network. Here we tackle this problem for the popular class of integrate-and-fire neurons and study a self-consistent statistics of input and output spectra of neural spike trains. Instead of actually using a large network, we use an iterative scheme, in which we simulate a single neuron over several generations. In each of these generations, the neuron is stimulated with surrogate stochastic input that has a similar statistics as the output of the previous generation. For the surrogate input, we employ two distinct approximations: (i) a superposition of renewal spike trains with the same interspike interval density as observed in the previous generation and (ii) a Gaussian current with a power spectrum proportional to that observed in the previous generation. For input parameters that correspond to balanced input in the network, both the renewal and the Gaussian iteration procedure converge quickly and yield comparable results for the self-consistent spike-train power spectrum. We compare our results to large-scale simulations of a random sparsely connected network of leaky integrate-and-fire neurons (Brunel, 2000) and show that in the asynchronous regime close to a state of balanced synaptic input from the network, our iterative schemes provide an excellent approximations to the autocorrelation of spike trains in the recurrent network.

A canonical circuit for generating phase-amplitude coupling

‘Phase amplitude coupling’ (PAC) in oscillatory neural activity describes a phenomenon whereby the amplitude of higher frequency activity is modulated by the phase of lower frequency activity. Such coupled oscillatory activity – also referred to as ‘cross-frequency coupling’ or ‘nested rhythms’ – has been shown to occur in a number of brain regions and at behaviorally relevant time points during cognitive tasks; this suggests functional relevance, but the circuit mechanisms of PAC generation remain unclear. In this paper we present a model of a canonical circuit for generating PAC activity, showing how interconnected excitatory and inhibitory neural populations can be periodically shifted in to and out of oscillatory firing patterns by afferent drive, hence generating higher frequency oscillations phase-locked to a lower frequency, oscillating input signal. Since many brain regions contain mutually connected excitatory-inhibitory populations receiving oscillatory input, the simplicity of the mechanism generating PAC in such networks may explain the ubiquity of PAC across diverse neural systems and behaviors. Analytic treatment of this circuit as a nonlinear dynamical system demonstrates how connection strengths and inputs to the populations can be varied in order to change the extent and nature of PAC activity, importantly which phase of the lower frequency rhythm the higher frequency activity is locked to. Consequently, this model can inform attempts to associate distinct types of PAC with different network topologies and physiologies in real data.

Mechanisms of zero-lag synchronization in cortical motifs

Zero-lag synchronization between distant cortical areas has been observed in a diversity of experimental data sets and between many different regions of the brain. Several computational mechanisms have been proposed to account for such isochronous synchronization in the presence of long conduction delays: Of these, the phenomenon of “dynamical relaying” – a mechanism that relies on a specific network motif – has proven to be the most robust with respect to parameter mismatch and system noise. Surprisingly, despite a contrary belief in the community, the common driving motif is an unreliable means of establishing zero-lag synchrony. Although dynamical relaying has been validated in empirical and computational studies, the deeper dynamical mechanisms and comparison to dynamics on other motifs is lacking. By systematically comparing synchronization on a variety of small motifs, we establish that the presence of a single reciprocally connected pair – a “resonance pair” – plays a crucial role in disambiguating those motifs that foster zero-lag synchrony in the presence of conduction delays (such as dynamical relaying) from those that do not (such as the common driving triad). Remarkably, minor structural changes to the common driving motif that incorporate a reciprocal pair recover robust zero-lag synchrony. The findings are observed in computational models of spiking neurons, populations of spiking neurons and neural mass models, and arise whether the oscillatory systems are periodic, chaotic, noise-free or driven by stochastic inputs. The influence of the resonance pair is also robust to parameter mismatch and asymmetrical time delays amongst the elements of the motif. We call this manner of facilitating zero-lag synchrony resonance-induced synchronization, outline the conditions for its occurrence, and propose that it may be a general mechanism to promote zero-lag synchrony in the brain.

Understanding large-scale neuronal dynamics – and how they relate to the cortical anatomy – is one of the key areas of neuroscience research. Despite a wealth of recent research, the key principles of this relationship have yet to be established. Here we employ computational modeling to study neuronal dynamics on small subgraphs – or motifs – across a hierarchy of spatial scales. We establish a novel organizing principle that we term a “resonance pair” (two mutually coupled nodes), which promotes stable, zero-lag synchrony amongst motif nodes. The bidirectional coupling between a resonance pair acts to mutually adjust their dynamics onto a common and relatively stable synchronized regime, which then propagates and stabilizes the synchronization of other nodes within the motif. Remarkably, we find that this effect can propagate along chains of coupled nodes and hence holds the potential to promote stable zero-lag synchrony in larger sub-networks of cortical systems. Our findings hence suggest a potential unifying account of the existence of zero-lag synchrony, an important phenomenon that may underlie crucial cognitive processes in the brain. Moreover, such pairs of mutually coupled oscillators are found in a wide variety of physical and biological systems suggesting a new, broadly relevant and unifying principle.

Small modifications to network topology can induce stochastic bistable spiking dynamics in a balanced cortical model

Directed random graph models frequently are used successfully in modeling the population dynamics of networks of cortical neurons connected by chemical synapses. Experimental results consistently reveal that neuronal network topology is complex, however, in the sense that it differs statistically from a random network, and differs for classes of neurons that are physiologically different. This suggests that complex network models whose subnetworks have distinct topological structure may be a useful, and more biologically realistic, alternative to random networks. Here we demonstrate that the balanced excitation and inhibition frequently observed in small cortical regions can transiently disappear in otherwise standard neuronal-scale models of fluctuation-driven dynamics, solely because the random network topology was replaced by a complex clustered one, whilst not changing the in-degree of any neurons. In this network, a small subset of cells whose inhibition comes only from outside their local cluster are the cause of bistable population dynamics, where different clusters of these cells irregularly switch back and forth from a sparsely firing state to a highly active state. Transitions to the highly active state occur when a cluster of these cells spikes sufficiently often to cause strong unbalanced positive feedback to each other. Transitions back to the sparsely firing state rely on occasional large fluctuations in the amount of non-local inhibition received. Neurons in the model are homogeneous in their intrinsic dynamics and in-degrees, but differ in the abundance of various directed feedback motifs in which they participate. Our findings suggest that (i) models and simulations should take into account complex structure that varies for neuron and synapse classes; (ii) differences in the dynamics of neurons with similar intrinsic properties may be caused by their membership in distinctive local networks; (iii) it is important to identify neurons that share physiological properties and location, but differ in their connectivity.

Single neuron dynamics and computation

At the single neuron level, information processing involves the transformation of input spike trains into an appropriate output spike train. Building upon the classical view of a neuron as a threshold device, models have been developed in recent years that take into account the diverse electrophysiological make-up of neurons and accurately describe their input-output relations. Here, we review these recent advances and survey the computational roles that they have uncovered for various electrophysiological properties, for dendritic arbor anatomy as well as for short-term synaptic plasticity.

Oscillatory multiplexing of population codes for selective communication in the mammalian brain

Mammalian brains exhibit population oscillations, the structures of which vary in time and space according to behavioural state. A proposed function of these oscillations is to control the flow of signals among anatomically connected networks. However, the nature of neural coding that may support selective communication that depends on oscillations has received relatively little attention. Here, we consider the role of multiplexing, whereby multiple information streams share a common neural substrate. We suggest that multiplexing implemented through periodic modulation of firing-rate population codes enables flexible reconfiguration of effective connectivity among brain areas.

Critical phenomena and noise-induced phase transitions in neuronal networks

We study numerically and analytically first- and second-order phase transitions in neuronal networks stimulated by shot noise (a flow of random spikes bombarding neurons). Using an exactly solvable cortical model of neuronal networks on classical random networks, we find critical phenomena accompanying the transitions and their dependence on the shot noise intensity. We show that a pattern of spontaneous neuronal activity near a critical point of a phase transition is a characteristic property that can be used to identify the bifurcation mechanism of the transition. We demonstrate that bursts and avalanches are precursors of a first-order phase transition, paroxysmal-like spikes of activity precede a second-order phase transition caused by a saddle-node bifurcation, while irregular spindle oscillations represent spontaneous activity near a second-order phase transition caused by a supercritical Hopf bifurcation. Our most interesting result is the observation of the paroxysmal-like spikes. We show that a paroxysmal-like spike is a single nonlinear event that appears instantly from a low background activity with a rapid onset, reaches a large amplitude, and ends up with an abrupt return to lower activity. These spikes are similar to single paroxysmal spikes and sharp waves observed in electroencephalographic (EEG) measurements. Our analysis shows that above the saddle-node bifurcation, sustained network oscillations appear with a large amplitude but a small frequency in contrast to network oscillations near the Hopf bifurcation that have a small amplitude but a large frequency. We discuss an amazing similarity between excitability of the cortical model stimulated by shot noise and excitability of the Morris-Lecar neuron stimulated by an applied current.

Efficient "communication through coherence" requires oscillations structured to minimize interference between signals

The ‘communication through coherence’ (CTC) hypothesis proposes that selective communication among neural networks is achieved by coherence between firing rate oscillation in a sending region and gain modulation in a receiving region. Although this hypothesis has stimulated extensive work, it remains unclear whether the mechanism can in principle allow reliable and selective information transfer. Here we use a simple mathematical model to investigate how accurately coherent gain modulation can filter a population-coded target signal from task-irrelevant distracting inputs. We show that selective communication can indeed be achieved, although the structure of oscillatory activity in the target and distracting networks must satisfy certain previously unrecognized constraints. Firstly, the target input must be differentiated from distractors by the amplitude, phase or frequency of its oscillatory modulation. When distracting inputs oscillate incoherently in the same frequency band as the target, communication accuracy is severely degraded because of varying overlap between the firing rate oscillations of distracting inputs and the gain modulation in the receiving region. Secondly, the oscillatory modulation of the target input must be strong in order to achieve a high signal-to-noise ratio relative to stochastic spiking of individual neurons. Thus, whilst providing a quantitative demonstration of the power of coherent oscillatory gain modulation to flexibly control information flow, our results identify constraints imposed by the need to avoid interference between signals, and reveal a likely organizing principle for the structure of neural oscillations in the brain.

Distributed regions of mammalian brains transiently engage in coherent oscillations, often at specific stages of behavioral or cognitive tasks. This activity may play a role in controlling information flow among connected regions, allowing the brain's connectivity structure to be flexibly reconfigured in response to changing task demands. We have used a computational model to investigate the conditions under which oscillations can generate selective communication through a mechanism in which the excitability of neurons in one region is modulated coherently with a firing rate oscillation in another region. Our results demonstrate that this mechanism is able to accurately and selectively control the flow of signals encoded as spatial patterns of firing rate. However, we found that the requirement to avoid interference between different signals imposes previously unrecognised constraints on the structures of oscillatory activity that can efficiently support this mechanism. These constraints may be an organizing principle for the structured oscillatory activity observed in vivo.

Stochastic amplification of fluctuations in cortical up-states

Cortical neurons are bistable; as a consequence their local field potentials can fluctuate between quiescent and active states, generating slow Hz oscillations which are widely known as transitions between Up and Down States. Despite a large number of studies on Up-Down transitions, deciphering its nature, mechanisms and function are still today challenging tasks. In this paper we focus on recent experimental evidence, showing that a class of spontaneous oscillations can emerge within the Up states. In particular, a non-trivial peak around Hz appears in their associated power-spectra, what produces an enhancement of the activity power for higher frequencies (in the Hz band). Moreover, this rhythm within Ups seems to be an emergent or collective phenomenon given that individual neurons do not lock to it as they remain mostly unsynchronized. Remarkably, similar oscillations (and the concomitant peak in the spectrum) do not appear in the Down states. Here we shed light on these findings by using different computational models for the dynamics of cortical networks in presence of different levels of physiological complexity. Our conclusion, supported by both theory and simulations, is that the collective phenomenon of “stochastic amplification of fluctuations” – previously described in other contexts such as Ecology and Epidemiology – explains in an elegant and parsimonious manner, beyond model-dependent details, this extra-rhythm emerging only in the Up states but not in the Downs.

Oscillatory dynamics in the hippocampus support dentate gyrus–CA3 coupling

Gamma oscillations in the dentate gyrus and hippocampal CA3 show variable coherence in vivo, but the mechanisms and relevance for information flow are unknown. We found that carbachol-induced oscillations in rat CA3 have biphasic phase-response curves, consistent with the ability to couple with oscillations in afferent projections. Differences in response to stimulation of either the intrinsic feedback circuit or the dentate gyrus were well described by varying an impulse vector in a two-dimensional dynamical system, representing the relative input to excitatory and inhibitory neurons. Responses to sinusoidally modulated optogenetic stimulation confirmed that the CA3 network oscillation can entrain to periodic inputs, with a steep dependence of entrainment phase on input frequency. CA3 oscillations are therefore suited to coupling with oscillations in the dentate gyrus over a broad range of frequencies.

Dynamic effective connectivity of inter-areal brain circuits

Anatomic connections between brain areas affect information flow between neuronal circuits and the synchronization of neuronal activity. However, such structural connectivity does not coincide with effective connectivity (or, more precisely, causal connectivity), related to the elusive question “Which areas cause the present activity of which others?”. Effective connectivity is directed and depends flexibly on contexts and tasks. Here we show that dynamic effective connectivity can emerge from transitions in the collective organization of coherent neural activity. Integrating simulation and semi-analytic approaches, we study mesoscale network motifs of interacting cortical areas, modeled as large random networks of spiking neurons or as simple rate units. Through a causal analysis of time-series of model neural activity, we show that different dynamical states generated by a same structural connectivity motif correspond to distinct effective connectivity motifs. Such effective motifs can display a dominant directionality, due to spontaneous symmetry breaking and effective entrainment between local brain rhythms, although all connections in the considered structural motifs are reciprocal. We show then that transitions between effective connectivity configurations (like, for instance, reversal in the direction of inter-areal interactions) can be triggered reliably by brief perturbation inputs, properly timed with respect to an ongoing local oscillation, without the need for plastic synaptic changes. Finally, we analyze how the information encoded in spiking patterns of a local neuronal population is propagated across a fixed structural connectivity motif, demonstrating that changes in the active effective connectivity regulate both the efficiency and the directionality of information transfer. Previous studies stressed the role played by coherent oscillations in establishing efficient communication between distant areas. Going beyond these early proposals, we advance here that dynamic interactions between brain rhythms provide as well the basis for the self-organized control of this “communication-through-coherence”, making thus possible a fast “on-demand” reconfiguration of global information routing modalities.

The circuits of the brain must perform a daunting amount of functions. But how can “brain states” be flexibly controlled, given that anatomic inter-areal connections can be considered as fixed, on timescales relevant for behavior? We hypothesize that, thanks to the nonlinear interaction between brain rhythms, even a simple circuit involving few brain areas can originate a multitude of effective circuits, associated with alternative functions selectable “on demand”. A distinction is usually made between structural connectivity, which describes actual synaptic connections, and effective connectivity, quantifying, beyond correlation, directed inter-areal causal influences. In our study, we measure effective connectivity based on time-series of neural activity generated by model inter-areal circuits. We find that “causality follows dynamics”. We show indeed that different effective networks correspond to different dynamical states associated to a same structural network (in particular, different phase-locking patterns between local neuronal oscillations). We then find that “information follows causality” (and thus, again, dynamics). We demonstrate that different effective networks give rise to alternative modalities of information routing between brain areas wired together in a fixed structural network. In particular, we show that the self-organization of interacting “analog” rate oscillations control the flow of “digital-like” information encoded in complex spiking patterns.

Synchronous chaos and broad band gamma rhythm in a minimal multi-layer model of primary visual cortex

Visually induced neuronal activity in V1 displays a marked gamma-band component which is modulated by stimulus properties. It has been argued that synchronized oscillations contribute to these gamma-band activity. However, analysis of Local Field Potentials (LFPs) across different experiments reveals considerable diversity in the degree of oscillatory behavior of this induced activity. Contrast-dependent power enhancements can indeed occur over a broad band in the gamma frequency range and spectral peaks may not arise at all. Furthermore, even when oscillations are observed, they undergo temporal decorrelation over very few cycles. This is not easily accounted for in previous network modeling of gamma oscillations. We argue here that interactions between cortical layers can be responsible for this fast decorrelation. We study a model of a V1 hypercolumn, embedding a simplified description of the multi-layered structure of the cortex. When the stimulus contrast is low, the induced activity is only weakly synchronous and the network resonates transiently without developing collective oscillations. When the contrast is high, on the other hand, the induced activity undergoes synchronous oscillations with an irregular spatiotemporal structure expressing a synchronous chaotic state. As a consequence the population activity undergoes fast temporal decorrelation, with concomitant rapid damping of the oscillations in LFPs autocorrelograms and peak broadening in LFPs power spectra. We show that the strength of the inter-layer coupling crucially affects this spatiotemporal structure. We predict that layer VI inactivation should induce global changes in the spectral properties of induced LFPs, reflecting their slower temporal decorrelation in the absence of inter-layer feedback. Finally, we argue that the mechanism underlying the emergence of synchronous chaos in our model is in fact very general. It stems from the fact that gamma oscillations induced by local delayed inhibition tend to develop chaos when coupled by sufficiently strong excitation.

Visual stimulation elicits neuronal responses in visual cortex. When the contrast of the used stimuli increases, the power of this induced activity is boosted over a broad frequency range (30–100 Hz), called the “gamma band.” It would be tempting to hypothesize that this phenomenon is due to the emergence of oscillations in which many neurons fire collectively in a rhythmic way. However, previous models trying to explain contrast-related power enhancements using synchronous oscillations failed to reproduce the observed spectra because they originated unrealistically sharp spectral peaks. The aim of our study is to reconcile synchronous oscillations with broad-band power spectra. We argue here that, thanks to the interaction between neuronal populations at different depths in the cortical tissue, the induced oscillatory responses are synchronous, but, at the same time, chaotic. The chaotic nature of the dynamics makes it possible to have broad-band power spectra together with synchrony. Our modeling study allows us formulating qualitative experimental predictions that provide a potential test for our theory. We predict that if the interactions between cortical layers are suppressed, for instance by inactivating neurons in deep layers, the induced responses might become more regular and narrow isolated peaks might develop in their power spectra.

Input-rate modulation of γ oscillations is sensitive to network topology, delays and short-term plasticity

Simulated networks of excitatory and inhibitory neurons have previously been shown to reproduce critical features of experimental data regarding neural coding in V1, such as a positive relationship between thalamic input spike rate and the power of gamma frequency oscillations. This effect, referred to as modulated gamma power, could represent a neural code in V1 for stimulus characteristics that affect thalamic spike rate such as contrast or intensity. The simulated network's assumptions included homogeneous random connectivity, equal synaptic delays after spike arrival, and constant synaptic efficacies. Plausible alternative assumptions include small world connectivity, a wide distribution of axonal propagation delays, and short term synaptic plasticity, and here we assess the individual impact of each of these on the model's success in reproducing modulated gamma power. First, we developed several alternative algorithms for simulating directed networks with clustered connectivity and balanced excitation and inhibition. We found that modulated gamma power was absent in all small-world networks that had a relatively low abundance of reciprocal connectivity, which suggests that such motifs are present in V1 cortical networks at levels at least equal to those found in random networks. We also found in a different network type that the balance of excitation and inhibition could be destroyed when the network was in the small-world regime. Given all neurons had identical in-degrees, this result suggests that balance relies on motif distributions as well as mean connectivity. Second, altering the distribution of axonal delays had little effect, but increasing the mean delay led to a secondary gamma modulation at harmonics of the main peak, and since this is not observed experimentally, it suggests a mean delay in V1 networks less than 2 ms. Finally, we compared two types of excitatory synaptic plasticity, and found that modulated beta power emerged in addition to gamma power for one type, in the presence of short term depression in interneurons. This article is part of a Special Issue entitled "Neural Coding".

Dynamics of networks of excitatory and inhibitory neurons in response to time-dependent inputs

We investigate the dynamics of recurrent networks of excitatory (E) and inhibitory (I) neurons in the presence of time-dependent inputs. The dynamics is characterized by the network dynamical transfer function, i.e., how the population firing rate is modulated by sinusoidal inputs at arbitrary frequencies. Two types of networks are studied and compared: (i) a Wilson-Cowan type firing rate model; and (ii) a fully connected network of leaky integrate-and-fire (LIF) neurons, in a strong noise regime. We first characterize the region of stability of the "asynchronous state" (a state in which population activity is constant in time when external inputs are constant) in the space of parameters characterizing the connectivity of the network. We then systematically characterize the qualitative behaviors of the dynamical transfer function, as a function of the connectivity. We find that the transfer function can be either low-pass, or with a single or double resonance, depending on the connection strengths and synaptic time constants. Resonances appear when the system is close to Hopf bifurcations, that can be induced by two separate mechanisms: the I-I connectivity and the E-I connectivity. Double resonances can appear when excitatory delays are larger than inhibitory delays, due to the fact that two distinct instabilities exist with a finite gap between the corresponding frequencies. In networks of LIF neurons, changes in external inputs and external noise are shown to be able to change qualitatively the network transfer function. Firing rate models are shown to exhibit the same diversity of transfer functions as the LIF network, provided delays are present. They can also exhibit input-dependent changes of the transfer function, provided a suitable static non-linearity is incorporated.

Emergent oscillations in networks of stochastic spiking neurons

Networks of neurons produce diverse patterns of oscillations, arising from the network's global properties, the propensity of individual neurons to oscillate, or a mixture of the two. Here we describe noisy limit cycles and quasi-cycles, two related mechanisms underlying emergent oscillations in neuronal networks whose individual components, stochastic spiking neurons, do not themselves oscillate. Both mechanisms are shown to produce gamma band oscillations at the population level while individual neurons fire at a rate much lower than the population frequency. Spike trains in a network undergoing noisy limit cycles display a preferred period which is not found in the case of quasi-cycles, due to the even faster decay of phase information in quasi-cycles. These oscillations persist in sparsely connected networks, and variation of the network's connectivity results in variation of the oscillation frequency. A network of such neurons behaves as a stochastic perturbation of the deterministic Wilson-Cowan equations, and the network undergoes noisy limit cycles or quasi-cycles depending on whether these have limit cycles or a weakly stable focus. These mechanisms provide a new perspective on the emergence of rhythmic firing in neural networks, showing the coexistence of population-level oscillations with very irregular individual spike trains in a simple and general framework.

The role of degree distribution in shaping the dynamics in networks of sparsely connected spiking neurons

Neuronal network models often assume a fixed probability of connection between neurons. This assumption leads to random networks with binomial in-degree and out-degree distributions which are relatively narrow. Here I study the effect of broad degree distributions on network dynamics by interpolating between a binomial and a truncated power-law distribution for the in-degree and out-degree independently. This is done both for an inhibitory network (I network) as well as for the recurrent excitatory connections in a network of excitatory and inhibitory neurons (EI network). In both cases increasing the width of the in-degree distribution affects the global state of the network by driving transitions between asynchronous behavior and oscillations. This effect is reproduced in a simplified rate model which includes the heterogeneity in neuronal input due to the in-degree of cells. On the other hand, broadening the out-degree distribution is shown to increase the fraction of common inputs to pairs of neurons. This leads to increases in the amplitude of the cross-correlation (CC) of synaptic currents. In the case of the I network, despite strong oscillatory CCs in the currents, CCs of the membrane potential are low due to filtering and reset effects, leading to very weak CCs of the spike-count. In the asynchronous regime of the EI network, broadening the out-degree increases the amplitude of CCs in the recurrent excitatory currents, while CC of the total current is essentially unaffected as are pairwise spiking correlations. This is due to a dynamic balance between excitatory and inhibitory synaptic currents. In the oscillatory regime, changes in the out-degree can have a large effect on spiking correlations and even on the qualitative dynamical state of the network.

Oscillations and filtering networks support flexible routing of information

The mammalian brain exhibits profuse interregional connectivity. How information flow is rapidly and flexibly switched among connected areas remains poorly understood. Task-dependent changes in the power and interregion coherence of network oscillations suggest that such oscillations play a role in signal routing. We show that switching one of several convergent pathways from an asynchronous to an oscillatory state allows accurate selective transmission of population-coded information, which can be extracted even when other convergent pathways fire asynchronously at comparable rates. We further show that the band-pass filtering required to perform this information extraction can be implemented in a simple spiking network model with a single feed-forward interneuron layer. This constitutes a mechanism for flexible signal routing in neural circuits, which exploits sparsely synchronized network oscillations and temporal filtering by feed-forward inhibition.

Neurophysiological and computational principles of cortical rhythms in cognition

Synchronous rhythms represent a core mechanism for sculpting temporal coordination of neural activity in the brain-wide network. This review focuses on oscillations in the cerebral cortex that occur during cognition, in alert behaving conditions. Over the last two decades, experimental and modeling work has made great strides in elucidating the detailed cellular and circuit basis of these rhythms, particularly gamma and theta rhythms. The underlying physiological mechanisms are diverse (ranging from resonance and pacemaker properties of single cells to multiple scenarios for population synchronization and wave propagation), but also exhibit unifying principles. A major conceptual advance was the realization that synaptic inhibition plays a fundamental role in rhythmogenesis, either in an interneuronal network or in a reciprocal excitatory-inhibitory loop. Computational functions of synchronous oscillations in cognition are still a matter of debate among systems neuroscientists, in part because the notion of regular oscillation seems to contradict the common observation that spiking discharges of individual neurons in the cortex are highly stochastic and far from being clocklike. However, recent findings have led to a framework that goes beyond the conventional theory of coupled oscillators and reconciles the apparent dichotomy between irregular single neuron activity and field potential oscillations. From this perspective, a plethora of studies will be reviewed on the involvement of long-distance neuronal coherence in cognitive functions such as multisensory integration, working memory, and selective attention. Finally, implications of abnormal neural synchronization are discussed as they relate to mental disorders like schizophrenia and autism.

Maximum likelihood decoding of neuronal inputs from an interspike interval distribution

An expression for the probability distribution of the interspike interval of a leaky integrate-and-fire (LIF) model neuron is rigorously derived, based on recent theoretical developments in the theory of stochastic processes. This enables us to find for the first time a way of developing maximum likelihood estimates (MLE) of the input information (e.g., afferent rate and variance) for an LIF neuron from a set of recorded spike trains. Dynamic inputs to pools of LIF neurons both with and without interactions are efficiently and reliably decoded by applying the MLE, even within time windows as short as 25 msec.

Cluster synchronization in community networks with nonidentical nodes

In this paper dynamical networks with community structure and nonidentical nodes and with identical local dynamics for all individual nodes in each community are considered. The cluster synchronization of these networks with or without time delay is studied by using some feedback control schemes. Several sufficient conditions for achieving cluster synchronization are obtained analytically and are further verified numerically by some examples with chaotic or nonchaotic nodes. In addition, an essential relation between synchronization dynamics and local dynamics is found by detailed analysis of dynamical networks without delay through the stage detection of cluster synchronization.

Information encoding in an oscillatory network

Information encoding in a globally coupled network is studied. When the network is in an oscillatory state, the network activities are dominated by the intrinsic oscillatory current and the stimulus is poorly encoded. However, when the amplitude of the input signal is large, the input can still be well read from the population rate and the temporal correlation between spike trains. The underlying reason is that there exists a competition between the intrinsic correlation caused by the oscillatory current and the external correlation caused by the input signal. With small input signal, the rate code performs better than the temporal correlation code. Our results provide insights into the effects of network dynamics on neuronal computations.

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.

Delays and oscillations in networks of spiking neurons: a two-timescale analysis

Oscillations are a ubiquitous feature of many neural systems, spanning many orders of magnitude in frequency. One of the most prominent oscillatory patterns, with possible functional implications, is that occurring in the mammalian thalamocortical system during sleep. This system is characterized by relatively long delays (reaching up to 40 msec) and gives rise to low-frequency oscillatory waves. Motivated by these phenomena, we study networks of excitatory and inhibitory integrate-and-fire neurons within a Fokker-Planck delay partial differential equation formalism and establish explicit conditions for the emergence of oscillatory solutions, and for the amplitude and period of the ensuing oscillations, for relatively large values of the delays. When a two-timescale analysis is employed, the full partial differential equation is replaced in this limit by a discrete time iterative map, leading to a relatively simple dynamic interpretation. This asymptotic result is shown numerically to hold, to a good approximation, over a wide range of parameter values, leading to an accurate characterization of the behavior in terms of the underlying physical parameters. Our results provide a simple mechanistic explanation for one type of slow oscillation based on delayed inhibition, which may play an important role in the slow spindle oscillations occurring during sleep. Moreover, they are consistent with experimental findings related to human motor behavior with visual feedback.

Introduction to focus issue: mixed mode oscillations: experiment, computation, and analysis

Mixed mode oscillations (MMOs) occur when a dynamical system switches between fast and slow motion and small and large amplitude. MMOs appear in a variety of systems in nature, and may be simple or complex. This focus issue presents a series of articles on theoretical, numerical, and experimental aspects of MMOs. The applications cover physical, chemical, and biological systems.