Fast and accurate spike sorting in vitro and in vivo for up to thousands of electrodes

Understanding how assemblies of neurons encode information requires recording large populations of cells in the brain. In recent years, multi-electrode arrays and large silicon probes have been developed to record simultaneously from hundreds or thousands of electrodes packed with a high density. However, these new devices challenge the classical way to do spike sorting. Here we developed a new method to solve these issues, based on a highly automated algorithm to extract spikes from extracellular data, and show that this algorithm reached near optimal performance both in vitro and in vivo. The algorithm is composed of two main steps: 1) a “template-finding” phase to extract the cell templates, i.e. the pattern of activity evoked over many electrodes when one neuron fires an action potential; 2) a “template-matching” phase where the templates were matched to the raw data to find the location of the spikes. The manual intervention by the user was reduced to the minimal, and the time spent on manual curation did not scale with the number of electrodes. We tested our algorithm with large-scale data from in vitro and in vivo recordings, from 32 to 4225 electrodes. We performed simultaneous extracellular and patch recordings to obtain “ground truth” data, i.e. cases where the solution to the sorting problem is at least partially known. The performance of our algorithm was always close to the best expected performance. We thus provide a general solution to sort spikes from large-scale extracellular recordings.

[1]  Lisandro Dalcin,et al.  Parallel distributed computing using Python , 2011 .

[2]  Markus Meister,et al.  Multi-neuronal signals from the retina: acquisition and analysis , 1994, Journal of Neuroscience Methods.

[3]  A.M. Litke,et al.  What does the eye tell the brain?: Development of a system for the large scale recording of retinal output activity , 2003, 2003 IEEE Nuclear Science Symposium. Conference Record (IEEE Cat. No.03CH37515).

[4]  Gaute T. Einevoll,et al.  ViSAPy: A Python tool for biophysics-based generation of virtual spiking activity for evaluation of spike-sorting algorithms , 2015, Journal of Neuroscience Methods.

[5]  Jan Müller,et al.  High-density microelectrode array recordings and real-time spike sorting for closed-loop experiments: an emerging technology to study neural plasticity , 2012, Front. Neural Circuits.

[6]  Stefan Rotter,et al.  Higher-Order Statistics of Input Ensembles and the Response of Simple Model Neurons , 2003, Neural Computation.

[7]  Andreas Hierlemann,et al.  Recording from defined populations of retinal ganglion cells using a high-density CMOS-integrated microelectrode array with real-time switchable electrode selection , 2012, Journal of Neuroscience Methods.

[8]  Sean Hughes,et al.  Clustering by Fast Search and Find of Density Peaks , 2016 .

[9]  Rongfang Bie,et al.  Clustering by fast search and find of density peaks via heat diffusion , 2016, Neurocomputing.

[10]  Luca Berdondini,et al.  Unsupervised spike sorting for large scale, high density multielectrode arrays , 2016 .

[11]  Alessandro Laio,et al.  Clustering by fast search and find of density peaks , 2014, Science.

[12]  Emad N. Eskandar,et al.  Neural signal processing and closed-loop control algorithm design for an implanted neural recording and stimulation system , 2015, 2015 37th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC).

[13]  Kenneth D Harris,et al.  Towards reliable spike-train recordings from thousands of neurons with multielectrodes , 2012, Current Opinion in Neurobiology.

[14]  Stéphane Mallat,et al.  Matching pursuits with time-frequency dictionaries , 1993, IEEE Trans. Signal Process..

[15]  Eero P. Simoncelli,et al.  A Model-Based Spike Sorting Algorithm for Removing Correlation Artifacts in Multi-Neuron Recordings , 2013, PloS one.

[16]  Jason S. Prentice,et al.  Fast, Scalable, Bayesian Spike Identification for Multi-Electrode Arrays , 2010, PloS one.

[17]  György Buzsáki,et al.  Neural Syntax: Cell Assemblies, Synapsembles, and Readers , 2010, Neuron.

[18]  A. Lambacher,et al.  Electrical imaging of neuronal activity by multi-transistor-array (MTA) recording at 7.8 μm resolution , 2004 .

[19]  Gilles Laurent,et al.  Using noise signature to optimize spike-sorting and to assess neuronal classification quality , 2002, Journal of Neuroscience Methods.

[20]  Nicholas V. Swindale,et al.  Spike sorting for polytrodes: a divide and conquer approach , 2014, Front. Syst. Neurosci..

[21]  Cyrille Rossant,et al.  Spike sorting for large, dense electrode arrays , 2015 .

[22]  D. Hubel Tungsten Microelectrode for Recording from Single Units. , 1957, Science.

[23]  Matteo Carandini,et al.  Kilosort: realtime spike-sorting for extracellular electrophysiology with hundreds of channels , 2016, bioRxiv.

[24]  Tim Gollisch,et al.  From response to stimulus: adaptive sampling in sensory physiology , 2007, Current Opinion in Neurobiology.

[25]  Daniel N Hill,et al.  Quality Metrics to Accompany Spike Sorting of Extracellular Signals , 2011, The Journal of Neuroscience.

[26]  Thomas Wachtler,et al.  Unsupervised neural spike sorting for high-density microelectrode arrays with convolutive independent component analysis , 2016, Journal of Neuroscience Methods.

[27]  Michael J. Berry,et al.  Mapping a Complete Neural Population in the Retina , 2012, The Journal of Neuroscience.

[28]  R. Quian Quiroga,et al.  Unsupervised Spike Detection and Sorting with Wavelets and Superparamagnetic Clustering , 2004, Neural Computation.

[29]  Armin Lambacher,et al.  Axonal Transmission in the Retina Introduces a Small Dispersion of Relative Timing in the Ganglion Cell Population Response , 2011, PloS one.

[30]  Pedro Barquinha,et al.  Validating silicon polytrodes with paired juxtacellular recordings: method and dataset , 2016, bioRxiv.

[31]  Andreas Moller,et al.  A CMOS-based sensor array for in-vitro neural tissue interfacing with 4225 recording sites and 1024 stimulation sites , 2014, 2014 IEEE Biomedical Circuits and Systems Conference (BioCAS) Proceedings.

[32]  Kenneth D. Harris,et al.  Hardware-accelerated interactive data visualization for neuroscience in Python , 2013, Front. Neuroinform..

[33]  Luca Berdondini,et al.  Unsupervised Spike Sorting for Large-Scale, High-Density Multielectrode Arrays. , 2017, Cell reports.

[34]  L. Berdondinia,et al.  High-density electrode array for imaging in vitro electrophysiological activity , 2005 .

[35]  Michael J. Berry,et al.  Recording spikes from a large fraction of the ganglion cells in a retinal patch , 2004, Nature Neuroscience.

[36]  J. Csicsvari,et al.  Accuracy of tetrode spike separation as determined by simultaneous intracellular and extracellular measurements. , 2000, Journal of neurophysiology.

[37]  Kenneth D. Harris,et al.  High-Dimensional Cluster Analysis with the Masked EM Algorithm , 2013, Neural Computation.

[38]  Michael S. Lewicki,et al.  Bayesian Modeling and Classification of Neural Signals , 1993, Neural Computation.