Unsupervised spike sorting for large scale, high density multielectrode arrays

A new method for automated spike sorting for recordings with high density, large scale multielectrode arrays is presented. Exploiting the dense sampling of single neurons by multiple electrodes, we obtain an efficient, low-dimensional representation of detected spikes consisting of estimated spatial spike locations and dominant spike shape features, which enables fast and reliable clustering into single units. Millions of events can be sorted in minutes, and the method is parallelized and scales better than quadratically with the number of detected spikes. We demonstrate this method using recordings with a 4,096 channel array, and present validation based on anatomical imaging, optogenetic stimulation and model-based quality control. A comparison with semi-automated, shape-based spike sorting exposes significant limitations of conventional methods. Our analysis shows that it is feasible to reliably isolate the activity of hundreds to thousands of neurons in a single recording, and that dense, multi-channel probes substantially aid reliable spike sorting.