A fully automated spike sorting algorithm using t-distributed neighbor embedding and density based clustering

In this study, a new spike sorting method was developed based on a combination of two methods, t-Distributed Stochastic Neighbor Embedding (t-SNE) and Density-Based Spatial Clustering of Applications with Noise (DBSCAN). Parameters of both methods were simultaneously optimized using a Genetic Algorithm (GA) using a simulated dataset containing 2 to 20 simultaneously recorded neurons. The performance of this method was evaluated using both a stimulated dataset as well as real multichannel electrophysiological data. The results indicated that our fully automated algorithm using t-SNE-DBSCAN outperforms other state-of-the-art algorithms and human experts in spike sorting especially when there are a large number of simultaneously recorded units. Our algorithm also determines the noise waveforms and has an overall high sensitivity, precision and accuracy for correctly classifying waveforms belonging to each neuron (all >90%) without the need for manual corrections afterwards. Our method can be a crucial part of the analysis pipeline in particular when manual sorting of units is becoming prohibitive due to the sheer number of recorded neurons per session.

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

[2]  Pierre Yger,et al.  Fast and accurate spike sorting in vitro and in vivo for up to thousands of electrodes , 2016, bioRxiv.

[3]  Hans-Peter Kriegel,et al.  A Density-Based Algorithm for Discovering Clusters in Large Spatial Databases with Noise , 1996, KDD.

[4]  Luciano Fadiga,et al.  Automatic online spike sorting with singular value decomposition and fuzzy C-mean clustering , 2012, BMC Neuroscience.

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

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

[7]  Rodrigo Quian Quiroga,et al.  Past, present and future of spike sorting techniques , 2015, Brain Research Bulletin.

[8]  Daniel Novak,et al.  Performance comparison of extracellular spike sorting algorithms for single-channel recordings , 2012, Journal of Neuroscience Methods.

[9]  Valérie Ventura,et al.  To sort or not to sort: the impact of spike-sorting on neural decoding performance , 2014, Journal of neural engineering.

[10]  Winnie Jensen,et al.  Introduction to Neural Engineering for Motor Rehabilitation , 2013 .

[11]  Geoffrey E. Hinton,et al.  Visualizing Data using t-SNE , 2008 .

[12]  Liam Paninski,et al.  Kalman Filter Mixture Model for Spike Sorting of Non-stationary Data , 2010 .

[13]  Ali Ghazizadeh,et al.  Prefrontal Cortex Mediates Extinction of Responding by Two Distinct Neural Mechanisms in Accumbens Shell , 2012, The Journal of Neuroscience.

[14]  Kenneth D Harris,et al.  Spike sorting for large, dense electrode arrays , 2015, Nature Neuroscience.

[15]  Sung June Kim,et al.  A wavelet-based method for action potential detection from extracellular neural signal recording with low signal-to-noise ratio , 2003, IEEE Transactions on Biomedical Engineering.

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

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

[18]  Yuzo Takayama,et al.  Efficient sequential Bayesian inference method for real-time detection and sorting of overlapped neural spikes , 2013, Journal of Neuroscience Methods.

[19]  Eero P. Simoncelli,et al.  Journal of Neuroscience Methods , 2022 .

[20]  Klaus Obermayer,et al.  An online spike detection and spike classification algorithm capable of instantaneous resolution of overlapping spikes , 2009, Journal of Computational Neuroscience.

[21]  Magdy A. Bayoumi,et al.  Implementable Spike Sorting techniques for VLSI wireless BCI/BMI implants: A survey , 2015, 5th International Conference on Energy Aware Computing Systems & Applications.

[22]  J. O’Keefe.  Place units in the hippocampus of the freely moving rat , 1976, Experimental Neurology.

[23]  I. Johnstone,et al.  Ideal spatial adaptation by wavelet shrinkage , 1994 .

[24]  David J. Field,et al.  Sparse coding with an overcomplete basis set: A strategy employed by V1? , 1997, Vision Research.

[25]  David B. Dunson,et al.  Multichannel Electrophysiological Spike Sorting via Joint Dictionary Learning and Mixture Modeling , 2013, IEEE Transactions on Biomedical Engineering.

[26]  Matias J. Ison,et al.  Realistic simulation of extracellular recordings , 2009, Journal of Neuroscience Methods.

[27]  Edvard I. Moser,et al.  Grid Cells and Neural Coding in High-End Cortices , 2013, Neuron.

[28]  Dimitrios A. Adamos,et al.  Performance evaluation of PCA-based spike sorting algorithms , 2008, Comput. Methods Programs Biomed..

[29]  Konrad P Kording,et al.  How advances in neural recording affect data analysis , 2011, Nature Neuroscience.

[30]  Matthew Fellows,et al.  On the variability of manual spike sorting , 2004, IEEE Transactions on Biomedical Engineering.

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

[32]  M S Lewicki,et al.  A review of methods for spike sorting: the detection and classification of neural action potentials. , 1998, Network.

[33]  Dejan Markovic,et al.  Spike Sorting: The First Step in Decoding the Brain: The first step in decoding the brain , 2012, IEEE Signal Processing Magazine.

[34]  John P. Donoghue,et al.  Automated spike sorting using density grid contour clustering and subtractive waveform decomposition , 2007, Journal of Neuroscience Methods.

[35]  C. Koch,et al.  The origin of extracellular fields and currents — EEG, ECoG, LFP and spikes , 2012, Nature Reviews Neuroscience.

[36]  Edward R. Dougherty,et al.  Genomic Signal Processing [Life Sciences] , 2012, IEEE Signal Process. Mag..

[37]  Rodrigo Quian Quiroga,et al.  How many neurons can we see with current spike sorting algorithms? , 2012, Journal of Neuroscience Methods.