Sparsely synchronized neuronal oscillations

Firing rate models for gamma oscillations in I-I and E-I networks

Firing rate models for describing the mean-field activities of neuronal ensembles can be used effectively to study network function and dynamics, including synchronization and rhythmicity of excitatory-inhibitory populations. However, traditional Wilson-Cowan-like models, even when extended to include an explicit dynamic synaptic activation variable, are found unable to capture some dynamics such as Interneuronal Network Gamma oscillations (ING). Use of an explicit delay is helpful in simulations at the expense of complicating mathematical analysis. We resolve this issue by introducing a dynamic variable, u, that acts as an effective delay in the negative feedback loop between firing rate (r) and synaptic gating of inhibition (s). In effect, u endows synaptic activation with second order dynamics. With linear stability analysis, numerical branch-tracking and simulations, we show that our r-u-s rate model captures some key qualitative features of spiking network models for ING. We also propose an alternative formulation, a v-u-s model, in which mean membrane potential v satisfies an averaged current-balance equation. Furthermore, we extend the framework to E-I networks. With our six-variable v-u-s model, we demonstrate in firing rate models the transition from Pyramidal-Interneuronal Network Gamma (PING) to ING by increasing the external drive to the inhibitory population without adjusting synaptic weights. Having PING and ING available in a single network, without invoking synaptic blockers, is plausible and natural for explaining the emergence and transition of two different types of gamma oscillations.

Intra-ripple frequency accommodation in an inhibitory network model for hippocampal ripple oscillations

Hippocampal ripple oscillations have been implicated in important cognitive functions such as memory consolidation and planning. Multiple computational models have been proposed to explain the emergence of ripple oscillations, relying either on excitation or inhibition as the main pacemaker. Nevertheless, the generating mechanism of ripples remains unclear. An interesting dynamical feature of experimentally measured ripples, which may advance model selection, is intra-ripple frequency accommodation (IFA): a decay of the instantaneous ripple frequency over the course of a ripple event. So far, only a feedback-based inhibition-first model, which relies on delayed inhibitory synaptic coupling, has been shown to reproduce IFA. Here we use an analytical mean-field approach and numerical simulations of a leaky integrate-and-fire spiking network to explain the mechanism of IFA. We develop a drift-based approximation for the oscillation dynamics of the population rate and the mean membrane potential of interneurons under strong excitatory drive and strong inhibitory coupling. For IFA, the speed at which the excitatory drive changes is critical. We demonstrate that IFA arises due to a speed-dependent hysteresis effect in the dynamics of the mean membrane potential, when the interneurons receive transient, sharp wave-associated excitation. We thus predict that the IFA asymmetry vanishes in the limit of slowly changing drive, but is otherwise a robust feature of the feedback-based inhibition-first ripple model.

Exact low-dimensional description for fast neural oscillations with low firing rates

Recently, low-dimensional models of neuronal activity have been exactly derived for large networks of deterministic, quadratic integrate-and-fire (QIF) neurons. Such firing rate models (FRM) describe the emergence of fast collective oscillations (>30 Hz) via the frequency locking of a subset of neurons to the global oscillation frequency. However, the suitability of such models to describe realistic neuronal states is seriously challenged by the fact that during episodes of fast collective oscillations, neuronal discharges are often very irregular and have low firing rates compared to the global oscillation frequency. Here we extend the theory to derive exact FRM for QIF neurons to include noise and show that networks of stochastic neurons displaying irregular discharges at low firing rates during episodes of fast oscillations are governed by exactly the same evolution equations as deterministic networks. Our results reconcile two traditionally confronted views on neuronal synchronization and upgrade the applicability of exact FRM to describe a broad range of biologically realistic neuronal states.

Comparison between an exact and a heuristic neural mass model with second-order synapses

Neural mass models (NMMs) are designed to reproduce the collective dynamics of neuronal populations. A common framework for NMMs assumes heuristically that the output firing rate of a neural population can be described by a static nonlinear transfer function (NMM1). However, a recent exact mean-field theory for quadratic integrate-and-fire (QIF) neurons challenges this view by showing that the mean firing rate is not a static function of the neuronal state but follows two coupled nonlinear differential equations (NMM2). Here we analyze and compare these two descriptions in the presence of second-order synaptic dynamics. First, we derive the mathematical equivalence between the two models in the infinitely slow synapse limit, i.e., we show that NMM1 is an approximation of NMM2 in this regime. Next, we evaluate the applicability of this limit in the context of realistic physiological parameter values by analyzing the dynamics of models with inhibitory or excitatory synapses. We show that NMM1 fails to reproduce important dynamical features of the exact model, such as the self-sustained oscillations of an inhibitory interneuron QIF network. Furthermore, in the exact model but not in the limit one, stimulation of a pyramidal cell population induces resonant oscillatory activity whose peak frequency and amplitude increase with the self-coupling gain and the external excitatory input. This may play a role in the enhanced response of densely connected networks to weak uniform inputs, such as the electric fields produced by noninvasive brain stimulation.

Scale free avalanches in excitatory-inhibitory populations of spiking neurons with conductance based synaptic currents

We investigate spontaneous critical dynamics of excitatory and inhibitory (EI) sparsely connected populations of spiking leaky integrate-and-fire neurons with conductance-based synapses. We use a bottom-up approach to derive a single neuron gain function and a linear Poisson neuron approximation which we use to study mean-field dynamics of the EI population and its bifurcations. In the low firing rate regime, the quiescent state loses stability due to saddle-node or Hopf bifurcations. In particular, at the Bogdanov-Takens (BT) bifurcation point which is the intersection of the Hopf bifurcation and the saddle-node bifurcation lines of the 2D dynamical system, the network shows avalanche dynamics with power-law avalanche size and duration distributions. This matches the characteristics of low firing spontaneous activity in the cortex. By linearizing gain functions and excitatory and inhibitory nullclines, we can approximate the location of the BT bifurcation point. This point in the control parameter phase space corresponds to the internal balance of excitation and inhibition and a slight excess of external excitatory input to the excitatory population. Due to the tight balance of average excitation and inhibition currents, the firing of the individual cells is fluctuation-driven. Around the BT point, the spiking of neurons is a Poisson process and the population average membrane potential of neurons is approximately at the middle of the operating interval [Formula: see text]. Moreover, the EI network is close to both oscillatory and active-inactive phase transition regimes.

Disynaptic effect of hilar cells on pattern separation in a spiking neural network of hippocampal dentate gyrus

We study the disynaptic effect of the hilar cells on pattern separation in a spiking neural network of the hippocampal dentate gyrus (DG). The principal granule cells (GCs) in the DG perform pattern separation, transforming similar input patterns into less-similar output patterns. In our DG network, the hilus consists of excitatory mossy cells (MCs) and inhibitory HIPP (hilar perforant path-associated) cells. Here, we consider the disynaptic effects of the MCs and the HIPP cells on the GCs, mediated by the inhibitory basket cells (BCs) in the granular layer; MC → BC → GC and HIPP → BC → GC. The MCs provide disynaptic inhibitory input (mediated by the intermediate BCs) to the GCs, which decreases the firing activity of the GCs. On the other hand, the HIPP cells disinhibit the intermediate BCs, which leads to increasing the firing activity of the GCs. In this way, the disynaptic effects of the MCs and the HIPP cells are opposite. We investigate change in the pattern separation efficacy by varying the synaptic strength K ( BC , X ) [from the pre-synaptic X (= MC or HIPP) to the post-synaptic BC]. Thus, sparsity for the firing activity of the GCs is found to improve the efficacy of pattern separation, and hence the disynaptic effects of the MCs and the HIPP cells on the pattern separation become opposite ones. In the combined case when simultaneously changing both K ( BC , MC ) and K ( BC , HIPP ) , as a result of balance between the two competing disynaptic effects of the MCs and the HIPP cells, the efficacy of pattern separation is found to become the highest at their original default values where the activation degree of the GCs is the lowest. We also note that, while the GCs perform pattern separation, sparsely synchronized rhythm is found to appear in the population of the GCs. Hence, we examine quantitative association between population and individual firing behaviors in the sparsely synchronized rhythm and pattern separation. They are found to be strongly correlated. Consequently, the better the population and individual firing behaviors in the sparsely synchronized rhythm are, the more pattern separation efficacy becomes enhanced.

Population and individual firing behaviors in sparsely synchronized rhythms in the hippocampal dentate gyrus

We investigate population and individual firing behaviors in sparsely synchronized rhythms (SSRs) in a spiking neural network of the hippocampal dentate gyrus (DG). The main encoding granule cells (GCs) are grouped into lamellar clusters. In each GC cluster, there is one inhibitory (I) basket cell (BC) along with excitatory (E) GCs, and they form the E-I loop. Winner-take-all competition, leading to sparse activation of the GCs, occurs in each GC cluster. Such sparsity has been thought to enhance pattern separation performed in the DG. During the winner-take-all competition, SSRs are found to appear in each population of the GCs and the BCs through interaction of excitation of the GCs with inhibition of the BCs. Sparsely synchronized spiking stripes appear successively with the population frequencyf p ( = 13.1 Hz) in the raster plots of spikes. We also note that excitatory hilar mossy cells (MCs) control the firing activity of the GC-BC loop by providing excitation to both the GCs and the BCs. SSR also appears in the population of MCs via interaction with the GCs (i.e., GC-MC loop). Population behaviors in the SSRs are quantitatively characterized in terms of the synchronization measures. In addition, we investigate individual firing activity of GCs, BCs, and MCs in the SSRs. Individual GCs exhibit random spike skipping, leading to a multi-peaked inter-spike-interval histogram, which is well characterized in terms of the random phase-locking degree. In this case, population-averaged mean-firing-rate (MFR) < f i ( GC ) > is less than the population frequency f p . On the other hand, both BCs and MCs show "intrastripe" burstings within stripes, together with random spike skipping. Thus, the population-averaged MFR ⟨ f i ( X ) ⟩ ( X = MC and BC) is larger than f p , in contrast to the case of the GCs. MC loss may occur during epileptogenesis. With decreasing the fraction of the MCs, changes in the population and individual firings in the SSRs are also studied. Finally, quantitative association between the population/individual firing behaviors in the SSRs and the winner-take-all competition is discussed.

Mean-field equations for neural populations with q-Gaussian heterogeneities

Describing the collective dynamics of large neural populations using low-dimensional models for averaged variables has long been an attractive task in theoretical neuroscience. Recently developed reduction methods make it possible to derive such models directly from the microscopic dynamics of individual neurons. To simplify the reduction, the Cauchy distribution is usually assumed for heterogeneous network parameters. Here we extend the reduction method for a wider class of heterogeneities defined by the q-Gaussian distribution. The shape of this distribution depends on the Tsallis index q and gradually changes from the Cauchy distribution to the normal Gaussian distribution as this index changes. We derive the mean-field equations for an inhibitory network of quadratic integrate-and-fire neurons with a q-Gaussian-distributed excitability parameter. It is shown that the dynamic modes of the network significantly depend on the form of the distribution determined by the Tsallis index. The results obtained from the mean-field equations are confirmed by numerical simulation of the microscopic model.

Criticality enhances the multilevel reliability of stimulus responses in cortical neural networks

Cortical neural networks exhibit high internal variability in spontaneous dynamic activities and they can robustly and reliably respond to external stimuli with multilevel features–from microscopic irregular spiking of neurons to macroscopic oscillatory local field potential. A comprehensive study integrating these multilevel features in spontaneous and stimulus–evoked dynamics with seemingly distinct mechanisms is still lacking. Here, we study the stimulus–response dynamics of biologically plausible excitation–inhibition (E–I) balanced networks. We confirm that networks around critical synchronous transition states can maintain strong internal variability but are sensitive to external stimuli. In this dynamical region, applying a stimulus to the network can reduce the trial-to-trial variability and shift the network oscillatory frequency while preserving the dynamical criticality. These multilevel features widely observed in different experiments cannot simultaneously occur in non-critical dynamical states. Furthermore, the dynamical mechanisms underlying these multilevel features are revealed using a semi-analytical mean-field theory that derives the macroscopic network field equations from the microscopic neuronal networks, enabling the analysis by nonlinear dynamics theory and linear noise approximation. The generic dynamical principle revealed here contributes to a more integrative understanding of neural systems and brain functions and incorporates multimodal and multilevel experimental observations. The E–I balanced neural network in combination with the effective mean-field theory can serve as a mechanistic modeling framework to study the multilevel neural dynamics underlying neural information and cognitive processes.

The complexity and variability of brain dynamical activity range from neuronal spiking and neural avalanches to oscillatory local field potentials of local neural circuits in both spontaneous and stimulus-evoked states. Such multilevel variable brain dynamics are functionally and behaviorally relevant and are principal components of the underlying circuit organization. To more comprehensively clarify their neural mechanisms, we use a bottom-up approach to study the stimulus–response dynamics of neural circuits. Our model assumes the following key biologically plausible components: excitation–inhibition (E–I) neuronal interaction and chemical synaptic coupling. We show that the circuits with E–I balance have a special dynamic sub-region, the critical region. Circuits around this region could account for the emergence of multilevel brain response patterns, both ongoing and stimulus-induced, observed in different experiments, including the reduction of trial-to-trial variability, effective modulation of gamma frequency, and preservation of criticality in the presence of a stimulus. We further analyze the corresponding nonlinear dynamical principles using a novel and highly generalizable semi-analytical mean-field theory. Our computational and theoretical studies explain the cross-level brain dynamical organization of spontaneous and evoked states in a more integrative manner.

Dynamical origin for winner-take-all competition in a biological network of the hippocampal dentate gyrus

We consider a biological network of the hippocampal dentate gyrus (DG). Computational models suggest that the DG would be a preprocessor for pattern separation (i.e., a process transforming a set of similar input patterns into distinct nonoverlapping output patterns) which could facilitate pattern storage and retrieval in the CA3 area of the hippocampus. The main encoding cells in the DG are the granule cells (GCs) which receive the input from the entorhinal cortex (EC) and send their output to the CA3. We note that the activation degree of GCs is very low (∼5%). This sparsity has been thought to enhance the pattern separation. We investigate the dynamical origin for winner-take-all (WTA) competition which leads to sparse activation of the GCs. The whole GCs are grouped into lamellar clusters. In each cluster, there is one inhibitory (I) basket cell (BC) along with excitatory (E) GCs. There are three kinds of external inputs into the GCs: the direct excitatory EC input; the indirect feedforward inhibitory EC input, mediated by the HIPP (hilar perforant path-associated) cells; and the excitatory input from the hilar mossy cells (MCs). The firing activities of the GCs are determined via competition between the external E and I inputs. The E-I conductance ratio R_{E-I}^{(con)}^{*} (given by the time average of the ratio of the external E to I conductances) may represent well the degree of such external E-I input competition. It is thus found that GCs become active when their R_{E-I}^{(con)}^{*} is larger than a threshold R_{th}^{*}, and then the mean firing rates of the active GCs are strongly correlated with R_{E-I}^{(con)}^{*}. In each cluster, the feedback inhibition from the BC may select the winner GCs. GCs with larger R_{E-I}^{(con)}^{*} than the threshold R_{th}^{*} survive, and they become winners; all the other GCs with smaller R_{E-I}^{(con)}^{*} become silent. In this way, WTA competition occurs via competition between the firing activity of the GCs and the feedback inhibition from the BC in each cluster. Finally, we also study the effects of MC death and adult-born immature GCs on the WTA competition.

Influence of various temporal recoding on pavlovian eyeblink conditioning in the cerebellum

We consider the Pavlovian eyeblink conditioning (EBC) via repeated presentation of paired conditioned stimulus (tone) and unconditioned stimulus (US; airpuff). In an effective cerebellar ring network, we change the connection probability p c from Golgi to granule (GR) cells, and make a dynamical classification of various firing patterns of the GR cells. Individual GR cells are thus found to show various well- and ill-matched firing patterns relative to the US timing signal. Then, these variously-recoded signals are fed into the Purkinje cells (PCs) through the parallel-fibers (PFs). Based on such unique dynamical classification of various firing patterns, we make intensive investigations on the influence of various temporal recoding (i.e., firing patterns) of the GR cells on the synaptic plasticity of the PF-PC synapses and the subsequent learning process for the EBC. We first note that the variously-recoded PF signals are effectively depressed by the (error-teaching) instructor climbing-fiber (CF) signals from the inferior olive neuron. In the case of well-matched PF signals, they are strongly depressed through strong long-term depression (LTD) by the instructor CF signals due to good association between the in-phase PF and the instructor CF signals. On the other hand, practically no LTD occurs for the ill-matched PF signals because most of them have no association with the instructor CF signals. This kind of "effective" depression at the PF-PC synapses coordinates firings of PCs effectively, which then makes effective inhibitory coordination on the cerebellar nucleus neuron [which elicits conditioned response (CR; eyeblink)]. When the learning trial passes a threshold, acquisition of CR begins. In this case, the timing degree T d of CR becomes good due to presence of the ill-matched firing group which plays a role of protection barrier for the timing. With further increase in the number of trials, strength S of CR (corresponding to the amplitude of eyelid closure) increases due to strong LTD in the well-matched firing group, while its timing degree T d decreases. In this way, the well- and the ill-matched firing groups play their own roles for the strength and the timing of CR, respectively. Thus, with increasing the number of learning trials, the (overall) learning efficiency degree L e (taking into consideration both timing and strength of CR) for the CR is increased, and eventually it becomes saturated. Finally, we also discuss dependence of the variety degree for firing patterns of the GR cells and the saturated learning efficiency degree L e of the CR on p c and their relations.

Robust neuromorphic coupled oscillators for adaptive pacemakers

Neural coupled oscillators are a useful building block in numerous models and applications. They were analyzed extensively in theoretical studies and more recently in biologically realistic simulations of spiking neural networks. The advent of mixed-signal analog/digital neuromorphic electronic circuits provides new means for implementing neural coupled oscillators on compact, low-power, spiking neural network hardware platforms. However, their implementation on this noisy, low-precision and inhomogeneous computing substrate raises new challenges with regards to stability and controllability. In this work, we present a robust, spiking neural network model of neural coupled oscillators and validate it with an implementation on a mixed-signal neuromorphic processor. We demonstrate its robustness showing how to reliably control and modulate the oscillator's frequency and phase shift, despite the variability of the silicon synapse and neuron properties. We show how this ultra-low power neural processing system can be used to build an adaptive cardiac pacemaker modulating the heart rate with respect to the respiration phases and compare it with surface ECG and respiratory signal recordings from dogs at rest. The implementation of our model in neuromorphic electronic hardware shows its robustness on a highly variable substrate and extends the toolbox for applications requiring rhythmic outputs such as pacemakers.

Transmission delays and frequency detuning can regulate information flow between brain regions

Brain networks exhibit very variable and dynamical functional connectivity and flexible configurations of information exchange despite their overall fixed structure. Brain oscillations are hypothesized to underlie time-dependent functional connectivity by periodically changing the excitability of neural populations. In this paper, we investigate the role of the connection delay and the detuning between the natural frequencies of neural populations in the transmission of signals. Based on numerical simulations and analytical arguments, we show that the amount of information transfer between two oscillating neural populations could be determined by their connection delay and the mismatch in their oscillation frequencies. Our results highlight the role of the collective phase response curve of the oscillating neural populations for the efficacy of signal transmission and the quality of the information transfer in brain networks.

The Refractory Period Matters: Unifying Mechanisms of Macroscopic Brain Waves

The relationship between complex brain oscillations and the dynamics of individual neurons is poorly understood. Here we utilize maximum caliber, a dynamical inference principle, to build a minimal yet general model of the collective (mean field) dynamics of large populations of neurons. In agreement with previous experimental observations, we describe a simple, testable mechanism, involving only a single type of neuron, by which many of these complex oscillatory patterns may emerge. Our model predicts that the refractory period of neurons, which has often been neglected, is essential for these behaviors.

Modulation of Coordinated Activity across Cortical Layers by Plasticity of Inhibitory Synapses

In the neocortex, synaptic inhibition shapes all forms of spontaneous and sensory evoked activity. Importantly, inhibitory transmission is highly plastic, but the functional role of inhibitory synaptic plasticity is unknown. In the mouse barrel cortex, activation of layer (L) 2/3 pyramidal neurons (PNs) elicits strong feedforward inhibition (FFI) onto L5 PNs. We find that FFI involving parvalbumin (PV)-expressing cells is strongly potentiated by postsynaptic PN burst firing. FFI plasticity modifies the PN excitation-to-inhibition (E/I) ratio, strongly modulates PN gain, and alters information transfer across cortical layers. Moreover, our LTPi-inducing protocol modifies firing of L5 PNs and alters the temporal association of PN spikes to γ-oscillations both in vitro and in vivo. All of these effects are captured by unbalancing the E/I ratio in a feedforward inhibition circuit model. Altogether, our results indicate that activity-dependent modulation of perisomatic inhibitory strength effectively influences the participation of single principal cortical neurons to cognition-relevant network activity.

•

Feedforward inhibition (FFI) of layer 5 pyramidal neurons (PNs) is highly plastic

•

Long-term potentiation of FFI modulates spiking activity of layer 5 PNs

•

LTPi affects information transfer across cortical layers

•

The strength of LTPi determines the phase locking of PN firing to γ-oscillations

Effect of diverse recoding of granule cells on optokinetic response in a cerebellar ring network with synaptic plasticity

We consider a cerebellar ring network for the optokinetic response (OKR), and investigate the effect of diverse recoding of granule (GR) cells on OKR by varying the connection probability p_c from Golgi to GR cells. For an optimal value of p_c^∗(=0.06), individual GR cells exhibit diverse spiking patterns which are in-phase, anti-phase, or complex out-of-phase with respect to their population-averaged firing activity. Then, these diversely-recoded signals via parallel fibers (PFs) from GR cells are effectively depressed by the error-teaching signals via climbing fibers from the inferior olive which are also in-phase ones. Synaptic weights at in-phase PF-Purkinje cell (PC) synapses of active GR cells are strongly depressed via strong long-term depression (LTD), while those at anti-phase and complex out-of-phase PF-PC synapses are weakly depressed through weak LTD. This kind of "effective" depression (i.e., strong/weak LTD) at the PF-PC synapses causes a big modulation in firings of PCs, which then exert effective inhibitory coordination on the vestibular nucleus (VN) neuron (which evokes OKR). For the firing of the VN neuron, the learning gain degree L_g, corresponding to the modulation gain ratio, increases with increasing the learning cycle, and it saturates at about the 300th cycle. By varying p_c from p_c^∗, we find that a plot of saturated learning gain degree L_g^∗ versus p_c forms a bell-shaped curve with a peak at p_c^∗ (where the diversity degree in spiking patterns of GR cells is also maximum). Consequently, the more diverse in recoding of GR cells, the more effective in motor learning for the OKR adaptation.

Hopf Bifurcation in Mean Field Explains Critical Avalanches in Excitation-Inhibition Balanced Neuronal Networks: A Mechanism for Multiscale Variability

Cortical neural circuits display highly irregular spiking in individual neurons but variably sized collective firing, oscillations and critical avalanches at the population level, all of which have functional importance for information processing. Theoretically, the balance of excitation and inhibition inputs is thought to account for spiking irregularity and critical avalanches may originate from an underlying phase transition. However, the theoretical reconciliation of these multilevel dynamic aspects in neural circuits remains an open question. Herein, we study excitation-inhibition (E-I) balanced neuronal network with biologically realistic synaptic kinetics. It can maintain irregular spiking dynamics with different levels of synchrony and critical avalanches emerge near the synchronous transition point. We propose a novel semi-analytical mean-field theory to derive the field equations governing the network macroscopic dynamics. It reveals that the E-I balanced state of the network manifesting irregular individual spiking is characterized by a macroscopic stable state, which can be either a fixed point or a periodic motion and the transition is predicted by a Hopf bifurcation in the macroscopic field. Furthermore, by analyzing public data, we find the coexistence of irregular spiking and critical avalanches in the spontaneous spiking activities of mouse cortical slice in vitro, indicating the universality of the observed phenomena. Our theory unveils the mechanism that permits complex neural activities in different spatiotemporal scales to coexist and elucidates a possible origin of the criticality of neural systems. It also provides a novel tool for analyzing the macroscopic dynamics of E-I balanced networks and its relationship to the microscopic counterparts, which can be useful for large-scale modeling and computation of cortical dynamics.

Synchronization, Stochasticity, and Phase Waves in Neuronal Networks With Spatially-Structured Connectivity

Oscillations in the beta/low gamma range (10–45 Hz) are recorded in diverse neural structures. They have successfully been modeled as sparsely synchronized oscillations arising from reciprocal interactions between randomly connected excitatory (E) pyramidal cells and local interneurons (I). The synchronization of spatially distant oscillatory spiking E–I modules has been well-studied in the rate model framework but less so for modules of spiking neurons. Here, we first show that previously proposed modifications of rate models provide a quantitative description of spiking E–I modules of Exponential Integrate-and-Fire (EIF) neurons. This allows us to analyze the dynamical regimes of sparsely synchronized oscillatory E–I modules connected by long-range excitatory interactions, for two modules, as well as for a chain of such modules. For modules with a large number of neurons (> 10⁵), we obtain results similar to previously obtained ones based on the classic deterministic Wilson-Cowan rate model, with the added bonus that the results quantitatively describe simulations of spiking EIF neurons. However, for modules with a moderate (~ 10⁴) number of neurons, stochastic variations in the spike emission of neurons are important and need to be taken into account. On the one hand, they modify the oscillations in a way that tends to promote synchronization between different modules. On the other hand, independent fluctuations on different modules tend to disrupt synchronization. The correlations between distant oscillatory modules can be described by stochastic equations for the oscillator phases that have been intensely studied in other contexts. On shorter distances, we develop a description that also takes into account amplitude modes and that quantitatively accounts for our simulation data. Stochastic dephasing of neighboring modules produces transient phase gradients and the transient appearance of phase waves. We propose that these stochastically-induced phase waves provide an explanative framework for the observations of traveling waves in the cortex during beta oscillations.

What is the true discharge rate and pattern of the striatal projection neurons in Parkinson's disease and Dystonia?

Dopamine and striatal dysfunctions play a key role in the pathophysiology of Parkinson's disease (PD) and Dystonia, but our understanding of the changes in the discharge rate and pattern of striatal projection neurons (SPNs) remains limited. Here, we recorded and examined multi-unit signals from the striatum of PD and dystonic patients undergoing deep brain stimulation surgeries. Contrary to earlier human findings, we found no drastic changes in the spontaneous discharge of the well-isolated and stationary SPNs of the PD patients compared to the dystonic patients or to the normal levels of striatal activity reported in healthy animals. Moreover, cluster analysis using SPN discharge properties did not characterize two well-separated SPN subpopulations, indicating no SPN subpopulation-specific (D1 or D2 SPNs) discharge alterations in the pathological state. Our results imply that small to moderate changes in spontaneous SPN discharge related to PD and Dystonia are likely amplified by basal ganglia downstream structures.

mPFC spindle cycles organize sparse thalamic activation and recently active CA1 cells during non-REM sleep

Sleep oscillations in the neocortex and hippocampus are critical for the integration of new memories into stable generalized representations in neocortex. However, the role of the thalamus in this process is poorly understood. To determine the thalamic contribution to non-REM oscillations (sharp-wave ripples, SWRs; slow/delta; spindles), we recorded units and local field potentials (LFPs) simultaneously in the limbic thalamus, mPFC, and CA1 in rats. We report that the cycles of neocortical spindles provide a key temporal window that coordinates CA1 SWRs with sparse but consistent activation of thalamic units. Thalamic units were phase-locked to delta and spindles in mPFC, and fired at consistent lags with other thalamic units within spindles, while CA1 units that were active during spatial exploration were engaged in SWR-coupled spindles after behavior. The sparse thalamic firing could promote an incremental integration of recently acquired memory traces into neocortical schemas through the interleaved activation of thalamocortical cells.

Synaptic Diversity Suppresses Complex Collective Behavior in Networks of Theta Neurons

Comprehending how the brain functions requires an understanding of the dynamics of neuronal assemblies. Previous work used a mean-field reduction method to determine the collective dynamics of a large heterogeneous network of uniformly and globally coupled theta neurons, which are a canonical formulation of Type I neurons. However, in modeling neuronal networks, it is unreasonable to assume that the coupling strength between every pair of neurons is identical. The goal in the present work is to analytically examine the collective macroscopic behavior of a network of theta neurons that is more realistic in that it includes heterogeneity in the coupling strength as well as in neuronal excitability. We consider the occurrence of dynamical structures that give rise to complicated dynamics via bifurcations of macroscopic collective quantities, concentrating on two biophysically relevant cases: (1) predominantly excitable neurons with mostly excitatory connections, and (2) predominantly spiking neurons with inhibitory connections. We find that increasing the synaptic diversity moves these dynamical structures to distant extremes of parameter space, leaving simple collective equilibrium states in the physiologically relevant region. We also study the node vs. focus nature of stable macroscopic equilibrium solutions and discuss our results in the context of recent literature.

Mechanisms of neural organization and rhythmogenesis during hippocampal and cortical ripples

Neural activity during ripples has attracted great theoretical and experimental attention over the last three decades. Perhaps one reason for such interest is that ripples occur during quiet waking moments and during sleep, times when we reflect and dream about what has just occurred and what we expect to happen next. The hope is that understanding such 'offline' activity may yield insights into reflection, planning, and the purposes of sleep. This review focuses on the mechanisms by which neurons organize during these high-frequency events. In studying ripples, broader principles have emerged that relate intrinsic neural properties, network topology and synaptic plasticity in controlling neural activity. Ripples, therefore, serve as an excellent model for studying how properties of a neural network relate to neural dynamics. This article is part of the Theo Murphy meeting issue 'Memory reactivation: replaying events past, present and future'.

Effect of interpopulation spike-timing-dependent plasticity on synchronized rhythms in neuronal networks with inhibitory and excitatory populations

We consider a two-population network consisting of both inhibitory (I) interneurons and excitatory (E) pyramidal cells. This I-E neuronal network has adaptive dynamic I to E and E to I interpopulation synaptic strengths, governed by interpopulation spike-timing-dependent plasticity (STDP). In previous works without STDPs, fast sparsely synchronized rhythms, related to diverse cognitive functions, were found to appear in a range of noise intensity D for static synaptic strengths. Here, by varying D, we investigate the effect of interpopulation STDPs on fast sparsely synchronized rhythms that emerge in both the I- and the E-populations. Depending on values of D, long-term potentiation (LTP) and long-term depression (LTD) for population-averaged values of saturated interpopulation synaptic strengths are found to occur. Then, the degree of fast sparse synchronization varies due to effects of LTP and LTD. In a broad region of intermediate D, the degree of good synchronization (with higher synchronization degree) becomes decreased, while in a region of large D, the degree of bad synchronization (with lower synchronization degree) gets increased. Consequently, in each I- or E-population, the synchronization degree becomes nearly the same in a wide range of D (including both the intermediate and the large D regions). This kind of "equalization effect" is found to occur via cooperative interplay between the average occupation and pacing degrees of spikes (i.e., the average fraction of firing neurons and the average degree of phase coherence between spikes in each synchronized stripe of spikes in the raster plot of spikes) in fast sparsely synchronized rhythms. Finally, emergences of LTP and LTD of interpopulation synaptic strengths (leading to occurrence of equalization effect) are intensively investigated via a microscopic method based on the distributions of time delays between the pre- and the post-synaptic spike times.

Jensen's force and the statistical mechanics of cortical asynchronous states

Cortical networks are shaped by the combined action of excitatory and inhibitory interactions. Among other important functions, inhibition solves the problem of the all-or-none type of response that comes about in purely excitatory networks, allowing the network to operate in regimes of moderate or low activity, between quiescent and saturated regimes. Here, we elucidate a noise-induced effect that we call “Jensen’s force” –stemming from the combined effect of excitation/inhibition balance and network sparsity– which is responsible for generating a phase of self-sustained low activity in excitation-inhibition networks. The uncovered phase reproduces the main empirically-observed features of cortical networks in the so-called asynchronous state, characterized by low, un-correlated and highly-irregular activity. The parsimonious model analyzed here allows us to resolve a number of long-standing issues, such as proving that activity can be self-sustained even in the complete absence of external stimuli or driving. The simplicity of our approach allows for a deep understanding of asynchronous states and of the phase transitions to other standard phases it exhibits, opening the door to reconcile, asynchronous-state and critical-state hypotheses, putting them within a unified framework. We argue that Jensen’s forces are measurable experimentally and might be relevant in contexts beyond neuroscience.

Firing rate models for gamma oscillations

Gamma oscillations are readily observed in a variety of brain regions during both waking and sleeping states. Computational models of gamma oscillations typically involve simulations of large networks of synaptically coupled spiking units. These networks can exhibit strongly synchronized gamma behavior, whereby neurons fire in near synchrony on every cycle, or weakly modulated gamma behavior, corresponding to stochastic, sparse firing of the individual units on each cycle of the population gamma rhythm. These spiking models offer valuable biophysical descriptions of gamma oscillations; however, because they involve many individual neuronal units they are limited in their ability to communicate general network-level dynamics. Here we demonstrate that few-variable firing rate models with established synaptic timescales can account for both strongly synchronized and weakly modulated gamma oscillations. These models go beyond the classical formulations of rate models by including at least two dynamic variables per population: firing rate and synaptic activation. The models' flexibility to capture the broad range of gamma behavior depends directly on the timescales that represent recruitment of the excitatory and inhibitory firing rates. In particular, we find that weakly modulated gamma oscillations occur robustly when the recruitment timescale of inhibition is faster than that of excitation. We present our findings by using an extended Wilson-Cowan model and a rate model derived from a network of quadratic integrate-and-fire neurons. These biophysical rate models capture the range of weakly modulated and coherent gamma oscillations observed in spiking network models, while additionally allowing for greater tractability and systems analysis. NEW & NOTEWORTHY Here we develop simple and tractable models of gamma oscillations, a dynamic feature observed throughout much of the brain with significant correlates to behavior and cognitive performance in a variety of experimental contexts. Our models depend on only a few dynamic variables per population, but despite this they qualitatively capture features observed in previous biophysical models of gamma oscillations that involve many individual spiking units.

Effect of inhibitory spike-timing-dependent plasticity on fast sparsely synchronized rhythms in a small-world neuronal network

We consider the Watts-Strogatz small-world network (SWN) consisting of inhibitory fast spiking Izhikevich interneurons. This inhibitory neuronal population has adaptive dynamic synaptic strengths governed by the inhibitory spike-timing-dependent plasticity (iSTDP). In previous works without iSTDP, fast sparsely synchronized rhythms, associated with diverse cognitive functions, were found to appear in a range of large noise intensities for fixed strong synaptic inhibition strengths. Here, we investigate the effect of iSTDP on fast sparse synchronization (FSS) by varying the noise intensity D. We employ an asymmetric anti-Hebbian time window for the iSTDP update rule [which is in contrast to the Hebbian time window for the excitatory STDP (eSTDP)]. Depending on values of D, population-averaged values of saturated synaptic inhibition strengths are potentiated [long-term potentiation (LTP)] or depressed [long-term depression (LTD)] in comparison with the initial mean value, and dispersions from the mean values of LTP/LTD are much increased when compared with the initial dispersion, independently of D. In most cases of LTD where the effect of mean LTD is dominant in comparison with the effect of dispersion, good synchronization (with higher spiking measure) is found to get better via LTD, while bad synchronization (with lower spiking measure) is found to get worse via LTP. This kind of Matthew effect in inhibitory synaptic plasticity is in contrast to that in excitatory synaptic plasticity where good (bad) synchronization gets better (worse) via LTP (LTD). Emergences of LTD and LTP of synaptic inhibition strengths are intensively investigated via a microscopic method based on the distributions of time delays between the pre- and the post-synaptic spike times. Furthermore, we also investigate the effects of network architecture on FSS by changing the rewiring probability p of the SWN in the presence of iSTDP.

Hippocampal Ripple Oscillations and Inhibition-First Network Models: Frequency Dynamics and Response to GABA Modulators

Hippocampal ripples are involved in memory consolidation, but the mechanisms underlying their generation remain unclear. Models relying on interneuron networks in the CA1 region disagree on the predominant source of excitation to interneurons: either "direct," via the Schaffer collaterals that provide feedforward input from CA3 to CA1, or "indirect," via the local pyramidal cells in CA1, which are embedded in a recurrent excitatory-inhibitory network. Here, we used physiologically constrained computational models of basket-cell networks to investigate how they respond to different conditions of transient, noisy excitation. We found that direct excitation of interneurons could evoke ripples (140-220 Hz) that exhibited intraripple frequency accommodation and were frequency-insensitive to GABA modulators, as previously shown in in vitro experiments. In addition, the indirect excitation of the basket-cell network enabled the expression of intraripple frequency accommodation in the fast-gamma range (90-140 Hz), as in vivo In our model, intraripple frequency accommodation results from a hysteresis phenomenon in which the frequency responds differentially to the rising and descending phases of the transient excitation. Such a phenomenon predicts a maximum oscillation frequency occurring several milliseconds before the peak of excitation. We confirmed this prediction for ripples in brain slices from male mice. These results suggest that ripple and fast-gamma episodes are produced by the same interneuron network that is recruited via different excitatory input pathways, which could be supported by the previously reported intralaminar connectivity bias between basket cells and functionally distinct subpopulations of pyramidal cells in CA1. Together, our findings unify competing inhibition-first models of rhythm generation in the hippocampus.SIGNIFICANCE STATEMENT The hippocampus is a part of the brain of humans and other mammals that is critical for the acquisition and consolidation of memories. During deep sleep and resting periods, the hippocampus generates high-frequency (∼200 Hz) oscillations called ripples, which are important for memory consolidation. The mechanisms underlying ripple generation are not well understood. A prominent hypothesis holds that the ripples are generated by local recurrent networks of inhibitory neurons. Using computational models and experiments in brain slices from rodents, we show that the dynamics of interneuron networks clarify several previously unexplained characteristics of ripple oscillations, which advances our understanding of hippocampus-dependent memory consolidation.

Landau-Ginzburg theory of cortex dynamics: Scale-free avalanches emerge at the edge of synchronization

Understanding the origin, nature, and functional significance of complex patterns of neural activity, as recorded by diverse electrophysiological and neuroimaging techniques, is a central challenge in neuroscience. Such patterns include collective oscillations emerging out of neural synchronization as well as highly heterogeneous outbursts of activity interspersed by periods of quiescence, called "neuronal avalanches." Much debate has been generated about the possible scale invariance or criticality of such avalanches and its relevance for brain function. Aimed at shedding light onto this, here we analyze the large-scale collective properties of the cortex by using a mesoscopic approach following the principle of parsimony of Landau-Ginzburg. Our model is similar to that of Wilson-Cowan for neural dynamics but crucially, includes stochasticity and space; synaptic plasticity and inhibition are considered as possible regulatory mechanisms. Detailed analyses uncover a phase diagram including down-state, synchronous, asynchronous, and up-state phases and reveal that empirical findings for neuronal avalanches are consistently reproduced by tuning our model to the edge of synchronization. This reveals that the putative criticality of cortical dynamics does not correspond to a quiescent-to-active phase transition as usually assumed in theoretical approaches but to a synchronization phase transition, at which incipient oscillations and scale-free avalanches coexist. Furthermore, our model also accounts for up and down states as they occur (e.g., during deep sleep). This approach constitutes a framework to rationalize the possible collective phases and phase transitions of cortical networks in simple terms, thus helping to shed light on basic aspects of brain functioning from a very broad perspective.

Firing rate equations require a spike synchrony mechanism to correctly describe fast oscillations in inhibitory networks

Recurrently coupled networks of inhibitory neurons robustly generate oscillations in the gamma band. Nonetheless, the corresponding Wilson-Cowan type firing rate equation for such an inhibitory population does not generate such oscillations without an explicit time delay. We show that this discrepancy is due to a voltage-dependent spike-synchronization mechanism inherent in networks of spiking neurons which is not captured by standard firing rate equations. Here we investigate an exact low-dimensional description for a network of heterogeneous canonical Class 1 inhibitory neurons which includes the sub-threshold dynamics crucial for generating synchronous states. In the limit of slow synaptic kinetics the spike-synchrony mechanism is suppressed and the standard Wilson-Cowan equations are formally recovered as long as external inputs are also slow. However, even in this limit synchronous spiking can be elicited by inputs which fluctuate on a time-scale of the membrane time-constant of the neurons. Our meanfield equations therefore represent an extension of the standard Wilson-Cowan equations in which spike synchrony is also correctly described.

Population models describing the average activity of large neuronal ensembles are a powerful mathematical tool to investigate the principles underlying cooperative function of large neuronal systems. However, these models do not properly describe the phenomenon of spike synchrony in networks of neurons. In particular, they fail to capture the onset of synchronous oscillations in networks of inhibitory neurons. We show that this limitation is due to a voltage-dependent synchronization mechanism which is naturally present in spiking neuron models but not captured by traditional firing rate equations. Here we investigate a novel set of macroscopic equations which incorporate both firing rate and membrane potential dynamics, and that correctly generate fast inhibition-based synchronous oscillations. In the limit of slow-synaptic processing oscillations are suppressed, and the model reduces to an equation formally equivalent to the Wilson-Cowan model.

Noisy Juxtacellular Stimulation In Vivo Leads to Reliable Spiking and Reveals High-Frequency Coding in Single Neurons

Single cells in the motor and somatosensory cortex of rats were stimulated in vivo with broadband fluctuating currents applied juxtacellularly. Unlike the DC current steps used previously, fluctuating stimulation currents reliably evoked spike trains with precise timing of individual spikes. Fluctuating currents resulted in strong cellular responses at stimulation frequencies beyond the inverse membrane time constant and the mean firing rate of the neuron. Neuronal firing was associated with high rates of information transmission, even for the high-frequency components of the stimulus. Such response characteristics were also revealed in additional experiments with sinusoidal juxtacellular stimulation. For selected cells, we could reproduce these statistics with compartmental models of varying complexity. We also developed a method to generate Gaussian stimuli that evoke spike trains with prescribed spike times (under the constraint of a certain rate and coefficient of variation) and exemplify its ability to achieve precise and reliable spiking in cortical neurons in vivo Our results demonstrate a novel method for precise control of spike timing by juxtacellular stimulation, confirm and extend earlier conclusions from ex vivo work about the capacity of cortical neurons to generate precise discharges, and contribute to the understanding of the biophysics of information transfer of single neurons in vivo at high frequencies.

SIGNIFICANCE STATEMENT

Nanostimulation of single identified neurons in vivo can control spike frequency parametrically and, surprisingly, can even bias the animal's behavioral response. Here, we extend this stimulation protocol to time-dependent broadband noise stimulation in sensory and motor cortices of rat. In response to such stimuli, we found increased temporal spike-time reliability. The information transmission properties reveal, both experimentally and theoretically, that the neurons support high-frequency stimulation beyond the inverse membrane time. Generating a stimulus using the neuron's response properties, we could evoke prescribed spike times with high precision. Our work helps to establish a novel method for precise temporal control of single-cell spiking and provides a simplified biophysical description of single-neuron spiking under time-dependent in vivo-like stimulation.

Dynamical responses to external stimuli for both cases of excitatory and inhibitory synchronization in a complex neuronal network

For studying how dynamical responses to external stimuli depend on the synaptic-coupling type, we consider two types of excitatory and inhibitory synchronization (i.e., synchronization via synaptic excitation and inhibition) in complex small-world networks of excitatory regular spiking (RS) pyramidal neurons and inhibitory fast spiking (FS) interneurons. For both cases of excitatory and inhibitory synchronization, effects of synaptic couplings on dynamical responses to external time-periodic stimuli S(t) (applied to a fraction of neurons) are investigated by varying the driving amplitude A of S(t). Stimulated neurons are phase-locked to external stimuli for both cases of excitatory and inhibitory couplings. On the other hand, the stimulation effect on non-stimulated neurons depends on the type of synaptic coupling. The external stimulus S(t) makes a constructive effect on excitatory non-stimulated RS neurons (i.e., it causes external phase lockings in the non-stimulated sub-population), while S(t) makes a destructive effect on inhibitory non-stimulated FS interneurons (i.e., it breaks up original inhibitory synchronization in the non-stimulated sub-population). As results of these different effects of S(t), the type and degree of dynamical response (e.g., synchronization enhancement or suppression), characterized by the dynamical response factor [Formula: see text] (given by the ratio of synchronization degree in the presence and absence of stimulus), are found to vary in a distinctly different way, depending on the synaptic-coupling type. Furthermore, we also measure the matching degree between the dynamics of the two sub-populations of stimulated and non-stimulated neurons in terms of a "cross-correlation" measure [Formula: see text]. With increasing A, based on [Formula: see text], we discuss the cross-correlations between the two sub-populations, affecting the dynamical responses to S(t).

Emergence of ultrafast sparsely synchronized rhythms and their responses to external stimuli in an inhomogeneous small-world complex neuronal network

We consider an inhomogeneous small-world network (SWN) composed of inhibitory short-range (SR) and long-range (LR) interneurons, and investigate the effect of network architecture on emergence of synchronized brain rhythms by varying the fraction of LR interneurons p_long. The betweenness centralities of the LR and SR interneurons (characterizing the potentiality in controlling communication between other interneurons) are distinctly different. Hence, in view of the betweenness, SWNs we consider are inhomogeneous, unlike the "canonical" Watts-Strogatz SWN with nearly the same betweenness centralities. For small p_long, the load of communication traffic is much concentrated on a few LR interneurons. However, as p_long is increased, the number of LR connections (coming from LR interneurons) increases, and then the load of communication traffic is less concentrated on LR interneurons, which leads to better efficiency of global communication between interneurons. Sparsely synchronized rhythms are thus found to emerge when passing a small critical value p_long^(c)(≃0.16). The population frequency of the sparsely synchronized rhythm is ultrafast (higher than 100 Hz), while the mean firing rate of individual interneurons is much lower (∼30 Hz) due to stochastic and intermittent neural discharges. These dynamical behaviors in the inhomogeneous SWN are also compared with those in the homogeneous Watts-Strogatz SWN, in connection with their network topologies. Particularly, we note that the main difference between the two types of SWNs lies in the distribution of betweenness centralities. Unlike the case of the Watts-Strogatz SWN, dynamical responses to external stimuli vary depending on the type of stimulated interneurons in the inhomogeneous SWN. We consider two cases of external time-periodic stimuli applied to sub-populations of the LR and SR interneurons, respectively. Dynamical responses (such as synchronization suppression and enhancement) to these two cases of stimuli are studied and discussed in relation to the betweenness centralities of stimulated interneurons, representing the effectiveness for transfer of stimulation effect in the whole network.

Modeling fast and slow gamma oscillations with interneurons of different subtype

Experimental and theoretical studies demonstrate that neuronal gamma oscillations crucially depend on interneurons, but current models do not consider the diversity of known interneuron subtypes. Moreover, in CA1 of the hippocampus, experimental evidence indicates the presence of multiple gamma oscillators, two of which may be coordinated by differing interneuron populations. In this article, we show that models of networks with competing interneuron populations with different postsynaptic effects are sufficient to generate, within CA1, distinct oscillatory regimes. We find that strong mutual inhibition between the interneuron populations permits distinct fast and slow gamma states, whereas weak mutual inhibition generates mixed gamma states. We develop idealized firing rate models to illuminate dynamic properties of these competitive gamma networks, and reinforce these concepts with basic spiking models. The models make several explicit predictions about gamma oscillators in CA1. Specifically, interneurons of different subtype phase-lock to different gamma states, and one population of interneurons is silenced and the other active during fast and slow gamma events. Finally, mutual inhibition between interneuron populations is necessary to generate distinct gamma states. Previous experimental studies indicate that fast and slow gamma oscillations reflect different information processing modes, although it is unclear whether these rhythms are intrinsic or imposed. The models outlined demonstrate that basic architectures can locally generate these oscillations, as well as capture other features of fast and slow gamma, including theta-phase preference and spontaneous transitions between gamma states. These models may extend to describe general dynamics in networks with diverse interneuron populations.NEW & NOTEWORTHY The oscillatory coordination of neural signals is crucial to healthy brain function. We have developed an idealized neuronal model that generates distinct fast and slow gamma oscillations, a known feature of the rodent hippocampus. Our work provides a mechanism of this phenomenon, as well as a theoretical framework for future experiments concerning hippocampal gamma. It moreover offers a tractable model of competitive gamma oscillations that is generalizable across the nervous system.

Mathematical Frameworks for Oscillatory Network Dynamics in Neuroscience

The tools of weakly coupled phase oscillator theory have had a profound impact on the neuroscience community, providing insight into a variety of network behaviours ranging from central pattern generation to synchronisation, as well as predicting novel network states such as chimeras. However, there are many instances where this theory is expected to break down, say in the presence of strong coupling, or must be carefully interpreted, as in the presence of stochastic forcing. There are also surprises in the dynamical complexity of the attractors that can robustly appear-for example, heteroclinic network attractors. In this review we present a set of mathematical tools that are suitable for addressing the dynamics of oscillatory neural networks, broadening from a standard phase oscillator perspective to provide a practical framework for further successful applications of mathematics to understanding network dynamics in neuroscience.

Subthalamic, not striatal, activity correlates with basal ganglia downstream activity in normal and parkinsonian monkeys

The striatum and the subthalamic nucleus (STN) constitute the input stage of the basal ganglia (BG) network and together innervate BG downstream structures using GABA and glutamate, respectively. Comparison of the neuronal activity in BG input and downstream structures reveals that subthalamic, not striatal, activity fluctuations correlate with modulations in the increase/decrease discharge balance of BG downstream neurons during temporal discounting classical condition task. After induction of parkinsonism with 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP), abnormal low beta (8-15 Hz) spiking and local field potential (LFP) oscillations resonate across the BG network. Nevertheless, LFP beta oscillations entrain spiking activity of STN, striatal cholinergic interneurons and BG downstream structures, but do not entrain spiking activity of striatal projection neurons. Our results highlight the pivotal role of STN divergent projections in BG physiology and pathophysiology and may explain why STN is such an effective site for invasive treatment of advanced Parkinson's disease and other BG-related disorders.

DOI:http://dx.doi.org/10.7554/eLife.16443.001

The symptoms of Parkinson’s disease include tremor and slow movement, as well as loss of balance, depression and problems with sleep and memory. The death of neurons in a region of the brain called the substantia nigra pars compacta is one of the major hallmarks of Parkinson’s disease. These neurons produce a chemical called dopamine, and their death reduces dopamine levels in another area of the brain called the striatum. This structure is one of five brain regions known collectively as the basal ganglia, which form a circuit that helps to control movement.

The most effective treatment currently available for advanced Parkinson’s disease entails lowering electrodes deep into the brain in order to shut down the activity of part of the basal ganglia. However, the target is not the striatum; instead it is a structure called the subthalamic nucleus. The striatum and the subthalamic nucleus are the two input regions of the basal ganglia: each sends signals to the other three structures downstream. So why does targeting the subthalamic nucleus, but not the striatum, reduce the symptoms of Parkinson’s disease?

To shed some light on this issue, Deffains et al. recorded the activity of neurons in the basal ganglia before and after injecting two monkeys with a drug called MPTP. Related to heroin, MPTP produces symptoms in animals that resemble those of Parkinson’s disease. Before the injections, spontaneous fluctuations in the activity of the subthalamic nucleus produced matching changes in the activity of the three downstream basal ganglia structures. Fluctuations in the activity of the striatum, by contrast, had no such effect. Moreover, injecting the monkeys with MPTP caused the basal ganglia to fire in an abnormal highly synchronized rhythm, similar to that seen in Parkinson’s disease. Crucially, the subthalamic nucleus contributed to this abnormal rhythm, whereas the striatum did not.

The results presented by Deffains et al. provide a concrete explanation for why inactivating the subthalamic nucleus, but not the striatum, reduces the symptoms of Parkinson’s disease. Further research is now needed to explore how the striatum controls the activity of downstream regions of the basal ganglia, both in healthy people and in those with Parkinson's disease.

DOI:http://dx.doi.org/10.7554/eLife.16443.002

Global oscillation regime change by gated inhibition

The role of sensory inputs in the modelling of synchrony regimes is exhibited by means of networks of spiking cells where the relative strength of the inhibitory interaction is controlled by the activation of a linear unit working as a gating variable. Adaptation to stimulus size is determined by the value of a changing length scale, modelled by the time-varying radius of a circular receptive field. In this set-up, 'consolidation' time intervals relevant to attentional effects are shown to depend on the dynamics governing the evolution of the introduced length scale.

Oscillations in the bistable regime of neuronal networks

Bistability between attracting fixed points in neuronal networks has been hypothesized to underlie persistent activity observed in several cortical areas during working memory tasks. In network models this kind of bistability arises due to strong recurrent excitation, sufficient to generate a state of high activity created in a saddle-node (SN) bifurcation. On the other hand, canonical network models of excitatory and inhibitory neurons (E-I networks) robustly produce oscillatory states via a Hopf (H) bifurcation due to the E-I loop. This mechanism for generating oscillations has been invoked to explain the emergence of brain rhythms in the β to γ bands. Although both bistability and oscillatory activity have been intensively studied in network models, there has not been much focus on the coincidence of the two. Here we show that when oscillations emerge in E-I networks in the bistable regime, their phenomenology can be explained to a large extent by considering coincident SN and H bifurcations, known as a codimension two Takens-Bogdanov bifurcation. In particular, we find that such oscillations are not composed of a stable limit cycle, but rather are due to noise-driven oscillatory fluctuations. Furthermore, oscillations in the bistable regime can, in principle, have arbitrarily low frequency.

From Quasiperiodic Partial Synchronization to Collective Chaos in Populations of Inhibitory Neurons with Delay

Collective chaos is shown to emerge, via a period-doubling cascade, from quasiperiodic partial synchronization in a population of identical inhibitory neurons with delayed global coupling. This system is thoroughly investigated by means of an exact model of the macroscopic dynamics, valid in the thermodynamic limit. The collective chaotic state is reproduced numerically with a finite population, and persists in the presence of weak heterogeneities. Finally, the relationship of the model's dynamics with fast neuronal oscillations is discussed.

Dynamical state of the network determines the efficacy of single neuron properties in shaping the network activity

Spike patterns are among the most common electrophysiological descriptors of neuron types. Surprisingly, it is not clear how the diversity in firing patterns of the neurons in a network affects its activity dynamics. Here, we introduce the state-dependent stochastic bursting neuron model allowing for a change in its firing patterns independent of changes in its input-output firing rate relationship. Using this model, we show that the effect of single neuron spiking on the network dynamics is contingent on the network activity state. While spike bursting can both generate and disrupt oscillations, these patterns are ineffective in large regions of the network state space in changing the network activity qualitatively. Finally, we show that when single-neuron properties are made dependent on the population activity, a hysteresis like dynamics emerges. This novel phenomenon has important implications for determining the network response to time-varying inputs and for the network sensitivity at different operating points.

Abnormal Gamma Oscillations in N-Methyl-D-Aspartate Receptor Hypofunction Models of Schizophrenia

N-methyl-D-aspartate receptor (NMDAR) hypofunction in parvalbumin-expressing (PV+) inhibitory neurons (INs) may contribute to symptoms in patients with schizophrenia (SZ). This hypothesis was inspired by studies in humans involving NMDAR antagonists that trigger SZ symptoms. Animal models of SZ using neuropharmacology and genetic knockouts have successfully replicated some of the key observations in human subjects involving alteration of gamma band oscillations (GBO) observed in electroencephalography and magnetoencephalography signals. However, it remains to be seen if NMDAR hypofunction in PV+ neurons is fundamental to the phenotype observed in these models. In this review, we discuss some of the key computational models of GBO and their predictions in the context of NMDAR hypofunction in INs. While PV+ INs have been the main focus of SZ studies in animal models, we also discuss the implications of NMDAR hypofunction in other types of INs using computational models for GBO modulation in the visual cortex.

Hippocampal CA1 Ripples as Inhibitory Transients

Memories are stored and consolidated as a result of a dialogue between the hippocampus and cortex during sleep. Neurons active during behavior reactivate in both structures during sleep, in conjunction with characteristic brain oscillations that may form the neural substrate of memory consolidation. In the hippocampus, replay occurs within sharp wave-ripples: short bouts of high-frequency activity in area CA1 caused by excitatory activation from area CA3. In this work, we develop a computational model of ripple generation, motivated by in vivo rat data showing that ripples have a broad frequency distribution, exponential inter-arrival times and yet highly non-variable durations. Our study predicts that ripples are not persistent oscillations but result from a transient network behavior, induced by input from CA3, in which the high frequency synchronous firing of perisomatic interneurons does not depend on the time scale of synaptic inhibition. We found that noise-induced loss of synchrony among CA1 interneurons dynamically constrains individual ripple duration. Our study proposes a novel mechanism of hippocampal ripple generation consistent with a broad range of experimental data, and highlights the role of noise in regulating the duration of input-driven oscillatory spiking in an inhibitory network.

Our memories are consolidated while we sleep through a bidirectional exchange of information between two brain areas called cortex and hippocampus. Neurons that were active in behavioral tasks reactivate again in both structures during sleep in a process of linking and modifying memories from the short term storage of the hippocampus to permanent storage in the neocortex. This process occurs mainly during short oscillatory hippocampal electrical events called sharp wave-ripples. We propose a novel mechanism of ripple generation consistent with a wide range of experimental data, to explain how hippocampal network properties shape ripple frequency and duration. Understanding the neuronal mechanism underlying ripples is crucial to explaining how the interaction between hippocampus and cortex during sleep enables memory consolidation.

Recovery of Dynamics and Function in Spiking Neural Networks with Closed-Loop Control

There is a growing interest in developing novel brain stimulation methods to control disease-related aberrant neural activity and to address basic neuroscience questions. Conventional methods for manipulating brain activity rely on open-loop approaches that usually lead to excessive stimulation and, crucially, do not restore the original computations performed by the network. Thus, they are often accompanied by undesired side-effects. Here, we introduce delayed feedback control (DFC), a conceptually simple but effective method, to control pathological oscillations in spiking neural networks (SNNs). Using mathematical analysis and numerical simulations we show that DFC can restore a wide range of aberrant network dynamics either by suppressing or enhancing synchronous irregular activity. Importantly, DFC, besides steering the system back to a healthy state, also recovers the computations performed by the underlying network. Finally, using our theory we identify the role of single neuron and synapse properties in determining the stability of the closed-loop system.

Brain stimulation is being used to ease symptoms in several neurological disorders in cases where pharmacological treatment is not effective (anymore). The most common way for stimulation so far has been to apply a fixed, predetermined stimulus irrespective of the actual state of the brain or the condition of the patient. Recently, alternative strategies such as event-triggered stimulation protocols have attracted the interest of researchers. In these protocols the state of the affected brain area is continuously monitored, but the stimulus is only applied if certain criteria are met. Here we go one step further and present a truly closed-loop stimulation protocol. That is, a stimulus is being continuously provided and the magnitude of the stimulus depends, at any point in time, on the ongoing neural activity dynamics of the affected brain area. This results not only in suppression of the pathological activity, but also in a partial recovery of the transfer function of the activity dynamics. Thus, the ability of the lesioned brain area to carry out relevant computations is restored up to a point as well.

Effect of intermodular connection on fast sparse synchronization in clustered small-world neural networks

We consider a clustered network with small-world subnetworks of inhibitory fast spiking interneurons and investigate the effect of intermodular connection on the emergence of fast sparsely synchronized rhythms by varying both the intermodular coupling strength J(inter) and the average number of intermodular links per interneuron M(syn)(inter). In contrast to the case of nonclustered networks, two kinds of sparsely synchronized states such as modular and global synchronization are found. For the case of modular sparse synchronization, the population behavior reveals the modular structure, because the intramodular dynamics of subnetworks make some mismatching. On the other hand, in the case of global sparse synchronization, the population behavior is globally identical, independently of the cluster structure, because the intramodular dynamics of subnetworks make perfect matching. We introduce a realistic cross-correlation modularity measure, representing the matching degree between the instantaneous subpopulation spike rates of the subnetworks, and examine whether the sparse synchronization is global or modular. Depending on its magnitude, the intermodular coupling strength J(inter) seems to play "dual" roles for the pacing between spikes in each subnetwork. For large J(inter), due to strong inhibition it plays a destructive role to "spoil" the pacing between spikes, while for small J(inter) it plays a constructive role to "favor" the pacing between spikes. Through competition between the constructive and the destructive roles of J(inter), there exists an intermediate optimal J(inter) at which the pacing degree between spikes becomes maximal. In contrast, the average number of intermodular links per interneuron M(syn)(inter) seems to play a role just to favor the pacing between spikes. With increasing M(syn)(inter), the pacing degree between spikes increases monotonically thanks to the increase in the degree of effectiveness of global communication between spikes. Furthermore, we employ the realistic sub- and whole-population order parameters, based on the instantaneous sub- and whole-population spike rates, to determine the threshold values for the synchronization-unsynchronization transition in the sub- and whole populations, and the degrees of global and modular sparse synchronization are also measured in terms of the realistic sub- and whole-population statistical-mechanical spiking measures defined by considering both the occupation and the pacing degrees of spikes. It is expected that our results could have implications for the role of the brain plasticity in some functional behaviors associated with population synchronization.

Fast sparsely synchronized brain rhythms in a scale-free neural network

We consider a directed version of the Barabási-Albert scale-free network model with symmetric preferential attachment with the same in- and out-degrees and study the emergence of sparsely synchronized rhythms for a fixed attachment degree in an inhibitory population of fast-spiking Izhikevich interneurons. Fast sparsely synchronized rhythms with stochastic and intermittent neuronal discharges are found to appear for large values of J (synaptic inhibition strength) and D (noise intensity). For an intensive study we fix J at a sufficiently large value and investigate the population states by increasing D. For small D, full synchronization with the same population-rhythm frequency fp and mean firing rate (MFR) fi of individual neurons occurs, while for large D partial synchronization with fp>〈fi〉 (〈fi〉: ensemble-averaged MFR) appears due to intermittent discharge of individual neurons; in particular, the case of fp>4〈fi〉 is referred to as sparse synchronization. For the case of partial and sparse synchronization, MFRs of individual neurons vary depending on their degrees. As D passes a critical value D* (which is determined by employing an order parameter), a transition to unsynchronization occurs due to the destructive role of noise to spoil the pacing between sparse spikes. For D<D*, population synchronization emerges in the whole population because the spatial correlation length between the neuronal pairs covers the whole system. Furthermore, the degree of population synchronization is also measured in terms of two types of realistic statistical-mechanical measures. Only for the partial and sparse synchronization do contributions of individual neuronal dynamics to population synchronization change depending on their degrees, unlike in the case of full synchronization. Consequently, dynamics of individual neurons reveal the inhomogeneous network structure for the case of partial and sparse synchronization, which is in contrast to the case of statistically homogeneous random graphs and small-world networks. Finally, we investigate the effect of network architecture on sparse synchronization for fixed values of J and D in the following three cases: (1) variation in the degree of symmetric attachment, (2) asymmetric preferential attachment of new nodes with different in- and out-degrees, and (3) preferential attachment between pre-existing nodes (without addition of new nodes). In these three cases, both relation between network topology (e.g., average path length and betweenness centralization) and sparse synchronization and contributions of individual dynamics to the sparse synchronization are discussed.

Neuronal morphology generates high-frequency firing resonance

The attenuation of neuronal voltage responses to high-frequency current inputs by the membrane capacitance is believed to limit single-cell bandwidth. However, neuronal populations subject to stochastic fluctuations can follow inputs beyond this limit. We investigated this apparent paradox theoretically and experimentally using Purkinje cells in the cerebellum, a motor structure that benefits from rapid information transfer. We analyzed the modulation of firing in response to the somatic injection of sinusoidal currents. Computational modeling suggested that, instead of decreasing with frequency, modulation amplitude can increase up to high frequencies because of cellular morphology. Electrophysiological measurements in adult rat slices confirmed this prediction and displayed a marked resonance at 200 Hz. We elucidated the underlying mechanism, showing that the two-compartment morphology of the Purkinje cell, interacting with a simple spiking mechanism and dendritic fluctuations, is sufficient to create high-frequency signal amplification. This mechanism, which we term morphology-induced resonance, is selective for somatic inputs, which in the Purkinje cell are exclusively inhibitory. The resonance sensitizes Purkinje cells in the frequency range of population oscillations observed in vivo.

Renewal Approach to the Analysis of the Asynchronous State for Coupled Noisy Oscillators

We develop a framework in which the activity of nonlinear pulse-coupled oscillators is posed within the renewal theory. In this approach, the evolution of the interevent density allows for a self-consistent calculation that determines the asynchronous state and its stability. This framework can readily be extended to the analysis of systems with more state variables and provides a population density treatment to evolve them in their thermodynamical limits. To demonstrate this we study a nonlinear pulse-coupled system, where couplings are dynamic and activity dependent. We investigate its stability and numerically study the nonequilibrium behavior of the system after the bifurcation. We show that this system undergoes a supercritical Hopf bifurcation to collective synchronization.