A transcriptomic axis predicts state modulation of cortical interneurons.

Transcriptomics has revealed that cortical inhibitory neurons exhibit a great diversity of fine molecular subtypes1-6, but it is not known whether these subtypes have correspondingly diverse patterns of activity in the living brain. Here we show that inhibitory subtypes in primary visual cortex (V1) have diverse correlates with brain state, which are organized by a single factor: position along the main axis of transcriptomic variation. We combined in vivo two-photon calcium imaging of mouse V1 with a transcriptomic method to identify mRNA for 72 selected genes in ex vivo slices. We classified inhibitory neurons imaged in layers 1-3 into a three-level hierarchy of 5 subclasses, 11 types and 35 subtypes using previously defined transcriptomic clusters3. Responses to visual stimuli differed significantly only between subclasses, with cells in the Sncg subclass uniformly suppressed, and cells in the other subclasses predominantly excited. Modulation by brain state differed at all hierarchical levels but could be largely predicted from the first transcriptomic principal component, which also predicted correlations with simultaneously recorded cells. Inhibitory subtypes that fired more in resting, oscillatory brain states had a smaller fraction of their axonal projections in layer 1, narrower spikes, lower input resistance and weaker adaptation as determined in vitro7, and expressed more inhibitory cholinergic receptors. Subtypes that fired more during arousal had the opposite properties. Thus, a simple principle may largely explain how diverse inhibitory V1 subtypes shape state-dependent cortical processing.

Striatal activity topographically reflects cortical activity.

The cortex projects to the dorsal striatum topographically1,2 to regulate behaviour3-5, but spiking activity in the two structures has previously been reported to have markedly different relations to sensorimotor events6-9. Here we show that the relationship between activity in the cortex and striatum is spatiotemporally precise, topographic, causal and invariant to behaviour. We simultaneously recorded activity across large regions of the cortex and across the width of the dorsal striatum in mice that performed a visually guided task. Striatal activity followed a mediolateral gradient in which behavioural correlates progressed from visual cue to response movement to reward licking. The summed activity in each part of the striatum closely and specifically mirrored activity in topographically associated cortical regions, regardless of task engagement. This relationship held for medium spiny neurons and fast-spiking interneurons, whereas the activity of tonically active neurons differed from cortical activity with stereotypical responses to sensory or reward events. Inactivation of the visual cortex abolished striatal responses to visual stimuli, supporting a causal role of cortical inputs in driving the striatum. Striatal visual responses were larger in trained mice than untrained mice, with no corresponding change in overall activity in the visual cortex. Striatal activity therefore reflects a consistent, causal and scalable topographical mapping of cortical activity.

Spatial connectivity matches direction selectivity in visual cortex.

The selectivity of neuronal responses arises from the architecture of excitatory and inhibitory connections. In the primary visual cortex, the selectivity of a neuron in layer 2/3 for stimulus orientation and direction is thought to arise from intracortical inputs that are similarly selective1-8. However, the excitatory inputs of a neuron can have diverse stimulus preferences1-4,6,7,9, and inhibitory inputs can be promiscuous10 and unselective11. Here we show that the excitatory and inhibitory intracortical connections to a layer 2/3 neuron accord with its selectivity by obeying precise spatial patterns. We used rabies tracing1,12 to label and functionally image the excitatory and inhibitory inputs to individual pyramidal neurons of layer 2/3 of the mouse visual cortex. Presynaptic excitatory neurons spanned layers 2/3 and 4 and were distributed coaxial to the preferred orientation of the postsynaptic neuron, favouring the region opposite to its preferred direction. By contrast, presynaptic inhibitory neurons resided within layer 2/3 and favoured locations near the postsynaptic neuron and ahead of its preferred direction. The direction selectivity of a postsynaptic neuron was unrelated to the selectivity of presynaptic neurons, but correlated with the spatial displacement between excitatory and inhibitory presynaptic ensembles. Similar asymmetric connectivity establishes direction selectivity in the retina13-17. This suggests that this circuit motif might be canonical in sensory processing.

Distributed coding of choice, action and engagement across the mouse brain.

Vision, choice, action and behavioural engagement arise from neuronal activity that may be distributed across brain regions. Here we delineate the spatial distribution of neurons underlying these processes. We used Neuropixels probes1,2 to record from approximately 30,000 neurons in 42 brain regions of mice performing a visual discrimination task3. Neurons in nearly all regions responded non-specifically when the mouse initiated an action. By contrast, neurons encoding visual stimuli and upcoming choices occupied restricted regions in the neocortex, basal ganglia and midbrain. Choice signals were rare and emerged with indistinguishable timing across regions. Midbrain neurons were activated before contralateral choices and were suppressed before ipsilateral choices, whereas forebrain neurons could prefer either side. Brain-wide pre-stimulus activity predicted engagement in individual trials and in the overall task, with enhanced subcortical but suppressed neocortical activity during engagement. These results reveal organizing principles for the distribution of neurons encoding behaviourally relevant variables across the mouse brain.

High-dimensional geometry of population responses in visual cortex.

A neuronal population encodes information most efficiently when its stimulus responses are high-dimensional and uncorrelated, and most robustly when they are lower-dimensional and correlated. Here we analysed the dimensionality of the encoding of natural images by large populations of neurons in the visual cortex of awake mice. The evoked population activity was high-dimensional, and correlations obeyed an unexpected power law: the nth principal component variance scaled as 1/n. This scaling was not inherited from the power law spectrum of natural images, because it persisted after stimulus whitening. We proved mathematically that if the variance spectrum was to decay more slowly then the population code could not be smooth, allowing small changes in input to dominate population activity. The theory also predicts larger power-law exponents for lower-dimensional stimulus ensembles, which we validated experimentally. These results suggest that coding smoothness may represent a fundamental constraint that determines correlations in neural population codes.

Coherent encoding of subjective spatial position in visual cortex and hippocampus.

A major role of vision is to guide navigation, and navigation is strongly driven by vision1-4. Indeed, the brain's visual and navigational systems are known to interact5,6, and signals related to position in the environment have been suggested to appear as early as in the visual cortex6,7. Here, to establish the nature of these signals, we recorded in the primary visual cortex (V1) and hippocampal area CA1 while mice traversed a corridor in virtual reality. The corridor contained identical visual landmarks in two positions, so that a purely visual neuron would respond similarly at those positions. Most V1 neurons, however, responded solely or more strongly to the landmarks in one position rather than the other. This modulation of visual responses by spatial location was not explained by factors such as running speed. To assess whether the modulation is related to navigational signals and to the animal's subjective estimate of position, we trained the mice to lick for a water reward upon reaching a reward zone in the corridor. Neuronal populations in both CA1 and V1 encoded the animal's position along the corridor, and the errors in their representations were correlated. Moreover, both representations reflected the animal's subjective estimate of position, inferred from the animal's licks, better than its actual position. When animals licked in a given location-whether correctly or incorrectly-neural populations in both V1 and CA1 placed the animal in the reward zone. We conclude that visual responses in V1 are controlled by navigational signals, which are coherent with those encoded in hippocampus and reflect the animal's subjective position. The presence of such navigational signals as early as a primary sensory area suggests that they permeate sensory processing in the cortex.

Fully integrated silicon probes for high-density recording of neural activity.

Sensory, motor and cognitive operations involve the coordinated action of large neuronal populations across multiple brain regions in both superficial and deep structures. Existing extracellular probes record neural activity with excellent spatial and temporal (sub-millisecond) resolution, but from only a few dozen neurons per shank. Optical Ca2+ imaging offers more coverage but lacks the temporal resolution needed to distinguish individual spikes reliably and does not measure local field potentials. Until now, no technology compatible with use in unrestrained animals has combined high spatiotemporal resolution with large volume coverage. Here we design, fabricate and test a new silicon probe known as Neuropixels to meet this need. Each probe has 384 recording channels that can programmably address 960 complementary metal-oxide-semiconductor (CMOS) processing-compatible low-impedance TiN sites that tile a single 10-mm long, 70 × 20-µm cross-section shank. The 6 × 9-mm probe base is fabricated with the shank on a single chip. Voltage signals are filtered, amplified, multiplexed and digitized on the base, allowing the direct transmission of noise-free digital data from the probe. The combination of dense recording sites and high channel count yielded well-isolated spiking activity from hundreds of neurons per probe implanted in mice and rats. Using two probes, more than 700 well-isolated single neurons were recorded simultaneously from five brain structures in an awake mouse. The fully integrated functionality and small size of Neuropixels probes allowed large populations of neurons from several brain structures to be recorded in freely moving animals. This combination of high-performance electrode technology and scalable chip fabrication methods opens a path towards recording of brain-wide neural activity during behaviour.

Dopamine Axons in Dorsal Striatum Encode Contralateral Visual Stimuli and Choices.

The striatum plays critical roles in visually-guided decision-making and receives dense axonal projections from midbrain dopamine neurons. However, the roles of striatal dopamine in visual decision-making are poorly understood. We trained male and female mice to perform a visual decision task with asymmetric reward payoff, and we recorded the activity of dopamine axons innervating striatum. Dopamine axons in the dorsomedial striatum (DMS) responded to contralateral visual stimuli and contralateral rewarded actions. Neural responses to contralateral stimuli could not be explained by orienting behavior such as eye movements. Moreover, these contralateral stimulus responses persisted in sessions where the animals were instructed to not move to obtain reward, further indicating that these signals are stimulus-related. Lastly, we show that DMS dopamine signals were qualitatively different from dopamine signals in the ventral striatum (VS), which responded to both ipsilateral and contralateral stimuli, conforming to canonical prediction error signaling under sensory uncertainty. Thus, during visual decisions, DMS dopamine encodes visual stimuli and rewarded actions in a lateralized fashion, and could facilitate associations between specific visual stimuli and actions.SIGNIFICANCE STATEMENT While the striatum is central to goal-directed behavior, the precise roles of its rich dopaminergic innervation in perceptual decision-making are poorly understood. We found that in a visual decision task, dopamine axons in the dorsomedial striatum (DMS) signaled stimuli presented contralaterally to the recorded hemisphere, as well as the onset of rewarded actions. Stimulus-evoked signals persisted in a no-movement task variant. We distinguish the patterns of these signals from those in the ventral striatum (VS). Our results contribute to the characterization of region-specific dopaminergic signaling in the striatum and highlight a role in stimulus-action association learning.

Diverse coupling of neurons to populations in sensory cortex.

A large population of neurons can, in principle, produce an astronomical number of distinct firing patterns. In cortex, however, these patterns lie in a space of lower dimension, as if individual neurons were "obedient members of a huge orchestra". Here we use recordings from the visual cortex of mouse (Mus musculus) and monkey (Macaca mulatta) to investigate the relationship between individual neurons and the population, and to establish the underlying circuit mechanisms. We show that neighbouring neurons can differ in their coupling to the overall firing of the population, ranging from strongly coupled 'choristers' to weakly coupled 'soloists'. Population coupling is largely independent of sensory preferences, and it is a fixed cellular attribute, invariant to stimulus conditions. Neurons with high population coupling are more strongly affected by non-sensory behavioural variables such as motor intention. Population coupling reflects a causal relationship, predicting the response of a neuron to optogenetically driven increases in local activity. Moreover, population coupling indicates synaptic connectivity; the population coupling of a neuron, measured in vivo, predicted subsequent in vitro estimates of the number of synapses received from its neighbours. Finally, population coupling provides a compact summary of population activity; knowledge of the population couplings of n neurons predicts a substantial portion of their n(2) pairwise correlations. Population coupling therefore represents a novel, simple measure that characterizes the relationship of each neuron to a larger population, explaining seemingly complex network firing patterns in terms of basic circuit variables.

Inhibition dominates sensory responses in the awake cortex.

The activity of the cerebral cortex is thought to depend on the precise relationship between synaptic excitation and inhibition. In the visual cortex, in particular, intracellular measurements have related response selectivity to coordinated increases in excitation and inhibition. These measurements, however, have all been made during anaesthesia, which strongly influences cortical state and therefore sensory processing. The synaptic activity that is evoked by visual stimulation during wakefulness is unknown. Here we measured visually evoked responses--and the underlying synaptic conductances--in the visual cortex of anaesthetized and awake mice. Under anaesthesia, responses could be elicited from a large region of visual space and were prolonged. During wakefulness, responses were more spatially selective and much briefer. Whole-cell patch-clamp recordings of synaptic conductances showed a difference in synaptic inhibition between the two conditions. Under anaesthesia, inhibition tracked excitation in amplitude and spatial selectivity. By contrast, during wakefulness, inhibition was much stronger than excitation and had extremely broad spatial selectivity. We conclude that during wakefulness, cortical responses to visual stimulation are dominated by synaptic inhibition, restricting the spatial spread and temporal persistence of neural activity. These results provide a direct glimpse of synaptic mechanisms that control sensory responses in the awake cortex.

Adaptable history biases in human perceptual decisions.

When making choices under conditions of perceptual uncertainty, past experience can play a vital role. However, it can also lead to biases that worsen decisions. Consistent with previous observations, we found that human choices are influenced by the success or failure of past choices even in a standard two-alternative detection task, where choice history is irrelevant. The typical bias was one that made the subject switch choices after a failure. These choice history biases led to poorer performance and were similar for observers in different countries. They were well captured by a simple logistic regression model that had been previously applied to describe psychophysical performance in mice. Such irrational biases seem at odds with the principles of reinforcement learning, which would predict exquisite adaptability to choice history. We therefore asked whether subjects could adapt their irrational biases following changes in trial order statistics. Adaptability was strong in the direction that confirmed a subject's default biases, but weaker in the opposite direction, so that existing biases could not be eradicated. We conclude that humans can adapt choice history biases, but cannot easily overcome existing biases even if irrational in the current context: adaptation is more sensitive to confirmatory than contradictory statistics.

Imaging the awake visual cortex with a genetically encoded voltage indicator.

Genetically encoded voltage indicators (GEVIs) promise to reveal the membrane potential of genetically targeted neuronal populations through noninvasive, chronic imaging of large portions of cortical space. Here we test a promising GEVI in mouse cortex during wakefulness, a challenging condition due to large hemodynamic activity, and we introduce a straightforward projection method to separate a signal dominated by membrane voltage from a signal dominated by hemodynamic activity. We expressed VSFP-Butterfly 1.2 plasmid in layer 2/3 pyramidal cells of visual cortex through electroporation in utero. We then used wide-field imaging with two cameras to measure both fluorophores of the indicator in response to visual stimuli. By taking weighted sums and differences of the two measurements, we obtained clear separation of hemodynamic and voltage signals. The hemodynamic signal showed strong heartbeat oscillations, superimposed on slow dynamics similar to blood oxygen level-dependent (BOLD) or "intrinsic" signals. The voltage signal had fast dynamics similar to neural responses measured electrically, and showed an orderly retinotopic mapping. We compared this voltage signal with calcium signals imaged in transgenic mice that express a calcium indicator (GCaMP3) throughout cortex. The voltage signal from VSFP had similar signal-to-noise ratios as the calcium signal, it was more immune to vascular artifacts, and it integrated over larger regions of visual space, which was consistent with its reporting mostly subthreshold activity rather than the spiking activity revealed by calcium signals. These results demonstrate that GEVIs provide a powerful tool to study the dynamics of neural populations at mesoscopic spatial scales in the awake cortex.

Cortical state determines global variability and correlations in visual cortex.

The response of neurons in sensory cortex to repeated stimulus presentations is highly variable. To investigate the nature of this variability, we compared the spike activity of neurons in the primary visual cortex (V1) of cats with that of their afferents from lateral geniculate nucleus (LGN), in response to similar stimuli. We found variability to be much higher in V1 than in LGN. To investigate the sources of the additional variability, we measured the spiking activity of large V1 populations and found that much of the variability was shared across neurons: the variable portion of the responses of one neuron could be well predicted from the summed activity of the rest of the neurons. Variability thus mostly reflected global fluctuations affecting all neurons. The size and prevalence of these fluctuations, both in responses to stimuli and in ongoing activity, depended on cortical state, being larger in synchronized states than in more desynchronized states. Contrary to previous reports, these fluctuations invested the overall population, regardless of preferred orientation. The global fluctuations substantially increased variability in single neurons and correlations among pairs of neurons. Once this effect was removed, pairwise correlations were reduced and were similar regardless of cortical state. These results highlight the importance of cortical state in controlling cortical operation and can help reconcile previous studies, which differed widely in their estimate of neuronal variability and pairwise correlations.

Local and global contributions to hemodynamic activity in mouse cortex.

Imaging techniques such as functional magnetic resonance imaging seek to estimate neural signals in local brain regions through measurements of hemodynamic activity. However, hemodynamic activity is accompanied by large vascular fluctuations of unclear significance. To characterize these fluctuations and their impact on estimates of neural signals, we used optical imaging in visual cortex of awake mice. We found that hemodynamic activity can be expressed as the sum of two components, one local and one global. The local component reflected presumed neural signals driven by visual stimuli in the appropriate retinotopic region. The global component constituted large fluctuations shared by larger cortical regions, which extend beyond visual cortex. These fluctuations varied from trial to trial, but they did not constitute noise; they correlated with pupil diameter, suggesting that they reflect variations in arousal or alertness. Distinguishing local and global contributions to hemodynamic activity may help understand neurovascular coupling and interpret measurements of hemodynamic responses.

Normalization as a canonical neural computation.

There is increasing evidence that the brain relies on a set of canonical neural computations, repeating them across brain regions and modalities to apply similar operations to different problems. A promising candidate for such a computation is normalization, in which the responses of neurons are divided by a common factor that typically includes the summed activity of a pool of neurons. Normalization was developed to explain responses in the primary visual cortex and is now thought to operate throughout the visual system, and in many other sensory modalities and brain regions. Normalization may underlie operations such as the representation of odours, the modulatory effects of visual attention, the encoding of value and the integration of multisensory information. Its presence in such a diversity of neural systems in multiple species, from invertebrates to mammals, suggests that it serves as a canonical neural computation.

Five key factors determining pairwise correlations in visual cortex.

The responses of cortical neurons to repeated presentation of a stimulus are highly variable, yet correlated. These "noise correlations" reflect a low-dimensional structure of population dynamics. Here, we examine noise correlations in 22,705 pairs of neurons in primary visual cortex (V1) of anesthetized cats, during ongoing activity and in response to artificial and natural visual stimuli. We measured how noise correlations depend on 11 factors. Because these factors are themselves not independent, we distinguished their influences using a nonlinear additive model. The model revealed that five key factors play a predominant role in determining pairwise correlations. Two of these are distance in cortex and difference in sensory tuning: these are known to decrease correlation. A third factor is firing rate: confirming most earlier observations, it markedly increased pairwise correlations. A fourth factor is spike width: cells with a broad spike were more strongly correlated amongst each other. A fifth factor is spike isolation: neurons with worse isolation were more correlated, even if they were recorded on different electrodes. For pairs of neurons with poor isolation, this last factor was the main determinant of correlations. These results were generally independent of stimulus type and timescale of analysis, but there were exceptions. For instance, pairwise correlations depended on difference in orientation tuning more during responses to gratings than to natural stimuli. These results consolidate disjoint observations in a vast literature on pairwise correlations and point towards regularities of population coding in sensory cortex.

Fast hemodynamic responses in the visual cortex of the awake mouse.

Hemodynamic responses in mice and other species are typically measured under anesthesia. However, anesthesia could influence their relationship to neural activity. To investigate this relationship, we used optical imaging in mouse primary visual cortex (V1). Hemodynamic responses yielded clear maps of retinotopy in both anesthetized and awake mice. However, during wakefulness, responses were four times larger and twice as fast. These differences held whether we induced anesthesia with urethane or isoflurane and whether awake mice were stationary or running on a treadmill. With electrode recordings, we established that the effects of wakefulness reflect changes in neurovascular coupling, not in neural activity. By activating V1 directly via optogenetics, we replicated the effects of wakefulness in terms of timing but not of amplitude. We conclude that neurovascular coupling depends critically on anesthesia and wakefulness: during wakefulness, neural activity is followed by much stronger and quicker hemodynamic responses.

Ultra-high density electrodes improve detection, yield, and cell type identification in neuronal recordings.

To understand the neural basis of behavior, it is essential to sensitively and accurately measure neural activity at single neuron and single spike resolution. Extracellular electrophysiology delivers this, but it has biases in the neurons it detects and it imperfectly resolves their action potentials. To minimize these limitations, we developed a silicon probe with much smaller and denser recording sites than previous designs, called Neuropixels Ultra (NP Ultra). This device samples neuronal activity at ultra-high spatial density (~10 times higher than previous probes) with low noise levels, while trading off recording span. NP Ultra is effectively an implantable voltage-sensing camera that captures a planar image of a neuron's electrical field. We use a spike sorting algorithm optimized for these probes to demonstrate that the yield of visually-responsive neurons in recordings from mouse visual cortex improves up to ~3-fold. We show that NP Ultra can record from small neuronal structures including axons and dendrites. Recordings across multiple brain regions and four species revealed a subset of extracellular action potentials with unexpectedly small spatial spread and axon-like features. We share a large-scale dataset of these brain-wide recordings in mice as a resource for studies of neuronal biophysics. Finally, using ground-truth identification of three major inhibitory cortical cell types, we found that these cell types were discriminable with approximately 75% success, a significant improvement over lower-resolution recordings. NP Ultra improves spike sorting performance, detection of subcellular compartments, and cell type classification to enable more powerful dissection of neural circuit activity during behavior.