Automatic online spike sorting with singular value decomposition and fuzzy C-mean clustering

BackgroundUnderstanding how neurons contribute to perception, motor functions and cognition requires the reliable detection of spiking activity of individual neurons during a number of different experimental conditions. An important problem in computational neuroscience is thus to develop algorithms to automatically detect and sort the spiking activity of individual neurons from extracellular recordings. While many algorithms for spike sorting exist, the problem of accurate and fast online sorting still remains a challenging issue.ResultsHere we present a novel software tool, called FSPS (Fuzzy SPike Sorting), which is designed to optimize: (i) fast and accurate detection, (ii) offline sorting and (iii) online classification of neuronal spikes with very limited or null human intervention. The method is based on a combination of Singular Value Decomposition for fast and highly accurate pre-processing of spike shapes, unsupervised Fuzzy C-mean, high-resolution alignment of extracted spike waveforms, optimal selection of the number of features to retain, automatic identification the number of clusters, and quantitative quality assessment of resulting clusters independent on their size. After being trained on a short testing data stream, the method can reliably perform supervised online classification and monitoring of single neuron activity. The generalized procedure has been implemented in our FSPS spike sorting software (available free for non-commercial academic applications at the address: http://www.spikesorting.com) using LabVIEW (National Instruments, USA). We evaluated the performance of our algorithm both on benchmark simulated datasets with different levels of background noise and on real extracellular recordings from premotor cortex of Macaque monkeys. The results of these tests showed an excellent accuracy in discriminating low-amplitude and overlapping spikes under strong background noise. The performance of our method is competitive with respect to other robust spike sorting algorithms.ConclusionsThis new software provides neuroscience laboratories with a new tool for fast and robust online classification of single neuron activity. This feature could become crucial in situations when online spike detection from multiple electrodes is paramount, such as in human clinical recordings or in brain-computer interfaces.

[1]  Kyung Hwan Kim.  Improved Algorithm for Fully-automated Neural Spike Sorting based on Projection Pursuit and Gaussian Mixture Model , 2006 .

[2]  Gordon Pipa,et al.  Single-Unit Recordings Revisited: Activity in Recurrent Microcircuits , 2005, ICANN.

[3]  William F. Christensen,et al.  Methods of Multivariate Analysis: Rencher/Methods , 2012 .

[4]  Susan T. Dumais,et al.  Using Linear Algebra for Intelligent Information Retrieval , 1995, SIAM Rev..

[5]  Steven J. Cox,et al.  Mathematics for Neuroscientists , 2010 .

[6]  Nicolas Fourcaud-Trocmé,et al.  OpenElectrophy: An Electrophysiological Data- and Analysis-Sharing Framework , 2008, Front. Neuroinform.

[7]  Azriel Rosenfeld,et al.  Computer Vision , 1988, Adv. Comput..

[8]  Ferial El-Hawary,et al.  The Electrical Engineering Handbook Series , 2004 .

[9]  R. Bro,et al.  Centering and scaling in component analysis , 2003 .

[10]  C Rosen,et al.  Supervisory control of wastewater treatment plants by combining principal component analysis and fuzzy c-means clustering. , 2001, Water science and technology : a journal of the International Association on Water Pollution Research.

[11]  Jason W. Osborne,et al.  Best practices in exploratory factor analysis: four recommendations for getting the most from your analysis. , 2005 .

[12]  Iyad Obeid,et al.  Fuzzy logic-based spike sorting system , 2011, Journal of Neuroscience Methods.

[13]  James C. Hayton,et al.  Factor Retention Decisions in Exploratory Factor Analysis: a Tutorial on Parallel Analysis , 2004 .

[14]  P. M. Horton,et al.  Spike sorting based upon machine learning algorithms (SOMA) , 2007, Journal of Neuroscience Methods.

[15]  G. Buzsáki.  Large-scale recording of neuronal ensembles , 2004, Nature Neuroscience.

[16]  R. Quiroga,et al.  Extracting information from neuronal populations : information theory and decoding approaches , 2022 .

[17]  Artur Luczak,et al.  Spectral representation—analyzing single-unit activity in extracellularly recorded neuronal data without spike sorting , 2005, Journal of Neuroscience Methods.

[18]  Daniela Calvetti,et al.  Matrix methods in data mining and pattern recognition , 2009, Math. Comput..

[19]  D. W. Scott.  On optimal and data based histograms , 1979 .

[20]  Werner Dubitzky,et al.  A Practical Approach to Microarray Data Analysis , 2003, Springer US.

[21]  Linda G. Shapiro,et al.  Computer Vision , 2001 .

[22]  Michael S. Dudzic,et al.  An industrial perspective on implementing on-line applications of multivariate statistics , 2004 .

[23]  I. Jolliffe.  Principal Component Analysis , 2002 .

[24]  K J Burchiel,et al.  Computer analysis of the tonic, phasic, and kinesthetic activity of pallidal discharges in Parkinson patients. , 1999, Surgical neurology.

[25]  R. K. Snider,et al.  Classification of non-stationary neural signals , 1998, Journal of Neuroscience Methods.

[26]  Shengrui Wang,et al.  FCM-Based Model Selection Algorithms for Determining the Number of Clusters , 2004, Pattern Recognit..

[27]  T. Suzuki,et al.  Fast Automatic Template Matching for Spike Sorting Based on Davies-Bouldin Validation Indices , 2007, 2007 29th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[28]  J. Bezdek.  Cluster Validity with Fuzzy Sets , 1973 .

[29]  Alex Zviagintsev,et al.  Algorithms and architectures for low power spike detection and alignment , 2006, Journal of neural engineering.

[30]  A. C. Rencher.  Methods of multivariate analysis , 1995 .

[31]  Matthew Brand,et al.  Fast Online SVD Revisions for Lightweight Recommender Systems , 2003, SDM.

[32]  J. C. Dunn,et al.  A Fuzzy Relative of the ISODATA Process and Its Use in Detecting Compact Well-Separated Clusters , 1973 .

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

[34]  Taejeong Kim,et al.  Solving alignment problems in neural spike sorting using frequency domain PCA , 2006, Neurocomputing.

[35]  James C. Bezdek,et al.  Efficient Implementation of the Fuzzy c-Means Clustering Algorithms , 1986, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[36]  F. Gabbiani,et al.  Comprar Mathematics For Neuroscientists | F. Gabbiani | 9780123748829 | Academic Press , 2010 .

[37]  Timothy J. Blanche,et al.  Python for Large-Scale Electrophysiology , 2009, Front. Neuroinform..

[38]  J. Letelier,et al.  Spike sorting based on discrete wavelet transform coefficients , 2000, Journal of Neuroscience Methods.

[39]  Pekka Teppola,et al.  Adaptive Fuzzy C-Means clustering in process monitoring , 1999 .

[40]  G. Rizzolatti,et al.  Somatotopic representation in inferior area 6 of the macaque monkey. , 1989, Brain, behavior and evolution.

[41]  Byoung Gu Min,et al.  Comparative Analysis of Texture Characteristics of Malignant and Benign Tumors in Breast Ultrasonograms , 2010, Journal of Digital Imaging.

[42]  Jadin C. Jackson,et al.  Quantitative measures of cluster quality for use in extracellular recordings , 2005, Neuroscience.

[43]  V. J. Rayward-Smith,et al.  Fuzzy Cluster Analysis: Methods for Classification, Data Analysis and Image Recognition , 1999 .

[44]  BMC Neuroscience , 2003 .

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

[46]  Y. Ritov,et al.  Failure in identification of overlapping spikes from multiple neuron activity causes artificial correlations , 2001, Journal of Neuroscience Methods.

[47]  J. Bezdek,et al.  FCM: The fuzzy c-means clustering algorithm , 1984 .

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

[49]  Vijay K. Madisetti.  Wireless, Networking, Radar, Sensor Array Processing, and Nonlinear Signal Processing , 2009 .

[50]  James C. Bezdek,et al.  Pattern Recognition with Fuzzy Objective Function Algorithms , 1981, Advanced Applications in Pattern Recognition.

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

[52]  John C. Nash,et al.  Simple Algorithms for the Partial Singular Value Decomposition , 1987, Comput. J..

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

[54]  In-Young Kim,et al.  Adaptable fuzzy C-Means for improved classification as a preprocessing procedure of brain parcellation , 2010, Journal of Digital Imaging.

[55]  R. Quiroga.  Spike sorting , 2012, Current Biology.

[56]  G H Ball,et al.  A clustering technique for summarizing multivariate data. , 1967, Behavioral science.

[57]  Theodore A. Walls,et al.  Non-Graphical Solutions for Cattell’s Scree Test , 2013 .

[58]  R. Kass,et al.  Multiple neural spike train data analysis: state-of-the-art and future challenges , 2004, Nature Neuroscience.

[59]  C. M. Stewart,et al.  Spike detection, characterization, and discrimination using feature analysis software written in LabVIEW , 2004, Comput. Methods Programs Biomed..

[60]  Ramon Gutiérrez-Castrejón,et al.  Using LabVIEWTM for advanced nonlinear optoelectronic device simulations in high-speed optical communications , 2006, Comput. Phys. Commun..