Turning on and off with excitation: the role of spike-timing asynchrony and synchrony in sustained neural activity

How synaptic function controls critical transitions in spiking neuron networks: insight from a Kuramoto model reduction

The dynamics of synaptic interactions within spiking neuron networks play a fundamental role in shaping emergent collective behavior. This paper studies a finite-size network of quadratic integrate-and-fire neurons interconnected via a general synaptic function that accounts for synaptic dynamics and time delays. Through asymptotic analysis, we transform this integrate-and-fire network into the Kuramoto-Sakaguchi model, whose parameters are explicitly expressed via synaptic function characteristics. This reduction yields analytical conditions on synaptic activation rates and time delays determining whether the synaptic coupling is attractive or repulsive. Our analysis reveals alternating stability regions for synchronous and partially synchronous firing, dependent on slow synaptic activation and time delay. We also demonstrate that the reduced microscopic model predicts the emergence of synchronization, weakly stable cyclops states, and non-stationary regimes remarkably well in the original integrate-and-fire network and its theta neuron counterpart. Our reduction approach promises to open the door to rigorous analysis of rhythmogenesis in networks with synaptic adaptation and plasticity.

Pulse Shape and Voltage-Dependent Synchronization in Spiking Neuron Networks

Pulse-coupled spiking neural networks are a powerful tool to gain mechanistic insights into how neurons self-organize to produce coherent collective behavior. These networks use simple spiking neuron models, such as the θ-neuron or the quadratic integrate-and-fire (QIF) neuron, that replicate the essential features of real neural dynamics. Interactions between neurons are modeled with infinitely narrow pulses, or spikes, rather than the more complex dynamics of real synapses. To make these networks biologically more plausible, it has been proposed that they must also account for the finite width of the pulses, which can have a significant impact on the network dynamics. However, the derivation and interpretation of these pulses are contradictory, and the impact of the pulse shape on the network dynamics is largely unexplored. Here, I take a comprehensive approach to pulse coupling in networks of QIF and θ-neurons. I argue that narrow pulses activate voltage-dependent synaptic conductances and show how to implement them in QIF neurons such that their effect can last through the phase after the spike. Using an exact low-dimensional description for networks of globally coupled spiking neurons, I prove for instantaneous interactions that collective oscillations emerge due to an effective coupling through the mean voltage. I analyze the impact of the pulse shape by means of a family of smooth pulse functions with arbitrary finite width and symmetric or asymmetric shapes. For symmetric pulses, the resulting voltage coupling is not very effective in synchronizing neurons, but pulses that are slightly skewed to the phase after the spike readily generate collective oscillations. The results unveil a voltage-dependent spike synchronization mechanism at the heart of emergent collective behavior, which is facilitated by pulses of finite width and complementary to traditional synaptic transmission in spiking neuron networks.

The h-current controls cortical recurrent network activity through modulation of dendrosomatic communication

A fundamental feature of the cerebral cortex is the ability to rapidly turn on and off maintained activity within ensembles of neurons through recurrent excitation balanced by inhibition. Here we demonstrate that reduction of the h-current, which is especially prominent in pyramidal cell dendrites, strongly increases the ability of local cortical networks to generate maintained recurrent activity. Reduction of the h-current resulted in hyperpolarization and increase in input resistance of both the somata and apical dendrites of layer 5 pyramidal cells, while strongly increasing the dendrosomatic transfer of low (<20 Hz) frequencies, causing an increased responsiveness to dynamic clamp-induced recurrent network-like activity injected into the dendrites and substantially increasing the duration of spontaneous Up states. We propose that modulation of the h-current may strongly control the ability of cortical networks to generate recurrent persistent activity and the formation and dissolution of neuronal ensembles.

Multilayer brain networks can identify the epileptogenic zone and seizure dynamics

Seizure generation, propagation, and termination occur through spatiotemporal brain networks. In this paper, we demonstrate the significance of large-scale brain interactions in high-frequency (80-200Hz) for the identification of the epileptogenic zone (EZ) and seizure evolution. To incorporate the continuity of neural dynamics, here we have modeled brain connectivity constructed from stereoelectroencephalography (SEEG) data during seizures using multilayer networks. After introducing a new measure of brain connectivity for temporal networks, named multilayer eigenvector centrality (mlEVC), we applied a consensus hierarchical clustering on the developed model to identify the EZ as a cluster of nodes with distinctive brain connectivity in the ictal period. Our algorithm could successfully predict electrodes inside the resected volume as EZ for 88% of participants, who all were seizure-free for at least 12 months after surgery. Our findings illustrated significant and unique desynchronization between EZ and the rest of the brain in the early to mid-seizure. We showed that aging and the duration of epilepsy intensify this desynchronization, which can be the outcome of abnormal neuroplasticity. Additionally, we illustrated that seizures evolve with various network topologies, confirming the existence of different epileptogenic networks in each patient. Our findings suggest not only the importance of early intervention in epilepsy but possible factors that correlate with disease severity. Moreover, by analyzing the propagation patterns of different seizures, we demonstrate the necessity of collecting sufficient data for identifying epileptogenic networks.

Collective states in a ring network of theta neurons

We consider a ring network of theta neurons with non-local homogeneous coupling. We analyse the corresponding continuum evolution equation, analytically describing all possible steady states and their stability. By considering a number of different parameter sets, we determine the typical bifurcation scenarios of the network, and put on a rigorous footing some previously observed numerical results.

Interpolating between bumps and chimeras

A "bump" refers to a group of active neurons surrounded by quiescent ones while a "chimera" refers to a pattern in a network in which some oscillators are synchronized while the remainder are asynchronous. Both types of patterns have been studied intensively but are sometimes conflated due to their similar appearance and existence in similar types of networks. Here, we numerically study a hybrid system that linearly interpolates between a network of theta neurons that supports a bump at one extreme and a network of phase oscillators that supports a chimera at the other extreme. Using the Ott/Antonsen ansatz, we derive the equation describing the hybrid network in the limit of an infinite number of oscillators and perform bifurcation analysis on this equation. We find that neither the bump nor chimera persists over the whole range of parameters, and the hybrid system shows a variety of other states such as spatiotemporal chaos, traveling waves, and modulated traveling waves.

Activity Stabilization in a Population Model of Working Memory by Sinusoidal and Noisy Inputs

According to mechanistic theories of working memory (WM), information is retained as stimulus-dependent persistent spiking activity of cortical neural networks. Yet, how this activity is related to changes in the oscillatory profile observed during WM tasks remains a largely open issue. We explore joint effects of input gamma-band oscillations and noise on the dynamics of several firing rate models of WM. The considered models have a metastable active regime, i.e., they demonstrate long-lasting transient post-stimulus firing rate elevation. We start from a single excitatory-inhibitory circuit and demonstrate that either gamma-band or noise input could stabilize the active regime, thus supporting WM retention. We then consider a system of two circuits with excitatory intercoupling. We find that fast coupling allows for better stabilization by common noise compared to independent noise and stronger amplification of this effect by in-phase gamma inputs compared to anti-phase inputs. Finally, we consider a multi-circuit system comprised of two clusters, each containing a group of circuits receiving a common noise input and a group of circuits receiving independent noise. Each cluster is associated with its own local gamma generator, so all its circuits receive gamma-band input in the same phase. We find that gamma-band input differentially stabilizes the activity of the "common-noise" groups compared to the "independent-noise" groups. If the inter-cluster connections are fast, this effect is more pronounced when the gamma-band input is delivered to the clusters in the same phase rather than in the anti-phase. Assuming that the common noise comes from a large-scale distributed WM representation, our results demonstrate that local gamma oscillations can stabilize the activity of the corresponding parts of this representation, with stronger effect for fast long-range connections and synchronized gamma oscillations.

Exact firing rate model reveals the differential effects of chemical versus electrical synapses in spiking networks

Chemical and electrical synapses shape the dynamics of neuronal networks. Numerous theoretical studies have investigated how each of these types of synapses contributes to the generation of neuronal oscillations, but their combined effect is less understood. This limitation is further magnified by the impossibility of traditional neuronal mean-field models-also known as firing rate models or firing rate equations-to account for electrical synapses. Here, we introduce a firing rate model that exactly describes the mean-field dynamics of heterogeneous populations of quadratic integrate-and-fire (QIF) neurons with both chemical and electrical synapses. The mathematical analysis of the firing rate model reveals a well-established bifurcation scenario for networks with chemical synapses, characterized by a codimension-2 cusp point and persistent states for strong recurrent excitatory coupling. The inclusion of electrical coupling generally implies neuronal synchrony by virtue of a supercritical Hopf bifurcation. This transforms the cusp scenario into a bifurcation scenario characterized by three codimension-2 points (cusp, Takens-Bogdanov, and saddle-node separatrix loop), which greatly reduces the possibility for persistent states. This is generic for heterogeneous QIF networks with both chemical and electrical couplings. Our results agree with several numerical studies on the dynamics of large networks of heterogeneous spiking neurons with electrical and chemical couplings.

Network Analysis of Murine Cortical Dynamics Implicates Untuned Neurons in Visual Stimulus Coding

Unbiased and dense sampling of large populations of layer 2/3 pyramidal neurons in mouse primary visual cortex (V1) reveals two functional sub-populations: neurons tuned and untuned to drifting gratings. Whether functional interactions between these two groups contribute to the representation of visual stimuli is unclear. To examine these interactions, we summarize the population partial pairwise correlation structure as a directed and weighted graph. We find that tuned and untuned neurons have distinct topological properties, with untuned neurons occupying central positions in functional networks (FNs). Implementation of a decoder that utilizes the topology of these FNs yields accurate decoding of visual stimuli. We further show that decoding performance degrades comparably following manipulations of either tuned or untuned neurons. Our results demonstrate that untuned neurons are an integral component of V1 FNs and suggest that network interactions contain information about the stimulus that is accessible to downstream elements.

Stochastic neural fields as gradient dynamical systems

Continuous attractor neural networks are used extensively to model a variety of experimentally observed coherent brain states, ranging from cortical waves of activity to stationary activity bumps. The latter are thought to play an important role in various forms of neural information processing, including population coding in primary visual cortex (V1) and working memory in prefrontal cortex. However, one limitation of continuous attractor networks is that the location of the peak of an activity bump (or wave) can diffuse due to intrinsic network noise. This reflects marginal stability of bump solutions with respect to the action of an underlying continuous symmetry group. Previous studies have used perturbation theory to derive an approximate stochastic differential equation for the location of the peak (phase) of the bump. Although this method captures the diffusive wandering of a bump solution, it ignores fluctuations in the amplitude of the bump. In this paper, we show how amplitude fluctuations can be analyzed by reducing the underlying stochastic neural field equation to a finite-dimensional stochastic gradient dynamical system that tracks the stochastic motion of both the amplitude and phase of bump solutions. This allows us to derive exact expressions for the steady-state probability density and its moments, which are then used to investigate two major issues: (i) the input-dependent suppression of neural variability and (ii) noise-induced transitions to bump extinction. We develop the theory by considering the particular example of a ring attractor network with SO(2) symmetry, which is the most common architecture used in attractor models of working memory and population tuning in V1. However, we also extend the analysis to a higher-dimensional spherical attractor network with SO(3) symmetry which has previously been proposed as a model of orientation and spatial frequency tuning in V1. We thus establish how a combination of stochastic analysis and group theoretic methods provides a powerful tool for investigating the effects of noise in continuous attractor networks.

A loop-based neural architecture for structured behavior encoding and decoding

We present a new type of artificial neural network that generalizes on anatomical and dynamical aspects of the mammal brain. Its main novelty lies in its topological structure which is built as an array of interacting elementary motifs shaped like loops. These loops come in various types and can implement functions such as gating, inhibitory or executive control, or encoding of task elements to name a few. Each loop features two sets of neurons and a control region, linked together by non-recurrent projections. The two neural sets do the bulk of the loop's computations while the control unit specifies the timing and the conditions under which the computations implemented by the loop are to be performed. By functionally linking many such loops together, a neural network is obtained that may perform complex cognitive computations. To demonstrate the potential offered by such a system, we present two neural network simulations. The first illustrates the structure and dynamics of a single loop implementing a simple gating mechanism. The second simulation shows how connecting four loops in series can produce neural activity patterns that are sufficient to pass a simplified delayed-response task. We also show that this network reproduces electrophysiological measurements gathered in various regions of the brain of monkeys performing similar tasks. We also demonstrate connections between this type of neural network and recurrent or long short-term memory network models, and suggest ways to generalize them for future artificial intelligence research.

On Rhythms in Neuronal Networks with Recurrent Excitation

We investigate rhythms in networks of neurons with recurrent excitation, that is, with excitatory cells exciting each other. Recurrent excitation can sustain activity even when the cells in the network are driven below threshold, too weak to fire on their own. This sort of "reverberating" activity is often thought to be the basis of working memory. Recurrent excitation can also lead to "runaway" transitions, sudden transitions to high-frequency firing; this may be related to epileptic seizures. Not all fundamental questions about these phenomena have been answered with clarity in the literature. We focus on three questions here: (1) How much recurrent excitation is needed to sustain reverberating activity? How does the answer depend on parameters? (2) Is there a positive minimum frequency of reverberating activity, a positive "onset frequency"? How does it depend on parameters? (3) When do runaway transitions occur? For reduced models, we give mathematical answers to these questions. We also examine computationally to which extent our findings are reflected in the behavior of biophysically more realistic model networks. Our main results can be summarized as follows. (1) Reverberating activity can be fueled by extremely weak slow recurrent excitation, but only by sufficiently strong fast recurrent excitation. (2) The onset of reverberating activity, as recurrent excitation is strengthened or external drive is raised, occurs at a positive frequency. It is faster when the external drive is weaker (and the recurrent excitation stronger). It is slower when the recurrent excitation has a longer decay time constant. (3) Runaway transitions occur only with fast, not with slow, recurrent excitation. We also demonstrate that the relation between reverberating activity fueled by recurrent excitation and runaway transitions can be visualized in an instructive way by a (generalized) cusp catastrophe surface.

Chimera states in a multilayer network of coupled and uncoupled neurons

We study the emergence of chimera states in a multilayer neuronal network, where one layer is composed of coupled and the other layer of uncoupled neurons. Through the multilayer structure, the layer with coupled neurons acts as the medium by means of which neurons in the uncoupled layer share information in spite of the absence of physical connections among them. Neurons in the coupled layer are connected with electrical synapses, while across the two layers, neurons are connected through chemical synapses. In both layers, the dynamics of each neuron is described by the Hindmarsh-Rose square wave bursting dynamics. We show that the presence of two different types of connecting synapses within and between the two layers, together with the multilayer network structure, plays a key role in the emergence of between-layer synchronous chimera states and patterns of synchronous clusters. In particular, we find that these chimera states can emerge in the coupled layer regardless of the range of electrical synapses. Even in all-to-all and nearest-neighbor coupling within the coupled layer, we observe qualitatively identical between-layer chimera states. Moreover, we show that the role of information transmission delay between the two layers must not be neglected, and we obtain precise parameter bounds at which chimera states can be observed. The expansion of the chimera region and annihilation of cluster and fully coherent states in the parameter plane for increasing values of inter-layer chemical synaptic time delay are illustrated using effective range measurements. These results are discussed in the light of neuronal evolution, where the coexistence of coherent and incoherent dynamics during the developmental stage is particularly likely.

Inverse stochastic resonance in networks of spiking neurons

Inverse Stochastic Resonance (ISR) is a phenomenon in which the average spiking rate of a neuron exhibits a minimum with respect to noise. ISR has been studied in individual neurons, but here, we investigate ISR in scale-free networks, where the average spiking rate is calculated over the neuronal population. We use Hodgkin-Huxley model neurons with channel noise (i.e., stochastic gating variable dynamics), and the network connectivity is implemented via electrical or chemical connections (i.e., gap junctions or excitatory/inhibitory synapses). We find that the emergence of ISR depends on the interplay between each neuron's intrinsic dynamical structure, channel noise, and network inputs, where the latter in turn depend on network structure parameters. We observe that with weak gap junction or excitatory synaptic coupling, network heterogeneity and sparseness tend to favor the emergence of ISR. With inhibitory coupling, ISR is quite robust. We also identify dynamical mechanisms that underlie various features of this ISR behavior. Our results suggest possible ways of experimentally observing ISR in actual neuronal systems.

Controlling Working Memory Operations by Selective Gating: The Roles of Oscillations and Synchrony

Working memory (WM) is a primary cognitive function that corresponds to the ability to update, stably maintain, and manipulate short-term memory (ST M) rapidly to perform ongoing cognitive tasks. A prevalent neural substrate of WM coding is persistent neural activity, the property of neurons to remain active after having been activated by a transient sensory stimulus. This persistent activity allows for online maintenance of memory as well as its active manipulation necessary for task performance. WM is tightly capacity limited. Therefore, selective gating of sensory and internally generated information is crucial for WM function. While the exact neural substrate of selective gating remains unclear, increasing evidence suggests that it might be controlled by modulating ongoing oscillatory brain activity. Here, we review experiments and models that linked selective gating, persistent activity, and brain oscillations, putting them in the more general mechanistic context of WM. We do so by defining several operations necessary for successful WM function and then discussing how such operations may be carried out by mechanisms suggested by computational models. We specifically show how oscillatory mechanisms may provide a rapid and flexible active gating mechanism for WM operations.

Robustness of persistent spiking to partial synchronization in a minimal model of synaptically driven self-sustained activity

We study the behavior of a minimal model of synaptically sustained persistent activity that consists of two quadratic integrate-and-fire neurons mutually coupled via excitatory synapses. Importantly, each of the neurons is excitable, as opposed to an oscillator; hence when uncoupled it sits at a subthreshold rest state. When the constituent neurons are mutually coupled via sufficiently strong fast excitatory synapses, the system demonstrates bistability between a fixed point (quiescent background state) and a limit cycle (memory state with synaptically driven spiking activity). Previous work showed that this persistent activity can be stopped by an excitatory input that synchronizes the network. Here we analyzed how this persistent state reacts to partial synchronization. We considered three types of progressively more complex excitatory synaptic kernels: delta pulse, square, and exponential. The first two cases were treated analytically, and the latter case numerically. Using phase-plane methods, we characterized the shape of the region, such that all orbits starting within it correspond to infinite spike trains; this constitutes the persistent activity region. In the case of instant coupling, all such active orbits were neutrally stable; in the case of noninstant coupling, the activity region contained a unique stable limit cycle (so the activity region was the basin of attraction for the limit cycle). This limit cycle corresponded to purely antiphase spiking of two neurons. Increasing synchronization shifted the system toward the border of the activity region, eventually terminating spiking activity. We calculated three measures of robustness of the active state: width of the activity region in the phase plane, critical level of synchronization that can be tolerated by the persistent spiking activity, and speed of reconvergence to the limit cycle. Our analysis revealed that the self-sustained activity is more robust to synchronization when each individual neuron is closer to SNIC bifurcation (closer to being an intrinsic oscillator), the recurrent synaptic excitation is stronger, and the synaptic decay is slower, which is in agreement with the existing data on local circuits in the cortex that show sustained activity.

Neural Computations in a Dynamical System with Multiple Time Scales

Neural systems display rich short-term dynamics at various levels, e.g., spike-frequency adaptation (SFA) at the single-neuron level, and short-term facilitation (STF) and depression (STD) at the synapse level. These dynamical features typically cover a broad range of time scales and exhibit large diversity in different brain regions. It remains unclear what is the computational benefit for the brain to have such variability in short-term dynamics. In this study, we propose that the brain can exploit such dynamical features to implement multiple seemingly contradictory computations in a single neural circuit. To demonstrate this idea, we use continuous attractor neural network (CANN) as a working model and include STF, SFA and STD with increasing time constants in its dynamics. Three computational tasks are considered, which are persistent activity, adaptation, and anticipative tracking. These tasks require conflicting neural mechanisms, and hence cannot be implemented by a single dynamical feature or any combination with similar time constants. However, with properly coordinated STF, SFA and STD, we show that the network is able to implement the three computational tasks concurrently. We hope this study will shed light on the understanding of how the brain orchestrates its rich dynamics at various levels to realize diverse cognitive functions.

Reaction-diffusion-like formalism for plastic neural networks reveals dissipative solitons at criticality

Self-organized structures in networks with spike-timing dependent synaptic plasticity (STDP) are likely to play a central role for information processing in the brain. In the present study we derive a reaction-diffusion-like formalism for plastic feed-forward networks of nonlinear rate-based model neurons with a correlation sensitive learning rule inspired by and being qualitatively similar to STDP. After obtaining equations that describe the change of the spatial shape of the signal from layer to layer, we derive a criterion for the nonlinearity necessary to obtain stable dynamics for arbitrary input. We classify the possible scenarios of signal evolution and find that close to the transition to the unstable regime metastable solutions appear. The form of these dissipative solitons is determined analytically and the evolution and interaction of several such coexistent objects is investigated.

Bumps in Small-World Networks

We consider a network of coupled excitatory and inhibitory theta neurons which is capable of supporting stable spatially-localized “bump” solutions. We randomly add long-range and simultaneously remove short-range connections within the network to form a small-world network and investigate the effects of this rewiring on the existence and stability of the bump solution. We consider two limits in which continuum equations can be derived; bump solutions are fixed points of these equations. We can thus use standard numerical bifurcation analysis to determine the stability of these bumps and to follow them as parameters (such as rewiring probabilities) are varied. We find that under some rewiring schemes bumps are quite robust, whereas in other schemes they can become unstable via Hopf bifurcation or even be destroyed in saddle-node bifurcations.

Multiscale modeling of brain dynamics: from single neurons and networks to mathematical tools

The extreme complexity of the brain naturally requires mathematical modeling approaches on a large variety of scales; the spectrum ranges from single neuron dynamics over the behavior of groups of neurons to neuronal network activity. Thus, the connection between the microscopic scale (single neuron activity) to macroscopic behavior (emergent behavior of the collective dynamics) and vice versa is a key to understand the brain in its complexity. In this work, we attempt a review of a wide range of approaches, ranging from the modeling of single neuron dynamics to machine learning. The models include biophysical as well as data-driven phenomenological models. The discussed models include Hodgkin-Huxley, FitzHugh-Nagumo, coupled oscillators (Kuramoto oscillators, Rössler oscillators, and the Hindmarsh-Rose neuron), Integrate and Fire, networks of neurons, and neural field equations. In addition to the mathematical models, important mathematical methods in multiscale modeling and reconstruction of the causal connectivity are sketched. The methods include linear and nonlinear tools from statistics, data analysis, and time series analysis up to differential equations, dynamical systems, and bifurcation theory, including Granger causal connectivity analysis, phase synchronization connectivity analysis, principal component analysis (PCA), independent component analysis (ICA), and manifold learning algorithms such as ISOMAP, and diffusion maps and equation-free techniques. WIREs Syst Biol Med 2016, 8:438-458. doi: 10.1002/wsbm.1348 For further resources related to this article, please visit the WIREs website.

Stimulus-dependent synchronization in delayed-coupled neuronal networks

Time delay is a general feature of all interactions. Although the effects of delayed interaction are often neglected when the intrinsic dynamics is much slower than the coupling delay, they can be crucial otherwise. We show that delayed coupled neuronal networks support transitions between synchronous and asynchronous states when the level of input to the network changes. The level of input determines the oscillation period of neurons and hence whether time-delayed connections are synchronizing or desynchronizing. We find that synchronizing connections lead to synchronous dynamics, whereas desynchronizing connections lead to out-of-phase oscillations in network motifs and to frustrated states with asynchronous dynamics in large networks. Since the impact of a neuronal network to downstream neurons increases when spikes are synchronous, networks with delayed connections can serve as gatekeeper layers mediating the firing transfer to other regions. This mechanism can regulate the opening and closing of communicating channels between cortical layers on demand.

Anions Govern Cell Volume: A Case Study of Relative Astrocytic and Neuronal Swelling in Spreading Depolarization

Cell volume changes are ubiquitous in normal and pathological activity of the brain. Nevertheless, we know little about the dynamics of cell and tissue swelling, and the differential changes in the volumes of neurons and glia during pathological states such as spreading depolarizations (SD) under ischemic and non–ischemic conditions, and epileptic seizures. By combining the Hodgkin–Huxley type spiking dynamics, dynamic ion concentrations, and simultaneous neuronal and astroglial volume changes into a comprehensive model, we elucidate why glial cells swell more than neurons in SD and the special case of anoxic depolarization (AD), and explore the relative contributions of the two cell types to tissue swelling. Our results demonstrate that anion channels, particularly Cl⁻, are intrinsically connected to cell swelling and blocking these currents prevents changes in cell volume. The model is based on a simple and physiologically realistic description. We introduce model extensions that are either derived purely from first physical principles of electroneutrality, osmosis, and conservation of particles, or by a phenomenological combination of these principles and known physiological facts. This work provides insights into numerous studies related to neuronal and glial volume changes in SD that otherwise seem contradictory, and is broadly applicable to swelling in other cell types and conditions.

Ketamine-Induced Changes in the Signal and Noise of Rule Representation in Working Memory by Lateral Prefrontal Neurons

Working memory dysfunction is an especially debilitating symptom in schizophrenia. The NMDA antagonist ketamine has been successfully used to model working memory deficits in both rodents and nonhuman primates, but how it affects the strength and the consistency of working memory representations remains unclear. Here we recorded single-neuron activity in the lateral prefrontal cortex of macaque monkeys before and after the administration of subanesthetic doses of ketamine in a rule-based working memory task. The rule was instructed with a color cue before each delay period and dictated the correct prosaccadic or antisaccadic response to a peripheral stimulus appearing after the delay. We found that acute ketamine injections both weakened the rule signal across all delay periods and amplified the trial-to-trial variance in neural activities (i.e., noise), both within individual neurons and at the ensemble level, resulting in impaired performance. In the minority of postinjection trials when the animals responded correctly, the preservation of the signal strength during the delay periods was predictive of their subsequent success. Our findings suggest that NMDA receptor function may be critical for establishing the optimal signal-to-noise ratio in information representation by ensembles of prefrontal cortex neurons.

SIGNIFICANCE STATEMENT

In schizophrenia patients, working memory deficit is highly debilitating and currently without any efficacious treatment. An improved understanding of the pathophysiology of this symptom may provide critical information to treatment development. The NMDA antagonist ketamine, when injected at a subanesthetic dose, produces working memory deficit and other schizophrenia-like symptoms in humans and other animals. Here we investigated the effects of ketamine on the representation of abstract rules by prefrontal neurons, while macaque monkeys held the rules in working memory before responding accordingly. We found that ketamine weakened the signal-to-noise ratio in rule representation by simultaneously weakening the signal and augmenting noise. Both processes may be relevant in an effective therapy for working memory impairment in schizophrenia.

Spatiotemporal Spike Coding of Behavioral Adaptation in the Dorsal Anterior Cingulate Cortex

The frontal cortex controls behavioral adaptation in environments governed by complex rules. Many studies have established the relevance of firing rate modulation after informative events signaling whether and how to update the behavioral policy. However, whether the spatiotemporal features of these neuronal activities contribute to encoding imminent behavioral updates remains unclear. We investigated this issue in the dorsal anterior cingulate cortex (dACC) of monkeys while they adapted their behavior based on their memory of feedback from past choices. We analyzed spike trains of both single units and pairs of simultaneously recorded neurons using an algorithm that emulates different biologically plausible decoding circuits. This method permits the assessment of the performance of both spike-count and spike-timing sensitive decoders. In response to the feedback, single neurons emitted stereotypical spike trains whose temporal structure identified informative events with higher accuracy than mere spike count. The optimal decoding time scale was in the range of 70-200 ms, which is significantly shorter than the memory time scale required by the behavioral task. Importantly, the temporal spiking patterns of single units were predictive of the monkeys' behavioral response time. Furthermore, some features of these spiking patterns often varied between jointly recorded neurons. All together, our results suggest that dACC drives behavioral adaptation through complex spatiotemporal spike coding. They also indicate that downstream networks, which decode dACC feedback signals, are unlikely to act as mere neural integrators.

The Role of Cell Volume in the Dynamics of Seizure, Spreading Depression, and Anoxic Depolarization

Cell volume changes are ubiquitous in normal and pathological activity of the brain. Nevertheless, we know little of how cell volume affects neuronal dynamics. We here performed the first detailed study of the effects of cell volume on neuronal dynamics. By incorporating cell swelling together with dynamic ion concentrations and oxygen supply into Hodgkin-Huxley type spiking dynamics, we demonstrate the spontaneous transition between epileptic seizure and spreading depression states as the cell swells and contracts in response to changes in osmotic pressure. Our use of volume as an order parameter further revealed a dynamical definition for the experimentally described physiological ceiling that separates seizure from spreading depression, as well as predicted a second ceiling that demarcates spreading depression from anoxic depolarization. Our model highlights the neuroprotective role of glial K buffering against seizures and spreading depression, and provides novel insights into anoxic depolarization and the relevant cell swelling during ischemia. We argue that the dynamics of seizures, spreading depression, and anoxic depolarization lie along a continuum of the repertoire of the neuron membrane that can be understood only when the dynamic ion concentrations, oxygen homeostasis,and cell swelling in response to osmotic pressure are taken into consideration. Our results demonstrate the feasibility of a unified framework for a wide range of neuronal behaviors that may be of substantial importance in the understanding of and potentially developing universal intervention strategies for these pathological states.

A robust in vivo-like persistent firing supported by a hybrid of intracellular and synaptic mechanisms

Persistent firing is believed to support short-term information retention in the brain. Established hypotheses make use of the recurrent synaptic connectivity to support persistent firing. However, this mechanism is known to suffer from a lack of robustness. On the other hand, persistent firing can be supported by an intrinsic cellular mechanism in multiple brain areas. However, the consequences of having both the intrinsic and the synaptic mechanisms (a hybrid model) on persistent firing remain largely unknown. The goal of this study is to investigate whether a hybrid neural network model with these two mechanisms has advantages over a conventional recurrent network based model. Our computer simulations were based on in vitro recordings obtained from hippocampal CA3 pyramidal cells under cholinergic receptor activation. Calcium activated non-specific cationic (CAN) current supported persistent firing in the Hodgkin-Huxley style cellular models. Our results suggest that the hybrid model supports persistent firing within a physiological frequency range over a wide range of different parameters, eliminating parameter sensitivity issues generally recognized in network based persistent firing. In addition, persistent firing in the hybrid model is substantially more robust against distracting inputs, can coexist with theta frequency oscillations, and supports pattern completion.

Chronic cocaine disrupts mesocortical learning mechanisms

The addictive power of drugs of abuse such as cocaine comes from their ability to hijack natural reward and plasticity mechanisms mediated by dopamine signaling in the brain. Reward learning involves burst firing of midbrain dopamine neurons in response to rewards and cues predictive of reward. The resulting release of dopamine in terminal regions is thought to act as a teaching signaling to areas such as the prefrontal cortex and striatum. In this review, we posit that a pool of extrasynaptic dopaminergic D1-like receptors activated in response to dopamine neuron burst firing serve to enable synaptic plasticity in the prefrontal cortex in response to rewards and their cues. We propose that disruptions in these mechanisms following chronic cocaine use contribute to addiction pathology, in part due to the unique architecture of the mesocortical pathway. By blocking dopamine reuptake in the cortex, cocaine elevates dopamine signaling at these extrasynaptic receptors, prolonging D1-receptor activation and the subsequent activation of intracellular signaling cascades, and thus inducing long-lasting maladaptive plasticity. These cellular adaptations may account for many of the changes in cortical function observed in drug addicts, including an enduring vulnerability to relapse. Therefore, understanding and targeting these neuroadaptations may provide cognitive benefits and help prevent relapse in human drug addicts.

Resolution revolution: epilepsy dynamics at the microscale

Our understanding of the neuronal mechanisms behind epilepsy dynamics has recently advanced due to the application of novel technologies, monitoring hundreds of neurons with single cell resolution. These developments have provided new theories on the relationship between physiological and pathological states, as well as common motifs for the propagation of paroxysmal activity. Although traditional electroencephalogram (EEG) recordings continue to describe normal network oscillations and abnormal epileptic events within and outside of the seizure focus, analysis of epilepsy dynamics at the microscale has found variability in the composition of macroscopically repetitive epileptiform events. These novel results point to heterogeneity in the underlying dynamics of the disorder, highlighting both the need and potential for more specific and targeted therapies.

Modulation of a Single Neuron Has State-Dependent Actions on Circuit Dynamics,

We use computational models to examine the circumstances under which modulation of a single neuron can alter the dynamics of an entire circuit. We show that under some circumstances, the circuit's behavior is robust to the neuromodulation of a single hub neuron, while under other circumstances the same neuromodulatory action of a single neuron can produce a large variety of circuit outcomes.

Cover Figure

A Tradeoff Between Accuracy and Flexibility in a Working Memory Circuit Endowed with Slow Feedback Mechanisms

Recent studies have shown that reverberation underlying mnemonic persistent activity must be slow, to ensure the stability of a working memory system and to give rise to long neural transients capable of accumulation of information over time. Is the slower the underlying process, the better? To address this question, we investigated 3 slow biophysical mechanisms that are activity-dependent and prominently present in the prefrontal cortex: Depolarization-induced suppression of inhibition (DSI), calcium-dependent nonspecific cationic current (ICAN), and short-term facilitation. Using a spiking network model for spatial working memory, we found that these processes enhance the memory accuracy by counteracting noise-induced drifts, heterogeneity-induced biases, and distractors. Furthermore, the incorporation of DSI and ICAN enlarges the range of network's parameter values required for working memory function. However, when a progressively slower process dominates the network, it becomes increasingly more difficult to erase a memory trace. We demonstrate this accuracy-flexibility tradeoff quantitatively and interpret it using a state-space analysis. Our results supports the scenario where N-methyl-d-aspartate receptor-dependent recurrent excitation is the workhorse for the maintenance of persistent activity, whereas slow synaptic or cellular processes contribute to the robustness of mnemonic function in a tradeoff that potentially can be adjusted according to behavioral demands.

Dendritic nonlinearities reduce network size requirements and mediate ON and OFF states of persistent activity in a PFC microcircuit model

Technological advances have unraveled the existence of small clusters of co-active neurons in the neocortex. The functional implications of these microcircuits are in large part unexplored. Using a heavily constrained biophysical model of a L5 PFC microcircuit, we recently showed that these structures act as tunable modules of persistent activity, the cellular correlate of working memory. Here, we investigate the mechanisms that underlie persistent activity emergence (ON) and termination (OFF) and search for the minimum network size required for expressing these states within physiological regimes. We show that (a) NMDA-mediated dendritic spikes gate the induction of persistent firing in the microcircuit. (b) The minimum network size required for persistent activity induction is inversely proportional to the synaptic drive of each excitatory neuron. (c) Relaxation of connectivity and synaptic delay constraints eliminates the gating effect of NMDA spikes, albeit at a cost of much larger networks. (d) Persistent activity termination by increased inhibition depends on the strength of the synaptic input and is negatively modulated by dADP. (e) Slow synaptic mechanisms and network activity contain predictive information regarding the ability of a given stimulus to turn ON and/or OFF persistent firing in the microcircuit model. Overall, this study zooms out from dendrites to cell assemblies and suggests a tight interaction between dendritic non-linearities and network properties (size/connectivity) that may facilitate the short-memory function of the PFC.

Working memory, the ability to retain information for a short period of time, is a fundamental cognitive function that shapes behavior. The cellular correlate of working memory is the prolonged spiking (persistent) activity of neurons in the prefrontal cortex. Impairments of prefrontal cortex functionalities and working memory have been associated with a variety of cognitive disorders, such as schizophrenia, the attention deficit hyperactivity disorder, and drug addiction. Hence, understanding how neurons embedded in the local circuitry support and maintain persistent activity is of outmost importance. Our work uses a multi-level integrative approach spanning from the dendritic, to the neuronal and network levels to identify the key biophysical and anatomical mechanisms contributing to persistent activity, leading to a number of high impact findings: it predicts a tradeoff between dendritic regenerative events and the size of a network expressing persistent activity. It also proposes when and how the persistent state can be stabilized, opening new avenues for pharmacological interventions. Finally, it describes decoding mechanisms for upcoming ON/OFF state transitions, furthering our understanding of information processing in the PFC and shedding new light on the emergence of anticipatory behaviors.

A Hamilton-Jacobi-Bellman approach for termination of seizure-like bursting

We use Hamilton-Jacobi-Bellman methods to find minimum-time and energy-optimal control strategies to terminate seizure-like bursting behavior in a conductance-based neural model. Averaging is used to eliminate fast variables from the model, and a target set is defined through bifurcation analysis of the slow variables of the model. This method is illustrated for a single neuron model and for a network model to illustrate its efficacy in terminating bursting once it begins. This work represents a numerical proof-of-concept that a new class of control strategies can be employed to mitigate bursting, and could ultimately be adapted to treat medically intractible epilepsy in patient-specific models.

EEG oscillatory states: universality, uniqueness and specificity across healthy-normal, altered and pathological brain conditions

For the first time the dynamic repertoires and oscillatory types of local EEG states in 13 diverse conditions (examined over 9 studies) that covered healthy-normal, altered and pathological brain states were quantified within the same methodological and conceptual framework. EEG oscillatory states were assessed by the probability-classification analysis of short-term EEG spectral patterns. The results demonstrated that brain activity consists of a limited repertoire of local EEG states in any of the examined conditions. The size of the state repertoires was associated with changes in cognition and vigilance or neuropsychopathologic conditions. Additionally universal, optional and unique EEG states across 13 diverse conditions were observed. It was demonstrated also that EEG oscillations which constituted EEG states were characteristic for different groups of conditions in accordance to oscillations' functional significance. The results suggested that (a) there is a limit in the number of local states available to the cortex and many ways in which these local states can rearrange themselves and still produce the same global state and (b) EEG individuality is determined by varying proportions of universal, optional and unique oscillatory states. The results enriched our understanding about dynamic microstructure of EEG-signal.

EEG oscillatory states: universality, uniqueness and specificity across healthy-normal, altered and pathological brain conditions

For the first time the dynamic repertoires and oscillatory types of local EEG states in 13 diverse conditions (examined over 9 studies) that covered healthy-normal, altered and pathological brain states were quantified within the same methodological and conceptual framework. EEG oscillatory states were assessed by the probability-classification analysis of short-term EEG spectral patterns. The results demonstrated that brain activity consists of a limited repertoire of local EEG states in any of the examined conditions. The size of the state repertoires was associated with changes in cognition and vigilance or neuropsychopathologic conditions. Additionally universal, optional and unique EEG states across 13 diverse conditions were observed. It was demonstrated also that EEG oscillations which constituted EEG states were characteristic for different groups of conditions in accordance to oscillations' functional significance. The results suggested that (a) there is a limit in the number of local states available to the cortex and many ways in which these local states can rearrange themselves and still produce the same global state and (b) EEG individuality is determined by varying proportions of universal, optional and unique oscillatory states. The results enriched our understanding about dynamic microstructure of EEG-signal.

From cognitive networks to seizures: stimulus evoked dynamics in a coupled cortical network

Epilepsy is one of the most common neuropathologies worldwide. Seizures arising in epilepsy or in seizure disorders are characterized generally by uncontrolled spread of excitation and electrical activity to a limited region or even over the entire cortex. While it is generally accepted that abnormal excessive firing and synchronization of neuron populations lead to seizures, little is known about the precise mechanisms underlying human epileptic seizures, the mechanisms of transitions from normal to paroxysmal activity, or about how seizures spread. Further complication arises in that seizures do not occur with a single type of dynamics but as many different phenotypes and genotypes with a range of patterns, synchronous oscillations, and time courses. The concept of preventing, terminating, or modulating seizures and/or paroxysmal activity through stimulation of brain has also received considerable attention. The ability of such stimulation to prevent or modulate such pathological activity may depend on identifiable parameters. In this work, firing rate networks with inhibitory and excitatory populations were modeled. Network parameters were chosen to model normal working memory behaviors. Two different models of cognitive activity were developed. The first model consists of a single network corresponding to a local area of the brain. The second incorporates two networks connected through sparser recurrent excitatory connectivity with transmission delays ranging from approximately 3 ms within local populations to 15 ms between populations residing in different cortical areas. The effect of excitatory stimulation to activate working memory behavior through selective persistent activation of populations is examined in the models, and the conditions and transition mechanisms through which that selective activation breaks down producing spreading paroxysmal activity and seizure states are characterized. Specifically, we determine critical parameters and architectural changes that produce the different seizure dynamics in the networks. This provides possible mechanisms for seizure generation. Because seizures arise as attractors in a multi-state system, the system may possibly be returned to its baseline state through some particular stimulation. The ability of stimulation to terminate seizure dynamics in the local and distributed models is studied. We systematically examine when this may occur and the form of the stimulation necessary for the range of seizure dynamics. In both the local and distributed network models, termination is possible for all seizure types observed by stimulation possessing some particular configuration of spatial and temporal characteristics.

Correlations in background activity control persistent state stability and allow execution of working memory tasks

Working memory (WM) requires selective information gating, active information maintenance, and rapid active updating. Hence performing a WM task needs rapid and controlled transitions between neural persistent activity and the resting state. We propose that changes in correlations in neural activity provides a mechanism for the required WM operations. As a proof of principle, we implement sustained activity and WM in recurrently coupled spiking networks with neurons receiving excitatory random background activity where background correlations are induced by a common noise source. We first characterize how the level of background correlations controls the stability of the persistent state. With sufficiently high correlations, the sustained state becomes practically unstable, so it cannot be initiated by a transient stimulus. We exploit this in WM models implementing the delay match to sample task by modulating flexibly in time the correlation level at different phases of the task. The modulation sets the network in different working regimes: more prompt to gate in a signal or clear the memory. We examine how the correlations affect the ability of the network to perform the task when distractors are present. We show that in a winner-take-all version of the model, where two populations cross-inhibit, correlations make the distractor blocking robust. In a version of the mode where no cross inhibition is present, we show that appropriate modulation of correlation levels is sufficient to also block the distractor access while leaving the relevant memory trace in tact. The findings presented in this manuscript can form the basis for a new paradigm about how correlations are flexibly controlled by the cortical circuits to execute WM operations.

Induction and modulation of persistent activity in a layer V PFC microcircuit model

Working memory refers to the temporary storage of information and is strongly associated with the prefrontal cortex (PFC). Persistent activity of cortical neurons, namely the activity that persists beyond the stimulus presentation, is considered the cellular correlate of working memory. Although past studies suggested that this type of activity is characteristic of large scale networks, recent experimental evidence imply that small, tightly interconnected clusters of neurons in the cortex may support similar functionalities. However, very little is known about the biophysical mechanisms giving rise to persistent activity in small-sized microcircuits in the PFC. Here, we present a detailed biophysically—yet morphologically simplified—microcircuit model of layer V PFC neurons that incorporates connectivity constraints and is validated against a multitude of experimental data. We show that (a) a small-sized network can exhibit persistent activity under realistic stimulus conditions. (b) Its emergence depends strongly on the interplay of dADP, NMDA, and GABA_B currents. (c) Although increases in stimulus duration increase the probability of persistent activity induction, variability in the stimulus firing frequency does not consistently influence it. (d) Modulation of ionic conductances (I_h, I_D, I_sAHP, I_caL, I_caN, I_caR) differentially controls persistent activity properties in a location dependent manner. These findings suggest that modulation of the microcircuit's firing characteristics is achieved primarily through changes in its intrinsic mechanism makeup, supporting the hypothesis of multiple bi-stable units in the PFC. Overall, the model generates a number of experimentally testable predictions that may lead to a better understanding of the biophysical mechanisms of persistent activity induction and modulation in the PFC.

Flexible frequency control of cortical oscillations enables computations required for working memory

Cognitive effort leads to a seeming cacophony of brain oscillations. For example, during tasks engaging working memory (WM), specific oscillatory frequency bands modulate in space and time. Despite ample data correlating such modulation to task performance, a mechanistic explanation remains elusive. We propose that flexible control of neural oscillations provides a unified mechanism for the rapid and controlled transitions between the computational operations required by WM. We show in a spiking network model that modulating the input oscillation frequency sets the network in different operating modes: rapid memory access and load is enabled by the beta-gamma oscillations, maintaining a memory while ignoring distractors by the theta, rapid memory clearance by the alpha. The various frequency bands determine the dynamic gating regimes enabling the necessary operations for WM, whose succession explains the need for the complex oscillatory brain dynamics during effortful cognition.

Dynamic excitatory and inhibitory gain modulation can produce flexible, robust and optimal decision-making

Behavioural and neurophysiological studies in primates have increasingly shown the involvement of urgency signals during the temporal integration of sensory evidence in perceptual decision-making. Neuronal correlates of such signals have been found in the parietal cortex, and in separate studies, demonstrated attention-induced gain modulation of both excitatory and inhibitory neurons. Although previous computational models of decision-making have incorporated gain modulation, their abstract forms do not permit an understanding of the contribution of inhibitory gain modulation. Thus, the effects of co-modulating both excitatory and inhibitory neuronal gains on decision-making dynamics and behavioural performance remain unclear. In this work, we incorporate time-dependent co-modulation of the gains of both excitatory and inhibitory neurons into our previous biologically based decision circuit model. We base our computational study in the context of two classic motion-discrimination tasks performed in animals. Our model shows that by simultaneously increasing the gains of both excitatory and inhibitory neurons, a variety of the observed dynamic neuronal firing activities can be replicated. In particular, the model can exhibit winner-take-all decision-making behaviour with higher firing rates and within a significantly more robust model parameter range. It also exhibits short-tailed reaction time distributions even when operating near a dynamical bifurcation point. The model further shows that neuronal gain modulation can compensate for weaker recurrent excitation in a decision neural circuit, and support decision formation and storage. Higher neuronal gain is also suggested in the more cognitively demanding reaction time than in the fixed delay version of the task. Using the exact temporal delays from the animal experiments, fast recruitment of gain co-modulation is shown to maximize reward rate, with a timescale that is surprisingly near the experimentally fitted value. Our work provides insights into the simultaneous and rapid modulation of excitatory and inhibitory neuronal gains, which enables flexible, robust, and optimal decision-making.

Perceptual decision-making involves not only simple transformation of sensory information to a motor decision, but can also be modulated by high-level cognition. For example, the latter may include strategic allocation of limited attentional resources over time in a decision task to improve performance. At the neurophysiological level, there is evidence supporting attention-induced neuronal gain modulation of both excitatory and inhibitory cortical neurons. In the context of perceptual discrimination tasks performed by animals, we make use of a biologically inspired computational model of decision-making to understand the computational capabilities of such co-modulation of neuronal gains. We find that dynamic co-modulation of both excitatory and inhibitory neurons is important for flexible, and cognitively demanding decision-making while also enhancing robustness in the decision circuit's functions. Our model captures the neuronal activity and behavioural data in the animal experiments remarkably well. Decision performance in a reaction time task can be optimized, maximizing the rate of receiving reward by using fast gain recruitment. Our experimentally fitted timescale is near the optimal one, suggesting that the animals performed almost optimally. By providing both computational simulations and theoretical analyses, our computational model sheds light into the multiple functions of rapid co-modulation of neuronal gains during decision-making.

Generation and annihilation of localized persistent-activity states in a two-population neural-field model

We investigate the generation and annihilation of persistent localized activity states, so-called bumps, in response to transient spatiotemporal external input in a two-population neural-field model of the Wilson-Cowan type. Such persistent cortical states have been implicated as a biological substrate for short-term working memory, that is, the ability to store stimulus-related information for a few seconds and discard it once it is no longer relevant. In previous studies of the same model it has been established that the stability of bump states hinges on the relative inhibitory constant τ, i.e., the ratio of the time constants governing the dynamics of the inhibitory and excitatory populations: persistent bump states are typically only stable for values of τ smaller than a critical value τcr. We find here that τ is also a key parameter determining whether a transient input can generate a persistent bump state (in the regime where τ<τcr) or not. For small values of τ generation of the persistent states is found to depend only on the overall strength of the transient input, i.e., as long as the magnitude and duration of the excitatory transient input are larger and/or long enough, the persistent state will be activated. For higher values of τ we find that only specific combinations of amplitude and duration leads to persistent activation. For the corresponding annihilation process, no such delicate selectivity on the transient input is observed.

Trading speed and accuracy by coding time: a coupled-circuit cortical model

Our actions take place in space and time, but despite the role of time in decision theory and the growing acknowledgement that the encoding of time is crucial to behaviour, few studies have considered the interactions between neural codes for objects in space and for elapsed time during perceptual decisions. The speed-accuracy trade-off (SAT) provides a window into spatiotemporal interactions. Our hypothesis is that temporal coding determines the rate at which spatial evidence is integrated, controlling the SAT by gain modulation. Here, we propose that local cortical circuits are inherently suited to the relevant spatial and temporal coding. In simulations of an interval estimation task, we use a generic local-circuit model to encode time by ‘climbing’ activity, seen in cortex during tasks with a timing requirement. The model is a network of simulated pyramidal cells and inhibitory interneurons, connected by conductance synapses. A simple learning rule enables the network to quickly produce new interval estimates, which show signature characteristics of estimates by experimental subjects. Analysis of network dynamics formally characterizes this generic, local-circuit timing mechanism. In simulations of a perceptual decision task, we couple two such networks. Network function is determined only by spatial selectivity and NMDA receptor conductance strength; all other parameters are identical. To trade speed and accuracy, the timing network simply learns longer or shorter intervals, driving the rate of downstream decision processing by spatially non-selective input, an established form of gain modulation. Like the timing network's interval estimates, decision times show signature characteristics of those by experimental subjects. Overall, we propose, demonstrate and analyse a generic mechanism for timing, a generic mechanism for modulation of decision processing by temporal codes, and we make predictions for experimental verification.

Studies in neuroscience have characterized how the brain represents objects in space and how these objects are selected for detailed perceptual processing. The selection process entails a decision about which object is favoured by the available evidence over time. This period of time is typically in the range of hundreds of milliseconds and is widely believed to be crucial for decisions, allowing neurons to filter noise in the evidence. Despite the widespread belief that time plays this role in decisions and the growing recognition that the brain estimates elapsed time during perceptual tasks, few studies have considered how the encoding of time effects decision making. We propose that neurons encode time in this range by the same general mechanisms used to select objects for detailed processing, and that these temporal representations determine how long evidence is filtered. To this end, we simulate a perceptual decision by coupling two instances of a neural network widely used to simulate localized regions of the cerebral cortex. One network encodes the passage of time and the other makes decisions based on noisy evidence. The former influences the performance of the latter, reproducing signature characteristics of temporal estimates and perceptual decisions.

Neural 'Bubble' Dynamics Revisited

In this paper, we revisit the work of John G Taylor on neural ‘bubble’ dynamics in two-dimensional neural field models. This builds on original work of Amari in a one-dimensional setting and makes use of the fact that mathematical treatments are much simpler when the firing rate function is chosen to be a Heaviside. In this case, the dynamics of an excited or active region, defining a ‘bubble’, reduce to the dynamics of the boundary. The focus of John’s work was on the properties of radially symmetric ‘bubbles’, including existence and radial stability, with applications to the theory of topographic map formation in self-organising neural networks. As well as reviewing John’s work in this area, we also include some recent results that treat more general classes of perturbations.

Impact of spike train autostructure on probability distribution of joint spike events

The discussion whether temporally coordinated spiking activity really exists and whether it is relevant has been heated over the past few years. To investigate this issue, several approaches have been taken to determine whether synchronized events occur significantly above chance, that is, whether they occur more often than expected if the neurons fire independently. Most investigations ignore or destroy the autostructure of the spiking activity of individual cells or assume Poissonian spiking as a model. Such methods that ignore the autostructure can significantly bias the coincidence statistics. Here, we study the influence of the autostructure on the probability distribution of coincident spiking events between tuples of mutually independent non-Poisson renewal processes. In particular, we consider two types of renewal processes that were suggested as appropriate models of experimental spike trains: a gamma and a log-normal process. For a gamma process, we characterize the shape of the distribution analytically with the Fano factor (FFc). In addition, we perform Monte Carlo estimations to derive the full shape of the distribution and the probability for false positives if a different process type is assumed as was actually present. We also determine how manipulations of such spike trains, here dithering, used for the generation of surrogate data change the distribution of coincident events and influence the significance estimation. We find, first, that the width of the coincidence count distribution and its FFc depend critically and in a nontrivial way on the detailed properties of the structure of the spike trains as characterized by the coefficient of variation CV. Second, the dependence of the FFc on the CV is complex and mostly nonmonotonic. Third, spike dithering, even if as small as a fraction of the interspike interval, can falsify the inference on coordinated firing.

Seizures as imbalanced up states: excitatory and inhibitory conductances during seizure-like events

Precisely timed and dynamically balanced excitatory (E) and inhibitory (I) conductances underlie the basis of neural network activity. Normal E/I balance is often shifted in epilepsy, resulting in neuronal network hyperexcitability and recurrent seizures. However, dynamics of the actual excitatory and inhibitory synaptic conductances (ge and gi, respectively) during seizures remain unknown. To study the dynamics of E and I network balance, we calculated ge and gi during the initiation, body, and termination of seizure-like events (SLEs) in the rat hippocampus in vitro. Repetitive emergent SLEs in 4-aminopyridine (100 μM) and reduced extracellular magnesium (0.6 mM) were recorded in the identified CA1 pyramidal cells (PC) and oriens-lacunosum moleculare (O-LM) interneurons. Calculated ge/gi ratio dynamics showed that the initiation stage of the SLEs was dominated by inhibition in the PCs and was more balanced in the O-LM cells. During the body of the SLEs, the balance shifted toward excitation, with ge and gi peaking in both cell types at nearly the same time. In the termination phase, PCs were again dominated by inhibition, whereas O-LM cells experienced persistent excitatory synaptic barrage. In this way, increased excitability of interneurons may play roles in both seizure initiation (žiburkus J, Cressman JR, Barreto E, Schiff SJ. J Neurophysiol 95: 3948-3954, 2006) and in their termination. Overall, SLE stages can be characterized in PC and O-LM cells by dynamically distinct changes in the balance of ge and gi, where a temporal sequence of imbalance shifts with the changing firing patterns of the cellular subtypes comprising the hyperexcitable microcircuits.

Synchronization and desynchronization in epilepsy: controversies and hypotheses

Epilepsy has been historically seen as a functional brain disorder associated with excessive synchronization of large neuronal populations leading to a hypersynchronous state. Recent evidence showed that epileptiform phenomena, particularly seizures, result from complex interactions between neuronal networks characterized by heterogeneity of neuronal firing and dynamical evolution of synchronization. Desynchronization is often observed preceding seizures or during their early stages; in contrast, high levels of synchronization observed towards the end of seizures may facilitate termination. In this review we discuss cellular and network mechanisms responsible for such complex changes in synchronization. Recent work has identified cell-type-specific inhibitory and excitatory interactions, the dichotomy between neuronal firing and the non-local measurement of local field potentials distant to that firing, and the reflection of the neuronal dark matter problem in non-firing neurons active in seizures. These recent advances have challenged long-established views and are leading to a more rigorous and realistic understanding of the pathophysiology of epilepsy.

Dynamics of feature categorization

In visual and auditory scenes, we are able to identify shared features among sensory objects and group them according to their similarity. This grouping is preattentive and fast and is thought of as an elementary form of categorization by which objects sharing similar features are clustered in some abstract perceptual space. It is unclear what neuronal mechanisms underlie this fast categorization. Here we propose a neuromechanistic model of fast feature categorization based on the framework of continuous attractor networks. The mechanism for category formation does not rely on learning and is based on biologically plausible assumptions, for example, the existence of populations of neurons tuned to feature values, feature-specific interactions, and subthreshold-evoked responses upon the presentation of single objects. When the network is presented with a sequence of stimuli characterized by some feature, the network sums the evoked responses and provides a running estimate of the distribution of features in the input stream. If the distribution of features is structured into different components or peaks (i.e., is multimodal), recurrent excitation amplifies the response of activated neurons, and categories are singled out as emerging localized patterns of elevated neuronal activity (bumps), centered at the centroid of each cluster. The emergence of bump states through sequential, subthreshold activation and the dependence on input statistics is a novel application of attractor networks. We show that the extraction and representation of multiple categories are facilitated by the rich attractor structure of the network, which can sustain multiple stable activity patterns for a robust range of connectivity parameters compatible with cortical physiology.

Saccades during object viewing modulate oscillatory phase in the superior temporal sulcus

Saccadic eye movements (SEMs) are the primary means of gating visual information in primates and strongly influence visual perception. The active exploration of the visual environment ("active vision") via SEMs produces suppression during saccades and enhancement afterward (i.e., during fixation) in occipital visual areas. In lateral temporal lobe visual areas, the influence, if any, of eye movements is less well understood, despite the necessity of these areas for forming coherent percepts of objects. The upper bank of the superior temporal sulcus (uSTS) is one such area whose sensitivity to SEMs is unknown. We therefore examined how saccades modulate local field potentials (LFPs) in the uSTS of macaque monkeys while they viewed face and nonface object stimuli. LFP phase concentration increased following fixation onset in the alpha (8-14 Hz), beta (14-30 Hz), and gamma (30-60 Hz) bands and was distinct from the image-evoked response. Furthermore, near-coincident onsets of fixation and image presentation--like those occurring in active vision--led to enhanced responses through greater phase concentration in the same frequency bands. Finally, single-unit activity was modulated by the phase of alpha, beta, and gamma oscillations, suggesting that the observed phase-locking influences spike timing in uSTS. Previous research implicates phase concentration in these frequency bands as a correlate of perceptual performance (Womelsdorf et al., 2006; Bosman et al., 2009). Together, these results demonstrate sensitivity to eye movements in an object-processing region of the brain and represent a plausible neural basis for the enhancement of object processing during active vision.