SpikeInterface, a unified framework for spike sorting
Abstract
Much development has been directed towards improving the performance and automation of spike sorting. This continuous development, while essential, has contributed to an over-saturation of new, incompatible tools that hinders rigorous benchmarking and complicates reproducible analysis. To address these limitations, we developed SpikeInterface, a Python framework designed to unify preexisting spike sorting technologies into a single codebase and to facilitate straightforward comparison and adoption of different approaches. With a few lines of code, researchers can reproducibly run, compare, and benchmark most modern spike sorting algorithms; pre-process, post-process, and visualize extracellular datasets; validate, curate, and export sorting outputs; and more. In this paper, we provide an overview of SpikeInterface and, with applications to real and simulated datasets, demonstrate how it can be utilized to reduce the burden of manual curation and to more comprehensively benchmark automated spike sorters.
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files. The datasets are uploaded to the DANDI archive, dataset 000034 (https://gui.dandiarchive.org/#/dandiset/000034). The source code for generating all figures is also publicly available at: https://spikeinterface.github.io/
Article and author information
Author details
Funding
Wellcome Trust (214431/Z/18/Z)
- Matthias H Hennig
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Copyright
© 2020, Buccino et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 16,046
- views
-
- 1,560
- downloads
-
- 191
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Neuroscience
At many vertebrate synapses, presynaptic functions are tuned by expression of different Cav2 channels. Most invertebrate genomes contain only one Cav2 gene. The Drosophila Cav2 homolog, cacophony (cac), induces synaptic vesicle release at presynaptic active zones (AZs). We hypothesize that Drosophila cac functional diversity is enhanced by two mutually exclusive exon pairs that are not conserved in vertebrates, one in the voltage sensor and one in the loop binding Caβ and Gβγ subunits. We find that alternative splicing in the voltage sensor affects channel activation voltage. Only the isoform with the higher activation voltage localizes to AZs at the glutamatergic Drosophila larval neuromuscular junction and is imperative for normal synapse function. By contrast, alternative splicing at the other alternative exon pair tunes multiple aspects of presynaptic function. While expression of one exon yields normal transmission, expression of the other reduces channel number in the AZ and thus release probability. This also abolishes presynaptic homeostatic plasticity. Moreover, reduced channel number affects short-term plasticity, which is rescued by increasing the external calcium concentration to match release probability to control. In sum, in Drosophila alternative splicing provides a mechanism to regulate different aspects of presynaptic functions with only one Cav2 gene.
-
- Neuroscience
The amygdala is a subcortical region in the mesiotemporal lobe that plays a key role in emotional and sensory functions. Conventional neuroimaging experiments treat this structure as a single, uniform entity, but there is ample histological evidence for subregional heterogeneity in microstructure and function. The current study characterized subregional structure-function coupling in the human amygdala, integrating post-mortem histology and in vivo MRI at ultra-high fields. Core to our work was a novel neuroinformatics approach that leveraged multiscale texture analysis as well as non-linear dimensionality reduction techniques to identify salient dimensions of microstructural variation in a 3D post-mortem histological reconstruction of the human amygdala. We observed two axes of subregional variation in this region, describing inferior-superior as well as mediolateral trends in microstructural differentiation that in part recapitulated established atlases of amygdala subnuclei. Translating our approach to in vivo MRI data acquired at 7 Tesla, we could demonstrate the generalizability of these spatial trends across 10 healthy adults. We then cross-referenced microstructural axes with functional blood-oxygen-level dependent (BOLD) signal analysis obtained during task-free conditions, and revealed a close association of structural axes with macroscale functional network embedding, notably the temporo-limbic, default mode, and sensory-motor networks. Our novel multiscale approach consolidates descriptions of amygdala anatomy and function obtained from histological and in vivo imaging techniques.