@article {67232, title = {Neuropixels 2.0: A miniaturized high-density probe for stable, long-term brain recordings.}, journal = {Science}, volume = {372}, year = {2021}, month = {2021 Apr 16}, abstract = {

Measuring the dynamics of neural processing across time scales requires following the spiking of thousands of individual neurons over milliseconds and months. To address this need, we introduce the Neuropixels 2.0 probe together with newly designed analysis algorithms. The probe has more than 5000 sites and is miniaturized to facilitate chronic implants in small mammals and recording during unrestrained behavior. High-quality recordings over long time scales were reliably obtained in mice and rats in six laboratories. Improved site density and arrangement combined with newly created data processing methods enable automatic post hoc correction for brain movements, allowing recording from the same neurons for more than 2 months. These probes and algorithms enable stable recordings from thousands of sites during free behavior, even in small animals such as mice.

}, issn = {1095-9203}, doi = {10.1126/science.abf4588}, author = {Steinmetz, Nicholas A and Ayd{\i}n, {\c C}a{\u g}atay and Lebedeva, Anna and Okun, Michael and Pachitariu, Marius and Bauza, Marius and Beau, Maxime and Bhagat, Jai and B{\"o}hm, Claudia and Broux, Martijn and Chen, Susu and Colonell, Jennifer and Gardner, Richard J and Karsh, Bill and Kloosterman, Fabian and Kostadinov, Dimitar and Mora-Lopez, Carolina and O{\textquoteright}Callaghan, John and Park, Junchol and Putzeys, Jan and Sauerbrei, Britton and van Daal, Rik J J and Vollan, Abraham Z and Wang, Shiwei and Welkenhuysen, Marleen and Ye, Zhiwen and Dudman, Joshua T and Dutta, Barundeb and Hantman, Adam W and Harris, Kenneth D and Lee, Albert K and Moser, Edvard I and O{\textquoteright}Keefe, John and Renart, Alfonso and Svoboda, Karel and H{\"a}usser, Michael and Haesler, Sebastian and Carandini, Matteo and Harris, Timothy D} } @article {65615, title = {Arousal modulates retinal output.}, journal = {Neuron}, volume = {107}, year = {2020}, month = {2020 Aug 01}, pages = {487}, abstract = {

At various stages of the visual system, visual responses are modulated by arousal. Here, we find that in mice this modulation operates as early as in the first synapse from the retina and even in retinal axons. To measure retinal activity in the awake, intact brain, we imaged the synaptic boutons of retinal axons in the superior colliculus. Their activity depended not only on vision but also on running speed and pupil size, regardless of retinal illumination. Arousal typically reduced their visual responses and selectivity for direction and orientation. Recordings from retinal axons in the optic tract revealed that arousal modulates the firing of some retinal ganglion cells. Arousal had similar effects postsynaptically in colliculus neurons, independent of activity in the other main source of visual inputs to the colliculus, the primary visual cortex. These results indicate that arousal modulates activity at every stage of the mouse visual system.

}, issn = {1097-4199}, doi = {10.1016/j.neuron.2020.04.026}, author = {Schr{\"o}der, Sylvia and Steinmetz, Nicholas A and Krumin, Michael and Pachitariu, Marius and Rizzi, Matteo and Lagnado, Leon and Harris, Kenneth D and Carandini, Matteo} } @article {49534, title = {High-dimensional geometry of population responses in visual cortex.}, journal = {Nature}, volume = {571}, year = {2019}, month = {2019 Jun 26}, pages = {361-65}, abstract = {

A neuronal population encodes information most efficiently when its activity is uncorrelated and high-dimensional, and most robustly when its activity is correlated and lower-dimensional. Here, we analyzed the correlation structure of natural image coding, in large visual cortical populations recorded from awake mice. Evoked population activity was high dimensional, with correlations obeying an unexpected power-law: the n-th principal component variance scaled as 1/n. This was not inherited from the 1/f spectrum of natural images, because it persisted after stimulus whitening. We proved mathematically that the variance spectrum must decay at least this fast if a population code is smooth, i.e. if small changes in input cannot 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 represents a fundamental constraint governing correlations in neural population codes.

}, doi = {10.1038/s41586-019-1346-5}, author = {Stringer, Carsen and Pachitariu, Marius and Steinmetz, Nicholas A and Carandini, Matteo and Harris, Kenneth D} } @article {49187, title = {Spontaneous behaviors drive multidimensional, brain-wide population activity.}, journal = {Science}, volume = {364}, year = {2019}, month = {2019 Apr 18}, pages = {255}, abstract = {

Sensory cortices are active in the absence of external sensory stimuli. To understand the nature of this ongoing activity, we used two-photon calcium imaging to record from over 10,000 neurons in the visual cortex of mice awake in darkness while monitoring their behavior videographically. Ongoing population activity was multidimensional, exhibiting at least 100 significant dimensions, some of which were related to the spontaneous behaviors of the mice. The largest single dimension was correlated with the running speed and pupil area, while a 16-dimensional summary of orofacial behaviors could predict ~45\% of the explainable neural variance. Electrophysiological recordings with 8 simultaneous Neuropixels probes revealed a similar encoding of high-dimensional orofacial behaviors across multiple forebrain regions. Representation of motor variables continued uninterrupted during visual stimulus presentation, occupying dimensions nearly orthogonal to the stimulus responses. Our results show that a multidimensional representation of motor state is encoded across the forebrain, and is integrated with visual input by neuronal populations in primary visual cortex.

}, doi = {10.1101/306019}, author = {Stringer, Carsen and Pachitariu, Marius and Steinmetz, Nicholas A and Reddy, Charu Bai and Carandini, Matteo and Harris, Kenneth D} } @article {49345, title = {Robustness of spike deconvolution for neuronal calcium imaging.}, journal = {The Journal of Neuroscience : the official journal of the Society for Neuroscience}, volume = {38}, year = {2018}, month = {2018 Aug 06}, pages = {7976-85}, abstract = {

Calcium imaging is a powerful method to record the activity of neural populations in many species, but inferring spike times from calcium signals is a challenging problem. We compared multiple approaches using multiple datasets with ground truth electrophysiology, and found that simple non-negative deconvolution (NND) outperformed all other algorithms on out-of-sample test data. We introduce a novel benchmark applicable to recordings without electrophysiological ground truth, based on the correlation of responses to two stimulus repeats, and used this to show that unconstrained NND also outperformed the other algorithms when run on "zoomed out" datasets of \~{}10,000 cell recordings from the visual cortex of mice of either sex. Finally, we show that NND-based methods match the performance of a supervised method based on convolutional neural networks, while avoiding some of the biases of such methods, and at much faster running times. We therefore recommend that spikes be inferred from calcium traces using simple NND, due to its simplicity, efficiency and accuracy.The experimental method that currently allows for recordings of the largest numbers of cells simultaneously is two-photon calcium imaging. However, use of this powerful method requires that neuronal firing times be inferred correctly from the large resulting datasets. Previous studies have claimed that complex supervised learning algorithms outperform simple deconvolution methods at this task. Unfortunately, these studies suffered from several problems and biases. When we repeated the analysis, using the same data and correcting these problems, we found that simpler spike inference methods perform better. Even more importantly, we found that supervised learning methods can introduce artifactual structure into spike trains, that can in turn lead to erroneous scientific conclusions. Of the algorithms we evaluated, we found that an extremely simple method performed best in all circumstances tested, was much faster to run, and was insensitive to parameter choices, making incorrect scientific conclusions much less likely.

}, issn = {1529-2401}, doi = {10.1523/JNEUROSCI.3339-17.2018}, author = {Pachitariu, Marius and Stringer, Carsen and Harris, Kenneth D} } @article {48839, title = {Vision and locomotion shape the interactions between neuron types in mouse visual cortex.}, journal = {Neuron}, volume = {98}, year = {2018}, month = {2018 May 2}, pages = {602-15}, abstract = {

Cortical computation arises from the interaction of multiple neuronal\ types, including pyramidal (Pyr) cells and interneurons expressing Sst, Vip, or Pvalb. To study the circuit underlying such\ interactions, we imaged these four\ types\ of cells in\ mouse\ primary\ visual\ cortex(V1). Our recordings in darkness were consistent with a "disinhibitory" model in which\ locomotion\ activates Vip cells, thus inhibiting Sst cells and disinhibiting Pyr cells. However, the disinhibitory model failed when\ visual\ stimuli were present:\ locomotion\ increased Sst cell responses to large stimuli and Vip cell responses to small stimuli. A recurrent network model successfully predicted each cell type{\textquoteright}s activity from the measured activity of other\ types. Capturing the effects of\ locomotion, however, required allowing it to increase feedforward synaptic weights and modulate recurrent weights. This network model summarizes interneuron\ interactions\ and suggests that\ locomotion\ may alter cortical computation by changing effective synaptic connectivity.

}, doi = {https://doi.org/10.1101/058396}, author = {Dipoppa, Mario and Ranson, Adam and Krumin, Michael and Pachitariu, Marius and Carandini, Matteo and Harris, Kenneth D} } @article {48836, title = {Community-based benchmarking improves spike inference from two-photon calcium imaging data.}, journal = {PLoS Computational Biology}, volume = {14}, year = {2017}, month = {2018 May 21}, pages = {e1006157}, abstract = {

In recent years, two-photon calcium imaging has become a standard tool to probe the function of neural circuits and to study computations in neuronal populations. However, the acquired signal is only an indirect measurement of neural activity due to the comparatively slow dynamics of fluorescent calcium indicators. Different algorithms for estimating spike rates from noisy calcium measurements have been proposed in the past, but it is an open question how far performance can be improved. Here, we report the results of the spikefinder challenge, launched to catalyze the development of new spike rate inference algorithms through crowd-sourcing. We present ten of the submitted algorithms which show improved performance compared to previously evaluated methods. Interestingly, the top-performing algorithms are based on a wide range of principles from deep neural networks to generative models, yet provide highly correlated estimates of the neural activity. The competition shows that benchmark challenges can drive algorithmic developments in neuroscience.

}, doi = {10.1371/journal.pcbi.1006157}, author = {Berens, Philipp and Freeman, Jeremy and Deneux, Thomas and Chenkov, Nicolay and McColgan, Thomas and Speiser, Artur and Macke, Jakob H and Turaga, Srinivas C and Mineault, Patrick and Rupprecht, Peter and Gerhard, Stephan and Friedrich, Rainer W and Friedrich, Johannes and Paninski, Liam and Pachitariu, Marius and Harris, Kenneth D and Bolte, Ben and Machado, Timothy A and Ringach, Dario and etal} } @article {48935, title = {Fully integrated silicon probes for high-density recording of neural activity.}, journal = {Nature}, volume = {551}, year = {2017}, month = {2017 Nov 08}, pages = {232-236}, abstract = {

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 Ca(2+) 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 {\texttimes} 20-μm cross-section shank. The 6 {\texttimes} 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.

}, issn = {1476-4687}, doi = {10.1038/nature24636}, author = {Jun, James J and Steinmetz, Nicholas A and Siegle, Joshua H and Denman, Daniel J and Bauza, Marius and Barbarits, Brian and Lee, Albert K and Anastassiou, Costas A and Andrei, Alexandru and Ayd{\i}n, {\c C}a{\u g}atay and Barbic, Mladen and Blanche, Timothy J and Bonin, Vincent and Couto, Jo{\~a}o and Dutta, Barundeb and Gratiy, Sergey L and Gutnisky, Diego A and H{\"a}usser, Michael and Karsh, Bill and Ledochowitsch, Peter and Lopez, Carolina Mora and Mitelut, Catalin and Musa, Silke and Okun, Michael and Pachitariu, Marius and Putzeys, Jan and Rich, P Dylan and Rossant, Cyrille and Sun, Wei-Lung and Svoboda, Karel and Carandini, Matteo and Harris, Kenneth D and Koch, Christof and O{\textquoteright}Keefe, John and Harris, Timothy D} } @article {48837, title = {Robustness of spike deconvolution for calcium imaging of neural spiking.}, journal = {bioRxiv}, year = {2017}, month = {2017 Jun 27}, pages = {156786}, abstract = {

Calcium imaging is a powerful method to record the activity of neural populations, but inferring spike times from calcium signals is a challenging problem. We compared multiple approaches using multiple datasets with ground truth electrophysiology, and found that simple non-negative deconvolution (NND) outperformed all other algorithms. We introduce a novel benchmark applicable to recordings without electrophysiological ground truth, based on the correlation of responses to two stimulus repeats, and used this to show that unconstrained NND also outperformed the other algorithms when run on {\textquoteright}zoomed out{\textquoteright} datasets of ~10,000 cell recordings. Finally, we show that NND-based methods match the performance of a supervised method based on convolutional neural networks, while avoiding some of the biases of such methods, and at much faster running times. We therefore recommend that spikes be inferred from calcium traces using simple NND, due to its simplicity, efficiency and accuracy.

}, doi = {https://doi.org/10.1101/156786}, author = {Pachitariu, Marius and Stringer, Carsen and Harris, Kenneth D} } @article {48840, title = {Suite2p: beyond 10,000 neurons with standard two-photon microscopy.}, journal = {bioRxiv}, year = {2017}, month = {2017 Jul 20}, pages = {061507}, abstract = {

Two-photon microscopy of calcium-dependent sensors has enabled unprecedented recordings from vast populations of neurons. While the sensors and microscopes have matured over several generations of development, computational methods to process the resulting movies remain inefficient and can give results that are hard to interpret. Here we introduce Suite2p: a fast, accurate and complete pipeline that registers raw movies, detects active cells, extracts their calcium traces and infers their spike times. Suite2p runs on standard workstations, operates faster than real time, and recovers ~2 times more cells than the previous state-of-the-art method. Its low computational load allows routine detection of ~10,000 cells simultaneously with standard two-photon resonant-scanning microscopes. Recordings at this scale promise to reveal the fine structure of activity in large populations of neurons or large populations of subcellular structures such as synaptic boutons.

}, doi = {https://doi.org/10.1101/061507}, author = {Pachitariu, Marius and Stringer, Carsen and Dipoppa, Mario and Schr{\"o}der, Sylvia and Rossi, L Federico and Dalgleish, Henry and Carandini, Matteo and Harris, Kenneth D} } @conference {48832, title = {Fast and accurate spike sorting of high-channel count probes with KiloSort.}, booktitle = {Neural Information Processing Systems (NIPS 2016)}, year = {2016}, month = {2016 Dec 05}, address = {Barcelona, Spain}, abstract = {

New silicon technology is enabling large-scale electrophysiological recordings in vivo from hundreds to thousands of channels. Interpreting these recordings requires scalable and accurate automated methods for spike sorting, which should minimize the time required for manual curation of the results. Here we introduce KiloSort, a new integrated spike sorting framework that uses template matching both during spike detection and during spike clustering. KiloSort models the electrical voltage as a sum of template waveforms triggered on the spike times, which allows overlapping spikes to be identified and resolved. Unlike previous algorithms that compress the data with PCA, KiloSort operates on the raw data which allows it to construct a more accurate model of the waveforms. Processing times are faster than in previous algorithms thanks to batch-based optimization on GPUs. We compare KiloSort to an established algorithm and show favorable performance, at much reduced processing times. A novel post-clustering merging step based on the continuity of the templates further reduced substantially the number of manual operations required on this data, for the neurons with near-zero error rates, paving the way for fully automated spike sorting of multichannel electrode recordings.

}, url = {http://papers.nips.cc/paper/6326-fast-and-accurate-spike-sorting-of-high-channel-count-probes-with-kilosort}, author = {Pachitariu, Marius and Steinmetz, Nicholas A and Kadir, Shabnam N and Carandini, Matteo and Harris, Kenneth D} } @article {47812, title = {Improving data quality in neuronal population recordings.}, journal = {Nature Neuroscience}, volume = {19}, year = {2016}, month = {2016 Aug 26}, pages = {1165-74}, abstract = {

Understanding how the brain operates requires understanding how large sets of neurons function together. Modern recording technology makes it possible to simultaneously record the activity of hundreds of neurons, and technological developments will soon allow recording of thousands or tens of thousands. As with all experimental techniques, these methods are subject to confounds that complicate the interpretation of such recordings, and could lead to erroneous scientific conclusions. Here we discuss methods for assessing and improving the quality of data from these techniques and outline likely future directions in this field.

}, issn = {1546-1726}, doi = {10.1038/nn.4365}, author = {Harris, Kenneth D and Quiroga, Rodrigo Quian and Freeman, Jeremy and Smith, Spencer L} } @article {48842, title = {Inhibitory control of correlated intrinsic variability in cortical networks.}, journal = {eLife}, volume = {5}, year = {2016}, month = {2016 Dec 07}, abstract = {

Cortical networks exhibit intrinsic dynamics that drive coordinated, large-scale fluctuations across neuronal populations and create noise correlations that impact sensory coding. To investigate the network-level mechanisms that underlie these dynamics, we developed novel computational techniques to fit a deterministic spiking network model directly to multi-neuron recordings from different rodent species, sensory modalities, and behavioral states. The model generated correlated variability without external noise and accurately reproduced the diverse activity patterns in our recordings. Analysis of the model parameters suggested that differences in noise correlations across recordings were due primarily to differences in the strength of feedback inhibition. Further analysis of our recordings confirmed that putative inhibitory neurons were indeed more active during desynchronized cortical states with weak noise correlations. Our results demonstrate that network models with intrinsically-generated variability can accurately reproduce the activity patterns observed in multi-neuron recordings and suggest that inhibition modulates the interactions between intrinsic dynamics and sensory inputs to control the strength of noise correlations.

}, issn = {2050-084X}, doi = {10.7554/eLife.19695}, author = {Stringer, Carsen and Pachitariu, Marius and Steinmetz, Nicholas A and Okun, Michael and Bartho, Peter and Harris, Kenneth D and Sahani, Maneesh and Lesica, Nicholas A} } @article {65751, title = {Inhibitory control of correlated intrinsic variability in cortical networks}, journal = {Elife}, volume = {5}, year = {2016}, month = {12/2016}, pages = {e19695}, abstract = {

Cortical networks exhibit intrinsic dynamics that drive coordinated, large-scale fluctuations across neuronal populations and create noise correlations that impact sensory coding. To investigate the network-level mechanisms that underlie these dynamics, we developed novel computational techniques to fit a deterministic spiking network model directly to multi-neuron recordings from different rodent species, sensory modalities, and behavioral states. The model generated correlated variability without external noise and accurately reproduced the diverse activity patterns in our recordings. Analysis of the model parameters suggested that differences in noise correlations across recordings were due primarily to differences in the strength of feedback inhibition. Further analysis of our recordings confirmed that putative inhibitory neurons were indeed more active during desynchronized cortical states with weak noise correlations. Our results demonstrate that network models with intrinsically-generated variability can accurately reproduce the activity patterns observed in multi-neuron recordings and suggest that inhibition modulates the interactions between intrinsic dynamics and sensory inputs to control the strength of noise correlations.

}, doi = {https://doi.org/10.7554/eLife.19695}, author = {Stringer, Carsen and Pachitariu, Marius and Steinmetz, Nicholas A and Okun, Michael and Bartho, Peter and Harris, Kenneth D and Sahani, Maneesh and Lesica, Nicholas A} } @article {2346, title = {Theta-mediated dynamics of spatial information in hippocampus.}, journal = {The Journal of Neuroscience: The Official Journal of the Society for Neuroscience}, volume = {28}, year = {2008}, month = {2008 Jun 4}, pages = {5959-64}, abstract = {

In rodent hippocampus, neuronal activity is organized by a 6-10 Hz theta oscillation. The spike timing of hippocampal pyramidal cells with respect to the theta rhythm correlates with an animal{\textquoteright}s position in space. This correlation has been suggested to indicate an explicit temporal code for position. Alternatively, it may be interpreted as a byproduct of theta-dependent dynamics of spatial information flow in hippocampus. Here we show that place cell activity on different phases of theta reflects positions shifted into the future or past along the animal{\textquoteright}s trajectory in a two-dimensional environment. The phases encoding future and past positions are consistent across recorded CA1 place cells, indicating a coherent representation at the network level. Consistent theta-dependent time offsets are not simply a consequence of phase-position correlation (phase precession), because they are no longer seen after data randomization that preserves the phase-position relationship. The scale of these time offsets, 100-300 ms, is similar to the latencies of hippocampal activity after sensory input and before motor output, suggesting that offset activity may maintain coherent brain activity in the face of information processing delays.

}, keywords = {Action Potentials, Animals, Hippocampus, Long-Evans, Male, Rats, Spatial Behavior, Theta Rhythm}, issn = {1529-2401}, doi = {10.1523/JNEUROSCI.3773-10.2011}, author = {Itskov, Vladimir and Pastalkova, Eva and Mizuseki, Kenji and Buzs{\'a}ki, Gy{\"o}rgy and Harris, Kenneth D} }