A Multimodal Fitting Approach to Construct Single-Neuron Models With Patch Clamp and High-Density Microelectrode Arrays.

In computational neuroscience, multicompartment models are among the most biophysically realistic representations of single neurons. Constructing such models usually involves the use of the patch-clamp technique to record somatic voltage signals under different experimental conditions. The experimental data are then used to fit the many parameters of the model. While patching of the soma is currently the gold-standard approach to build multicompartment models, several studies have also evidenced a richness of dynamics in dendritic and axonal sections. Recording from the soma alone makes it hard to observe and correctly parameterize the activity of nonsomatic compartments. In order to provide a richer set of data as input to multicompartment models, we here investigate the combination of somatic patch-clamp recordings with recordings of high-density microelectrode arrays (HD-MEAs). HD-MEAs enable the observation of extracellular potentials and neural activity of neuronal compartments at subcellular resolution. In this work, we introduce a novel framework to combine patch-clamp and HD-MEA data to construct multicompartment models. We first validate our method on a ground-truth model with known parameters and show that the use of features extracted from extracellular signals, in addition to intracellular ones, yields models enabling better fits than using intracellular features alone. We also demonstrate our procedure using experimental data by constructing cell models from in vitro cell cultures. The proposed multimodal fitting procedure has the potential to augment the modeling efforts of the computational neuroscience community and provide the field with neuronal models that are more realistic and can be better validated.

Combining biophysical modeling and deep learning for multielectrode array neuron localization and classification.

Neural circuits typically consist of many different types of neurons, and one faces a challenge in disentangling their individual contributions in measured neural activity. Classification of cells into inhibitory and excitatory neurons and localization of neurons on the basis of extracellular recordings are frequently employed procedures. Current approaches, however, need a lot of human intervention, which makes them slow, biased, and unreliable. In light of recent advances in deep learning techniques and exploiting the availability of neuron models with quasi-realistic three-dimensional morphology and physiological properties, we present a framework for automatized and objective classification and localization of cells based on the spatiotemporal profiles of the extracellular action potentials recorded by multielectrode arrays. We train convolutional neural networks on simulated signals from a large set of cell models and show that our framework can predict the position of neurons with high accuracy, more precisely than current state-of-the-art methods. Our method is also able to classify whether a neuron is excitatory or inhibitory with very high accuracy, substantially improving on commonly used clustering techniques. Furthermore, our new method seems to have the potential to separate certain subtypes of excitatory and inhibitory neurons. The possibility of automatically localizing and classifying all neurons recorded with large high-density extracellular electrodes contributes to a more accurate and more reliable mapping of neural circuits. NEW & NOTEWORTHY We propose a novel approach to localize and classify neurons from their extracellularly recorded action potentials with a combination of biophysically detailed neuron models and deep learning techniques. Applied to simulated data, this new combination of forward modeling and machine learning yields higher performance compared with state-of-the-art localization and classification methods.

A Modular Implementation to Handle and Benchmark Drift Correction for High-Density Extracellular Recordings.

High-density neural devices are now offering the possibility to record from neuronal populations in vivo at unprecedented scale. However, the mechanical drifts often observed in these recordings are currently a major issue for "spike sorting," an essential analysis step to identify the activity of single neurons from extracellular signals. Although several strategies have been proposed to compensate for such drifts, the lack of proper benchmarks makes it hard to assess the quality and effectiveness of motion correction. In this paper, we present a benchmark study to precisely and quantitatively evaluate the performance of several state-of-the-art motion correction algorithms introduced in the literature. Using simulated recordings with induced drifts, we dissect the origins of the errors performed while applying a motion correction algorithm as a preprocessing step in the spike sorting pipeline. We show how important it is to properly estimate the positions of the neurons from extracellular traces in order to correctly estimate the probe motion, compare several interpolation procedures, and highlight what are the current limits for motion correction approaches.

Mechanical stimulation and electrophysiological monitoring at subcellular resolution reveals differential mechanosensation of neurons within networks.

A growing consensus that the brain is a mechanosensitive organ is driving the need for tools that mechanically stimulate and simultaneously record the electrophysiological response of neurons within neuronal networks. Here we introduce a synchronized combination of atomic force microscopy, high-density microelectrode array and fluorescence microscopy to monitor neuronal networks and to mechanically characterize and stimulate individual neurons at piconewton force sensitivity and nanometre precision while monitoring their electrophysiological activity at subcellular spatial and millisecond temporal resolution. No correlation is found between mechanical stiffness and electrophysiological activity of neuronal compartments. Furthermore, spontaneously active neurons show exceptional functional resilience to static mechanical compression of their soma. However, application of fast transient (∼500 ms) mechanical stimuli to the neuronal soma can evoke action potentials, which depend on the anchoring of neuronal membrane and actin cytoskeleton. Neurons show higher responsivity, including bursts of action potentials, to slower transient mechanical stimuli (∼60 s). Moreover, transient and repetitive application of the same compression modulates the neuronal firing rate. Seemingly, neuronal networks can differentiate and respond to specific characteristics of mechanical stimulation. Ultimately, the developed multiparametric tool opens the door to explore manifold nanomechanobiological responses of neuronal systems and new ways of mechanical control.

Low-Power Fully Integrated 256-Channel Nanowire Electrode-on-Chip Neural Interface for Intracellular Electrophysiology.

Intracellular electrophysiology, a vital and versatile technique in cellular neuroscience, is typically conducted using the patch-clamp method. Despite its effectiveness, this method poses challenges due to its complexity and low throughput. The pursuit of multi-channel parallel neural intracellular recording has been a long-standing goal, yet achieving reliable and consistent scaling has been elusive because of several technological barriers. In this work, we introduce a micropower integrated circuit, optimized for scalable, high-throughput in vitro intrinsically intracellular electrophysiology. This system is capable of simultaneous recording and stimulation, implementing all essential functions such as signal amplification, acquisition, and control, with a direct interface to electrodes integrated on the chip. The electrophysiology system-on-chip (eSoC), fabricated in 180nm CMOS, measures 2.236 mm × 2.236 mm. It contains four 8 × 8 arrays of nanowire electrodes, each with a 50 µm pitch, placed over the top-metal layer on the chip surface, totaling 256 channels. Each channel has a power consumption of 0.47 µW, suitable for current stimulation and voltage recording, and covers 80 dB adjustable range at a sampling rate of 25 kHz. Experimental recordings with the eSoC from cultured neurons in vitro validate its functionality in accurately resolving chemically induced multi-unit intracellular electrical activity.

Spyglass: a framework for reproducible and shareable neuroscience research.

Scientific progress depends on reliable and reproducible results. Progress can also be accelerated when data are shared and re-analyzed to address new questions. Current approaches to storing and analyzing neural data typically involve bespoke formats and software that make replication, as well as the subsequent reuse of data, difficult if not impossible. To address these challenges, we created Spyglass, an open-source software framework that enables reproducible analyses and sharing of data and both intermediate and final results within and across labs. Spyglass uses the Neurodata Without Borders (NWB) standard and includes pipelines for several core analyses in neuroscience, including spectral filtering, spike sorting, pose tracking, and neural decoding. It can be easily extended to apply both existing and newly developed pipelines to datasets from multiple sources. We demonstrate these features in the context of a cross-laboratory replication by applying advanced state space decoding algorithms to publicly available data. New users can try out Spyglass on a Jupyter Hub hosted by HHMI and 2i2c: https://spyglass.hhmi.2i2c.cloud/.

Compression strategies for large-scale electrophysiology data.

Objective.With the rapid adoption of high-density electrode arrays for recording neural activity, electrophysiology data volumes within labs and across the field are growing at unprecedented rates. For example, a one-hour recording with a 384-channel Neuropixels probe generates over 80 GB of raw data. These large data volumes carry a high cost, especially if researchers plan to store and analyze their data in the cloud. Thus, there is a pressing need for strategies that can reduce the data footprint of each experiment.Approach.Here, we establish a set of benchmarks for comparing the performance of various compression algorithms on experimental and simulated recordings from Neuropixels 1.0 (NP1) and 2.0 (NP2) probes.Main results.For lossless compression, audio codecs (FLACandWavPack) achieve compression ratios (CRs) 6% higher for NP1 and 10% higher for NP2 than the best general-purpose codecs, at the expense of decompression speed. For lossy compression, theWavPackalgorithm in 'hybrid mode' increases the CR from 3.59 to 7.08 for NP1 and from 2.27 to 7.04 for NP2 (compressed file size of ∼14% for both types of probes), without adverse effects on spike sorting accuracy or spike waveforms.Significance.Along with the tools we have developed to make compression easier to deploy, these results should encourage all electrophysiologists to apply compression as part of their standard analysis workflows.

ProbeInterface: A Unified Framework for Probe Handling in Extracellular Electrophysiology.

Recording neuronal activity with penetrating extracellular multi-channel electrode arrays, more commonly known as neural probes, is one of the most widespread approaches to probe neuronal activity. Despite a plethora of available extracellular probe designs, the time-consuming process of mapping of electrode channel order and relative geometries, as required by spike-sorting software is invariably left to the end-user. Consequently, this manual process is prone to mis-mapping mistakes, which in turn lead to undesirable spike-sorting errors and inefficiencies. Here, we introduce ProbeInterface, an open-source project that aims to unify neural probe metadata descriptions by removing the manual step of probe mapping prior to spike-sorting for the analysis of extracellular neural recordings. ProbeInterface is first of all a Python API, which enables users to create and visualize probes and probe groups at any required complexity level. Second, ProbeInterface facilitates the generation of comprehensive wiring description in a reproducible fashion for any specific data-acquisition setup, which usually involves the use of a recording probe, a headstage, adapters, and an acquisition system. Third, we collaborate with probe manufacturers to compile an open library of available probes, which can be downloaded at run time using our Python API. Finally, with ProbeInterface we define a file format for probe handling which includes all necessary information for a FAIR probe description and is compatible with and complementary to other open standards in neuroscience.

How Do Spike Collisions Affect Spike Sorting Performance?

Recently, a new generation of devices have been developed to record neural activity simultaneously from hundreds of electrodes with a very high spatial density, both for in vitro and in vivo applications. While these advances enable to record from many more cells, they also challenge the already complicated process of spike sorting (i.e., extracting isolated single-neuron activity from extracellular signals). In this work, we used synthetic ground-truth recordings with controlled levels of correlations among neurons to quantitatively benchmark the performance of state-of-the-art spike sorters focusing specifically on spike collisions. Our results show that while modern template-matching-based algorithms are more accurate than density-based approaches, all methods, to some extent, failed to detect synchronous spike events of neurons with similar extracellular signals. Interestingly, the performance of the sorters is not largely affected by the spiking activity in the recordings, with respect to average firing rates and spike-train correlation levels. Since the performances of all modern spike sorting algorithms can be affected as function of the activity of the recorded neurons, scientific claims on correlations and synchrony should be carefully assessed based on the analysis provided in this paper.

An automated method for precise axon reconstruction from recordings of high-density micro-electrode arrays.

Objective:Neurons communicate with each other by sending action potentials (APs) through their axons. The velocity of axonal signal propagation describes how fast electrical APs can travel. This velocity can be affected in a human brain by several pathologies, including multiple sclerosis, traumatic brain injury and channelopathies. High-density microelectrode arrays (HD-MEAs) provide unprecedented spatio-temporal resolution to extracellularly record neural electrical activity. The high density of the recording electrodes enables to image the activity of individual neurons down to subcellular resolution, which includes the propagation of axonal signals. However, axon reconstruction, to date, mainly relies on manual approaches to select the electrodes and channels that seemingly record the signals along a specific axon, while an automated approach to track multiple axonal branches in extracellular action-potential recordings is still missing.Approach:In this article, we propose a fully automated approach to reconstruct axons from extracellular electrical-potential landscapes, so-called 'electrical footprints' of neurons. After an initial electrode and channel selection, the proposed method first constructs a graph based on the voltage signal amplitudes and latencies. Then, the graph is interrogated to extract possible axonal branches. Finally, the axonal branches are pruned, and axonal action-potential propagation velocities are computed.Main results:We first validate our method using simulated data from detailed reconstructions of neurons, showing that our approach is capable of accurately reconstructing axonal branches. We then apply the reconstruction algorithm to experimental recordings of HD-MEAs and show that it can be used to determine axonal morphologies and signal-propagation velocities at high throughput.Significance:We introduce a fully automated method to reconstruct axonal branches and estimate axonal action-potential propagation velocities using HD-MEA recordings. Our method yields highly reliable and reproducible velocity estimations, which constitute an important electrophysiological feature of neuronal preparations.

Tracking axon initial segment plasticity using high-density microelectrode arrays: A computational study.

Despite being composed of highly plastic neurons with extensive positive feedback, the nervous system maintains stable overall function. To keep activity within bounds, it relies on a set of negative feedback mechanisms that can induce stabilizing adjustments and that are collectively termed "homeostatic plasticity." Recently, a highly excitable microdomain, located at the proximal end of the axon-the axon initial segment (AIS)-was found to exhibit structural modifications in response to activity perturbations. Though AIS plasticity appears to serve a homeostatic purpose, many aspects governing its expression and its functional role in regulating neuronal excitability remain elusive. A central challenge in studying the phenomenon is the rich heterogeneity of its expression (distal/proximal relocation, shortening, lengthening) and the variability of its functional role. A potential solution is to track AISs of a large number of neurons over time and attempt to induce structural plasticity in them. To this end, a promising approach is to use extracellular electrophysiological readouts to track a large number of neurons at high spatiotemporal resolution by means of high-density microelectrode arrays (HD-MEAs). However, an analysis framework that reliably identifies specific activity signatures that uniquely map on to underlying microstructural changes is missing. In this study, we assessed the feasibility of such a task and used the distal relocation of the AIS as an exemplary problem. We used sophisticated computational models to systematically explore the relationship between incremental changes in AIS positions and the specific consequences observed in simulated extracellular field potentials. An ensemble of feature changes in the extracellular fields that reliably characterize AIS plasticity was identified. We trained models that could detect these signatures with remarkable accuracy. Based on these findings, we propose a hybrid analysis framework that could potentially enable high-throughput experimental studies of activity-dependent AIS plasticity using HD-MEAs.

Inferring monosynaptic connections from paired dendritic spine Ca2+imaging and large-scale recording of extracellular spiking.

Objective: Techniques to identify monosynaptic connections between neurons have been vital for neuroscience research, facilitating important advancements concerning network topology, synaptic plasticity, and synaptic integration, among others.Approach: Here, we introduce a novel approach to identify and monitor monosynaptic connections using high-resolution dendritic spine Ca2+imaging combined with simultaneous large-scale recording of extracellular electrical activity by means of high-density microelectrode arrays.Main results: We introduce an easily adoptable analysis pipeline that associates the imaged spine with its presynaptic unit and test it onin vitrorecordings. The method is further validated and optimized by simulating synaptically-evoked spine Ca2+transients based on measured spike trains in order to obtain simulated ground-truth connections.Significance: The proposed approach offers unique advantages as (a) it can be used to identify monosynaptic connections with an accurate localization of the synapse within the dendritic tree, (b) it provides precise information of presynaptic spiking, and (c) postsynaptic spine Ca2+signals and, finally, (d) the non-invasive nature of the proposed method allows for long-term measurements. The analysis toolkit together with the rich data sets that were acquired are made publicly available for further exploration by the research community.

Functional neuronal circuitry and oscillatory dynamics in human brain organoids.

Human brain organoids replicate much of the cellular diversity and developmental anatomy of the human brain. However, the physiology of neuronal circuits within organoids remains under-explored. With high-density CMOS microelectrode arrays and shank electrodes, we captured spontaneous extracellular activity from brain organoids derived from human induced pluripotent stem cells. We inferred functional connectivity from spike timing, revealing a large number of weak connections within a skeleton of significantly fewer strong connections. A benzodiazepine increased the uniformity of firing patterns and decreased the relative fraction of weakly connected edges. Our analysis of the local field potential demonstrate that brain organoids contain neuronal assemblies of sufficient size and functional connectivity to co-activate and generate field potentials from their collective transmembrane currents that phase-lock to spiking activity. These results point to the potential of brain organoids for the study of neuropsychiatric diseases, drug action, and the effects of external stimuli upon neuronal networks.

MEArec: A Fast and Customizable Testbench Simulator for Ground-truth Extracellular Spiking Activity.

When recording neural activity from extracellular electrodes, both in vivo and in vitro, spike sorting is a required and very important processing step that allows for identification of single neurons' activity. Spike sorting is a complex algorithmic procedure, and in recent years many groups have attempted to tackle this problem, resulting in numerous methods and software packages. However, validation of spike sorting techniques is complicated. It is an inherently unsupervised problem and it is hard to find universal metrics to evaluate performance. Simultaneous recordings that combine extracellular and patch-clamp or juxtacellular techniques can provide ground-truth data to evaluate spike sorting methods. However, their utility is limited by the fact that only a few cells can be measured at the same time. Simulated ground-truth recordings can provide a powerful alternative mean to rank the performance of spike sorters. We present here MEArec, a Python-based software which permits flexible and fast simulation of extracellular recordings. MEArec allows users to generate extracellular signals on various customizable electrode designs and can replicate various problematic aspects for spike sorting, such as bursting, spatio-temporal overlapping events, and drifts. We expect MEArec will provide a common testbench for spike sorting development and evaluation, in which spike sorting developers can rapidly generate and evaluate the performance of their algorithms.

A modulated template-matching approach to improve spike sorting of bursting neurons.

In extracellular neural electrophysiology, individual spikes have to be assigned to their cell of origin in a procedure called "spike sorting". Spike sorting is an unsupervised problem, since no ground-truth information is generally available. Here, we focus on improving spike sorting performance, particularly during periods of high synchronous activity or so-called "bursting". Bursting entails systematic changes in spike shapes and amplitudes and remains a challenge for current spike sorting schemes. We use realistic simulated bursting recordings of high-density micro-electrode arrays (HD-MEAs) and we present a fully automated algorithm based on template matching with a focus on recovering missed spikes during bursts. To compare and benchmark spike-sorting performance after applying our method, we used ground-truth information of simulated recordings. We show that our approach can be effective in improving spike sorting performance during bursting. Further validation with experimental recordings is necessary.

Optogenetic pacing of medial septum parvalbumin-positive cells disrupts temporal but not spatial firing in grid cells.

Grid cells in the medial entorhinal cortex (MEC) exhibit remarkable spatial activity patterns with spikes coordinated by theta oscillations driven by the medial septal area (MSA). Spikes from grid cells progress relative to the theta phase in a phenomenon called phase precession, which is suggested as essential to create the spatial periodicity of grid cells. Here, we show that optogenetic activation of parvalbumin-positive (PV+) cells in the MSA enabled selective pacing of local field potential (LFP) oscillations in MEC. During optogenetic stimulation, the grid cells were locked to the imposed pacing frequency but kept their spatial patterns. Phase precession was abolished, and speed information was no longer reflected in the LFP oscillations but was still carried by rate coding of individual MEC neurons. Together, these results support that theta oscillations are not critical to the spatial pattern of grid cells and do not carry a crucial velocity signal.

Electrophysiological Phenotype Characterization of Human iPSC-Derived Neuronal Cell Lines by Means of High-Density Microelectrode Arrays.

Recent advances in the field of cellular reprogramming have opened a route to studying the fundamental mechanisms underlying common neurological disorders. High-density microelectrode-arrays (HD-MEAs) provide unprecedented means to study neuronal physiology at different scales, ranging from network through single-neuron to subcellular features. In this work, HD-MEAs are used in vitro to characterize and compare human induced-pluripotent-stem-cell-derived dopaminergic and motor neurons, including isogenic neuronal lines modeling Parkinson's disease and amyotrophic lateral sclerosis. Reproducible electrophysiological network, single-cell and subcellular metrics are used for phenotype characterization and drug testing. Metrics, such as burst shape and axonal velocity, enable the distinction of healthy and diseased neurons. The HD-MEA metrics can also be used to detect the effects of dosing the drug retigabine to human motor neurons. Finally, it is shown that the ability to detect drug effects and the observed culture-to-culture variability critically depend on the number of available recording electrodes.

Experimental Pipeline (Expipe): A Lightweight Data Management Platform to Simplify the Steps From Experiment to Data Analysis.

As experimental neuroscience is moving toward more integrative approaches, with a variety of acquisition techniques covering multiple spatiotemporal scales, data management is becoming increasingly challenging for neuroscience laboratories. Often, datasets are too large to practically be stored on a laptop or a workstation. The ability to query metadata collections without retrieving complete datasets is therefore critical to efficiently perform new analyses and explore the data. At the same time, new experimental paradigms lead to constantly changing specifications for the metadata to be stored. Despite this, there is currently a serious lack of agile software tools for data management in neuroscience laboratories. To meet this need, we have developed Expipe, a lightweight data management framework that simplifies the steps from experiment to data analysis. Expipe provides the functionality to store and organize experimental data and metadata for easy retrieval in exploration and analysis throughout the experimental pipeline. It is flexible in terms of defining the metadata to store and aims to solve the storage and retrieval challenges of data/metadata due to ever changing experimental pipelines. Due to its simplicity and lightweight design, we envision Expipe as an easy-to-use data management solution for experimental laboratories, that can improve provenance, reproducibility, and sharing of scientific projects.

SpikeInterface, a unified framework for spike sorting.

Much development has been directed toward 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.

SpikeForest, reproducible web-facing ground-truth validation of automated neural spike sorters.

Spike sorting is a crucial step in electrophysiological studies of neuronal activity. While many spike sorting packages are available, there is little consensus about which are most accurate under different experimental conditions. SpikeForest is an open-source and reproducible software suite that benchmarks the performance of automated spike sorting algorithms across an extensive, curated database of ground-truth electrophysiological recordings, displaying results interactively on a continuously-updating website. With contributions from eleven laboratories, our database currently comprises 650 recordings (1.3 TB total size) with around 35,000 ground-truth units. These data include paired intracellular/extracellular recordings and state-of-the-art simulated recordings. Ten of the most popular spike sorting codes are wrapped in a Python package and evaluated on a compute cluster using an automated pipeline. SpikeForest documents community progress in automated spike sorting, and guides neuroscientists to an optimal choice of sorter and parameters for a wide range of probes and brain regions.