Differentially synchronized spiking enables multiplexed neural coding

Freezing of movements and its correspondence with MLG1 neuron response to looming stimuli in the crab Neohelice

Upon visually detecting a moving predator, animals often freeze, i.e. stop moving, to minimize being uncovered and to gather detailed information of the object's movements and properties. In certain conditions, the freezing behavior can be enough to avoid a predatory menace but, when the risk is high or increases to a higher level, animals switch strategy and engage in an escape response. The neural bases underlying escape responses to visual stimuli have been extensively investigated both in vertebrates and arthropods. However, those involved in freezing behaviors are much less studied. Here, we investigated the freezing behavior displayed by the crab Neohelice granulata when confronted with a variety of looming stimuli simulating objects of distinct sizes approaching on a collision course at different speeds. The experiments were performed in a treadmill-like device. Animals engaged in exploratory walks responded to the looming stimulus with freezing followed by escaping. The analysis of the stimulus optical variables shows that regardless of the looming dynamic, the freezing decision is made when the angular size of the object increases by 1.4 deg. In vivo intracellular recording responses of monostratified lobula giant neurons (MLG1) to the same looming stimuli show that the freezing times correlate with the times predicted by a hypothetical spike counter of this neuron.

Encoding of Vibrotactile Stimuli by Mechanoreceptors in Rodent Glabrous Skin

Somatosensory coding in rodents has been mostly studied in the whisker system and hairy skin, whereas the function of low-threshold mechanoreceptors (LTMRs) in the rodent glabrous skin has received scant attention, unlike in primates where the glabrous skin has been the focus. The relative activation of different LTMR subtypes carries information about vibrotactile stimuli, as does the rate and temporal patterning of LTMR spikes. Rate coding depends on the probability of a spike occurring on each stimulus cycle (reliability), whereas temporal coding depends on the timing of spikes relative to the stimulus cycle (precision). Using in vivo extracellular recordings in male rats and mice of either sex, we measured the reliability and precision of LTMR responses to tactile stimuli including sustained pressure and vibration. Similar to other species, rodent LTMRs were separated into rapid-adapting (RA) or slow-adapting based on their response to sustained pressure. However, unlike the dichotomous frequency preference characteristic of RA1 and RA2/Pacinian afferents in other species, rodent RAs fell along a continuum. Fitting generalized linear models to experimental data reproduced the reliability and precision of rodent RAs. The resulting model parameters highlight key mechanistic differences across the RA spectrum; specifically, the integration window of different RAs transitions from wide to narrow as tuning preferences across the population move from low to high frequencies. Our results show that rodent RAs can support both rate and temporal coding, but their heterogeneity suggests that coactivation patterns play a greater role in population coding than for dichotomously tuned primate RAs.

Multimodal neuroimaging of hierarchical cognitive control

Cognitive control enables us to translate our knowledge into actions, allowing us to flexibly adjust our behavior, according to environmental contexts, our internal goals, and future plans. Multimodal neuroimaging and neurostimulation techniques have proven essential for advancing our understanding of how cognitive control emerges from the coordination of distributed neuronal activities in the brain. In this review, we examine the literature on multimodal studies of cognitive control. We explore how these studies provide converging evidence for a novel, multiplexed model of cognitive control, in which neural oscillations support different levels of control processing along a functionally hierarchical organization of distinct frontoparietal networks.

Hippocampal synchronization in a realistic CA1 neuron model

This work delves into studying the synchronization in two realistic neuron models using Hodgkin-Huxley dynamics. Unlike simplistic pointlike models, excitatory synapses are here randomly distributed along the dendrites, introducing strong stochastic contributions into their signal propagation. To focus on the role of different synaptic locations, we use two copies of the same neuron whose synapses are located at different distances from the soma and activated by identical Poisson distributed pulses. The synchronization is investigated through a specifically defined spiking correlation function, and its behavior is analyzed as a function of several parameters: inhibition weight, distance from the soma of one synaptic group, weight and decay time constants of the excitatory synapses.

Disentangling Temporal and Rate Codes in the Primate Somatosensory Cortex

Millisecond-scale temporal spiking patterns encode sensory information in the periphery, but their role in the neocortex remains controversial. The sense of touch provides a window into temporal coding because tactile neurons often exhibit precise, repeatable, and informative temporal spiking patterns. In the somatosensory cortex (S1), responses to skin vibrations exhibit phase locking that faithfully carries information about vibratory frequency. However, the respective roles of spike timing and rate in frequency coding are confounded because vibratory frequency shapes both the timing and rates of responses. To disentangle the contributions of these two neural features, we measured S1 responses as rhesus macaques performed frequency discrimination tasks in which differences in frequency were accompanied by orthogonal variations in amplitude. We assessed the degree to which the strength and timing of responses could account for animal performance. First, we showed that animals can discriminate frequency, but their performance is biased by amplitude variations. Second, rate-based representations of frequency are susceptible to changes in amplitude but in ways that are inconsistent with the animals' behavioral biases, calling into question a rate-based neural code for frequency. In contrast, timing-based representations are highly informative about frequency but impervious to changes in amplitude, which is also inconsistent with the animals' behavior. We account for the animals' behavior with a model wherein frequency coding relies on a temporal code, but frequency judgments are biased by perceived magnitude. We conclude that information about vibratory frequency is not encoded in S1 firing rates but primarily in temporal patterning on millisecond timescales.

Clinico-physiological correlates of Parkinson's disease from multi-resolution basal ganglia recordings

Parkinson's disease (PD) has been associated with pathological neural activity within the basal ganglia. Herein, we analyzed resting-state single-neuron and local field potential (LFP) activities from people with PD who underwent awake deep brain stimulation surgery of the subthalamic nucleus (STN; n = 125) or globus pallidus internus (GPi; n = 44), and correlated rate-based and oscillatory features with UPDRSIII off-medication subscores. Rate-based single-neuron features did not correlate with PD symptoms. STN single-neuron and LFP low-beta (12-21 Hz) power and burst dynamics showed modest correlations with bradykinesia and rigidity severity, while STN spiketrain theta (4-8 Hz) power correlated modestly with tremor severity. GPi low- and high-beta (21-30 Hz) power and burst dynamics correlated moderately with bradykinesia and axial symptom severity. These findings suggest that elevated single-neuron and LFP oscillations may be linked to symptoms, though modest correlations imply that the pathophysiology of PD may extend beyond resting-state beta oscillations.

Multiplexing intensity and frequency sensations for artificial touch by modulating temporal features of electrical pulse trains

In neural prostheses, intensity modulation of a single channel (i.e., through a single stimulating electrode) has been achieved by increasing the magnitude or width of each stimulation pulse, which risks eliciting pain or paraesthesia; and by changing the stimulation rate, which leads to concurrent changes in perceived frequency. In this study, we sought to render a perception of tactile intensity and frequency independently, by means of temporal pulse train patterns of fixed magnitude, delivered non-invasively. Our psychophysical study exploits a previously discovered frequency coding mechanism, where the perceived frequency of stimulus pulses grouped into periodic bursts depends on the duration of the inter-burst interval, rather than the mean pulse rate or periodicity. When electrical stimulus pulses were organised into bursts, perceived intensity was influenced by the number of pulses within a burst, while perceived frequency was determined by the time between the end of one burst envelope and the start of the next. The perceived amplitude was modulated by 1.6× while perceived frequency was varied independently by 2× within the tested range (20-40 Hz). Thus, the sensation of intensity might be controlled independently from frequency through a single stimulation channel without having to vary the injected electrical current. This can form the basis for improving strategies in delivering more complex and natural sensations for prosthetic hand users.

Basic Properties of Coordinated Neuronal Ensembles in the Auditory Thalamus

Coordinated neuronal activity has been identified to play an important role in information processing and transmission in the brain. However, current research predominantly focuses on understanding the properties and functions of neuronal coordination in hippocampal and cortical areas, leaving subcortical regions relatively unexplored. In this study, we use single-unit recordings in female Sprague Dawley rats to investigate the properties and functions of groups of neurons exhibiting coordinated activity in the auditory thalamus-the medial geniculate body (MGB). We reliably identify coordinated neuronal ensembles (cNEs), which are groups of neurons that fire synchronously, in the MGB. cNEs are shown not to be the result of false-positive detections or by-products of slow-state oscillations in anesthetized animals. We demonstrate that cNEs in the MGB have enhanced information-encoding properties over individual neurons. Their neuronal composition is stable between spontaneous and evoked activity, suggesting limited stimulus-induced ensemble dynamics. These MGB cNE properties are similar to what is observed in cNEs in the primary auditory cortex (A1), suggesting that ensembles serve as a ubiquitous mechanism for organizing local networks and play a fundamental role in sensory processing within the brain.

Absence of paresthesia during high-rate spinal cord stimulation reveals importance of synchrony for sensations evoked by electrical stimulation

Electrically activating mechanoreceptive afferents inhibits pain. However, paresthesia evoked by spinal cord stimulation (SCS) at 40-60 Hz becomes uncomfortable at high pulse amplitudes, limiting SCS "dosage." Kilohertz-frequency SCS produces analgesia without paresthesia and is thought, therefore, not to activate afferent axons. We show that paresthesia is absent not because axons do not spike but because they spike asynchronously. In a pain patient, selectively increasing SCS frequency abolished paresthesia and epidurally recorded evoked compound action potentials (ECAPs). Dependence of ECAP amplitude on SCS frequency was reproduced in pigs, rats, and computer simulations and is explained by overdrive desynchronization: spikes desychronize when axons are stimulated faster than their refractory period. Unlike synchronous spikes, asynchronous spikes fail to produce paresthesia because their transmission to somatosensory cortex is blocked by feedforward inhibition. Our results demonstrate how stimulation frequency impacts synchrony based on axon properties and how synchrony impacts sensation based on circuit properties.

Robust and consistent measures of pattern separation based on information theory and demonstrated in the dentate gyrus

Pattern separation is a valuable computational function performed by neuronal circuits, such as the dentate gyrus, where dissimilarity between inputs is increased, reducing noise and increasing the storage capacity of downstream networks. Pattern separation is studied from both in vivo experimental and computational perspectives and, a number of different measures (such as orthogonalisation, decorrelation, or spike train distance) have been applied to quantify the process of pattern separation. However, these are known to give conclusions that can differ qualitatively depending on the choice of measure and the parameters used to calculate it. We here demonstrate that arbitrarily increasing sparsity, a noticeable feature of dentate granule cell firing and one that is believed to be key to pattern separation, typically leads to improved classical measures for pattern separation even, inappropriately, up to the point where almost all information about the inputs is lost. Standard measures therefore both cannot differentiate between pattern separation and pattern destruction, and give results that may depend on arbitrary parameter choices. We propose that techniques from information theory, in particular mutual information, transfer entropy, and redundancy, should be applied to penalise the potential for lost information (often due to increased sparsity) that is neglected by existing measures. We compare five commonly-used measures of pattern separation with three novel techniques based on information theory, showing that the latter can be applied in a principled way and provide a robust and reliable measure for comparing the pattern separation performance of different neurons and networks. We demonstrate our new measures on detailed compartmental models of individual dentate granule cells and a dentate microcircuit, and show how structural changes associated with epilepsy affect pattern separation performance. We also demonstrate how our measures of pattern separation can predict pattern completion accuracy. Overall, our measures solve a widely acknowledged problem in assessing the pattern separation of neural circuits such as the dentate gyrus, as well as the cerebellum and mushroom body. Finally we provide a publicly available toolbox allowing for easy analysis of pattern separation in spike train ensembles.

Global and local neuronal coding of tactile information in the barrel cortex

During tactile sensation in rodents, the whisker movements across surfaces give rise to intricate whisker motions that encompass discrete and transient stick-slip events, effectively conveying valuable information regarding surface properties. These surface characteristics are transformed into cortical neuronal responses. This study examined the coding strategies underlying these transformations in rat whiskers. We found that changes in surface coarseness modified the number and magnitude of stick-slip events, which in turn both modulated properties of neuronal responses. Global changes in the number of stick-slip events primarily affected neuronal discharge rates and the degree of neuronal synchronization. In contrast, local changes in the magnitude of stick-slip events affected the transformation of these kinematic and kinetic characteristics into neuronal discharges. Most cortical neurons exhibited surface coarseness selectivity through global and local stick-slip event properties. However, this selectivity varied across coding strategies in the same neurons, given that each coding strategy reflected different aspects of changes in whisker-surface interactions. The degree of spatial similarity in surface coarseness preference in adjacently recorded neurons differed among these coding strategies. Adjacently recorded neurons exhibited the same surface coarseness preference in their firing rates but not through other coding strategies. Through these results, we were able to show that local stick-slip event properties contribute to texture discrimination, complementing and surpassing global coding in this context. These findings suggest that the representation of surface coarseness in the cortex may rely on concurrent coding strategies that integrate tactile information across different spatiotemporal scales.

Inter-areal transmission of multiple neural signals through frequency-division-multiplexing communication

Inter-areal information transmission in the brain cortex relates to cognitive functions. Researches used to pay attention to activity pattern transmission, signals gating, or routing in neuronal networks. However, the underlying mechanism of simultaneous transmission of multiple neural signals in the same channel across networks remains unclear. In this work, we construct a two-layer feedforward neuronal network (sender-receiver) with each layer's intrinsic rhythms consisting of slow- (low-frequency) and fast- gamma rhythms (high-frequency), investigating how to realize simultaneous transmission of multiple signals in neuronal systems. With the aid of resonance and frequency analysis, it is shown that low- and high-frequency signals can be transmitted simultaneously in such a feedforward network through frequency division multiplexing (FDM) communication. The transmission performance is related to the local resonance, connectivity, as well as background noise. Moreover, low- and high-frequency signals can also be gated or selected with appropriate adjustments of recurrent connection strength and delay, and background noise. Our model might provide a novel insight into the underlying mechanism of complex signals communication between different cortex areas.

Distinctive properties of biological neural networks and recent advances in bottom-up approaches toward a better biologically plausible neural network

Although it may appear infeasible and impractical, building artificial intelligence (AI) using a bottom-up approach based on the understanding of neuroscience is straightforward. The lack of a generalized governing principle for biological neural networks (BNNs) forces us to address this problem by converting piecemeal information on the diverse features of neurons, synapses, and neural circuits into AI. In this review, we described recent attempts to build a biologically plausible neural network by following neuroscientifically similar strategies of neural network optimization or by implanting the outcome of the optimization, such as the properties of single computational units and the characteristics of the network architecture. In addition, we proposed a formalism of the relationship between the set of objectives that neural networks attempt to achieve, and neural network classes categorized by how closely their architectural features resemble those of BNN. This formalism is expected to define the potential roles of top-down and bottom-up approaches for building a biologically plausible neural network and offer a map helping the navigation of the gap between neuroscience and AI engineering.

Synchrony-Division Neural Multiplexing: An Encoding Model

Cortical neurons receive mixed information from the collective spiking activities of primary sensory neurons in response to a sensory stimulus. A recent study demonstrated an abrupt increase or decrease in stimulus intensity and the stimulus intensity itself can be respectively represented by the synchronous and asynchronous spikes of S1 neurons in rats. This evidence capitalized on the ability of an ensemble of homogeneous neurons to multiplex, a coding strategy that was referred to as synchrony-division multiplexing (SDM). Although neural multiplexing can be conceived by distinct functions of individual neurons in a heterogeneous neural ensemble, the extent to which nearly identical neurons in a homogeneous neural ensemble encode multiple features of a mixed stimulus remains unknown. Here, we present a computational framework to provide a system-level understanding on how an ensemble of homogeneous neurons enable SDM. First, we simulate SDM with an ensemble of homogeneous conductance-based model neurons receiving a mixed stimulus comprising slow and fast features. Using feature-estimation techniques, we show that both features of the stimulus can be inferred from the generated spikes. Second, we utilize linear nonlinear (LNL) cascade models and calculate temporal filters and static nonlinearities of differentially synchronized spikes. We demonstrate that these filters and nonlinearities are distinct for synchronous and asynchronous spikes. Finally, we develop an augmented LNL cascade model as an encoding model for the SDM by combining individual LNLs calculated for each type of spike. The augmented LNL model reveals that a homogeneous neural ensemble model can perform two different functions, namely, temporal- and rate-coding, simultaneously.

Coexisting neuronal coding strategies in the barrel cortex

During tactile sensation by rodents, whisker movements across surfaces generate complex whisker motions, including discrete, transient stick-slip events, which carry information about surface properties. The characteristics of these events and how the brain encodes this tactile information remain enigmatic. We found that cortical neurons show a mixture of synchronized and nontemporally correlated spikes in their tactile responses. Synchronous spikes convey the magnitude of stick-slip events by numerous aspects of temporal coding. These spikes show preferential selectivity for kinetic and kinematic whisker motion. By contrast, asynchronous spikes in each neuron convey the magnitude of stick-slip events by their discharge rates, response probability, and interspike intervals. We further show that the differentiation between these two types of activity is highly dependent on the magnitude of stick-slip events and stimulus and response history. These results suggest that cortical neurons transmit multiple components of tactile information through numerous coding strategies.

Precise and stable edge orientation signaling by human first-order tactile neurons

Fast-adapting type 1 (FA-1) and slow-adapting type 1 (SA-1) first-order neurons in the human tactile system have distal axons that branch in the skin and form many transduction sites, yielding receptive fields with many highly sensitive zones or 'subfields.' We previously demonstrated that this arrangement allows FA-1 and SA-1 neurons to signal the geometric features of touched objects, specifically the orientation of raised edges scanned with the fingertips. Here, we show that such signaling operates for fine edge orientation differences (5-20°) and is stable across a broad range of scanning speeds (15-180 mm/s); that is, under conditions relevant for real-world hand use. We found that both FA-1 and SA-1 neurons weakly signal fine edge orientation differences via the intensity of their spiking responses and only when considering a single scanning speed. Both neuron types showed much stronger edge orientation signaling in the sequential structure of the evoked spike trains, and FA-1 neurons performed better than SA-1 neurons. Represented in the spatial domain, the sequential structure was strikingly invariant across scanning speeds, especially those naturally used in tactile spatial discrimination tasks. This speed invariance suggests that neurons' responses are structured via sequential stimulation of their subfields and thus links this capacity to their terminal organization in the skin. Indeed, the spatial precision of elicited action potentials rationally matched spatial acuity of subfield arrangements, which corresponds to a spatial period similar to the dimensions of individual fingertip ridges.

Physiological noise facilitates multiplexed coding of vibrotactile-like signals in somatosensory cortex

Neurons can use different aspects of their spiking to simultaneously represent (multiplex) different features of a stimulus. For example, some pyramidal neurons in primary somatosensory cortex (S1) use the rate and timing of their spikes to, respectively, encode the intensity and frequency of vibrotactile stimuli. Doing so has several requirements. Because they fire at low rates, pyramidal neurons cannot entrain 1:1 with high-frequency (100 to 600 Hz) inputs and, instead, must skip (i.e., not respond to) some stimulus cycles. The proportion of skipped cycles must vary inversely with stimulus intensity for firing rate to encode stimulus intensity. Spikes must phase-lock to the stimulus for spike times (intervals) to encode stimulus frequency, but, in addition, skipping must occur irregularly to avoid aliasing. Using simulations and in vitro experiments in which mouse S1 pyramidal neurons were stimulated with inputs emulating those induced by vibrotactile stimuli, we show that fewer cycles are skipped as stimulus intensity increases, as required for rate coding, and that intrinsic or synaptic noise can induce irregular skipping without disrupting phase locking, as required for temporal coding. This occurs because noise can modulate the reliability without disrupting the precision of spikes evoked by small-amplitude, fast-onset signals. Specifically, in the fluctuation-driven regime associated with sparse spiking, rate and temporal coding are both paradoxically improved by the strong synaptic noise characteristic of the intact cortex. Our results demonstrate that multiplexed coding by S1 pyramidal neurons is not only feasible under in vivo conditions, but that background synaptic noise is actually beneficial.

Single-neuron bursts encode pathological oscillations in subcortical nuclei of patients with Parkinson's disease and essential tremor

Deep brain stimulation procedures offer an invaluable opportunity to study disease through intracranial recordings from awake patients. Here, we address the relationship between single-neuron and aggregate-level (local field potential; LFP) activities in the subthalamic nucleus (STN) and thalamic ventral intermediate nucleus (Vim) of patients with Parkinson's disease (n = 19) and essential tremor (n = 16), respectively. Both disorders have been characterized by pathologically elevated LFP oscillations, as well as an increased tendency for neuronal bursting. Our findings suggest that periodic single-neuron bursts encode both pathophysiological beta (13 to 33 Hz; STN) and tremor (4 to 10 Hz; Vim) LFP oscillations, evidenced by strong time-frequency and phase-coupling relationships between the bursting and LFP signals. Spiking activity occurring outside of bursts had no relationship to the LFP. In STN, bursting activity most commonly preceded the LFP oscillation, suggesting that neuronal bursting generated within STN may give rise to an aggregate-level LFP oscillation. In Vim, LFP oscillations most commonly preceded bursting activity, suggesting that neuronal firing may be entrained by periodic afferent inputs. In both STN and Vim, the phase-coupling relationship between LFP and high-frequency oscillation (HFO) signals closely resembled the relationships between the LFP and single-neuron bursting. This suggests that periodic single-neuron bursting is likely representative of a higher spatial and temporal resolution readout of periodic increases in the amplitude of HFOs, which themselves may be a higher resolution readout of aggregate-level LFP oscillations. Overall, our results may reconcile "rate" and "oscillation" models of Parkinson's disease and shed light on the single-neuron basis and origin of pathophysiological oscillations in movement disorders.

Are Grid-Like Representations a Component of All Perception and Cognition?

Grid cells or grid-like responses have been reported in the rodent, bat and human brains during various spatial and non-spatial tasks. However, the functions of grid-like representations beyond the classical hippocampal formation remain elusive. Based on accumulating evidence from recent rodent recordings and human fMRI data, we make speculative accounts regarding the mechanisms and functional significance of the sensory cortical grid cells and further make theory-driven predictions. We argue and reason the rationale why grid responses may be universal in the brain for a wide range of perceptual and cognitive tasks that involve locomotion and mental navigation. Computational modeling may provide an alternative and complementary means to investigate the grid code or grid-like map. We hope that the new discussion will lead to experimentally testable hypotheses and drive future experimental data collection.

Synchronous spiking of cerebellar Purkinje cells during control of movements

SignificanceThe information that one region of the brain transmits to another is usually viewed through the lens of firing rates. However, if the output neurons could vary the timing of their spikes, then, through synchronization, they would spotlight information that may be critical for control of behavior. Here we report that, in the cerebellum, Purkinje cell populations that share a preference for error convey, to the nucleus, when to decelerate the movement, by reducing their firing rates and temporally synchronizing the remaining spikes.

Multiple gamma rhythms carry distinct spatial frequency information in primary visual cortex

Gamma rhythms in many brain regions, including the primary visual cortex (V1), are thought to play a role in information processing. Here, we report a surprising finding of 3 narrowband gamma rhythms in V1 that processed distinct spatial frequency (SF) signals and had different neural origins. The low gamma (LG; 25 to 40 Hz) rhythm was generated at the V1 superficial layer and preferred a higher SF compared with spike activity, whereas both the medium gamma (MG; 40 to 65 Hz), generated at the cortical level, and the high gamma HG; (65 to 85 Hz), originated precortically, preferred lower SF information. Furthermore, compared with the rates of spike activity, the powers of the 3 gammas had better performance in discriminating the edge and surface of simple objects. These findings suggest that gamma rhythms reflect the neural dynamics of neural circuitries that process different SF information in the visual system, which may be crucial for multiplexing SF information and synchronizing different features of an object.

Corticothalamic Pathways From Layer 5: Emerging Roles in Computation and Pathology

Large portions of the thalamus receive strong driving input from cortical layer 5 (L5) neurons but the role of this important pathway in cortical and thalamic computations is not well understood. L5-recipient "higher-order" thalamic regions participate in cortico-thalamo-cortical (CTC) circuits that are increasingly recognized to be (1) anatomically and functionally distinct from better-studied "first-order" CTC networks, and (2) integral to cortical activity related to learning and perception. Additionally, studies are beginning to elucidate the clinical relevance of these networks, as dysfunction across these pathways have been implicated in several pathological states. In this review, we highlight recent advances in understanding L5 CTC networks across sensory modalities and brain regions, particularly studies leveraging cell-type-specific tools that allow precise experimental access to L5 CTC circuits. We aim to provide a focused and accessible summary of the anatomical, physiological, and computational properties of L5-originating CTC networks, and outline their underappreciated contribution in pathology. We particularly seek to connect single-neuron and synaptic properties to network (dys)function and emerging theories of cortical computation, and highlight information processing in L5 CTC networks as a promising focus for computational studies.

Information diversity in individual auditory cortical neurons is associated with functionally distinct coordinated neuronal ensembles

Neuronal activity in auditory cortex is often highly synchronous between neighboring neurons. Such coordinated activity is thought to be crucial for information processing. We determined the functional properties of coordinated neuronal ensembles (cNEs) within primary auditory cortical (AI) columns relative to the contributing neurons. Nearly half of AI cNEs showed robust spectro-temporal receptive fields whereas the remaining cNEs showed little or no acoustic feature selectivity. cNEs can therefore capture either specific, time-locked information of spectro-temporal stimulus features or reflect stimulus-unspecific, less-time specific processing aspects. By contrast, we show that individual neurons can represent both of those aspects through membership in multiple cNEs with either high or absent feature selectivity. These associations produce functionally heterogeneous spikes identifiable by instantaneous association with different cNEs. This demonstrates that single neuron spike trains can sequentially convey multiple aspects that contribute to cortical processing, including stimulus-specific and unspecific information.

Orientation processing by synaptic integration across first-order tactile neurons

Our ability to manipulate objects relies on tactile inputs from first-order tactile neurons that innervate the glabrous skin of the hand. The distal axon of these neurons branches in the skin and innervates many mechanoreceptors, yielding spatially-complex receptive fields. Here we show that synaptic integration across the complex signals from the first-order neuronal population could underlie human ability to accurately (< 3°) and rapidly process the orientation of edges moving across the fingertip. We first derive spiking models of human first-order tactile neurons that fit and predict responses to moving edges with high accuracy. We then use the model neurons in simulating the peripheral neuronal population that innervates a fingertip. We train classifiers performing synaptic integration across the neuronal population activity, and show that synaptic integration across first-order neurons can process edge orientations with high acuity and speed. In particular, our models suggest that integration of fast-decaying (AMPA-like) synaptic inputs within short timescales is critical for discriminating fine orientations, whereas integration of slow-decaying (NMDA-like) synaptic inputs supports discrimination of coarser orientations and maintains robustness over longer timescales. Taken together, our results provide new insight into the computations occurring in the earliest stages of the human tactile processing pathway and how they may be critical for supporting hand function.

Development of Natural Scene Representation in Primary Visual Cortex Requires Early Postnatal Experience

The development of the visual system is known to be shaped by early-life experience. To identify response properties that contribute to enhanced natural scene representation, we performed calcium imaging of excitatory neurons in the primary visual cortex (V1) of awake mice raised in three different conditions (standard-reared, dark-reared, and delayed-visual experience) and compared neuronal responses to natural scene features in relation to simpler grating stimuli that varied in orientation and spatial frequency. We assessed population selectivity in the V1 by using decoding methods and found that natural scene discriminability increased by 75% between the ages of 4 and 6 weeks. Both natural scene and grating discriminability were higher in standard-reared animals than in those raised in the dark. This increase in discriminability was accompanied by a reduction in the number of neurons that responded to low-spatial-frequency gratings. At the same time, there was an increase in neuronal preference for natural scenes. Light exposure restricted to a 2- to 4-week window during adulthood did not induce improvements in natural scene or in grating stimulus discriminability. Our results demonstrate that experience reduces the number of neurons needed to effectively encode grating stimuli and that early visual experience enhances natural scene discriminability by directly increasing responsiveness to natural scene features.

Context-Dependent Multiplexing by Individual VTA Dopamine Neurons

Dopamine (DA) neurons of the VTA track cues and rewards to generate a reward prediction error signal during Pavlovian conditioning. Here we explored how these neurons respond to a self-paced, operant task in freely moving mice. The animal could trigger a reward-predicting cue by remaining in a specific location of an operant box for a brief time before moving to a spout for reward collection. VTA DA neurons were identified using DAT-Cre male mice that carried an optrode with minimal impact on the behavioral task. In vivo single-unit recordings revealed transient fast spiking responses to the cue and reward in correct trials, while for incorrect ones the activity paused, reflecting positive and negative error signals of a reward prediction. In parallel, a majority of VTA DA neurons simultaneously encoded multiple actions (e.g., movement velocity, acceleration, distance to goal, and licking) in sustained slow firing modulation. Applying a GLM, we show that such multiplexed encoding of rewarding and motor variables by individual DA neurons was only apparent while the mouse was engaged in the task. Downstream targets may exploit such goal-directed multiplexing of VTA DA neurons to adjust actions to optimize the task's outcome.SIGNIFICANCE STATEMENT VTA DA neurons code for multiple functions, including the reward prediction error but also motivation and locomotion. Here we show that about half of the recorded VTA DA neurons perform multiplexing: they exploit the phasic and tonic activity modes to encode, respectively, the cue/reward responses and motor parameters, most prominently when the mouse engages in a self-paced operand task. VTA non-DA neurons, by contrast, encode motor parameters regardless of task engagement.

A Time-Varying Information Measure for Tracking Dynamics of Neural Codes in a Neural Ensemble

The amount of information that differentially correlated spikes in a neural ensemble carry is not the same; the information of different types of spikes is associated with different features of the stimulus. By calculating a neural ensemble’s information in response to a mixed stimulus comprising slow and fast signals, we show that the entropy of synchronous and asynchronous spikes are different, and their probability distributions are distinctively separable. We further show that these spikes carry a different amount of information. We propose a time-varying entropy (TVE) measure to track the dynamics of a neural code in an ensemble of neurons at each time bin. By applying the TVE to a multiplexed code, we show that synchronous and asynchronous spikes carry information in different time scales. Finally, a decoder based on the Kalman filtering approach is developed to reconstruct the stimulus from the spikes. We demonstrate that slow and fast features of the stimulus can be entirely reconstructed when this decoder is applied to asynchronous and synchronous spikes, respectively. The significance of this work is that the TVE can identify different types of information (for example, corresponding to synchronous and asynchronous spikes) that might simultaneously exist in a neural code.

Necessary Conditions for Reliable Propagation of Slowly Time-Varying Firing Rate

Reliable propagation of slow-modulations of the firing rate across multiple layers of a feedforward network (FFN) has proven difficult to capture in spiking neural models. In this paper, we explore necessary conditions for reliable and stable propagation of time-varying asynchronous spikes whose instantaneous rate of changes-in fairly short time windows [20-100] msec-represents information of slow fluctuations of the stimulus. Specifically, we study the effect of network size, level of background synaptic noise, and the variability of synaptic delays in an FFN with all-to-all connectivity. We show that network size and the level of background synaptic noise, together with the strength of synapses, are substantial factors enabling the propagation of asynchronous spikes in deep layers of an FFN. In contrast, the variability of synaptic delays has a minor effect on signal propagation.

Multiplexing rhythmic information by spike timing dependent plasticity

Rhythmic activity has been associated with a wide range of cognitive processes including the encoding of sensory information, navigation, the transfer of information and others. Rhythmic activity in the brain has also been suggested to be used for multiplexing information. Multiplexing is the ability to transmit more than one signal via the same channel. Here we focus on frequency division multiplexing, in which different signals are transmitted in different frequency bands. Recent work showed that spike-timing-dependent plasticity (STDP) can facilitate the transfer of rhythmic activity downstream the information processing pathway. However, STDP has also been known to generate strong winner-take-all like competition between subgroups of correlated synaptic inputs. This competition between different rhythmicity channels, induced by STDP, may prevent the multiplexing of information. Thus, raising doubts whether STDP is consistent with the idea of multiplexing. This study explores whether STDP can facilitate the multiplexing of information across multiple frequency channels, and if so, under what conditions. We address this question in a modelling study, investigating the STDP dynamics of two populations synapsing downstream onto the same neuron in a feed-forward manner. Each population was assumed to exhibit rhythmic activity, albeit in a different frequency band. Our theory reveals that the winner-take-all like competitions between the two populations is limited, in the sense that different rhythmic populations will not necessarily fully suppress each other. Furthermore, we found that for a wide range of parameters, the network converged to a solution in which the downstream neuron responded to both rhythms. Yet, the synaptic weights themselves did not converge to a fixed point, rather remained dynamic. These findings imply that STDP can support the multiplexing of rhythmic information, and demonstrate how functionality (multiplexing of information) can be retained in the face of continuous remodeling of all the synaptic weights. The constraints on the types of STDP rules that can support multiplexing provide a natural test for our theory.

Spike timing dependent plasticity (STDP) quantifies the change in the synaptic efficacy as a function of the temporal relationship between pre- and post-synaptic firing. STDP can be viewed as a microscopic unsupervised learning rule, and a wide range of such microscopic learning rules have been described empirically. Since there is no supervisor in unsupervised learning (which would provide with the system its goal), theoreticians have struggled with the question of the possible computational roles of the various STDP rules. Previous studies have focused on the possible contribution of STDP to the spontaneous development of spatial structure. However, the rich temporal repertoire of reported STDP rules has largely been ignored. Here we studied the contribution of STDP to the development of temporal structure. We show how STDP can shape synaptic efficacies to facilitate the transfer of rhythmic information downstream and to enable the multiplexing of information across different frequency channels. Our work emphasizes the relationship between the temporal structure of the STDP rule and the rhythmic activity it can support.

Understanding the computation of time using neural network models

Perceiving, maintaining, and using time intervals in working memory are crucial for animals to anticipate or act correctly at the right time in the ever-changing world. Here, we systematically study the underlying neural mechanisms by training recurrent neural networks to perform temporal tasks or complex tasks in combination with spatial information processing and decision making. We found that neural networks perceive time through state evolution along stereotypical trajectories and produce time intervals by scaling evolution speed. Temporal and nontemporal information is jointly coded in a way that facilitates decoding generalizability. We also provided potential sources for the temporal signals observed in nontiming tasks. Our study revealed the computational principles of a number of experimental phenomena and provided several predictions.

To maximize future rewards in this ever-changing world, animals must be able to discover the temporal structure of stimuli and then anticipate or act correctly at the right time. How do animals perceive, maintain, and use time intervals ranging from hundreds of milliseconds to multiseconds in working memory? How is temporal information processed concurrently with spatial information and decision making? Why are there strong neuronal temporal signals in tasks in which temporal information is not required? A systematic understanding of the underlying neural mechanisms is still lacking. Here, we addressed these problems using supervised training of recurrent neural network models. We revealed that neural networks perceive elapsed time through state evolution along stereotypical trajectory, maintain time intervals in working memory in the monotonic increase or decrease of the firing rates of interval-tuned neurons, and compare or produce time intervals by scaling state evolution speed. Temporal and nontemporal information is coded in subspaces orthogonal with each other, and the state trajectories with time at different nontemporal information are quasiparallel and isomorphic. Such coding geometry facilitates the decoding generalizability of temporal and nontemporal information across each other. The network structure exhibits multiple feedforward sequences that mutually excite or inhibit depending on whether their preferences of nontemporal information are similar or not. We identified four factors that facilitate strong temporal signals in nontiming tasks, including the anticipation of coming events. Our work discloses fundamental computational principles of temporal processing, and it is supported by and gives predictions to a number of experimental phenomena.

Clustering of Neural Activity: A Design Principle for Population Codes

We propose that correlations among neurons are generically strong enough to organize neural activity patterns into a discrete set of clusters, which can each be viewed as a population codeword. Our reasoning starts with the analysis of retinal ganglion cell data using maximum entropy models, showing that the population is robustly in a frustrated, marginally sub-critical, or glassy, state. This leads to an argument that neural populations in many other brain areas might share this structure. Next, we use latent variable models to show that this glassy state possesses well-defined clusters of neural activity. Clusters have three appealing properties: (i) clusters exhibit error correction, i.e., they are reproducibly elicited by the same stimulus despite variability at the level of constituent neurons; (ii) clusters encode qualitatively different visual features than their constituent neurons; and (iii) clusters can be learned by downstream neural circuits in an unsupervised fashion. We hypothesize that these properties give rise to a "learnable" neural code which the cortical hierarchy uses to extract increasingly complex features without supervision or reinforcement.

Excitatory neurons are more disinhibited than inhibitory neurons by chloride dysregulation in the spinal dorsal horn

Neuropathic pain is a debilitating condition caused by the abnormal processing of somatosensory input. Synaptic inhibition in the spinal dorsal horn plays a key role in that processing. Mechanical allodynia - the misperception of light touch as painful - occurs when inhibition is compromised. Disinhibition is due primarily to chloride dysregulation caused by hypofunction of the potassium-chloride co-transporter KCC2. Here we show, in rats, that excitatory neurons are disproportionately affected. This is not because chloride is differentially dysregulated in excitatory and inhibitory neurons, but, rather, because excitatory neurons rely more heavily on inhibition to counterbalance strong excitation. Receptive fields in both cell types have a center-surround organization but disinhibition unmasks more excitatory input to excitatory neurons. Differences in intrinsic excitability also affect how chloride dysregulation affects spiking. These results deepen understanding of how excitation and inhibition are normally balanced in the spinal dorsal horn, and how their imbalance disrupts somatosensory processing.