Despite advances in experimental techniques and accumulation of large datasets concerning the composition and properties of the cortex, quantitative modeling of cortical circuits under in-vivo-like conditions remains challenging. Here we report and publicly release a biophysically detailed circuit model of layer 4 in the mouse primary visual cortex, receiving thalamo-cortical visual inputs. The 45,000-neuron model was subjected to a battery of visual stimuli, and results were compared to published work and new in vivo experiments. Simulations reproduced a variety of observations, including effects of optogenetic perturbations. Critical to the agreement between responses in silico and in vivo were the rules of functional synaptic connectivity between neurons. Interestingly, after extreme simplification the model still performed satisfactorily on many measurements, although quantitative agreement with experiments suffered. These results emphasize the importance of functional rules of cortical wiring and enable a next generation of data-driven models of in vivo neural activity and computations.
Spiking dataset from McBride et al. 2022: Influence of claustrum on cortex varies by area, layer, and cell type Contains spiking data collected with Neuropixels probes from mouse cortical brain areas during optogenetic stimulation of claustrum. More details can be found in the preprint here: https://www.biorxiv.org/content/10.1101/2022.02.22.481532v2 Final publication: (link coming soon) Spiking data has been spike sorted with Kilosort. Other data includes wheel movement, optogenetic and visual stimulus information, trial timing, and metadata about each experiment and animal. To access data, unzip "McBride_et_al_2022_spiking_data.zip" into a folder. One file "metadata.csv" contains metadata, while each ".hdf5" file contains the data from a single probe recording from a single experiment. In most experiments, 2 or more probes recorded simultaneously. For more detail in accessing the data and examples in python, see the jupyter notebook provided "McBride et al. 2022 dataset tutorial.ipynb" -- also provided in .html.
In complex environments, animals can adopt diverse strategies to find rewards. How distinct strategies differentially engage brain circuits is not well understood. Here, we investigate this question, focusing on the cortical Vip-Sst disinhibitory circuit between vasoactive intestinal peptide-postive (Vip) interneurons and somatostatin-positive (Sst) interneurons. We characterize the behavioral strategies used by mice during a visual change detection task. Using a dynamic logistic regression model, we find that individual mice use mixtures of a visual comparison strategy and a statistical timing strategy. Separately, mice also have periods of task engagement and disengagement. Two-photon calcium imaging shows large strategy-dependent differences in neural activity in excitatory, Sst inhibitory, and Vip inhibitory cells in response to both image changes and image omissions. In contrast, task engagement has limited effects on neural population activity. We find that the diversity of neural correlates of strategy can be understood parsimoniously as the increased activation of the Vip-Sst disinhibitory circuit during the visual comparison strategy, which facilitates task-appropriate responses.
Abstract When sensory information is incomplete or ambiguous, the brain relies on prior expectations to infer perceptual objects. Despite the centrality of this process to perception, the neural mechanism of sensory inference is not known. Illusory contours (ICs) are key tools to study sensory inference because they contain edges or objects that are implied only by their spatial context. Using cellular resolution, mesoscale two-photon calcium imaging and multi-Neuropixels recordings in the mouse visual cortex, we identified a sparse subset of neurons in the primary visual cortex (V1) and higher visual areas that respond emergently to ICs. We found that these highly selective ‘IC-encoders’ mediate the neural representation of IC inference. Strikingly, selective activation of these neurons using two-photon holographic optogenetics was sufficient to recreate IC representation in the rest of the V1 network, in the absence of any visual stimulus. This outlines a model in which primary sensory cortex facilitates sensory inference by selectively strengthening input patterns that match prior expectations through local, recurrent circuitry. Our data thus suggest a clear computational purpose for recurrence in the generation of holistic percepts under sensory ambiguity. More generally, selective reinforcement of top-down predictions by pattern-completing recurrent circuits in lower sensory cortices may constitute a key step in sensory inference.
Abstract The visual cortex is organized hierarchically, but the presence of extensive recurrent and parallel pathways make it challenging to decipher how signals flow between neuronal populations. Here, we tracked the flow of spiking activity recorded from six interconnected levels of the mouse visual hierarchy. By analyzing leading and lagging spike-timing relationships among pairs of simultaneously recorded neurons, we created a cellular-scale directed network graph. Using a module-detection algorithm to cluster neurons based on shared connectivity patterns, we uncovered two multi-regional communication modules distributed across the hierarchy. The direction of signal flow between and within these modules, differences in layer and area distributions, and distinct temporal dynamics suggest that one module is positioned to transmit feedforward sensory signals, whereas the other integrates inputs for recurrent processing. These results suggest that multi-regional functional modules may be a fundamental feature of organization beyond cortical areas that supports signal propagation across hierarchical recurrent networks.
The maintenance of short-term memories is critical for survival in a dynamically changing world. Previous studies suggest that this memory can be stored in the form of persistent neural activity or using a synaptic mechanism, such as with short-term plasticity. Here, we compare the predictions of these two mechanisms to neural and behavioral measurements in a visual change detection task. Mice were trained to respond to changes in a repeated sequence of natural images while neural activity was recorded using two-photon calcium imaging. We also trained two types of artificial neural networks on the same change detection task as the mice. Following fixed pre-processing using a pretrained convolutional neural network, either a recurrent neural network (RNN) or a feedforward neural network with short-term synaptic depression (STPNet) was trained to the same level of performance as the mice. While both networks are able to learn the task, the STPNet model contains units whose activity are more similar to the in vivo data and produces errors which are more similar to the mice. When images are omitted, an unexpected perturbation which was absent during training, mice often do not respond to the omission but are more likely to respond to the subsequent image. Unlike the RNN model, STPNet produces a similar pattern of behavior. These results suggest that simple neural adaptation mechanisms may serve as an important bottom-up memory signal in this task, which can be used by downstream areas in the decision-making process.
Abstract Neurophysiological differentiation (ND), a metric that quantifies the number of distinct activity states that the brain or its part visits over a period of time, has been used as a correlate of meaningfulness or subjective perception of visual stimuli. ND has largely been studied in non-invasive human whole-brain recordings where spatial resolution is limited. However, it is likely that perception is supported by discrete populations of spiking neurons rather than the whole brain. Therefore, in this study, we use Neuropixels recordings from the mouse brain to characterize the ND metric within neural populations recorded at single-cell resolution in localized regions. Using the spiking activity of thousands of simultaneously recorded neurons spanning 6 visual cortical areas as well as the visual thalamus, we show that the ND of stimulus-evoked activity of the entire visual cortex is higher for naturalistic stimuli relative to artificial ones. This finding holds in most individual areas throughout the visual hierarchy as well. For animals performing an image change detection task, ND of the entire visual cortex (though not individual areas) is higher for successful detection compared to failed trials, consistent with the assumed perception of the stimulus. Analysis of spiking activity allows us to characterize the ND metric across a wide range of timescales from 10s of milliseconds to a few seconds. This analysis reveals that although ND of activity of single neurons is often maximized at an optimal timescale around 100 ms, the optimum shifts to under 5 ms for ND of neuronal ensembles. Finally, we find that the ND of activations in convolutional neural networks (CNNs) trained on an image classification task shows distinct trends relative to the mouse visual system: ND is often higher for less naturalistic stimuli and varies by orders of magnitude across the hierarchy, compared to modest variation in the mouse brain. Together, these results suggest that ND computed on cellular-level neural recordings can be a useful tool highlighting cell populations that may be involved in subjective perception. Summary Advances in our understanding on neural coding has revealed that information about visual stimuli is represented across several brain regions. However, availability of information does not imply that it is necessarily utilized by the brain, much less that it is subjectively perceived. Since percepts originate in neural activity, distinct percepts must be associated with distinct ‘states’ of neural activity, at least within the brain region that supports the percepts. Thus, one approach developed in this direction is to quantify the number of distinct ‘states’ that the activity of the brain goes through, called neurophysiological differentiation (ND). ND of the entire brain has been shown to reflect subjective reports of visual stimulus meaningfulness. But what specific subpopulations within the brain could be supporting conscious perception, and what is the correct timescale on which states should be quantified? In this study, we analyze ND of spiking neural activity in the mouse visual cortex recorded using Neuropixels probes, allowing us to characterize the ND metric across a wide range of timescales all the way down from 5 ms to a few seconds. It also allows us to understand the ND of neural activity of different ensembles of neurons, from individual thalamic or cortical ensembles to those spanning across multiple visual areas in the mouse brain.
Perturbational complexity analysis predicts the presence of consciousness in volunteers and patients by stimulating the brain with brief pulses, recording electroencephalographic (EEG) responses, and computing their spatiotemporal complexity. We examined the underlying neural circuits in mice by directly stimulating cortex while recording with EEG and Neuropixels probes during wakefulness and isoflurane anesthesia. When mice are awake, stimulation of deep cortical layers reliably evokes locally a brief pulse of excitation, followed by a bi-phasic sequence of 120 ms profound off period and a rebound excitation. A similar pattern, partially attributed to burst spiking, is seen in thalamic nuclei, and is associated with a pronounced late component in the evoked EEG. We infer that cortico-thalamo-cortical interactions drive the long-lasting evoked EEG signals elicited by deep cortical stimulation during the awake state. The cortical and thalamic off period and rebound excitation, and the late component in the EEG, are reduced during running and absent during anesthesia.
Local field potential (LFP) recordings reflect the dynamics of the current source density (CSD) in brain tissue. The synaptic, cellular, and circuit contributions to current sinks and sources are ill-understood. We investigated these in mouse primary visual cortex using public Neuropixels recordings and a detailed circuit model based on simulating the Hodgkin–Huxley dynamics of >50,000 neurons belonging to 17 cell types. The model simultaneously captured spiking and CSD responses and demonstrated a two-way dissociation: firing rates are altered with minor effects on the CSD pattern by adjusting synaptic weights, and CSD is altered with minor effects on firing rates by adjusting synaptic placement on the dendrites. We describe how thalamocortical inputs and recurrent connections sculpt specific sinks and sources early in the visual response, whereas cortical feedback crucially alters them in later stages. These results establish quantitative links between macroscopic brain measurements (LFP/CSD) and microscopic biophysics-based understanding of neuron dynamics and show that CSD analysis provides powerful constraints for modeling beyond those from considering spikes.
