title: Applicability of independent component analysis on high-density microelectrode array recordings.
abstract: Emerging complementary metal oxide semiconductor (CMOS)-based, high-density microelectrode array (HD-MEA) devices provide high spatial resolution at subcellular level and a large number of readout channels. These devices allow for simultaneous recording of extracellular activity of a large number of neurons with every neuron being detected by multiple electrodes. To analyze the recorded signals, spiking events have to be assigned to individual neurons, a process referred to as "spike sorting." For a set of observed signals, which constitute a linear mixture of a set of source signals, independent component (IC) analysis (ICA) can be used to demix blindly the data and extract the individual source signals. This technique offers great potential to alleviate the problem of spike sorting in HD-MEA recordings, as it represents an unsupervised method to separate the neuronal sources. The separated sources or ICs then constitute estimates of single-neuron signals, and threshold detection on the ICs yields the sorted spike times. However, it is unknown to what extent extracellular neuronal recordings meet the requirements of ICA. In this paper, we evaluate the applicability of ICA to spike sorting of HD-MEA recordings. The analysis of extracellular neuronal signals, recorded at high spatiotemporal resolution, reveals that the recorded data cannot be modeled as a purely linear mixture. As a consequence, ICA fails to separate completely the neuronal signals and cannot be used as a stand-alone method for spike sorting in HD-MEA recordings. We assessed the demixing performance of ICA using simulated data sets and found that the performance strongly depends on neuronal density and spike amplitude. Furthermore, we show how postprocessing techniques can be used to overcome the most severe limitations of ICA. In combination with these postprocessing techniques, ICA represents a viable method to facilitate rapid spike sorting of multidimensional neuronal recordings.
authors.affiliation: Department of Biosystems Science and Engineering, ETH Zurich, Basel, Switzerland. david.jaeckel@bsse.ethz.ch
authors.forename: David
authors.initials: D
authors.lastname: Jäckel
authors.forename: Urs
authors.initials: U
authors.lastname: Frey
authors.forename: Michele
authors.initials: M
authors.lastname: Fiscella
authors.forename: Felix
authors.initials: F
authors.lastname: Franke
authors.forename: Andreas
authors.initials: A
authors.lastname: Hierlemann
mesh_terms: D000200:Action Potentials / Q000502:physiology*; D000465:Algorithms; D000818:Animals; D003198:Computer Simulation; D054707:Cyclic Nucleotide Phosphodiesterases, Type 6 / Q000235:genetics; D066298:In Vitro Techniques; D016014:Linear Models; D051379:Mice; D008817:Mice, Mutant Strains; D046228:Microarray Analysis; D008839:Microelectrodes; D008959:Models, Neurological; D017711:Nonlinear Dynamics; D025341:Principal Component Analysis; D012165:Retinal Ganglion Cells / Q000502:physiology*; D012815:Signal Processing, Computer-Assisted
chemical_list: D054707:Cyclic Nucleotide Phosphodiesterases, Type 6; C517438:Pde6b protein, mouse