The Role of Striatal Feedforward Inhibition in the Maintenance of Absence Seizures.

UNLABELLED: Absence seizures are characterized by brief interruptions of conscious experience accompanied by oscillations of activity synchronized across many brain areas. Although the dynamics of the thalamocortical circuits are traditionally thought to underlie absence seizures, converging experimental evidence supports the key involvement of the basal ganglia (BG). In this theoretical work, we argue that the BG are essential for the maintenance of absence seizures. To this end, we combine analytical calculations with numerical simulations to investigate a computational model of the BG-thalamo-cortical network. We demonstrate that abnormally strong striatal feedforward inhibition can promote synchronous oscillatory activity that persists in the network over several tens of seconds as observed during seizures. We show that these maintained oscillations result from an interplay between the negative feedback through the cortico-subthalamo-nigral pathway and the striatal feedforward inhibition. The negative feedback promotes epileptic oscillations whereas the striatal feedforward inhibition suppresses the positive feedback provided by the cortico-striato-nigral pathway. Our theory is consistent with experimental evidence regarding the influence of BG on seizures (e.g., with the fact that a pharmacological blockade of the subthalamo-nigral pathway suppresses seizures). It also accounts for the observed strong suppression of the striatal output during seizures. Our theory predicts that well-timed transient excitatory inputs to the cortex advance the termination of absence seizures. In contrast with the thalamocortical theory, it also predicts that reducing the synaptic transmission along the cortico-subthalamo-nigral pathway while keeping constant the average firing rate of substantia nigra pars reticulata reduces the incidence of seizures. SIGNIFICANCE STATEMENT: Absence seizures are characterized by brief interruptions of consciousness accompanied by abnormal brain oscillations persisting tens of seconds. Thalamocortical circuits are traditionally thought to underlie absence seizures. However, recent experiments have highlighted the key role of the basal ganglia (BG). This work argues for a novel theory according to which the BG drive the oscillatory patterns of activity occurring during the seizures. It demonstrates that abnormally strong striatal feedforward inhibition promotes synchronous oscillatory activity in the BG-thalamo-cortical network and relate this property to the observed strong suppression of the striatal output during seizures. The theory is compatible with virtually all known experimental results, and it predicts that well-timed transient excitatory inputs to the cortex advance the termination of absence seizures.

Asynchronous Rate Chaos in Spiking Neuronal Circuits.

The brain exhibits temporally complex patterns of activity with features similar to those of chaotic systems. Theoretical studies over the last twenty years have described various computational advantages for such regimes in neuronal systems. Nevertheless, it still remains unclear whether chaos requires specific cellular properties or network architectures, or whether it is a generic property of neuronal circuits. We investigate the dynamics of networks of excitatory-inhibitory (EI) spiking neurons with random sparse connectivity operating in the regime of balance of excitation and inhibition. Combining Dynamical Mean-Field Theory with numerical simulations, we show that chaotic, asynchronous firing rate fluctuations emerge generically for sufficiently strong synapses. Two different mechanisms can lead to these chaotic fluctuations. One mechanism relies on slow I-I inhibition which gives rise to slow subthreshold voltage and rate fluctuations. The decorrelation time of these fluctuations is proportional to the time constant of the inhibition. The second mechanism relies on the recurrent E-I-E feedback loop. It requires slow excitation but the inhibition can be fast. In the corresponding dynamical regime all neurons exhibit rate fluctuations on the time scale of the excitation. Another feature of this regime is that the population-averaged firing rate is substantially smaller in the excitatory population than in the inhibitory population. This is not necessarily the case in the I-I mechanism. Finally, we discuss the neurophysiological and computational significance of our results.

Interference and shaping in sensorimotor adaptations with rewards.

When a perturbation is applied in a sensorimotor transformation task, subjects can adapt and maintain performance by either relying on sensory feedback, or, in the absence of such feedback, on information provided by rewards. For example, in a classical rotation task where movement endpoints must be rotated to reach a fixed target, human subjects can successfully adapt their reaching movements solely on the basis of binary rewards, although this proves much more difficult than with visual feedback. Here, we investigate such a reward-driven sensorimotor adaptation process in a minimal computational model of the task. The key assumption of the model is that synaptic plasticity is gated by the reward. We study how the learning dynamics depend on the target size, the movement variability, the rotation angle and the number of targets. We show that when the movement is perturbed for multiple targets, the adaptation process for the different targets can interfere destructively or constructively depending on the similarities between the sensory stimuli (the targets) and the overlap in their neuronal representations. Destructive interferences can result in a drastic slowdown of the adaptation. As a result of interference, the time to adapt varies non-linearly with the number of targets. Our analysis shows that these interferences are weaker if the reward varies smoothly with the subject's performance instead of being binary. We demonstrate how shaping the reward or shaping the task can accelerate the adaptation dramatically by reducing the destructive interferences. We argue that experimentally investigating the dynamics of reward-driven sensorimotor adaptation for more than one sensory stimulus can shed light on the underlying learning rules.

Short-term plasticity explains irregular persistent activity in working memory tasks.

Persistent activity in cortex is the neural correlate of working memory (WM). In persistent activity, spike trains are highly irregular, even more than in baseline. This seemingly innocuous feature challenges our current understanding of the synaptic mechanisms underlying WM. Here we argue that in WM the prefrontal cortex (PFC) operates in a regime of balanced excitation and inhibition and that the observed temporal irregularity reflects this regime. We show that this requires that nonlinearities underlying the persistent activity are primarily in the neuronal interactions between PFC neurons. We also show that short-term synaptic facilitation can be the physiological substrate of these nonlinearities and that the resulting mechanism of balanced persistent activity is robust, in particular with respect to changes in the connectivity. As an example, we put forward a computational model of the PFC circuit involved in oculomotor delayed response task. The novelty of this model is that recurrent excitatory synapses are facilitating. We demonstrate that this model displays direction-selective persistent activity. We find that, even though the memory eventually degrades because of the heterogeneities, it can be stored for several seconds for plausible network size and connectivity. This model accounts for a large number of experimental findings, such as the findings that have shown that firing is more irregular during the persistent state than during baseline, that the neuronal responses are very diverse, and that the preferred directions during cue and delay periods are strongly correlated but tuning widths are not.

Polybrominated diphenyl ethers in indoor dust from Brazil: assessing demographic differences and human health exposure.

Indoor house dust is considered an important human exposure route to polybrominated diphenyl ethers (PBDEs), which has raised concern about their environmental persistence and toxicity properties. In this study, eight PBDEs (BDE-28, -47, -99, -100, -153, -154, -183, and -209) were determined in house dust from two cities with different socio-demographic characteristics from Brazil, examining possible relationships with factors that potentially influence contamination (population density, economic activities, presence of electronic equipment, and so on) and also estimating the risk of human exposure through oral ingestion and dermal contact. The Σ8PBDE concentration in Sorocaba city ranged between 380 and 4269 ng/g dw, while in Itapetininga city ranged from 106 to 1000 ng/g dw. In both regions, BDE-209 was the most abundantly found congener, followed by BDE-99. House dust from Sorocaba presented significantly greater concentrations of BDE-183 and BDE-209 than Itapetininga. Regarding risk exposure assessment, the estimated daily intake (EDI) of PBDEs was much lower than their respective reference doses (RfDs) in all pathways estimated (non-dietary ingestion and dermal contact). This study provided valuable data to improve the knowledge about the presence and exposure to PBDEs in Brazilian house dust in comparison to other developing countries and the need to control environmental pollution and protect human health.

The mechanism of orientation selectivity in primary visual cortex without a functional map.

Neurons in primary visual cortex (V1) display substantial orientation selectivity even in species where V1 lacks an orientation map, such as in mice and rats. The mechanism underlying orientation selectivity in V1 with such a salt-and-pepper organization is unknown; it is unclear whether a connectivity that depends on feature similarity is required, or a random connectivity suffices. Here we argue for the latter. We study the response to a drifting grating of a network model of layer 2/3 with random recurrent connectivity and feedforward input from layer 4 neurons with random preferred orientations. We show that even though the total feedforward and total recurrent excitatory and inhibitory inputs all have a very weak orientation selectivity, strong selectivity emerges in the neuronal spike responses if the network operates in the balanced excitation/inhibition regime. This is because in this regime the (large) untuned components in the excitatory and inhibitory contributions approximately cancel. As a result the untuned part of the input into a neuron as well as its modulation with orientation and time all have a size comparable to the neuronal threshold. However, the tuning of the F0 and F1 components of the input are uncorrelated and the high-frequency fluctuations are not tuned. This is reflected in the subthreshold voltage response. Remarkably, due to the nonlinear voltage-firing rate transfer function, the preferred orientation of the F0 and F1 components of the spike response are highly correlated.

On the distribution of firing rates in networks of cortical neurons.

The distribution of in vivo average firing rates within local cortical networks has been reported to be highly skewed and long tailed. The distribution of average single-cell inputs, conversely, is expected to be Gaussian by the central limit theorem. This raises the issue of how a skewed distribution of firing rates might result from a symmetric distribution of inputs. We argue that skewed rate distributions are a signature of the nonlinearity of the in vivo f-I curve. During in vivo conditions, ongoing synaptic activity produces significant fluctuations in the membrane potential of neurons, resulting in an expansive nonlinearity of the f-I curve for low and moderate inputs. Here, we investigate the effects of single-cell and network parameters on the shape of the f-I curve and, by extension, on the distribution of firing rates in randomly connected networks.

Bistability and spatiotemporal irregularity in neuronal networks with nonlinear synaptic transmission.

We present a mean-field theory for spiking networks operating in the balanced excitation-inhibition regime, with synapses displaying short-term plasticity. The theory reveals a novel mechanism for bistability which relies on the nonlinearity of the synaptic interactions. As synaptic nonlinearity is mainly controlled by the spiking rates, the different states are stabilized by dynamically generated changes in the noise level. Thus, in both states, the network operates in the fluctuation-driven regime, producing activity patterns characterized by strong spatiotemporal irregularity.

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.

Power-law input-output transfer functions explain the contrast-response and tuning properties of neurons in visual cortex.

We develop a unified model accounting simultaneously for the contrast invariance of the width of the orientation tuning curves (OT) and for the sigmoidal shape of the contrast response function (CRF) of neurons in the primary visual cortex (V1). We determine analytically the conditions for the structure of the afferent LGN and recurrent V1 inputs that lead to these properties for a hypercolumn composed of rate based neurons with a power-law transfer function. We investigate what are the relative contributions of single neuron and network properties in shaping the OT and the CRF. We test these results with numerical simulations of a network of conductance-based model (CBM) neurons and we demonstrate that they are valid and more robust here than in the rate model. The results indicate that because of the acceleration in the transfer function, described here by a power-law, the orientation tuning curves of V1 neurons are more tuned, and their CRF is steeper than those of their inputs. Last, we show that it is possible to account for the diversity in the measured CRFs by introducing heterogeneities either in single neuron properties or in the input to the neurons. We show how correlations among the parameters that characterize the CRF depend on these sources of heterogeneities. Comparison with experimental data suggests that both sources contribute nearly equally to the diversity of CRF shapes observed in V1 neurons.

Mechanisms underlying the response of mouse cortical networks to optogenetic manipulation.

GABAergic interneurons can be subdivided into three subclasses: parvalbumin positive (PV), somatostatin positive (SOM) and serotonin positive neurons. With principal cells (PCs) they form complex networks. We examine PCs and PV responses in mouse anterior lateral motor cortex (ALM) and barrel cortex (S1) upon PV photostimulation in vivo. In ALM layer five and S1, the PV response is paradoxical: photoexcitation reduces their activity. This is not the case in ALM layer 2/3. We combine analytical calculations and numerical simulations to investigate how these results constrain the architecture. Two-population models cannot explain the results. Four-population networks with V1-like architecture account for the data in ALM layer 2/3 and layer 5. Our data in S1 can be explained if SOM neurons receive inputs only from PCs and PV neurons. In both four-population models, the paradoxical effect implies not too strong recurrent excitation. It is not evidence for stabilization by inhibition.

Very long transients, irregular firing, and chaotic dynamics in networks of randomly connected inhibitory integrate-and-fire neurons.

We present results of an extensive numerical study of the dynamics of networks of integrate-and-fire neurons connected randomly through inhibitory interactions. We first consider delayed interactions with infinitely fast rise and decay. Depending on the parameters, the network displays transients which are short or exponentially long in the network size. At the end of these transients, the dynamics settle on a periodic attractor. If the number of connections per neuron is large ( approximately 1000) , this attractor is a cluster state with a short period. In contrast, if the number of connections per neuron is small ( approximately 100) , the attractor has complex dynamics and very long period. During the long transients the neurons fire in a highly irregular manner. They can be viewed as quasistationary states in which, depending on the coupling strength, the pattern of activity is asynchronous or displays population oscillations. In the first case, the average firing rates and the variability of the single-neuron activity are well described by a mean-field theory valid in the thermodynamic limit. Bifurcations of the long transient dynamics from asynchronous to synchronous activity are also well predicted by this theory. The transient dynamics display features reminiscent of stable chaos. In particular, despite being linearly stable, the trajectories of the transient dynamics are destabilized by finite perturbations as small as O(1/N) . We further show that stable chaos is also observed for postsynaptic currents with finite decay time. However, we report in this type of network that chaotic dynamics characterized by positive Lyapunov exponents can also be observed. We show in fact that chaos occurs when the decay time of the synaptic currents is long compared to the synaptic delay, provided that the network is sufficiently large.

Dynamics and orientation selectivity in a cortical model of rodent V1 with excess bidirectional connections.

Recent experiments have revealed fine structure in cortical microcircuitry. In particular, bidirectional connections are more prevalent than expected by chance. Whether this fine structure affects cortical dynamics and function has not yet been studied. Here we investigate the effects of excess bidirectionality in a strongly recurrent network model of rodent V1. We show that reciprocal connections have only a very weak effect on orientation selectivity. We find that excess reciprocity between inhibitory neurons slows down the dynamics and strongly increases the Fano factor, while for reciprocal connections between excitatory and inhibitory neurons it has the opposite effect. In contrast, excess bidirectionality within the excitatory population has a minor effect on the neuronal dynamics. These results can be explained by an effective delayed neuronal self-coupling which stems from the fine structure. Our work suggests that excess bidirectionality between inhibitory neurons decreases the efficiency of feature encoding in cortex by reducing the signal to noise ratio. On the other hand it implies that the experimentally observed strong reciprocity between excitatory and inhibitory neurons improves the feature encoding.

Publisher Correction: Idiosyncratic choice bias naturally emerges from intrinsic stochasticity in neuronal dynamics.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Idiosyncratic choice bias naturally emerges from intrinsic stochasticity in neuronal dynamics.

Idiosyncratic tendency to choose one alternative over others in the absence of an identified reason is a common observation in two-alternative forced-choice experiments. Here we quantify idiosyncratic choice biases in a perceptual discrimination task and a motor task. We report substantial and significant biases in both cases that cannot be accounted for by the experimental context. Then, we present theoretical evidence that even in an idealized experiment, in which the settings are symmetric, idiosyncratic choice bias is expected to emerge from the dynamics of competing neuronal networks. We thus argue that idiosyncratic choice bias reflects the microscopic dynamics of choice and therefore is virtually inevitable in any comparison or decision task.

Spatiotemporal constraints on optogenetic inactivation in cortical circuits.

Optogenetics allows manipulations of genetically and spatially defined neuronal populations with excellent temporal control. However, neurons are coupled with other neurons over multiple length scales, and the effects of localized manipulations thus spread beyond the targeted neurons. We benchmarked several optogenetic methods to inactivate small regions of neocortex. Optogenetic excitation of GABAergic neurons produced more effective inactivation than light-gated ion pumps. Transgenic mice expressing the light-dependent chloride channel GtACR1 produced the most potent inactivation. Generally, inactivation spread substantially beyond the photostimulation light, caused by strong coupling between cortical neurons. Over some range of light intensity, optogenetic excitation of inhibitory neurons reduced activity in these neurons, together with pyramidal neurons, a signature of inhibition-stabilized neural networks ('paradoxical effect'). The offset of optogenetic inactivation was followed by rebound excitation in a light dose-dependent manner, limiting temporal resolution. Our data offer guidance for the design of in vivo optogenetics experiments.

Mechanisms of firing patterns in fast-spiking cortical interneurons.

Cortical fast-spiking (FS) interneurons display highly variable electrophysiological properties. Their spike responses to step currents occur almost immediately following the step onset or after a substantial delay, during which subthreshold oscillations are frequently observed. Their firing patterns include high-frequency tonic firing and rhythmic or irregular bursting (stuttering). What is the origin of this variability? In the present paper, we hypothesize that it emerges naturally if one assumes a continuous distribution of properties in a small set of active channels. To test this hypothesis, we construct a minimal, single-compartment conductance-based model of FS cells that includes transient Na(+), delayed-rectifier K(+), and slowly inactivating d-type K(+) conductances. The model is analyzed using nonlinear dynamical system theory. For small Na(+) window current, the neuron exhibits high-frequency tonic firing. At current threshold, the spike response is almost instantaneous for small d-current conductance, gd, and it is delayed for larger gd. As gd further increases, the neuron stutters. Noise substantially reduces the delay duration and induces subthreshold oscillations. In contrast, when the Na(+) window current is large, the neuron always fires tonically. Near threshold, the firing rates are low, and the delay to firing is only weakly sensitive to noise; subthreshold oscillations are not observed. We propose that the variability in the response of cortical FS neurons is a consequence of heterogeneities in their gd and in the strength of their Na(+) window current. We predict the existence of two types of firing patterns in FS neurons, differing in the sensitivity of the delay duration to noise, in the minimal firing rate of the tonic discharge, and in the existence of subthreshold oscillations. We report experimental results from intracellular recordings supporting this prediction.

Emergent Orientation Selectivity from Random Networks in Mouse Visual Cortex.

The connectivity principles underlying the emergence of orientation selectivity in primary visual cortex (V1) of mammals lacking an orientation map (such as rodents and lagomorphs) are poorly understood. We present a computational model in which random connectivity gives rise to orientation selectivity that matches experimental observations. The model predicts that mouse V1 neurons should exhibit intricate receptive fields in the two-dimensional frequency domain, causing a shift in orientation preferences with spatial frequency. We find evidence for these features in mouse V1 using calcium imaging and intracellular whole-cell recordings.

Competition between feedback loops underlies normal and pathological dynamics in the basal ganglia.

Experiments performed in normal animals suggest that the basal ganglia (BG) are crucial in motor program selection. BG are also involved in movement disorders. In particular, BG neuronal activity in parkinsonian animals and patients is more oscillatory and more synchronous than in normal individuals. We propose a new model for the function and dysfunction of the motor part of BG. We hypothesize that the striatum, the subthalamic nucleus, the internal pallidum (GPi), the thalamus, and the cortex are involved in closed feedback loops. The direct (cortex-striatum-GPi-thalamus-cortex) and the hyperdirect loops (cortex-subthalamic nucleus-GPi-thalamus-cortex), which have different polarities, play a key role in the model. We show that the competition between these two loops provides the BG-cortex system with the ability to perform motor program selection. Under the assumption that dopamine potentiates corticostriatal synaptic transmission, we demonstrate that, in our model, moderate dopamine depletion leads to a complete loss of action selection ability. High depletion can lead to synchronous oscillations. These modifications of the network dynamical state stem from an imbalance between the feedback in the direct and hyperdirect loops when dopamine is depleted. Our model predicts that the loss of selection ability occurs before the appearance of oscillations, suggesting that Parkinson's disease motor impairments are not directly related to abnormal oscillatory activity. Another major prediction of our model is that synchronous oscillations driven by the hyperdirect loop appear in BG after inactivation of the striatum.

A canonical neural mechanism for behavioral variability.

The ability to generate variable movements is essential for learning and adjusting complex behaviours. This variability has been linked to the temporal irregularity of neuronal activity in the central nervous system. However, how neuronal irregularity actually translates into behavioural variability is unclear. Here we combine modelling, electrophysiological and behavioural studies to address this issue. We demonstrate that a model circuit comprising topographically organized and strongly recurrent neural networks can autonomously generate irregular motor behaviours. Simultaneous recordings of neurons in singing finches reveal that neural correlations increase across the circuit driving song variability, in agreement with the model predictions. Analysing behavioural data, we find remarkable similarities in the babbling statistics of 5-6-month-old human infants and juveniles from three songbird species and show that our model naturally accounts for these 'universal' statistics.