Survey of spiking in the mouse visual system reveals functional hierarchy.

The anatomy of the mammalian visual system, from the retina to the neocortex, is organized hierarchically1. However, direct observation of cellular-level functional interactions across this hierarchy is lacking due to the challenge of simultaneously recording activity across numerous regions. Here we describe a large, open dataset-part of the Allen Brain Observatory2-that surveys spiking from tens of thousands of units in six cortical and two thalamic regions in the brains of mice responding to a battery of visual stimuli. Using cross-correlation analysis, we reveal that the organization of inter-area functional connectivity during visual stimulation mirrors the anatomical hierarchy from the Allen Mouse Brain Connectivity Atlas3. We find that four classical hierarchical measures-response latency, receptive-field size, phase-locking to drifting gratings and response decay timescale-are all correlated with the hierarchy. Moreover, recordings obtained during a visual task reveal that the correlation between neural activity and behavioural choice also increases along the hierarchy. Our study provides a foundation for understanding coding and signal propagation across hierarchically organized cortical and thalamic visual areas.

Removing independent noise in systems neuroscience data using DeepInterpolation.

Progress in many scientific disciplines is hindered by the presence of independent noise. Technologies for measuring neural activity (calcium imaging, extracellular electrophysiology and functional magnetic resonance imaging (fMRI)) operate in domains in which independent noise (shot noise and/or thermal noise) can overwhelm physiological signals. Here, we introduce DeepInterpolation, a general-purpose denoising algorithm that trains a spatiotemporal nonlinear interpolation model using only raw noisy samples. Applying DeepInterpolation to two-photon calcium imaging data yielded up to six times more neuronal segments than those computed from raw data with a 15-fold increase in the single-pixel signal-to-noise ratio (SNR), uncovering single-trial network dynamics that were previously obscured by noise. Extracellular electrophysiology recordings processed with DeepInterpolation yielded 25% more high-quality spiking units than those computed from raw data, while DeepInterpolation produced a 1.6-fold increase in the SNR of individual voxels in fMRI datasets. Denoising was attained without sacrificing spatial or temporal resolution and without access to ground truth training data. We anticipate that DeepInterpolation will provide similar benefits in other domains in which independent noise contaminates spatiotemporally structured datasets.

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.

Population burst propagation across interacting areas of the brain.

For many perceptual and behavioral tasks, a prominent feature of neural spike trains involves high firing rates across relatively short intervals of time. We call these events "population bursts." Because during a population burst information is, presumably, transmitted from one part of the brain to another, burst timing should reveal activity related to the flow of information across neural circuits. We developed a statistical method (based on a point process model) of determining, accurately, the time of the maximum (peak) population firing rate on a trial-by-trial basis and used it to characterize burst propagation across areas. We then examined the tendency of peak firing rates in distinct brain areas to shift earlier or later in time, together, across repeated trials, and found this trial-to-trial coupling of peak times to be a sensitive indicator of interaction across populations. In the data we examined, from the Allen Brain Observatory, we found many very strong correlations (95% confidence intervals above 0.75) in cases where standard methods were unable to demonstrate cross-area correlation. The statistical model introduced cross-area covariation only through population-level trial-dependent time shifts and gain constants (values of which were learned from the data), yet it provided very good fits to data histograms, including histograms of spike count correlations within and across visual areas. Our results demonstrate the utility of carefully assessing timing and propagation, across brain regions, of transient bursts in neural population activity, based on multiple spike train recordings.NEW & NOTEWORTHY We developed a novel statistical method for identifying coordinated propagation of activity across populations of spiking neurons, with high temporal accuracy. Using simultaneous recordings from three visual areas we document precise timing relationships on a trial-by-trial basis, and we show how previously existing techniques can fail to discover coordinated activity in cases where the new approach finds very strong cross-area correlation.

Cross-population coupling of neural activity based on Gaussian process current source densities.

Because local field potentials (LFPs) arise from multiple sources in different spatial locations, they do not easily reveal coordinated activity across neural populations on a trial-to-trial basis. As we show here, however, once disparate source signals are decoupled, their trial-to-trial fluctuations become more accessible, and cross-population correlations become more apparent. To decouple sources we introduce a general framework for estimation of current source densities (CSDs). In this framework, the set of LFPs result from noise being added to the transform of the CSD by a biophysical forward model, while the CSD is considered to be the sum of a zero-mean, stationary, spatiotemporal Gaussian process, having fast and slow components, and a mean function, which is the sum of multiple time-varying functions distributed across space, each varying across trials. We derived biophysical forward models relevant to the data we analyzed. In simulation studies this approach improved identification of source signals compared to existing CSD estimation methods. Using data recorded from primate auditory cortex, we analyzed trial-to-trial fluctuations in both steady-state and task-evoked signals. We found cortical layer-specific phase coupling between two probes and showed that the same analysis applied directly to LFPs did not recover these patterns. We also found task-evoked CSDs to be correlated across probes, at specific cortical depths. Using data from Neuropixels probes in mouse visual areas, we again found evidence for depth-specific phase coupling of primary visual cortex and lateromedial area based on the CSDs.

High-density extracellular probes reveal dendritic backpropagation and facilitate neuron classification.

Different neuron types serve distinct roles in neural processing. Extracellular electrical recordings are extensively used to study brain function but are typically blind to cell identity. Morphoelectrical properties of neurons measured on spatially dense electrode arrays have the potential to distinguish neuron types. We used high-density silicon probes to record from cortical and subcortical regions of the mouse brain. Extracellular waveforms of each neuron were detected across many channels and showed distinct spatiotemporal profiles among brain regions. Classification of neurons by brain region was improved with multichannel compared with single-channel waveforms. In visual cortex, unsupervised clustering identified the canonical regular-spiking (RS) and fast-spiking (FS) classes but also indicated a subclass of RS units with unidirectional backpropagating action potentials (BAPs). Moreover, BAPs were observed in many hippocampal RS cells. Overall, waveform analysis of spikes from high-density probes aids neuron identification and can reveal dendritic backpropagation. NEW & NOTEWORTHY It is challenging to identify neuron types with extracellular electrophysiology in vivo. We show that spatiotemporal action potentials measured on high-density electrode arrays can capture cell type-specific morphoelectrical properties, allowing classification of neurons across brain structures and within the cortex. Moreover, backpropagating action potentials are reliably detected in vivo from subpopulations of cortical and hippocampal neurons. Together, these results enhance the utility of dense extracellular electrophysiology for cell-type interrogation of brain network function.

Spatial Organization of Chromatic Pathways in the Mouse Dorsal Lateral Geniculate Nucleus.

In both dichromats and trichromats, cone opsin signals are maintained independently in cones and combined at the bipolar and retinal ganglion cell level, creating parallel color opponent pathways to the central visual system. Like other dichromats, the mouse retina expresses a short-wavelength (S) and a medium-wavelength (M) opsin, with the S-opsin shifted to peak sensitivity in the ultraviolet (UV) range. Unlike in primates, nonuniform opsin expression across the retina and coexpression in single cones creates a mostly mixed chromatic signal. Here, we describe the visuotopic and chromatic organization of spiking responses in the dorsal lateral geniculate and of the local field potentials in their recipient zone in primary visual cortex (V1). We used an immersive visual stimulus dome that allowed us to present spatiotemporally modulated UV and green luminance in any region of the visual field of an awake, head-fixed mouse. Consistent with retinal expression of opsins, we observed graded UV-to-green dominated responses from the upper to lower visual fields, with a smaller difference across azimuth. In addition, we identified a subpopulation of cells (<10%) that exhibited spectrally opponent responses along the S-M axis. Luminance signals of each wavelength and color signals project to the middle layers of V1. SIGNIFICANCE STATEMENT: In natural environments, color information is useful for guiding behavior. How small terrestrial mammals such as mice use graded expression of cone opsins to extract visual information from their environments is not clear, even as the use of mice for studying visually guided behavior grows. In this study, we examined the color signals that the retina sends to the visual cortex via the lateral geniculate nucleus of the thalamus. We found that green dominated responses in the lower and nasal visual field and ultraviolet dominated responses in the upper visual field. We describe a subset of cells that exhibit color opponent responses.

Closed-loop, ultraprecise, automated craniotomies.

A large array of neuroscientific techniques, including in vivo electrophysiology, two-photon imaging, optogenetics, lesions, and microdialysis, require access to the brain through the skull. Ideally, the necessary craniotomies could be performed in a repeatable and automated fashion, without damaging the underlying brain tissue. Here we report that when drilling through the skull a stereotypical increase in conductance can be observed when the drill bit passes through the skull base. We present an architecture for a robotic device that can perform this algorithm, along with two implementations--one based on homebuilt hardware and one based on commercially available hardware--that can automatically detect such changes and create large numbers of precise craniotomies, even in a single skull. We also show that this technique can be adapted to automatically drill cranial windows several millimeters in diameter. Such robots will not only be useful for helping neuroscientists perform both small and large craniotomies more reliably but can also be used to create precisely aligned arrays of craniotomies with stereotaxic registration to standard brain atlases that would be difficult to drill by hand.

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.

SHIELD: Skull-shaped hemispheric implants enabling large-scale electrophysiology datasets in the mouse brain.

To understand the neural basis of behavior, it is essential to measure spiking dynamics across many interacting brain regions. Although new technologies, such as Neuropixels probes, facilitate multi-regional recordings, significant surgical and procedural hurdles remain for these experiments to achieve their full potential. Here, we describe skull-shaped hemispheric implants enabling large-scale electrophysiology datasets (SHIELD). These 3D-printed skull-replacement implants feature customizable insertion holes, allowing dozens of cortical and subcortical structures to be recorded in a single mouse using repeated multi-probe insertions over many days. We demonstrate the procedure's high success rate, biocompatibility, lack of adverse effects on behavior, and compatibility with imaging and optogenetics. To showcase SHIELD's scientific utility, we use multi-probe recordings to reveal novel insights into how alpha rhythms organize spiking activity across visual and sensorimotor networks. Overall, this method enables powerful, large-scale electrophysiological experiments for the study of distributed neural computation.

Sharing neurophysiology data from the Allen Brain Observatory.

Nullius in verba ('trust no one'), chosen as the motto of the Royal Society in 1660, implies that independently verifiable observations-rather than authoritative claims-are a defining feature of empirical science. As the complexity of modern scientific instrumentation has made exact replications prohibitive, sharing data is now essential for ensuring the trustworthiness of one's findings. While embraced in spirit by many, in practice open data sharing remains the exception in contemporary systems neuroscience. Here, we take stock of the Allen Brain Observatory, an effort to share data and metadata associated with surveys of neuronal activity in the visual system of laboratory mice. Data from these surveys have been used to produce new discoveries, to validate computational algorithms, and as a benchmark for comparison with other data, resulting in over 100 publications and preprints to date. We distill some of the lessons learned about open surveys and data reuse, including remaining barriers to data sharing and what might be done to address these.

Compression strategies for large-scale electrophysiology data.

Objective.With the rapid adoption of high-density electrode arrays for recording neural activity, electrophysiology data volumes within labs and across the field are growing at unprecedented rates. For example, a one-hour recording with a 384-channel Neuropixels probe generates over 80 GB of raw data. These large data volumes carry a high cost, especially if researchers plan to store and analyze their data in the cloud. Thus, there is a pressing need for strategies that can reduce the data footprint of each experiment.Approach.Here, we establish a set of benchmarks for comparing the performance of various compression algorithms on experimental and simulated recordings from Neuropixels 1.0 (NP1) and 2.0 (NP2) probes.Main results.For lossless compression, audio codecs (FLACandWavPack) achieve compression ratios (CRs) 6% higher for NP1 and 10% higher for NP2 than the best general-purpose codecs, at the expense of decompression speed. For lossy compression, theWavPackalgorithm in 'hybrid mode' increases the CR from 3.59 to 7.08 for NP1 and from 2.27 to 7.04 for NP2 (compressed file size of ∼14% for both types of probes), without adverse effects on spike sorting accuracy or spike waveforms.Significance.Along with the tools we have developed to make compression easier to deploy, these results should encourage all electrophysiologists to apply compression as part of their standard analysis workflows.

Pinpoint: trajectory planning for multi-probe electrophysiology and injections in an interactive web-based 3D environment.

Targeting deep brain structures during electrophysiology and injections requires intensive training and expertise. Even with experience, researchers often can't be certain that a probe is placed precisely in a target location and this complexity scales with the number of simultaneous probes used in an experiment. Here, we present Pinpoint, open-source software that allows for interactive exploration of stereotaxic insertion plans. Once an insertion plan is created, Pinpoint allows users to save these online and share them with collaborators. 3D modeling tools allow users to explore their insertions alongside rig and implant hardware and ensure plans are physically possible. Probes in Pinpoint can be linked to electronic micro-manipulators allowing real-time visualization of current brain region targets alongside neural data. In addition, Pinpoint can control manipulators to automate and parallelize the insertion process. Compared to previously available software, Pinpoint's easy access through web browsers, extensive features, and real-time experiment integration enable more efficient and reproducible recordings.

Acute head-fixed recordings in awake mice with multiple Neuropixels probes.

Multi-electrode arrays such as Neuropixels probes enable electrophysiological recordings from large populations of single neurons with high temporal resolution. By using such probes, the activity from functionally interacting, yet distinct, brain regions can be measured simultaneously by inserting multiple probes into the same subject. However, the use of multiple probes in small animals such as mice requires the removal of a sizable fraction of the skull, while also minimizing tissue damage and keeping the brain stable during the recordings. Here, we describe a step-by-step process designed to facilitate reliable recordings from up to six Neuropixels probes simultaneously in awake, head-fixed mice. The procedure involves four stages: the implantation of a headframe and a removable glass coverslip, the precise positioning of the Neuropixels probes at targeted points on the brain surface, the placement of a perforated plastic imaging window and the insertion of the probes into the brain of an awake mouse. The approach provides access to multiple brain regions and has been successfully applied across hundreds of mice. The procedure has been optimized for dense recordings from the mouse visual system, but it can be adapted for alternative recording configurations to target multiple probes in other brain areas. The protocol is suitable for users with experience in stereotaxic surgery in mice.

A survey of neurophysiological differentiation across mouse visual brain areas and timescales.

Neurophysiological differentiation (ND), a measure of the number of distinct activity states that a neural population visits over a time interval, 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 neuronal populations rather than the whole brain. Therefore, here we use Neuropixels recordings from the mouse brain to characterize the ND metric across a wide range of temporal scales, 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 and 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. Moreover, 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. Together, these results suggest that ND computed on cellular-level neural recordings is a useful tool highlighting cell populations that may be involved in subjective perception.

Associations between in vitro, in vivo and in silico cell classes in mouse primary visual cortex.

The brain consists of many cell classes yet in vivo electrophysiology recordings are typically unable to identify and monitor their activity in the behaving animal. Here, we employed a systematic approach to link cellular, multi-modal in vitro properties from experiments with in vivo recorded units via computational modeling and optotagging experiments. We found two one-channel and six multi-channel clusters in mouse visual cortex with distinct in vivo properties in terms of activity, cortical depth, and behavior. We used biophysical models to map the two one- and the six multi-channel clusters to specific in vitro classes with unique morphology, excitability and conductance properties that explain their distinct extracellular signatures and functional characteristics. These concepts were tested in ground-truth optotagging experiments with two inhibitory classes unveiling distinct in vivo properties. This multi-modal approach presents a powerful way to separate in vivo clusters and infer their cellular properties from first principles.

Associations between in vitro , in vivo and in silico cell classes in mouse primary visual cortex.

The brain consists of many cell classes yet in vivo electrophysiology recordings are typically unable to identify and monitor their activity in the behaving animal. Here, we employed a systematic approach to link cellular, multi-modal in vitro properties from experiments with in vivo recorded units via computational modeling and optotagging experiments. We found two one-channel and six multi-channel clusters in mouse visual cortex with distinct in vivo properties in terms of activity, cortical depth, and behavior. We used biophysical models to map the two one- and the six multi-channel clusters to specific in vitro classes with unique morphology, excitability and conductance properties that explain their distinct extracellular signatures and functional characteristics. These concepts were tested in ground-truth optotagging experiments with two inhibitory classes unveiling distinct in vivo properties. This multi-modal approach presents a powerful way to separate in vivo clusters and infer their cellular properties from first principles.

A unified open-source platform for multimodal neural recording and perturbation during naturalistic behavior.

Behavioral neuroscience faces two conflicting demands: long-duration recordings from large neural populations and unimpeded animal behavior. To meet this challenge, we developed ONIX, an open-source data acquisition system with high data throughput (2GB/sec) and low closed-loop latencies (<1ms) that uses a novel 0.3 mm thin tether to minimize behavioral impact. Head position and rotation are tracked in 3D and used to drive active commutation without torque measurements. ONIX can acquire from combinations of passive electrodes, Neuropixels probes, head-mounted microscopes, cameras, 3D-trackers, and other data sources. We used ONIX to perform uninterrupted, long (~7 hours) neural recordings in mice as they traversed complex 3-dimensional terrain. ONIX allowed exploration with similar mobility as non-implanted animals, in contrast to conventional tethered systems which restricted movement. By combining long recordings with full mobility, our technology will enable new progress on questions that require high-quality neural recordings during ethologically grounded behaviors.

Multi-regional module-based signal transmission in mouse visual cortex.

The visual cortex is hierarchically organized, yet 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 both between and within these modules, differences in layer and area distributions, and distinct temporal dynamics suggest that one module transmits 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.

CellExplorer: A framework for visualizing and characterizing single neurons.

The large diversity of neuron types provides the means by which cortical circuits perform complex operations. Neuron can be described by biophysical and molecular characteristics, afferent inputs, and neuron targets. To quantify, visualize, and standardize those features, we developed the open-source, MATLAB-based framework CellExplorer. It consists of three components: a processing module, a flexible data structure, and a powerful graphical interface. The processing module calculates standardized physiological metrics, performs neuron-type classification, finds putative monosynaptic connections, and saves them to a standardized, yet flexible, machine-readable format. The graphical interface makes it possible to explore the computed features at the speed of a mouse click. The framework allows users to process, curate, and relate their data to a growing public collection of neurons. CellExplorer can link genetically identified cell types to physiological properties of neurons collected across laboratories and potentially lead to interlaboratory standards of single-cell metrics.