ViSAPy: A Python tool for biophysics-based generation of virtual spiking activity for evaluation of spike-sorting algorithms

BACKGROUND New, silicon-based multielectrodes comprising hundreds or more electrode contacts offer the possibility to record spike trains from thousands of neurons simultaneously. This potential cannot be realized unless accurate, reliable automated methods for spike sorting are developed, in turn requiring benchmarking data sets with known ground-truth spike times. NEW METHOD We here present a general simulation tool for computing benchmarking data for evaluation of spike-sorting algorithms entitled ViSAPy (Virtual Spiking Activity in Python). The tool is based on a well-established biophysical forward-modeling scheme and is implemented as a Python package built on top of the neuronal simulator NEURON and the Python tool LFPy. RESULTS ViSAPy allows for arbitrary combinations of multicompartmental neuron models and geometries of recording multielectrodes. Three example benchmarking data sets are generated, i.e., tetrode and polytrode data mimicking in vivo cortical recordings and microelectrode array (MEA) recordings of in vitro activity in salamander retinas. The synthesized example benchmarking data mimics salient features of typical experimental recordings, for example, spike waveforms depending on interspike interval. COMPARISON WITH EXISTING METHODS ViSAPy goes beyond existing methods as it includes biologically realistic model noise, synaptic activation by recurrent spiking networks, finite-sized electrode contacts, and allows for inhomogeneous electrical conductivities. ViSAPy is optimized to allow for generation of long time series of benchmarking data, spanning minutes of biological time, by parallel execution on multi-core computers. CONCLUSION ViSAPy is an open-ended tool as it can be generalized to produce benchmarking data or arbitrary recording-electrode geometries and with various levels of complexity.

[1]  D. Robinson,et al.  The electrical properties of metal microelectrodes , 1968 .

[2]  Janet Wiles,et al.  Action Potential Waveform Variability Limits Multi-Unit Separation in Freely Behaving Rats , 2012, PloS one.

[3]  Gaute T. Einevoll,et al.  Estimation of population firing rates and current source densities from laminar electrode recordings , 2008, Journal of Computational Neuroscience.

[4]  Z. Somogyvári,et al.  Model-based source localization of extracellular action potentials , 2005, Journal of Neuroscience Methods.

[5]  Partha P. Mitra,et al.  Automatic sorting of multiple unit neuronal signals in the presence of anisotropic and non-Gaussian variability , 1996, Journal of Neuroscience Methods.

[6]  J. Fermaglich.  Electric Fields of the Brain: The Neurophysics of EEG , 1982 .

[7]  Leslie S. Smith,et al.  A tool for synthesizing spike trains with realistic interference , 2007, Journal of Neuroscience Methods.

[8]  Sergei Vassilvitskii,et al.  k-means++: the advantages of careful seeding , 2007, SODA '07.

[9]  W. Rall.  Electrophysiology of a dendritic neuron model. , 1962, Biophysical journal.

[10]  D H Brainard,et al.  The Psychophysics Toolbox. , 1997, Spatial vision.

[11]  Klaus Obermayer,et al.  An automated online positioning system and simulation environment for multi-electrodes in extracellular recordings , 2010, 2010 Annual International Conference of the IEEE Engineering in Medicine and Biology.

[12]  H. Markram,et al.  Morphological, Electrophysiological, and Synaptic Properties of Corticocallosal Pyramidal Cells in the Neonatal Rat Neocortex , 2006 .

[13]  Anders M. Dale,et al.  Handbook of Neural Activity Measurement: Extracellular spikes and CSD , 2012 .

[14]  Moritz Helias,et al.  Correlations in spiking neuronal networks with distance dependent connections , 2009, Journal of Computational Neuroscience.

[15]  István Ulbert,et al.  Multiple microelectrode-recording system for human intracortical applications , 2001, Journal of Neuroscience Methods.

[16]  Marc-Oliver Gewaltig,et al.  NEST (NEural Simulation Tool) , 2007, Scholarpedia.

[17]  D. Kleinfeld,et al.  Variability of extracellular spike waveforms of cortical neurons. , 1996, Journal of neurophysiology.

[18]  Mats G. Larson,et al.  The Finite Element , 2013 .

[19]  Henry C. Tuckwell,et al.  A spatial stochastic neuronal model with Ornstein–Uhlenbeck input current , 2002, Biological Cybernetics.

[20]  Nikos K. Logothetis,et al.  Physiological Foundations of Neural Signals , 2013 .

[21]  P. T. Thorbergsson,et al.  Statistical modelling of spike libraries for simulation of extracellular recordings in the cerebellum , 2010, 2010 Annual International Conference of the IEEE Engineering in Medicine and Biology.

[22]  Daryl R Kipke,et al.  Theoretical analysis of intracortical microelectrode recordings , 2011, Journal of neural engineering.

[23]  R. Mazo.  On the theory of brownian motion , 1973 .

[24]  J. B. Ranck,et al.  Studies on single neurons in dorsal hippocampal formation and septum in unrestrained rats. I. Behavioral correlates and firing repertoires. , 1973, Experimental neurology.

[25]  Klas H. Pettersen,et al.  Modeling the Spatial Reach of the LFP , 2011, Neuron.

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

[27]  Klas H. Pettersen,et al.  Amplitude variability and extracellular low-pass filtering of neuronal spikes. , 2008, Biophysical journal.

[28]  J. Csicsvari,et al.  Intracellular features predicted by extracellular recordings in the hippocampus in vivo. , 2000, Journal of neurophysiology.

[29]  C. Koch,et al.  On the origin of the extracellular action potential waveform: A modeling study. , 2006, Journal of neurophysiology.

[30]  Marc-Oliver Gewaltig,et al.  NEST (NEural Simulation Tool) , 2007, Scholarpedia.

[31]  Michael W. Reimann,et al.  A Biophysically Detailed Model of Neocortical Local Field Potentials Predicts the Critical Role of Active Membrane Currents , 2013, Neuron.

[32]  Craig S. Henriquez,et al.  Confounded spikes generated by synchrony within neural tissue models , 2005, Neurocomputing.

[33]  Eero P. Simoncelli,et al.  Spatio-temporal correlations and visual signalling in a complete neuronal population , 2008, Nature.

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

[35]  Martin Garwicz,et al.  Computationally efficient simulation of extracellular recordings with multielectrode arrays , 2012, Journal of Neuroscience Methods.

[36]  T. Sejnowski,et al.  Fluctuating synaptic conductances recreate in vivo-like activity in neocortical neurons , 2001, Neuroscience.

[37]  A. Lambacher,et al.  Identifying firing mammalian neurons in networks with high-resolution multi-transistor array (MTA) , 2011 .

[38]  Nicolas Brunel,et al.  Dynamics of Sparsely Connected Networks of Excitatory and Inhibitory Spiking Neurons , 2000, Journal of Computational Neuroscience.

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

[40]  Markus Diesmann,et al.  Dependence of Spike-Count Correlations on Spike-Train Statistics and Observation Time Scale , 2010 .

[41]  Roger W. Hockney.  The science of computer benchmarking , 1995, Software, environments, tools.

[42]  Christof Koch,et al.  Using extracellular action potential recordings to constrain compartmental models , 2007, Journal of Computational Neuroscience.

[43]  Gaute T. Einevoll,et al.  LFPy: a tool for biophysical simulation of extracellular potentials generated by detailed model neurons , 2014, Front. Neuroinform..

[44]  Stefano Panzeri,et al.  Modelling and analysis of local field potentials for studying the function of cortical circuits , 2013, Nature Reviews Neuroscience.

[45]  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.

[46]  Arnd Roth,et al.  Automated optimization of a reduced layer 5 pyramidal cell model based on experimental data , 2012, Journal of Neuroscience Methods.

[47]  Andreas Hierlemann,et al.  Switch-Matrix-Based High-Density Microelectrode Array in CMOS Technology , 2010, IEEE Journal of Solid-State Circuits.

[48]  Gaute T. Einevoll,et al.  Intrinsic dendritic filtering gives low-pass power spectra of local field potentials , 2010, Journal of Computational Neuroscience.

[49]  R. Barrera,et al.  Point charge in a three‐dielectric medium with planar interfaces , 1978 .

[50]  Christophe Pouzat,et al.  Efficient spike-sorting of multi-state neurons using inter-spike intervals information , 2005, Journal of Neuroscience Methods.

[51]  Matthew J Nelson,et al.  Do electrode properties create a problem in interpreting local field potential recordings? , 2010, Journal of neurophysiology.

[52]  S. Palmer.  Vision Science : Photons to Phenomenology , 1999 .

[53]  Henry Markram,et al.  A Novel Multiple Objective Optimization Framework for Constraining Conductance-Based Neuron Models by Experimental Data , 2007, Front. Neurosci..

[54]  Michael L. Hines,et al.  Neuroinformatics Original Research Article Neuron and Python , 2022 .

[55]  Michael J. Berry,et al.  Functional organization of ganglion cells in the salamander retina. , 2006, Journal of neurophysiology.

[56]  R. K. Ursem.  Multi-objective Optimization using Evolutionary Algorithms , 2009 .

[57]  N. Logothetis,et al.  In Vivo Measurement of Cortical Impedance Spectrum in Monkeys: Implications for Signal Propagation , 2007, Neuron.

[58]  Naftali Tishby,et al.  Multi-Electrode Spike Sorting by Clustering Transfer Functions , 1998, NIPS.

[59]  T. Insel,et al.  The NIH BRAIN Initiative , 2013, Science.

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

[61]  Nicholas T. Carnevale,et al.  ModelDB: A Database to Support Computational Neuroscience , 2004, Journal of Computational Neuroscience.

[62]  B. McNaughton,et al.  Tetrodes markedly improve the reliability and yield of multiple single-unit isolation from multi-unit recordings in cat striate cortex , 1995, Journal of Neuroscience Methods.

[63]  Jonathan D. Victor,et al.  Dipole characterization of single neurons from their extracellular action potentials , 2011, Journal of Computational Neuroscience.

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

[65]  Ulrich G. Hofmann,et al.  Test of spike-sorting algorithms on the basis of simulated network data , 2002, Neurocomputing.

[66]  G. Buzsáki,et al.  Temporal Interaction between Single Spikes and Complex Spike Bursts in Hippocampal Pyramidal Cells , 2001, Neuron.

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

[68]  Andreas Hierlemann,et al.  Growing Cells Atop Microelectronic Chips: Interfacing Electrogenic Cells In Vitro With CMOS-Based Microelectrode Arrays , 2011, Proceedings of the IEEE.

[69]  A. Destexhe,et al.  Impact of network activity on the integrative properties of neocortical pyramidal neurons in vivo. , 1999, Journal of neurophysiology.

[70]  T. Blanche,et al.  Polytrodes: high-density silicon electrode arrays for large-scale multiunit recording. , 2005, Journal of neurophysiology.

[71]  Jeffrey D. Schall,et al.  Review of signal distortion through metal microelectrode recording circuits and filters , 2008, Journal of Neuroscience Methods.

[72]  D J Field,et al.  Relations between the statistics of natural images and the response properties of cortical cells. , 1987, Journal of the Optical Society of America. A, Optics and image science.

[73]  Luis A. Camuñas-Mesa,et al.  A Detailed and Fast Model of Extracellular Recordings , 2013, Neural Computation.

[74]  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.

[75]  Moritz Helias,et al.  How pattern formation in ring networks of excitatory and inhibitory spiking neurons depends on the input current regime , 2013, Front. Comput. Neurosci..

[76]  Christof Koch,et al.  High-amplitude positive spikes recorded extracellularly in cat visual cortex. , 2009, Journal of neurophysiology.

[77]  Henry Markram,et al.  Models of Neocortical Layer 5b Pyramidal Cells Capturing a Wide Range of Dendritic and Perisomatic Active Properties , 2011, PLoS Comput. Biol..

[78]  Gaute T. Einevoll,et al.  Frequency Dependence of Signal Power and Spatial Reach of the Local Field Potential , 2013, PLoS Comput. Biol..

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

[80]  U. Frey,et al.  Microelectronic system for high-resolution mapping of extracellular electric fields applied to brain slices. , 2009, Biosensors & bioelectronics.

[81]  Dieter Jaeger,et al.  The capabilities and limitations of conductance-based compartmental neuron models with reduced branched or unbranched morphologies and active dendrites , 2010, Journal of Computational Neuroscience.

[82]  M. Quirk,et al.  Interaction between spike waveform classification and temporal sequence detection , 1999, Journal of Neuroscience Methods.

[83]  Emery N. Brown,et al.  The Time-Rescaling Theorem and Its Application to Neural Spike Train Data Analysis , 2002, Neural Computation.

[84]  Christof Koch,et al.  Electrical Interactions via the Extracellular Potential Near Cell Bodies , 1999, Journal of Computational Neuroscience.

[85]  Henry Markram,et al.  Preserving axosomatic spiking features despite diverse dendritic morphology. , 2013, Journal of neurophysiology.

[86]  R. Ilmoniemi,et al.  Magnetoencephalography-theory, instrumentation, and applications to noninvasive studies of the working human brain , 1993 .

[87]  Hans Petter Langtangen,et al.  A Primer on Scientific Programming with Python , 2009 .

[88]  Gaute T. Einevoll,et al.  Power Laws from Linear Neuronal Cable Theory: Power Spectral Densities of the Soma Potential, Soma Membrane Current and Single-Neuron Contribution to the EEG , 2013, PLoS Comput. Biol..

[89]  Gaute T. Einevoll,et al.  Modelling and Analysis of Electrical Potentials Recorded in Microelectrode Arrays (MEAs) , 2015, Neuroinformatics.

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

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

[92]  J. Fohlmeister,et al.  Impulse encoding across the dendritic morphologies of retinal ganglion cells. , 1999, Journal of neurophysiology.

[93]  Anders M. Dale,et al.  The Challenge of Connecting the Dots in the B.R.A.I.N. , 2013, Neuron.

[94]  Erik De Schutter,et al.  Automated neuron model optimization techniques: a review , 2008, Biological Cybernetics.

[95]  R. Kawashima,et al.  An evaluation of the conductivity profile in the somatosensory barrel cortex of Wistar rats. , 2010, Journal of neurophysiology.

[96]  Moritz Helias,et al.  Neuroinformatics Original Research Article Pynest: a Convenient Interface to the Nest Simulator , 2022 .