The impact of input fluctuations on the frequency-current relationships of layer 5 pyramidal neurons in the rat medial prefrontal cortex

Orbitofrontal cortex contributes to the comparison of values underlying economic choices

Economic choices between goods entail the computation and comparison of subjective values. Previous studies examined neuronal activity in the orbitofrontal cortex (OFC) of monkeys choosing between different types of juices. Three groups of neurons were identified: offer value cells encoding the value of individual offers, chosen juice cells encoding the identity of the chosen juice, and chosen value cells encoding the value of the chosen offer. The encoded variables capture both the input (offer value) and the output (chosen juice, chosen value) of the decision process, suggesting that values are compared within OFC. Recent work demonstrates that choices are causally linked to the activity of offer value cells. Conversely, the hypothesis that OFC contributes to value comparison has not been confirmed. Here we show that weak electrical stimulation of OFC specifically disrupts value comparison without altering offer values. This result implies that neuronal populations in OFC participate in value comparison.

Values encoded in orbitofrontal cortex are causally related to economic choices

In the eighteenth century, Daniel Bernoulli, Adam Smith and Jeremy Bentham proposed that economic choices rely on the computation and comparison of subjective values1. This hypothesis continues to inform modern economic theory2 and research in behavioural economics3, but behavioural measures are ultimately not sufficient to verify the proposal4. Consistent with the hypothesis, when agents make choices, neurons in the orbitofrontal cortex (OFC) encode the subjective value of offered and chosen goods5. Value-encoding cells integrate multiple dimensions6-9, variability in the activity of each cell group correlates with variability in choices10,11 and the population dynamics suggests the formation of a decision12. However, it is unclear whether these neural processes are causally related to choices. More generally, the evidence linking economic choices to value signals in the brain13-15 remains correlational16. Here we show that neuronal activity in the OFC is causal to economic choices. We conducted two experiments using electrical stimulation in rhesus monkeys (Macaca mulatta). Low-current stimulation increased the subjective value of individual offers and thus predictably biased choices. Conversely, high-current stimulation disrupted both the computation and the comparison of subjective values, and thus increased choice variability. These results demonstrate a causal chain linking subjective values encoded in OFC to valuation and choice.

Homogeneous and Narrow Bandwidth of Spike Initiation in Rat L1 Cortical Interneurons

The cortical layer 1 (L1) contains a population of GABAergic interneurons, considered a key component of information integration, processing, and relaying in neocortical networks. In fact, L1 interneurons combine top–down information with feed-forward sensory inputs in layer 2/3 and 5 pyramidal cells (PCs), while filtering their incoming signals. Despite the importance of L1 for network emerging phenomena, little is known on the dynamics of the spike initiation and the encoding properties of its neurons. Using acute brain tissue slices from the rat neocortex, combined with the analysis of an existing database of model neurons, we investigated the dynamical transfer properties of these cells by sampling an entire population of known “electrical classes” and comparing experiments and model predictions. We found the bandwidth of spike initiation to be significantly narrower than in L2/3 and 5 PCs, with values below 100 cycle/s, but without significant heterogeneity in the cell response properties across distinct electrical types. The upper limit of the neuronal bandwidth was significantly correlated to the mean firing rate, as anticipated from theoretical studies but not reported for PCs. At high spectral frequencies, the magnitude of the neuronal response attenuated as a power-law, with an exponent significantly smaller than what was reported for pyramidal neurons and reminiscent of the dynamics of a “leaky” integrate-and-fire model of spike initiation. Finally, most of our in vitro results matched quantitatively the numerical simulations of the models as a further contribution to independently validate the models against novel experimental data.

Mechanisms underlying gain modulation in the cortex

Cortical gain regulation allows neurons to respond adaptively to changing inputs. Neural gain is modulated by internal and external influences, including attentional and arousal states, motor activity and neuromodulatory input. These influences converge to a common set of mechanisms for gain modulation, including GABAergic inhibition, synaptically driven fluctuations in membrane potential, changes in cellular conductance and changes in other biophysical neural properties. Recent work has identified GABAergic interneurons as targets of neuromodulatory input and mediators of state-dependent gain modulation. Here, we review the engagement and effects of gain modulation in the cortex. We highlight key recent findings that link phenomenological observations of gain modulation to underlying cellular and circuit-level mechanisms. Finally, we place these cellular and circuit interactions in the larger context of their impact on perception and cognition.

Spike threshold adaptation diversifies neuronal operating modes in the auditory brain stem

Single neurons function along a spectrum of neuronal operating modes whose properties determine how the output firing activity is generated from synaptic input. The auditory brain stem contains a diversity of neurons, from pure coincidence detectors to pure integrators and those with intermediate properties. We investigated how intrinsic spike initiation mechanisms regulate neuronal operating mode in the avian cochlear nucleus. Although the neurons in one division of the avian cochlear nucleus, nucleus magnocellularis, have been studied in depth, the spike threshold dynamics of the tonically firing neurons of a second division of cochlear nucleus, nucleus angularis (NA), remained unexplained. The input-output functions of tonically firing NA neurons were interrogated with directly injected in vivo-like current stimuli during whole cell patch-clamp recordings in vitro. Increasing the amplitude of the noise fluctuations in the current stimulus enhanced the firing rates in one subset of tonically firing neurons ("differentiators") but not another ("integrators"). We found that spike thresholds showed significantly greater adaptation and variability in the differentiator neurons. A leaky integrate-and-fire neuronal model with an adaptive spike initiation process derived from sodium channel dynamics was fit to the firing responses and could recapitulate >80% of the precise temporal firing across a range of fluctuation and mean current levels. Greater threshold adaptation explained the frequency-current curve changes due to a hyperpolarized shift in the effective adaptation voltage range and longer-lasting threshold adaptation in differentiators. The fine-tuning of the intrinsic properties of different NA neurons suggests they may have specialized roles in spectrotemporal processing.NEW & NOTEWORTHY Avian cochlear nucleus angularis (NA) neurons are responsible for encoding sound intensity for sound localization and spectrotemporal processing. An adaptive spike threshold mechanism fine-tunes a subset of repetitive-spiking neurons in NA to confer coincidence detector-like properties. A model based on sodium channel inactivation properties reproduced the activity via a hyperpolarized shift in adaptation conferring fluctuation sensitivity.

Information transmission of mean and variance coding in integrate-and-fire neurons

Neurons process information by translating continuous signals into patterns of discrete spike times. An open question is how much information these spike times contain about signals which modulate either the mean or the variance of the somatic currents in neurons, as is observed experimentally. Here we calculate the exact information contained in discrete spike times about a continuous signal in both encoding strategies. We show that the information content about mean modulating signals is generally substantially larger than about variance modulating signals for biological parameters. Our analysis further reveals that higher information transmission is associated with a larger proportion of nonlinear signal encoding. Our study measures the complete information content of mean and variance coding and provides a method to determine what fraction of the total information is linearly decodable.

Large and fast human pyramidal neurons associate with intelligence

It is generally assumed that human intelligence relies on efficient processing by neurons in our brain. Although grey matter thickness and activity of temporal and frontal cortical areas correlate with IQ scores, no direct evidence exists that links structural and physiological properties of neurons to human intelligence. Here, we find that high IQ scores and large temporal cortical thickness associate with larger, more complex dendrites of human pyramidal neurons. We show in silico that larger dendritic trees enable pyramidal neurons to track activity of synaptic inputs with higher temporal precision, due to fast action potential kinetics. Indeed, we find that human pyramidal neurons of individuals with higher IQ scores sustain fast action potential kinetics during repeated firing. These findings provide the first evidence that human intelligence is associated with neuronal complexity, action potential kinetics and efficient information transfer from inputs to output within cortical neurons.

Our brains are made up of almost 100 billion brain cells. Each of them acts like a small chip: they collect, process and pass on information in the form of electrical signals. In brain areas that integrate different types of information, such as frontal and temporal lobes, brain cells have larger dendrites – long projections specialized to collect signals. Theoretical studies predict that larger dendrites help cells to initiate electrical signals faster.

Because of difficulty in accessing human neurons, it has been unknown whether any of these features also relate to human intelligence. Previous studies have revealed that people with a higher IQ have a thicker outer layer (the cortex) in areas such as the frontal and temporal lobes. But does a thicker cortex also contain cells with larger dendrites and is their role different?

To test whether smarter brains are equipped with faster and larger cells, Goriounova et al. studied 46 people who needed surgery for brain tumors or epilepsy. Each took an IQ test before the operation. To access the diseased tissue deep in the brain, the surgeon also removed small, undamaged samples of temporal lobe. These samples still contained living cells and their electrical signals were measured in the lab. The experiments showed that cells from people with a higher IQ had larger dendrites that transported information more quickly, especially when they are very active. Computer models were then used to understand how these findings can lead to more efficient information transfer in human neurons.

Traditionally, research on human intelligence has focused on three main strategies: to study brain structure and function, to find genes associated with intelligence and to study the connection between our mind and behavior. Goriounova et al. are the first to take the single-cell perspective and link cell properties to human intelligence. The findings could help connect these separate approaches, and explain how genes for intelligence lead to thicker cortices and faster reaction times in people with higher IQ.

The Changes of Intrinsic Excitability of Pyramidal Neurons in Anterior Cingulate Cortex in Neuropathic Pain

To find satisfactory treatment strategies for neuropathic pain syndromes, the cellular mechanisms should be illuminated. Central sensitization is a generator of pain hypersensitivity, and is mainly reflected in neuronal hyperexcitability in pain pathway. Neuronal excitability depends on two components, the synaptic inputs and the intrinsic excitability. Previous studies have focused on the synaptic plasticity in different forms of pain. But little is known about the changes of neuronal intrinsic excitability in neuropathic pain. To address this question, whole-cell patch clamp recordings were performed to study the synaptic transmission and neuronal intrinsic excitability 1 week after spared nerve injury (SNI) or sham operation in male C57BL/6J mice. We found increased spontaneous excitatory postsynaptic currents (sEPSC) frequency in layer II/III pyramidal neurons of anterior cingulate cortex (ACC) from mice with neuropathic pain. Elevated intrinsic excitability of these neurons after nerve injury was also picked up, which was reflected in gain of input-output curve, inter-spike interval (ISI), spike threshold and Refractory period (RP). Besides firing rate related to neuronal intrinsic excitability, spike timing also plays an important role in neural information processing. The precision of spike timing measured by standard deviation of spike timing (SDST) was decreased in neuropathic pain state. The electrophysiological studies revealed the elevated intrinsic excitation in layer II/III pyramidal neurons of ACC in mice with neuropathic pain, which might contribute to central excitation.

Dynamical response properties of neocortical neurons to conductance-driven time-varying inputs

Ensembles of cortical neurons can track fast-varying inputs and relay them in their spike trains, far beyond the cut-off imposed by membrane passive electrical properties and mean firing rates. Initially explored in silico and later demonstrated experimentally, investigating how neurons respond to sinusoidally modulated stimuli provides a deeper insight into spike initiation mechanisms and information processing than conventional F-I curve methodologies. Besides net membrane currents, physiological synaptic inputs can also induce a stimulus-dependent modulation of the total membrane conductance, which is not reproduced by standard current-clamp protocols. Here, we investigated whether rat cortical neurons can track fast temporal modulations over a noisy conductance background. We also determined input-output transfer properties over a range of conditions, including: distinct presynaptic activation rates, postsynaptic firing rates and variability and type of temporal modulations. We found a very broad signal transfer bandwidth across all conditions, similar large cut-off frequencies and power-law attenuations of fast-varying inputs. At slow and intermediate input modulations, the response gain decreased for increasing output mean firing rates. The gain also decreased significantly for increasing intensities of background synaptic activity, thus generalising earlier studies on F-I curves. We also found a direct correlation between the action potentials' onset rapidness and the neuronal bandwidth. Our novel results extend previous investigations of dynamical response properties to non-stationary and conductance-driven conditions, and provide computational neuroscientists with a novel set of observations that models must capture when aiming to replicate cortical cellular excitability.

The basis of sharp spike onset in standard biophysical models

In most vertebrate neurons, spikes initiate in the axonal initial segment (AIS). When recorded in the soma, they have a surprisingly sharp onset, as if sodium (Na) channels opened abruptly. The main view stipulates that spikes initiate in a conventional manner at the distal end of the AIS, then progressively sharpen as they backpropagate to the soma. We examined the biophysical models used to substantiate this view, and we found that spikes do not initiate through a local axonal current loop that propagates along the axon, but through a global current loop encompassing the AIS and soma, which forms an electrical dipole. Therefore, the phenomenon is not adequately modeled as the backpropagation of an electrical wave along the axon, since the wavelength would be as large as the entire system. Instead, in these models, we found that spike initiation rather follows the critical resistive coupling model proposed recently, where the Na current entering the AIS is matched by the axial resistive current flowing to the soma. Besides demonstrating it by examining the balance of currents at spike initiation, we show that the observed increase in spike sharpness along the axon is artifactual and disappears when an appropriate measure of rapidness is used; instead, somatic onset rapidness can be predicted from spike shape at initiation site. Finally, we reproduce the phenomenon in a two-compartment model, showing that it does not rely on propagation. In these models, the sharp onset of somatic spikes is therefore not an artifact of observing spikes at the incorrect location, but rather the signature that spikes are initiated through a global soma-AIS current loop forming an electrical dipole.

Enhanced Sensitivity to Rapid Input Fluctuations by Nonlinear Threshold Dynamics in Neocortical Pyramidal Neurons

The way in which single neurons transform input into output spike trains has fundamental consequences for network coding. Theories and modeling studies based on standard Integrate-and-Fire models implicitly assume that, in response to increasingly strong inputs, neurons modify their coding strategy by progressively reducing their selective sensitivity to rapid input fluctuations. Combining mathematical modeling with in vitro experiments, we demonstrate that, in L5 pyramidal neurons, the firing threshold dynamics adaptively adjust the effective timescale of somatic integration in order to preserve sensitivity to rapid signals over a broad range of input statistics. For that, a new Generalized Integrate-and-Fire model featuring nonlinear firing threshold dynamics and conductance-based adaptation is introduced that outperforms state-of-the-art neuron models in predicting the spiking activity of neurons responding to a variety of in vivo-like fluctuating currents. Our model allows for efficient parameter extraction and can be analytically mapped to a Generalized Linear Model in which both the input filter—describing somatic integration—and the spike-history filter—accounting for spike-frequency adaptation—dynamically adapt to the input statistics, as experimentally observed. Overall, our results provide new insights on the computational role of different biophysical processes known to underlie adaptive coding in single neurons and support previous theoretical findings indicating that the nonlinear dynamics of the firing threshold due to Na⁺-channel inactivation regulate the sensitivity to rapid input fluctuations.

Over the last decades, a variety of simplified spiking models have been shown to achieve a surprisingly high performance in predicting the neuronal responses to in vitro somatic current injections. Because of the complex adaptive behavior featured by cortical neurons, this success is however restricted to limited stimulus ranges: model parameters optimized for a specific input regime are often inappropriate to describe the response to input currents with different statistical properties. In the present study, a new spiking neuron model is introduced that captures single-neuron computation over a wide range of input statistics and explains different aspects of the neuronal dynamics within a single framework. Our results indicate that complex forms of single neuron adaptation are mediated by the nonlinear dynamics of the firing threshold and that the input-output transformation performed by cortical pyramidal neurons can be intuitively understood in terms of an enhanced Generalized Linear Model in which both the input filter and the spike-history filter adapt to the input statistics.

The domain of neuronal firing on a plane of input current and conductance

The activation of neurotransmitter receptors increases the current flow and membrane conductance and thus controls the firing rate of a neuron. In the present work, we justified the two-dimensional representation of a neuronal input by voltage-independent current and conductance and obtained experimentally and numerically a complete input-output (I/O) function. The dependence of the steady-state firing rate on the input current and conductance was studied as a two-parameter I/O function. We employed the dynamic patch clamp technique in slices to get this dependence for the whole domain of two input signals that evoke stationary spike trains in a single neuron (Ω-domain). As found, the Ω-domain is finite and an additional conductance decreases the range of spike-evoking currents. The I/O function has been reproduced in a Hodgkin-Huxley-like model. Among the simulated effects of different factors on the I/O function, including passive and active membrane properties, external conditions and input signal properties, the most interesting were: the shift of the right boundary of the Ω-domain (corresponding to the exCitation block) leftwards due to the decrease of the maximal potassium conductance; and the reduction of the Ω-domain by the decrease of the maximal sodium concentration. As found in experiments and simulations, the Ω-domain is reduced by the decrease of extracellular sodium concentration, by cooling, and by adding slow potassium currents providing interspike interval adaptation; the Ω-domain height is increased by adding color noise. Our modeling data provided a generalization of I/O dependencies that is consistent with previous studies and our experiments. Our results suggest that both current flow and membrane conductance should be taken into account when determining neuronal firing activity.

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.

A reconfigurable visual-programming library for real-time closed-loop cellular electrophysiology

Most of the software platforms for cellular electrophysiology are limited in terms of flexibility, hardware support, ease of use, or re-configuration and adaptation for non-expert users. Moreover, advanced experimental protocols requiring real-time closed-loop operation to investigate excitability, plasticity, dynamics, are largely inaccessible to users without moderate to substantial computer proficiency. Here we present an approach based on MATLAB/Simulink, exploiting the benefits of LEGO-like visual programming and configuration, combined to a small, but easily extendible library of functional software components. We provide and validate several examples, implementing conventional and more sophisticated experimental protocols such as dynamic-clamp or the combined use of intracellular and extracellular methods, involving closed-loop real-time control. The functionality of each of these examples is demonstrated with relevant experiments. These can be used as a starting point to create and support a larger variety of electrophysiological tools and methods, hopefully extending the range of default techniques and protocols currently employed in experimental labs across the world.

Non-linear Membrane Properties in Entorhinal Cortical Stellate Cells Reduce Modulation of Input-Output Responses by Voltage Fluctuations

The presence of voltage fluctuations arising from synaptic activity is a critical component in models of gain control, neuronal output gating, and spike rate coding. The degree to which individual neuronal input-output functions are modulated by voltage fluctuations, however, is not well established across different cortical areas. Additionally, the extent and mechanisms of input-output modulation through fluctuations have been explored largely in simplified models of spike generation, and with limited consideration for the role of non-linear and voltage-dependent membrane properties. To address these issues, we studied fluctuation-based modulation of input-output responses in medial entorhinal cortical (MEC) stellate cells of rats, which express strong sub-threshold non-linear membrane properties. Using in vitro recordings, dynamic clamp and modeling, we show that the modulation of input-output responses by random voltage fluctuations in stellate cells is significantly limited. In stellate cells, a voltage-dependent increase in membrane resistance at sub-threshold voltages mediated by Na⁺ conductance activation limits the ability of fluctuations to elicit spikes. Similarly, in exponential leaky integrate-and-fire models using a shallow voltage-dependence for the exponential term that matches stellate cell membrane properties, a low degree of fluctuation-based modulation of input-output responses can be attained. These results demonstrate that fluctuation-based modulation of input-output responses is not a universal feature of neurons and can be significantly limited by subthreshold voltage-gated conductances.

The membrane voltage of neurons in vivo is dominated by noisy “background” fluctuations generated by network-based synaptic activity from nearby cells. It has been speculated that membrane voltage fluctuations in neurons play an important role in scaling the relationship between input amplitude and spike rate response. For this to be true, neuronal spike input-output behavior must be sensitive to physiological membrane voltage fluctuations. Using a combination of single cell recordings and modeling, we investigated the mechanisms through which voltage fluctuations modulate neuronal input-output responses. We find that neurons that express an increase in membrane input resistance with depolarization show low levels of noise-mediated modulation of input-output responses due, in part, to voltage trajectories that suppress the likelihood of generating a spike in response to random current input fluctuations. Hence, non-linear membrane properties arising from certain types of voltage-gated conductances limit noise-based modulation of neuronal input-output responses.

High bandwidth synaptic communication and frequency tracking in human neocortex

Neuronal firing, synaptic transmission, and its plasticity form the building blocks for processing and storage of information in the brain. It is unknown whether adult human synapses are more efficient in transferring information between neurons than rodent synapses. To test this, we recorded from connected pairs of pyramidal neurons in acute brain slices of adult human and mouse temporal cortex and probed the dynamical properties of use-dependent plasticity. We found that human synaptic connections were purely depressing and that they recovered three to four times more swiftly from depression than synapses in rodent neocortex. Thereby, during realistic spike trains, the temporal resolution of synaptic information exchange in human synapses substantially surpasses that in mice. Using information theory, we calculate that information transfer between human pyramidal neurons exceeds that of mouse pyramidal neurons by four to nine times, well into the beta and gamma frequency range. In addition, we found that human principal cells tracked fine temporal features, conveyed in received synaptic inputs, at a wider bandwidth than for rodents. Action potential firing probability was reliably phase-locked to input transients up to 1,000 cycles/s because of a steep onset of action potentials in human pyramidal neurons during spike trains, unlike in rodent neurons. Our data show that, in contrast to the widely held views of limited information transfer in rodent depressing synapses, fast recovering synapses of human neurons can actually transfer substantial amounts of information during spike trains. In addition, human pyramidal neurons are equipped to encode high synaptic information content. Thus, adult human cortical microcircuits relay information at a wider bandwidth than rodent microcircuits.

Virtual Electrode Recording Tool for EXtracellular potentials (VERTEX): comparing multi-electrode recordings from simulated and biological mammalian cortical tissue

Local field potentials (LFPs) sampled with extracellular electrodes are frequently used as a measure of population neuronal activity. However, relating such measurements to underlying neuronal behaviour and connectivity is non-trivial. To help study this link, we developed the Virtual Electrode Recording Tool for EXtracellular potentials (VERTEX). We first identified a reduced neuron model that retained the spatial and frequency filtering characteristics of extracellular potentials from neocortical neurons. We then developed VERTEX as an easy-to-use Matlab tool for simulating LFPs from large populations (>100,000 neurons). A VERTEX-based simulation successfully reproduced features of the LFPs from an in vitro multi-electrode array recording of macaque neocortical tissue. Our model, with virtual electrodes placed anywhere in 3D, allows direct comparisons with the in vitro recording setup. We envisage that VERTEX will stimulate experimentalists, clinicians, and computational neuroscientists to use models to understand the mechanisms underlying measured brain dynamics in health and disease.

Emergence of band-pass filtering through adaptive spiking in the owl's cochlear nucleus

In the visual, auditory, and electrosensory modalities, stimuli are defined by first- and second-order attributes. The fast time-pressure signal of a sound, a first-order attribute, is important, for instance, in sound localization and pitch perception, while its slow amplitude-modulated envelope, a second-order attribute, can be used for sound recognition. Ascending the auditory pathway from ear to midbrain, neurons increasingly show a preference for the envelope and are most sensitive to particular envelope modulation frequencies, a tuning considered important for encoding sound identity. The level at which this tuning property emerges along the pathway varies across species, and the mechanism of how this occurs is a matter of debate. In this paper, we target the transition between auditory nerve fibers and the cochlear nucleus angularis (NA). While the owl's auditory nerve fibers simultaneously encode the fast and slow attributes of a sound, one synapse further, NA neurons encode the envelope more efficiently than the auditory nerve. Using in vivo and in vitro electrophysiology and computational analysis, we show that a single-cell mechanism inducing spike threshold adaptation can explain the difference in neural filtering between the two areas. We show that spike threshold adaptation can explain the increased selectivity to modulation frequency, as input level increases in NA. These results demonstrate that a spike generation nonlinearity can modulate the tuning to second-order stimulus features, without invoking network or synaptic mechanisms.

Command-line cellular electrophysiology for conventional and real-time closed-loop experiments

BACKGROUND

Current software tools for electrophysiological experiments are limited in flexibility and rarely offer adequate support for advanced techniques such as dynamic clamp and hybrid experiments, which are therefore limited to laboratories with a significant expertise in neuroinformatics.

NEW METHOD

We have developed lcg, a software suite based on a command-line interface (CLI) that allows performing both standard and advanced electrophysiological experiments. Stimulation protocols for classical voltage and current clamp experiments are defined by a concise and flexible meta description that allows representing complex waveforms as a piece-wise parametric decomposition of elementary sub-waveforms, abstracting the stimulation hardware. To perform complex experiments lcg provides a set of elementary building blocks that can be interconnected to yield a large variety of experimental paradigms.

RESULTS

We present various cellular electrophysiological experiments in which lcg has been employed, ranging from the automated application of current clamp protocols for characterizing basic electrophysiological properties of neurons, to dynamic clamp, response clamp, and hybrid experiments. We finally show how the scripting capabilities behind a CLI are suited for integrating experimental trials into complex workflows, where actual experiment, online data analysis and computational modeling seamlessly integrate.

COMPARISON WITH EXISTING METHODS

We compare lcg with two open source toolboxes, RTXI and RELACS.

CONCLUSIONS

We believe that lcg will greatly contribute to the standardization and reproducibility of both simple and complex experiments. Additionally, on the long run the increased efficiency due to a CLI will prove a great benefit for the experimental community.

Fastest strategy to achieve given number of neuronal firing in theta model

We investigate the constrained optimization of excitatory synaptic input patterns to fastest generate given number of spikes in theta neuron model. Optimal input timings and strengths are identified by using phase plane arguments for discrete input kicks with a given total magnitude. Furthermore, analytical results are conducted to estimate the firing time of given number of spikes resulting from a given input train. We obtain the fastest strategy as the total input size increases. In particular, when the parameter -b is large and total input size G is not so large, there are two candidate strategies to fastest achieve given number of spikes, which depend on the considered parameters. The fastest strategy for some cases of G≫-b to fire m spikes should partition m spikes into m-n+1 spikes for the highest band, with largest g, and one spike for each subsequent n-1 band. When G is sufficiently large, big kick is the fastest strategy. In addition, we establish an optimal value for the dependent variable, θ, where each input should be delivered in a non-threshold-based strategy to fastest achieve given output of subsequent spikes. Moreover, we find that reset and kick strategy is the fastest when G is small and G≫-b. The obtained results can lead to a better understanding of how the period of nonlinear oscillators are affected by different input timings and strengths.

Impact of neuronal properties on network coding: roles of spike initiation dynamics and robust synchrony transfer

Neural networks are more than the sum of their parts, but the properties of those parts are nonetheless important. For instance, neuronal properties affect the degree to which neurons receiving common input will spike synchronously, and whether that synchrony will propagate through the network. Stimulus-evoked synchrony can help or hinder network coding depending on the type of code. In this Perspective, we describe how spike initiation dynamics influence neuronal input-output properties, how those properties affect synchronization, and how synchronization affects network coding. We propose that synchronous and asynchronous spiking can be used to multiplex temporal (synchrony) and rate coding and discuss how pyramidal neurons would be well suited for that task.

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.

Intrinsic firing properties in the avian auditory brain stem allow both integration and encoding of temporally modulated noisy inputs in vitro

The intrinsic properties of tonically firing neurons in the cochlear nucleus contribute to representing average sound intensity by favoring synaptic integration across auditory nerve inputs, reducing phase locking to fine temporal acoustic structure and enhancing envelope locking. To determine whether tonically firing neurons of the avian cochlear nucleus angularis (NA) resemble ideal integrators, we investigated their firing responses to noisy current injections during whole cell patch-clamp recordings in brain slices. One subclass of neurons (36% of tonically firing neurons, mainly subtype tonic III) showed no significant changes in firing rate with noise fluctuations, acting like pure integrators. In contrast, many tonically firing neurons (>60%, mainly subtype tonic I or II) showed a robust sensitivity to noisy current fluctuations, increasing their firing rates with increased fluctuation amplitudes. For noise-sensitive tonic neurons, the firing rate vs. average current curves with noise had larger maximal firing rates, lower gains, and wider dynamic ranges compared with FI curves for current steps without noise. All NA neurons showed fluctuation-driven patterning of spikes with a high degree of temporal reliability and millisecond spike time precision. Single-spiking neurons in NA also responded to noisy currents with higher firing rates and reliable spike trains, although less precisely than nucleus magnocellularis neurons. Thus some NA neurons function as integrators by encoding average input levels over wide dynamic ranges regardless of current fluctuations, others detect the degree of coherence in the inputs, and most encode the temporal patterns contained in their inputs with a high degree of precision.

Predictive features of persistent activity emergence in regular spiking and intrinsic bursting model neurons

Proper functioning of working memory involves the expression of stimulus-selective persistent activity in pyramidal neurons of the prefrontal cortex (PFC), which refers to neural activity that persists for seconds beyond the end of the stimulus. The mechanisms which PFC pyramidal neurons use to discriminate between preferred vs. neutral inputs at the cellular level are largely unknown. Moreover, the presence of pyramidal cell subtypes with different firing patterns, such as regular spiking and intrinsic bursting, raises the question as to what their distinct role might be in persistent firing in the PFC. Here, we use a compartmental modeling approach to search for discriminatory features in the properties of incoming stimuli to a PFC pyramidal neuron and/or its response that signal which of these stimuli will result in persistent activity emergence. Furthermore, we use our modeling approach to study cell-type specific differences in persistent activity properties, via implementing a regular spiking (RS) and an intrinsic bursting (IB) model neuron. We identify synaptic location within the basal dendrites as a feature of stimulus selectivity. Specifically, persistent activity-inducing stimuli consist of activated synapses that are located more distally from the soma compared to non-inducing stimuli, in both model cells. In addition, the action potential (AP) latency and the first few inter-spike-intervals of the neuronal response can be used to reliably detect inducing vs. non-inducing inputs, suggesting a potential mechanism by which downstream neurons can rapidly decode the upcoming emergence of persistent activity. While the two model neurons did not differ in the coding features of persistent activity emergence, the properties of persistent activity, such as the firing pattern and the duration of temporally-restricted persistent activity were distinct. Collectively, our results pinpoint to specific features of the neuronal response to a given stimulus that code for its ability to induce persistent activity and predict differential roles of RS and IB neurons in persistent activity expression.

Memory, referred to as the ability to retain, store and recall information, represents one of the most fundamental cognitive functions in daily life. A significant feature of memory processes is selectivity to particular events or items that are important to our survival and relevant to specific situations. For long-term memory, the selectivity to a specific stimulus is seen both at the behavioral as well as the cellular level. For working memory, a type of short-term memory involved in decision making and attention processes, stimulus selectivity has been observed in vivo using spatial working memory tasks. In addition, persistent activity, which is the cellular correlate of working memory, is also selective to specific stimuli for each neuron, suggesting that each neuron has a ‘memory field’. Our study proposes that both the location of incoming inputs onto the neuronal dendritic tree and specific temporal features of the neuronal response can be used to predict the emergence of persistent activity in two neuron models with different firing patterns, revealing possible mechanisms for generating and propagating stimulus-selectivity in working memory processes. The study also reveals that neurons with different firing patterns may have different roles in persistent activity expression.

Single neuron firing properties impact correlation-based population coding

Correlated spiking has been widely observed, but its impact on neural coding remains controversial. Correlation arising from comodulation of rates across neurons has been shown to vary with the firing rates of individual neurons. This translates into rate and correlation being equivalently tuned to the stimulus; under those conditions, correlated spiking does not provide information beyond that already available from individual neuron firing rates. Such correlations are irrelevant and can reduce coding efficiency by introducing redundancy. Using simulations and experiments in rat hippocampal neurons, we show here that pairs of neurons receiving correlated input also exhibit correlations arising from precise spike-time synchronization. Contrary to rate comodulation, spike-time synchronization is unaffected by firing rate, thus enabling synchrony- and rate-based coding to operate independently. The type of output correlation depends on whether intrinsic neuron properties promote integration or coincidence detection: "ideal" integrators (with spike generation sensitive to stimulus mean) exhibit rate comodulation, whereas ideal coincidence detectors (with spike generation sensitive to stimulus variance) exhibit precise spike-time synchronization. Pyramidal neurons are sensitive to both stimulus mean and variance, and thus exhibit both types of output correlation proportioned according to which operating mode is dominant. Our results explain how different types of correlations arise based on how individual neurons generate spikes, and why spike-time synchronization and rate comodulation can encode different stimulus properties. Our results also highlight the importance of neuronal properties for population-level coding insofar as neural networks can employ different coding schemes depending on the dominant operating mode of their constituent neurons.

Inferring network dynamics and neuron properties from population recordings

Understanding the computational capabilities of the nervous system means to “identify” its emergent multiscale dynamics. For this purpose, we propose a novel model-driven identification procedure and apply it to sparsely connected populations of excitatory integrate-and-fire neurons with spike frequency adaptation (SFA). Our method does not characterize the system from its microscopic elements in a bottom-up fashion, and does not resort to any linearization. We investigate networks as a whole, inferring their properties from the response dynamics of the instantaneous discharge rate to brief and aspecific supra-threshold stimulations. While several available methods assume generic expressions for the system as a black box, we adopt a mean-field theory for the evolution of the network transparently parameterized by identified elements (such as dynamic timescales), which are in turn non-trivially related to single-neuron properties. In particular, from the elicited transient responses, the input–output gain function of the neurons in the network is extracted and direct links to the microscopic level are made available: indeed, we show how to extract the decay time constant of the SFA, the absolute refractory period and the average synaptic efficacy. In addition and contrary to previous attempts, our method captures the system dynamics across bifurcations separating qualitatively different dynamical regimes. The robustness and the generality of the methodology is tested on controlled simulations, reporting a good agreement between theoretically expected and identified values. The assumptions behind the underlying theoretical framework make the method readily applicable to biological preparations like cultured neuron networks and in vitro brain slices.

Cellular and network contributions to vestibular signal processing: impact of ion conductances, synaptic inhibition, and noise

Head motion-related sensory signals are transformed by second-order vestibular neurons (2°VNs) into appropriate commands for retinal image stabilization during body motion. In frogs, these 2°VNs form two distinct subpopulations that have either tonic or highly phasic intrinsic properties, essentially compatible with low-pass and bandpass filter characteristics, respectively. In the present study, physiological data on cellular properties of 2°VNs of the grass frog (Rana temporaria) have been used to construct conductance-based spiking cellular models that were fine-tuned by fitting to recorded spike-frequency data. The results of this approach suggest that low-threshold, voltage-dependent potassium channels in phasic and spike-dependent potassium channels in tonic 2°VNs are important contributors to the differential, yet complementary response characteristics of the two vestibular subtypes. Extension of the cellular model with conductance-based synapses allowed simulation of afferent excitation and evaluation of the emerging properties of local feedforward inhibitory circuits. This approach revealed the relative contributions of intrinsic and synaptic factors on afferent signal processing in phasic 2°VNs. Additional extension of the single-cell model to a population model allowed testing under more natural conditions including asynchronous afferent labyrinthine input and synaptic noise. This latter approach indicated that the feedforward inhibition from the local inhibitory network acts as a high-pass filter, which reinforces the impact of the intrinsic membrane properties of phasic 2°VNs on peak response amplitude and timing. Thus, the combination of cellular and network properties enables phasic 2°VNs to work as a noise-resistant detector, suitable for central processing of short-duration vestibular signals.

Effective stimuli for constructing reliable neuron models

The rich dynamical nature of neurons poses major conceptual and technical challenges for unraveling their nonlinear membrane properties. Traditionally, various current waveforms have been injected at the soma to probe neuron dynamics, but the rationale for selecting specific stimuli has never been rigorously justified. The present experimental and theoretical study proposes a novel framework, inspired by learning theory, for objectively selecting the stimuli that best unravel the neuron's dynamics. The efficacy of stimuli is assessed in terms of their ability to constrain the parameter space of biophysically detailed conductance-based models that faithfully replicate the neuron's dynamics as attested by their ability to generalize well to the neuron's response to novel experimental stimuli. We used this framework to evaluate a variety of stimuli in different types of cortical neurons, ages and animals. Despite their simplicity, a set of stimuli consisting of step and ramp current pulses outperforms synaptic-like noisy stimuli in revealing the dynamics of these neurons. The general framework that we propose paves a new way for defining, evaluating and standardizing effective electrical probing of neurons and will thus lay the foundation for a much deeper understanding of the electrical nature of these highly sophisticated and non-linear devices and of the neuronal networks that they compose.

Membrane voltage fluctuations reduce spike frequency adaptation and preserve output gain in CA1 pyramidal neurons in a high-conductance state

Modulating the gain of the input-output function of neurons is critical for processing of stimuli and network dynamics. Previous gain control mechanisms have suggested that voltage fluctuations play a key role in determining neuronal gain in vivo. Here we show that, under increased membrane conductance, voltage fluctuations restore Na(+) current and reduce spike frequency adaptation in rat hippocampal CA1 pyramidal neurons in vitro. As a consequence, membrane voltage fluctuations produce a leftward shift in the frequency-current relationship without a change in gain, relative to an increase in conductance alone. Furthermore, we show that these changes have important implications for the integration of inhibitory inputs. Due to the ability to restore Na(+) current, hyperpolarizing membrane voltage fluctuations mediated by GABA(A)-like inputs can increase firing rate in a high-conductance state. Finally, our data show that the effects on gain and synaptic integration are mediated by voltage fluctuations within a physiologically relevant range of frequencies (10-40 Hz).

Impact of fast sodium channel inactivation on spike threshold dynamics and synaptic integration

Neurons spike when their membrane potential exceeds a threshold value. In central neurons, the spike threshold is not constant but depends on the stimulation. Thus, input-output properties of neurons depend both on the effect of presynaptic spikes on the membrane potential and on the dynamics of the spike threshold. Among the possible mechanisms that may modulate the threshold, one strong candidate is Na channel inactivation, because it specifically impacts spike initiation without affecting the membrane potential. We collected voltage-clamp data from the literature and we found, based on a theoretical criterion, that the properties of Na inactivation could indeed cause substantial threshold variability by itself. By analyzing simple neuron models with fast Na inactivation (one channel subtype), we found that the spike threshold is correlated with the mean membrane potential and negatively correlated with the preceding depolarization slope, consistent with experiments. We then analyzed the impact of threshold dynamics on synaptic integration. The difference between the postsynaptic potential (PSP) and the dynamic threshold in response to a presynaptic spike defines an effective PSP. When the neuron is sufficiently depolarized, this effective PSP is briefer than the PSP. This mechanism regulates the temporal window of synaptic integration in an adaptive way. Finally, we discuss the role of other potential mechanisms. Distal spike initiation, channel noise and Na activation dynamics cannot account for the observed negative slope-threshold relationship, while adaptive conductances (e.g. K+) and Na inactivation can. We conclude that Na inactivation is a metabolically efficient mechanism to control the temporal resolution of synaptic integration.

Neurons spike when their combined inputs exceed a threshold value, but recent experimental findings have shown that this value also depends on the inputs. Thus, to understand how neurons respond to input spikes, it is important to know how inputs modify the spike threshold. Spikes are generated by sodium channels, which inactivate when the neuron is depolarized, raising the threshold for spike initiation. We found that inactivation properties of sodium channels could indeed cause substantial threshold variability in central neurons. We then analyzed in models the implications of this form of threshold modulation on neuronal function. We found that this mechanism makes neurons more sensitive to coincident spikes and provides them with an energetically efficient form of gain control.

Tailoring inputs to achieve maximal neuronal firing

We consider the constrained optimization of excitatory synaptic input patterns to maximize spike generation in leaky integrate-and-fire (LIF) and theta model neurons. In the case of discrete input kicks with a fixed total magnitude, optimal input timings and strengths are identified for each model using phase plane arguments. In both cases, optimal features relate to finding an input level at which the drop in input between successive spikes is minimized. A bounded minimizing level always exists in the theta model and may or may not exist in the LIF model, depending on parameter tuning. We also provide analytical formulas to estimate the number of spikes resulting from a given input train. In a second case of continuous inputs of fixed total magnitude, we analyze the tuning of an input shape parameter to maximize the number of spikes occurring in a fixed time interval. Results are obtained using numerical solution of a variational boundary value problem that we derive, as well as analysis, for the theta model and using a combination of simulation and analysis for the LIF model. In particular, consistent with the discrete case, the number of spikes in the theta model rises and then falls again as the input becomes more tightly peaked. Under a similar variation in the LIF case, we numerically show that the number of spikes increases monotonically up to some bound and we analytically constrain the times at which spikes can occur and estimate the bound on the number of spikes fired.

The dynamics of integrate-and-fire: mean versus variance modulations and dependence on baseline parameters

The leaky integrate-and-fire (LIF) is the simplest neuron model that captures the essential properties of neuronal signaling. Yet common intuitions are inadequate to explain basic properties of LIF responses to sinusoidal modulations of the input. Here we examine responses to low and moderate frequency modulations of both the mean and variance of the input current and quantify how these responses depend on baseline parameters. Across parameters, responses to modulations in the mean current are low pass, approaching zero in the limit of high frequencies. For very low baseline firing rates, the response cutoff frequency matches that expected from membrane integration. However, the cutoff shows a rapid, supralinear increase with firing rate, with a steeper increase in the case of lower noise. For modulations of the input variance, the gain at high frequency remains finite. Here, we show that the low-frequency responses depend strongly on baseline parameters and derive an analytic condition specifying the parameters at which responses switch from being dominated by low versus high frequencies. Additionally, we show that the resonant responses for variance modulations have properties not expected for common oscillatory resonances: they peak at frequencies higher than the baseline firing rate and persist when oscillatory spiking is disrupted by high noise. Finally, the responses to mean and variance modulations are shown to have a complementary dependence on baseline parameters at higher frequencies, resulting in responses to modulations of Poisson input rates that are independent of baseline input statistics.

Rate Dynamics of Leaky Integrate-and-Fire Neurons with Strong Synapses

Firing-rate models provide a practical tool for studying the dynamics of trial- or population-averaged neuronal signals. A wealth of theoretical and experimental studies has been dedicated to the derivation or extraction of such models by investigating the firing-rate response characteristics of ensembles of neurons. The majority of these studies assumes that neurons receive input spikes at a high rate through weak synapses (diffusion approximation). For many biological neural systems, however, this assumption cannot be justified. So far, it is unclear how time-varying presynaptic firing rates are transmitted by a population of neurons if the diffusion assumption is dropped. Here, we numerically investigate the stationary and non-stationary firing-rate response properties of leaky integrate-and-fire neurons receiving input spikes through excitatory synapses with alpha-function shaped postsynaptic currents for strong synaptic weights. Input spike trains are modeled by inhomogeneous Poisson point processes with sinusoidal rate. Average rates, modulation amplitudes, and phases of the period-averaged spike responses are measured for a broad range of stimulus, synapse, and neuron parameters. Across wide parameter regions, the resulting transfer functions can be approximated by a linear first-order low-pass filter. Below a critical synaptic weight, the cutoff frequencies are approximately constant and determined by the synaptic time constants. Only for synapses with unrealistically strong weights are the cutoff frequencies significantly increased. To account for stimuli with larger modulation depths, we combine the measured linear transfer function with the nonlinear response characteristics obtained for stationary inputs. The resulting linear-nonlinear model accurately predicts the population response for a variety of non-sinusoidal stimuli.

GABAB receptor-dependent modulation of network activity in the rat prefrontal cortex in vitro

GABA (gamma-aminobutyric acid) can mediate inhibition via pre- and post/extrasynaptic GABA receptors. In this paper we demonstrate potentially post/extrasynaptic GABA(B) receptor-dependent tonic inhibition in L2/3 pyramidal cells of rat medial prefrontal cortex (mPFC) in vitro. First, we show via voltage-clamp experiments the presence of a tonic GABA(B) receptor-dependent outward current in these neurons. This GABA(B)ergic current could be induced by ambient GABA when present at sufficient concentrations. To increase ambient GABA levels in the usually silent slice preparation, we amplified network activity and hence synaptic GABA release with a modified artificial cerebrospinal fluid. The amplitude of tonic GABA(B) current was similar at different temperatures. In addition to the tonic GABA(B) current, we found presynaptic GABA(B) effects, GABA(B)-mediated inhibitory postsynaptic currents and tonic GABA(A) currents. Second, we performed current-clamp experiments to evaluate the functional impact of GABA(B) receptor-mediated inhibition in the mPFC. Activating or inactivating GABA(B) receptors led to rightward (reduction of excitability) or leftward (increase of excitability) shifts, respectively, of the input-output function of mPFC L2/3 pyramidal cells without effects on the slope. Finally, we showed in electrophysiological recordings and epifluorescence Ca(2+)-imaging that GABA(B) receptor-mediated tonic inhibition is capable of regulating network activity. Blocking GABA(B) receptors increased the frequency of excitatory postsynaptic currents impinging on a neuron and prolonged network upstates. These results show that ambient GABA via GABA(B) receptors is powerful enough to modulate neuronal excitability and the activity of neural networks.

Estimating network parameters from combined dynamics of firing rate and irregularity of single neurons

High firing irregularity is a hallmark of cortical neurons in vivo, and modeling studies suggest a balance of excitation and inhibition is necessary to explain this high irregularity. Such a balance must be generated, at least partly, from local interconnected networks of excitatory and inhibitory neurons, but the details of the local network structure are largely unknown. The dynamics of the neural activity depends on the local network structure; this in turn suggests the possibility of estimating network structure from the dynamics of the firing statistics. Here we report a new method to estimate properties of the local cortical network from the instantaneous firing rate and irregularity (CV(2)) under the assumption that recorded neurons are a part of a randomly connected sparse network. The firing irregularity, measured in monkey motor cortex, exhibits two features; many neurons show relatively stable firing irregularity in time and across different task conditions; the time-averaged CV(2) is widely distributed from quasi-regular to irregular (CV(2) = 0.3-1.0). For each recorded neuron, we estimate the three parameters of a local network [balance of local excitation-inhibition, number of recurrent connections per neuron, and excitatory postsynaptic potential (EPSP) size] that best describe the dynamics of the measured firing rates and irregularities. Our analysis shows that optimal parameter sets form a two-dimensional manifold in the three-dimensional parameter space that is confined for most of the neurons to the inhibition-dominated region. High irregularity neurons tend to be more strongly connected to the local network, either in terms of larger EPSP and inhibitory PSP size or larger number of recurrent connections, compared with the low irregularity neurons, for a given excitatory/inhibitory balance. Incorporating either synaptic short-term depression or conductance-based synapses leads many low CV(2) neurons to move to the excitation-dominated region as well as to an increase of EPSP size.

Multiple timescale encoding of slowly varying whisker stimulus envelope in cortical and thalamic neurons in vivo

Adaptive processes over many timescales endow neurons with sensitivity to stimulus changes over a similarly wide range of scales. Although spike timing of single neurons can precisely signal rapid fluctuations in their inputs, the mean firing rate can convey information about slower-varying properties of the stimulus. Here, we investigate the firing rate response to a slowly varying envelope of whisker motion in two processing stages of the rat vibrissa pathway. The whiskers of anesthetized rats were moved through a noise trajectory with an amplitude that was sinusoidally modulated at one of several frequencies. In thalamic neurons, we found that the rate response to the stimulus envelope was also sinusoidal, with an approximately frequency-independent phase advance with respect to the input. Responses in cortex were similar but with a phase shift that was about three times larger, consistent with a larger amount of rate adaptation. These response properties can be described as a linear transformation of the input for which a single parameter quantifies the phase shift as well as the degree of adaptation. These results are reproduced by a model of adapting neurons connected by synapses with short-term plasticity, showing that the observed linear response and phase lead can be built up from a network that includes a sequence of nonlinear adapting elements. Our study elucidates how slowly varying envelope information under passive stimulation is preserved and transformed through the vibrissa processing pathway.

Frequency-dependent disynaptic inhibition in the pyramidal network: a ubiquitous pathway in the developing rat neocortex

The general structure of the mammalian neocortex is remarkably similar across different cortical areas. Despite certain cytoarchitectural specializations and deviations from the general blueprint, the principal organization of the neocortex is relatively uniform. It is not known, however, to what extent stereotypic synaptic pathways resemble each other between cortical areas, and how far they might reflect possible functional uniformity or specialization. Here, we show that frequency-dependent disynaptic inhibition (FDDI) is a generic circuit motif that is present in all neocortical areas we investigated (primary somatosensory, auditory and motor cortex, secondary visual cortex and medial prefrontal cortex of the developing rat). We did find, however, area-specific differences in occurrence and kinetics of FDDI and the short-term dynamics of monosynaptic connections between pyramidal cells (PCs). Connectivity between PCs, both monosynaptic and via FDDI, is higher in primary cortices. The long-term effectiveness of FDDI is likely to be limited by an activity-dependent attenuation of the PC-interneuron synaptic transmission. Our results suggest that the basic construction of neocortical synaptic pathways follows principles that are independent of modality or hierarchical order within the neocortex.

Feedback-induced gain control in stochastic spiking networks

The joint influence of recurrent feedback and noise on gain control in a network of globally coupled spiking leaky integrate-and-fire neurons is studied theoretically and numerically. The context of our work is the origin of divisive versus subtractive gain control, as mixtures of these effects are seen in a variety of experimental systems. We focus on changes in the slope of the mean firing frequency-versus-input bias (f-I) curve when the gain control signal to the cells comes from the cells' output spikes. Feedback spikes are modeled as alpha functions that produce an additive current in the current balance equation. For generality, they occur after a fixed minimum delay. We show that purely divisive gain control, i.e. changes in the slope of the f-I curve, arises naturally with this additive negative or positive feedback, due to a linearizing actions of feedback. Negative feedback alone lowers the gain, accounting in particular for gain changes in weakly electric fish upon pharmacological opening of the feedback loop as reported by Bastian (J Neurosci 6:553-562, 1986). When negative feedback is sufficiently strong it further causes oscillatory firing patterns which produce irregularities in the f-I curve. Small positive feedback alone increases the gain, but larger amounts cause abrupt jumps to higher firing frequencies. On the other hand, noise alone in open loop linearizes the f-I curve around threshold, and produces mixtures of divisive and subtractive gain control. With both noise and feedback, the combined gain control schemes produce a primarily divisive gain control shift, indicating the robustness of feedback gain control in stochastic networks. Similar results are found when the "input" parameter is the contrast of a time-varying signal rather than the bias current. Theoretical results are derived relating the slope of the f-I curve to feedback gain and noise strength. Good agreement with simulation results are found for inhibitory and excitatory feedback. Finally, divisive feedback is also found for conductance-based feedback (shunting or excitatory) with and without noise.

Divisive gain modulation with dynamic stimuli in integrate-and-fire neurons

The modulation of the sensitivity, or gain, of neural responses to input is an important component of neural computation. It has been shown that divisive gain modulation of neural responses can result from a stochastic shunting from balanced (mixed excitation and inhibition) background activity. This gain control scheme was developed and explored with static inputs, where the membrane and spike train statistics were stationary in time. However, input statistics, such as the firing rates of pre-synaptic neurons, are often dynamic, varying on timescales comparable to typical membrane time constants. Using a population density approach for integrate-and-fire neurons with dynamic and temporally rich inputs, we find that the same fluctuation-induced divisive gain modulation is operative for dynamic inputs driving nonequilibrium responses. Moreover, the degree of divisive scaling of the dynamic response is quantitatively the same as the steady-state responses—thus, gain modulation via balanced conductance fluctuations generalizes in a straight-forward way to a dynamic setting.

Many neural computations, including sensory and motor processing, require neurons to control their sensitivity (often termed ‘gain’) to stimuli. One common form of gain manipulation is divisive gain control, where the neural response to a specific stimulus is simply scaled by a constant. Most previous theoretical and experimental work on divisive gain control have assumed input statistics to be constant in time. However, realistic inputs can be highly time-varying, often with time-varying statistics, and divisive gain control remains to be extended to these cases. A widespread mechanism for divisive gain control for static inputs is through an increase in stimulus independent membrane fluctuations. We address the question of whether this divisive gain control scheme is indeed operative for time-varying inputs. Using simplified spiking neuron models, we employ accurate theoretical methods to estimate the dynamic neural response. We find that gain control via membrane fluctuations does indeed extend to the time-varying regime, and moreover, the degree of divisive scaling does not depend on the timescales of the driving input. This significantly increases the relevance of this form of divisive gain control for neural computations where input statistics change in time, as expected during normal sensory and motor behavior.

Sensitivity of firing rate to input fluctuations depends on time scale separation between fast and slow variables in single neurons

Neuronal responses are often characterized by the firing rate as a function of the stimulus mean, or the f-I curve. We introduce a novel classification of neurons into Types A, B-, and B+ according to how f-I curves are modulated by input fluctuations. In Type A neurons, the f-I curves display little sensitivity to input fluctuations when the mean current is large. In contrast, Type B neurons display sensitivity to fluctuations throughout the entire range of input means. Type B- neurons do not fire repetitively for any constant input, whereas Type B+ neurons do. We show that Type B+ behavior results from a separation of time scales between a slow and fast variable. A voltage-dependent time constant for the recovery variable can facilitate sensitivity to input fluctuations. Type B+ firing rates can be approximated using a simple "energy barrier" model.

Conditional bursting enhances resonant firing in neocortical layer 2-3 pyramidal neurons

The frequency response properties of neurons are critical for signal transmission and control of network oscillations. At subthreshold membrane potential, some neurons show resonance caused by voltage-gated channels. During action potential firing, resonance of the spike output may arise from subthreshold mechanisms and/or spike-dependent currents that cause afterhyperpolarizations (AHPs) and afterdepolarizations (ADPs). Layer 2-3 pyramidal neurons (L2-3 PNs) have a fast ADP that can trigger bursts. The present study investigated what stimuli elicit bursting in these cells and whether bursts transmit specific frequency components of the synaptic input, leading to resonance at particular frequencies. We found that two-spike bursts are triggered by step onsets, sine waves in two frequency bands, and noise. Using noise adjusted to elicit firing at approximately 10 Hz, we measured the gain for modulation of the time-varying firing rate as a function of stimulus frequency, finding a primary peak (7-16 Hz) and a high-frequency resonance (250-450 Hz). Gain was also measured separately for single and burst spikes. For a given spike rate, bursts provided higher gain at the primary peak and lower gain at intermediate frequencies, sharpening the high-frequency resonance. Suppression of bursting using automated current feedback weakened the primary and high-frequency resonances. The primary resonance was also influenced by the SK channel-mediated medium AHP (mAHP), because the SK blocker apamin reduced the sharpness of the primary peak. Our results suggest that resonance in L2-3 PNs depends on burst firing and the mAHP. Bursting enhances resonance in two distinct frequency bands.

The response of cortical neurons to in vivo-like input current: theory and experiment : I. Noisy inputs with stationary statistics

The study of several aspects of the collective dynamics of interacting neurons can be highly simplified if one assumes that the statistics of the synaptic input is the same for a large population of similarly behaving neurons (mean field approach). In particular, under such an assumption, it is possible to determine and study all the equilibrium points of the network dynamics when the neuronal response to noisy, in vivo-like, synaptic currents is known. The response function can be computed analytically for simple integrate-and-fire neuron models and it can be measured directly in experiments in vitro. Here we review theoretical and experimental results about the neural response to noisy inputs with stationary statistics. These response functions are important to characterize the collective neural dynamics that are proposed to be the neural substrate of working memory, decision making and other cognitive functions. Applications to the case of time-varying inputs are reviewed in a companion paper (Giugliano et al. in Biol Cybern, 2008). We conclude that modified integrate-and-fire neuron models are good enough to reproduce faithfully many of the relevant dynamical aspects of the neuronal response measured in experiments on real neurons in vitro.

Spike-train spectra and network response functions for non-linear integrate-and-fire neurons

Reduced models have long been used as a tool for the analysis of the complex activity taking place in neurons and their coupled networks. Recent advances in experimental and theoretical techniques have further demonstrated the usefulness of this approach. Despite the often gross simplification of the underlying biophysical properties, reduced models can still present significant difficulties in their analysis, with the majority of exact and perturbative results available only for the leaky integrate-and-fire model. Here an elementary numerical scheme is demonstrated which can be used to calculate a number of biologically important properties of the general class of non-linear integrate-and-fire models. Exact results for the first-passage-time density and spike-train spectrum are derived, as well as the linear response properties and emergent states of recurrent networks. Given that the exponential integrate-fire model has recently been shown to agree closely with the experimentally measured response of pyramidal cells, the methodology presented here promises to provide a convenient tool to facilitate the analysis of cortical-network dynamics.

The response of cortical neurons to in vivo-like input current: theory and experiment: II. Time-varying and spatially distributed inputs

The response of a population of neurons to time-varying synaptic inputs can show a rich phenomenology, hardly predictable from the dynamical properties of the membrane's inherent time constants. For example, a network of neurons in a state of spontaneous activity can respond significantly more rapidly than each single neuron taken individually. Under the assumption that the statistics of the synaptic input is the same for a population of similarly behaving neurons (mean field approximation), it is possible to greatly simplify the study of neural circuits, both in the case in which the statistics of the input are stationary (reviewed in La Camera et al. in Biol Cybern, 2008) and in the case in which they are time varying and unevenly distributed over the dendritic tree. Here, we review theoretical and experimental results on the single-neuron properties that are relevant for the dynamical collective behavior of a population of neurons. We focus on the response of integrate-and-fire neurons and real cortical neurons to long-lasting, noisy, in vivo-like stationary inputs and show how the theory can predict the observed rhythmic activity of cultures of neurons. We then show how cortical neurons adapt on multiple time scales in response to input with stationary statistics in vitro. Next, we review how it is possible to study the general response properties of a neural circuit to time-varying inputs by estimating the response of single neurons to noisy sinusoidal currents. Finally, we address the dendrite-soma interactions in cortical neurons leading to gain modulation and spike bursts, and show how these effects can be captured by a two-compartment integrate-and-fire neuron. Most of the experimental results reviewed in this article have been successfully reproduced by simple integrate-and-fire model neurons.

Intrinsic gain modulation and adaptive neural coding

In many cases, the computation of a neural system can be reduced to a receptive field, or a set of linear filters, and a thresholding function, or gain curve, which determines the firing probability; this is known as a linear/nonlinear model. In some forms of sensory adaptation, these linear filters and gain curve adjust very rapidly to changes in the variance of a randomly varying driving input. An apparently similar but previously unrelated issue is the observation of gain control by background noise in cortical neurons: the slope of the firing rate versus current (f-I) curve changes with the variance of background random input. Here, we show a direct correspondence between these two observations by relating variance-dependent changes in the gain of f-I curves to characteristics of the changing empirical linear/nonlinear model obtained by sampling. In the case that the underlying system is fixed, we derive relationships relating the change of the gain with respect to both mean and variance with the receptive fields derived from reverse correlation on a white noise stimulus. Using two conductance-based model neurons that display distinct gain modulation properties through a simple change in parameters, we show that coding properties of both these models quantitatively satisfy the predicted relationships. Our results describe how both variance-dependent gain modulation and adaptive neural computation result from intrinsic nonlinearity.

Many neurons are known to achieve a wide dynamic range by adaptively changing their computational input/output function according to the input statistics. These adaptive changes can be very rapid, and it has been suggested that a component of this adaptation could be purely input-driven: even a fixed neural system can show apparent adaptive behavior since inputs with different statistics interact with the nonlinearity of the system in different ways. In this paper, we show how a single neuron's intrinsic computational function can dictate such input-driven changes in its response to varying input statistics, which begets a relationship between two different characterizations of neural function—in terms of mean firing rate and in terms of generating precise spike timing. We then apply our results to two biophysically defined model neurons, which have significantly different response patterns to inputs with various statistics. Our model of intrinsic adaptation explains their behaviors well. Contrary to the picture that neurons carry out a stereotyped computation on their inputs, our results show that even in the simplest cases they have simple yet effective mechanisms by which they can adapt to their input. Adaptation to stimulus statistics, therefore, is built into the most basic single neuron computations.

Two computational regimes of a single-compartment neuron separated by a planar boundary in conductance space

Recent in vitro data show that neurons respond to input variance with varying sensitivities. Here we demonstrate that Hodgkin-Huxley (HH) neurons can operate in two computational regimes: one that is more sensitive to input variance (differentiating) and one that is less sensitive (integrating). A boundary plane in the 3D conductance space separates these two regimes. For a reduced HH model, this plane can be derived analytically from the V nullcline, thus suggesting a means of relating biophysical parameters to neural computation by analyzing the neuron's dynamical system.

Dopamine Increases the Gain of the Input-Output Response of Rat Prefrontal Pyramidal Neurons

Dopaminergic modulation of prefrontal cortical activity is known to affect cognitive functions like working memory. Little consensus on the role of dopamine modulation has been achieved, however, in part because quantities directly relating to the neuronal substrate of working memory are difficult to measure. Here we show that dopamine increases the gain of the frequency-current relationship of layer 5 pyramidal neurons in vitro in response to noisy input currents. The gain increase could be attributed to a reduction of the slow afterhyperpolarization by dopamine. Dopamine also increases neuronal excitability by shifting the input-output functions to lower inputs. The modulation of these response properties is mainly mediated by D1 receptors. Integrate-and-fire neurons were fitted to the experimentally recorded input-output functions and recurrently connected in a model network. The gain increase induced by dopamine application facilitated and stabilized persistent activity in this network. The results support the hypothesis that catecholamines increase the neuronal gain and suggest that dopamine improves working memory via gain modulation.