SpikeInterface, a unified framework for spike sorting

Localization of neurons from extracellular footprints

BACKGROUND

High density microelectrode arrays (HD-MEAs) are now widely used for both in-vitro and in-vivo recordings, as they allow spikes from hundreds of neurons to be recorded simultaneously. Since extracellular recordings do not allow visualization of the recorded neurons, algorithms are needed to estimate their physical positions, especially to track their movements when the are drifting away from recording devices.

NEW METHOD

The objective of this study was to evaluate the performance of multiple algorithms for neuron localization solely from extracellular traces (MEA recordings), either artificial or obtained from mouse retina. The algorithms compared included center-of-mass, monopolar, and grid-based algorithms. The first method is a barycenter calculation. The second algorithm infers the position of the cell using triangulation with the assumption that the neuron behaves as a monopole. Finally, grid-based methods rely on comparing the recorded spike with a projection of spikes of hypothetical neurons with different positions.

RESULTS

The Grid-Based algorithm yielded the most satisfactory outcomes. The center-of-mass exhibited a minimal computational cost, yet its average localization was suboptimal. Monopolar algorithms gave cell localizations with an average error of less than 10μm, but they had considerable variability and a high computational cost. For the grid-based method, the variability was smaller, with satisfactory performance and low computational cost.

COMPARISON WITH EXISTING METHOD(S)

The accuracy of the different localization methods benchmarked in this article had not been properly tested with ground-truth recordings before.

CONCLUSION

The objective of this article is to provide guidance to researchers on the selection of optimal methods for localizing neurons based on MEA recordings.

Pallidal neuronal activity in Gilles de la Tourette syndrome and dystonic patients: A comparative study

Gilles de la Tourette syndrome (GTS) and dystonia (DYS) are both hyperkinetic movement disorders effectively treated by deep brain stimulation (DBS) of the internal part of the globus pallidus (GPi). In this study, we compared single-neuron activity in the GPi between 18 GTS patients (with an average of 41 cells per patient) and 17 DYS patients (with an average of 54 cells per patient), all of whom underwent bilateral pallidal stimulation surgery, under general anesthesia or while awake at rest. We found no significant differences in GPi neuronal activity characteristics between patients operated on under general anesthesia versus those who were awake, irrespective of their diagnosis (GTS or DYS). We found higher firing rates, firing rate in bursts, pause duration and interspike interval coefficient of variation in GTS patients compared to DYS patients. On the opposite, we found higher number of pauses and bursts frequency in DYS patients. Lastly, we found a higher proportion of GPi oscillatory activities in DYS compared to GTS patients, with predominant activity within the low-frequency band (theta/alpha) in both patient groups. These findings underscore the complex relationship between the different neuronal discharge characteristic such as oscillatory or bursting activity within the GPi in shaping the clinical phenotypes of hyperkinetic disorders. Further research is warranted to deepen our understanding of how neuronal patterns are transmitted within deep brain structures and to develop strategies aimed at normalizing these pathological activities, by refining DBS techniques to enhance treatment efficacy and individual outcomes.

Graphene-based microelectrodes with bidirectional functionality for next-generation retinal electronic interfaces

Neuroelectronic prostheses are being developed for restoring vision at the retinal level in patients who have lost their sight due to photoreceptor loss. The core component of these devices is the electrode array, which enables interfacing with retinal neurons. Generating the perception of meaningful images requires high-density microelectrode arrays (MEAs) capable of precisely activating targeted retinal neurons. Achieving this precision necessitates the downscaling of electrodes to micrometer dimensions. However, miniaturization increases electrode impedance, which poses challenges by limiting the amount of current that can be delivered, thereby impairing the electrode's capability for effective neural modulation. Additionally, it elevates noise levels, reducing the signal quality of the recorded neural activity. This report focuses on evaluating reduced graphene oxide (rGO) based devices for interfacing with the retina, showcasing their potential in vision restoration. Our findings reveal low impedance and high charge injection limit for microscale rGO electrodes, confirming their suitability for developing next-generation high-density retinal devices. We successfully demonstrated bidirectional interfacing with cell cultures and explanted retinal tissue, enabling the identification and modulation of multiple cells' activity. Additionally, calcium imaging allowed real-time monitoring of retinal cell dynamics, demonstrating a significant reduction in activated areas with small-sized electrodes. Overall, this study lays the groundwork for developing advanced rGO-based MEAs for high-acuity visual prostheses.

Massive perturbation of sound representations by anesthesia in the auditory brainstem

Anesthesia modifies sensory representations in the thalamo-cortical circuit but is considered to have a milder impact on peripheral sensory processing. Here, tracking the same neurons across wakefulness and isoflurane or ketamine medetomidine anesthesia, we show that the amplitude and sign of single neuron responses to sounds are massively modified by anesthesia in the cochlear nucleus of the brainstem, the first relay of the auditory system. The reorganization of activity is so profound that decoding of sound representation under anesthesia is not possible based on awake activity. However, population-level parameters, such as average tuning strength and population decoding accuracy, are weakly affected by anesthesia, explaining why its effect has previously gone unnoticed when comparing independently sampled neurons. Together, our results indicate that the functional organization of the auditory brainstem largely depends on the network state and is ill-defined under anesthesia. This demonstrates a remarkable sensitivity of an early sensory stage to anesthesia, which is bound to disrupt downstream processing.

Human assembloids reveal the consequences of CACNA1G gene variants in the thalamocortical pathway

Abnormalities in thalamocortical crosstalk can lead to neuropsychiatric disorders. Variants in CACNA1G, which encodes the α1G subunit of the thalamus-enriched T-type calcium channel, are associated with absence seizures, intellectual disability, and schizophrenia, but the cellular and circuit consequences of these genetic variants in humans remain unknown. Here, we developed a human assembloid model of the thalamocortical pathway to dissect the contribution of genetic variants in T-type calcium channels. We discovered that the M1531V CACNA1G variant associated with seizures led to changes in T-type currents in thalamic neurons, as well as correlated hyperactivity of thalamic and cortical neurons in assembloids. By contrast, CACNA1G loss, which has been associated with risk of schizophrenia, resulted in abnormal thalamocortical connectivity that was related to both increased spontaneous thalamic activity and aberrant axonal projections. These results illustrate the utility of multi-cellular systems for interrogating human genetic disease risk variants at both cellular and circuit level.

Assembloid model to study loop circuits of the human nervous system

Neural circuits connecting the cerebral cortex, the basal ganglia and the thalamus are fundamental networks for sensorimotor processing and their dysfunction has been consistently implicated in neuropsychiatric disorders 1-9 . These recursive, loop circuits have been investigated in animal models and by clinical neuroimaging, however, direct functional access to developing human neurons forming these networks has been limited. Here, we use human pluripotent stem cells to reconstruct an in vitro cortico-striatal-thalamic-cortical circuit by creating a four-part loop assembloid. More specifically, we generate regionalized neural organoids that resemble the key elements of the cortico-striatal-thalamic-cortical circuit, and functionally integrate them into loop assembloids using custom 3D-printed biocompatible wells. Volumetric and mesoscale calcium imaging, as well as extracellular recordings from individual parts of these assembloids reveal the emergence of synchronized patterns of neuronal activity. In addition, a multi-step rabies retrograde tracing approach demonstrate the formation of neuronal connectivity across the network in loop assembloids. Lastly, we apply this system to study heterozygous loss of ASH1L gene associated with autism spectrum disorder and Tourette syndrome and discover aberrant synchronized activity in disease model assembloids. Taken together, this human multi-cellular platform will facilitate functional investigations of the cortico-striatal-thalamic-cortical circuit in the context of early human development and in disease conditions.

I-Spin live, an open-source software based on blind-source separation for real-time decoding of motor unit activity in humans

Decoding the activity of individual neural cells during natural behaviours allows neuroscientists to study how the nervous system generates and controls movements. Contrary to other neural cells, the activity of spinal motor neurons can be determined non-invasively (or minimally invasively) from the decomposition of electromyographic (EMG) signals into motor unit firing activities. For some interfacing and neuro-feedback investigations, EMG decomposition needs to be performed in real time. Here, we introduce an open-source software that performs real-time decoding of motor neurons using a blind-source separation approach for multichannel EMG signal processing. Separation vectors (motor unit filters) are optimised for each motor unit from baseline contractions and then re-applied in real time during test contractions. In this way, the firing activity of multiple motor neurons can be provided through different forms of visual feedback. We provide a complete framework with guidelines and examples of recordings to guide researchers who aim to study movement control at the motor neuron level. We first validated the software with synthetic EMG signals generated during a range of isometric contraction patterns. We then tested the software on data collected using either surface or intramuscular electrode arrays from five lower limb muscles (gastrocnemius lateralis and medialis, vastus lateralis and medialis, and tibialis anterior). We assessed how the muscle or variation of contraction intensity between the baseline contraction and the test contraction impacted the accuracy of the real-time decomposition. This open-source software provides a set of tools for neuroscientists to design experimental paradigms where participants can receive real-time feedback on the output of the spinal cord circuits.

Engineered modular neuronal networks-on-chip represent structure-function relationship

Brain function is substantially linked to the highly organized modular structure of neuronal networks. However, the structure of in vitro assembled neuronal circuits often exhibits variability, complicating the consistent recording of network functional output and its correlation to network structure. Therefore, engineering neuronal structures with predefined geometry and reproducible functional features is essential to precisely model in vivo neuronal circuits. Here, we engineered microchannel devices to assemble 2D and 3D modular networks. The microchannel devices were coupled with a multi-electrode array (MEA) electrophysiology system to enable recordings from circuits. Each network consisted of 64 modules connected to their adjacent modules by micron-sized channels. Modular circuits within microchannel devices showed enhanced activity and functional connectivity traits. This includes metrics such as connection weights, clustering coefficient, global efficiency, and the number of hub neurons with higher betweenness centrality. In addition, modular networks demonstrated an increased functional modularity score compared to the randomly formed circuits. Neurons within individual modules displayed uniform network characteristics and predominantly participated in their respective functional communities within the same or neighboring physical modules. These observations highlight that the modular network structure promotes the development of segregated functional connectivity traits while simultaneously enhancing the efficiency of overall network connectivity. Our findings emphasize the significant impact of physical constraints on the activity patterns and functional organization within engineered modular networks. These circuits, characterized by stable modular architecture and intricate functional dynamics-key features of the brain networks-offer a robust in vitro model for advancing neuroscience research.

Prefrontal synaptic regulation of homeostatic sleep pressure revealed through synaptic chemogenetics

Sleep is regulated by homeostatic processes, yet the biological basis of sleep pressure that accumulates during wakefulness, triggers sleep, and dissipates during sleep remains elusive. We explored a causal relationship between cellular synaptic strength and electroencephalography delta power indicating macro-level sleep pressure by developing a theoretical framework and a molecular tool to manipulate synaptic strength. The mathematical model predicted that increased synaptic strength promotes the neuronal "down state" and raises the delta power. Our molecular tool (synapse-targeted chemically induced translocation of Kalirin-7, SYNCit-K), which induces dendritic spine enlargement and synaptic potentiation through chemically induced translocation of protein Kalirin-7, demonstrated that synaptic potentiation of excitatory neurons in the prefrontal cortex (PFC) increases nonrapid eye movement sleep amounts and delta power. Thus, synaptic strength of PFC excitatory neurons dictates sleep pressure in mammals.

Tracking neurons across days with high-density probes

Neural activity spans multiple time scales, from milliseconds to months. Its evolution can be recorded with chronic high-density arrays such as Neuropixels probes, which can measure each spike at tens of sites and record hundreds of neurons. These probes produce vast amounts of data that require different approaches for tracking neurons across recordings. Here, to meet this need, we developed UnitMatch, a pipeline that operates after spike sorting, based only on each unit's average spike waveform. We tested UnitMatch in Neuropixels recordings from the mouse brain, where it tracked neurons across weeks. Across the brain, neurons had distinctive inter-spike interval distributions. Their correlations with other neurons remained stable over weeks. In the visual cortex, the neurons' selectivity for visual stimuli remained similarly stable. In the striatum, however, neuronal responses changed across days during learning of a task. UnitMatch is thus a promising tool to reveal both invariance and plasticity in neural activity across days.

Human Neural Organoid Microphysiological Systems Show the Building Blocks Necessary for Basic Learning and Memory

Brain Microphysiological Systems including neural organoids derived from human induced pluripotent stem cells offer a unique lens to study the intricate workings of the human brain. This paper investigates the foundational elements of learning and memory in neural organoids, also known as Organoid Intelligence by quantifying immediate early gene expression, synaptic plasticity, neuronal network dynamics, and criticality to demonstrate the utility of these organoids in basic science research. Neural organoids showed synapse formation, glutamatergic and GABAergic receptor expression, immediate early gene expression basally and evoked, functional connectivity, criticality, and synaptic plasticity in response to theta-burst stimulation. In addition, pharmacological interventions on GABAergic and glutamatergic receptors, and input specific theta-burst stimulation further shed light on the capacity of neural organoids to mirror synaptic modulation and short-term potentiation, demonstrating their potential as tools for studying neurophysiological and neurological processes and informing therapeutic strategies for diseases.

Unsupervised spike sorting for multielectrode arrays based on spike shape features and location methods

Microelectrode arrays (MEAs) enable simultaneous measurement of spike trains from numerous neurons, owing to advancements in microfabrication technology. These probes are highly valuable for comprehending the intricate dynamics of neuronal networks. Spike sorting is a pivotal step in comprehensively analyzing the activity of neuronal networks from extracellular recordings. However, the accuracy of spike sorting is relatively low due to the dense sampling of spikes in MEAs. Here, we propose an unsupervised pipeline named UMAP-COM method, which utilizes combined features to address this problem. These combined features comprise dominant spike shape features extracted by the uniform manifold approximation and projection (UMAP), as well as spike locations estimated by the center of mass (COM). We validate the UMAP-COM method on publicly available datasets from different kinds of probes, demonstrating that it is more accurate than other spike sorting methods. Furthermore, we conduct separate evaluations of spike shape feature extraction methods and spike localization methods. In this comparison, UMAP emerges as the superior feature extraction method, demonstrating its effectiveness in accurately representing spike shapes. Additionally, we find that the COM method outperforms other spike localization methods, highlighting its ability to enhance the accuracy of spike sorting.

A nonoscillatory, millisecond-scale embedding of brain state provides insight into behavior

The most robust and reliable signatures of brain states are enriched in rhythms between 0.1 and 20 Hz. Here we address the possibility that the fundamental unit of brain state could be at the scale of milliseconds and micrometers. By analyzing high-resolution neural activity recorded in ten mouse brain regions over 24 h, we reveal that brain states are reliably identifiable (embedded) in fast, nonoscillatory activity. Sleep and wake states could be classified from 100 to 101 ms of neuronal activity sampled from 100 µm of brain tissue. In contrast to canonical rhythms, this embedding persists above 1,000 Hz. This high-frequency embedding is robust to substates, sharp-wave ripples and cortical on/off states. Individual regions intermittently switched states independently of the rest of the brain, and such brief state discontinuities coincided with brief behavioral discontinuities. Our results suggest that the fundamental unit of state in the brain is consistent with the spatial and temporal scale of neuronal computation.

Failure in a population: Tauopathy disrupts homeostatic set-points in emergent dynamics despite stability in the constituent neurons

Homeostatic regulation of neuronal activity is essential for robust computation; set-points, such as firing rate, are actively stabilized to compensate for perturbations. The disruption of brain function central to neurodegenerative disease likely arises from impairments of computationally essential set-points. Here, we systematically investigated the effects of tau-mediated neurodegeneration on all known set-points in neuronal activity. We continuously tracked hippocampal neuronal activity across the lifetime of a mouse model of tauopathy. We were unable to detect effects of disease in measures of single-neuron firing activity. By contrast, as tauopathy progressed, there was disruption of network-level neuronal activity, quantified by measuring neuronal pairwise interactions and criticality, a homeostatically controlled, ideal computational regime. Deviations in criticality correlated with symptoms, predicted underlying anatomical pathology, occurred in a sleep-wake-dependent manner, and could be used to reliably classify an animal's genotype. This work illustrates how neurodegeneration may disrupt the computational capacity of neurobiological systems.

Cortical Coding of Gustatory and Thermal Signals in Active Licking Mice

Eating behaviors are influenced by the integration of gustatory, olfactory, and somatosensory signals, which all contribute to the perception of flavor. Although extensive research has explored the neural correlates of taste in the gustatory cortex (GC), less is known about its role in encoding thermal information. This study investigates the encoding of oral thermal and chemosensory signals by GC neurons compared to the oral somatosensory cortex. In this study, we recorded the spiking activity of more than 900 GC neurons and 500 neurons from the oral somatosensory cortex in mice allowed to freely lick small drops of gustatory stimuli or deionized water at varying non-nociceptive temperatures. We then developed and used a Bayesian-based analysis technique to assess neural classification scores based on spike rate and phase timing within the lick cycle. Our results indicate that GC neurons rely predominantly on rate information, although phase information is needed to achieve maximum accuracy, to effectively encode both chemosensory and thermosensory signals. GC neurons can effectively differentiate between thermal stimuli, excelling in distinguishing both large contrasts (14°C vs. 36°C) and, although less effectively, more subtle temperature differences. Finally, a direct comparison of the decoding accuracy of thermosensory signals between the two cortices reveals that while the somatosensory cortex showed higher overall accuracy, the GC still encodes significant thermosensory information. These findings highlight the GC's dual role in processing taste and temperature, emphasizing the importance of considering temperature in future studies of taste processing.

Neural Processing of Taste-Related Signals in the Mediodorsal Thalamus of Mice

Our consummatory decisions depend on the taste of food and the reward experienced while eating, which are processed through neural computations in interconnected brain areas. Although many gustatory regions of rodents have been explored, the mediodorsal nucleus of the thalamus (MD) remains understudied. The MD, a multimodal brain area connected with gustatory centers, is often studied for its role in processing associative and cognitive information and has been shown to represent intraorally-delivered chemosensory stimuli after strong retronasal odor-taste associations. Key questions remain about whether MD neurons can process taste quality independently of odor-taste associations and how they represent extraoral signals predicting rewarding and aversive gustatory outcomes. Here, we present electrophysiological evidence demonstrating how mouse MD neurons represent and encode 1) the identity and concentrations of basic taste qualities during active licking, and 2) auditory signals anticipating rewarding and aversive taste outcomes. Our data reveal that MD neurons can reliably and dynamically encode taste identity in a broadly tuned manner and taste concentrations with spiking activity positively and negatively correlated with stimulus intensity. Our data also show that MD can represent information related to predictive cues and their associated outcomes, regardless of whether the cue predicts a rewarding or aversive outcome. In summary, our findings suggest that the mediodorsal thalamus is integral to the taste pathway, as it can encode sensory-discriminative dimensions of tastants and participate in processing associative information essential for ingestive behaviors.

A fully automatic multichannel neural spike sorting algorithm with spike reduction and positional feature

Objective: The sorting of neural spike data recorded by multichannel and high channel neural probes such as Neuropixels, especially in real-time, remains a significant technical challenge. Most neural spike sorting algorithms focus on sorting neural spikes post-hoc for high sorting accuracy-but reducing the processing delay for fast sorting, potentially even live sorting, is generally not possible with these algorithms.Approach: Here we report our Graph nEtwork Multichannel sorting (GEMsort) algorithm, which is largely based on graph network, to allow rapid neural spike sorting for multiple neural recording channels. This was accomplished by two innovations: In GEMsort, duplicated neural spikes recorded from multiple channels were eliminated from duplicate channels by only selecting the highest amplitude neural spike in any channel for subsequent processing. In addition, the channel from which the representative neural spike was recorded was used as an additional feature to differentiate between neural spikes recorded from different neurons having similar temporal features.Main results: Synthetic and experimentally recorded multichannel neural recordings were used to evaluate the sorting performance of GEMsort. The sorting results of GEMsort were also compared with two other state-of-the-art sorting algorithms (Kilosort and Mountainsort) in sorting time and sorting agreements.Significance: GEMsort allows rapidly sort neural spikes and is highly suitable to be implemented with digital circuitry for high processing speed and channel scalability.

Assessing Cross-Contamination in Spike-Sorted Electrophysiology Data

Recent advances in extracellular electrophysiology now facilitate the recording of spikes from hundreds or thousands of neurons simultaneously. This has necessitated both the development of new computational methods for spike sorting and better methods to determine spike-sorting accuracy. One long-standing method of assessing the false discovery rate (FDR) of spike sorting-the rate at which spikes are assigned to the wrong cluster-has been the rate of interspike interval (ISI) violations. Despite their near ubiquitous usage in spike sorting, our understanding of how exactly ISI violations relate to FDR, as well as best practices for using ISI violations as a quality metric, remains limited. Here, we describe an analytical solution that can be used to predict FDR from the ISI violation rate (ISIv). We test this model in silico through Monte Carlo simulation and apply it to publicly available spike-sorted electrophysiology datasets. We find that the relationship between ISIv and FDR is highly nonlinear, with additional dependencies on firing frequency, the correlation in activity between neurons, and contaminant neuron count. Predicted median FDRs in public datasets recorded in mice were found to range from 3.1 to 50.0%. We found that stochasticity in the occurrence of ISI violations as well as uncertainty in cluster-specific parameters make it difficult to predict FDR for single clusters with high confidence but that FDR can be estimated accurately across a population of clusters. Our findings will help the growing community of researchers using extracellular electrophysiology assess spike-sorting accuracy in a principled manner.

Stimulus type shapes the topology of cellular functional networks in mouse visual cortex

On the timescale of sensory processing, neuronal networks have relatively fixed anatomical connectivity, while functional interactions between neurons can vary depending on the ongoing activity of the neurons within the network. We thus hypothesized that different types of stimuli could lead those networks to display stimulus-dependent functional connectivity patterns. To test this hypothesis, we analyzed single-cell resolution electrophysiological data from the Allen Institute, with simultaneous recordings of stimulus-evoked activity from neurons across 6 different regions of mouse visual cortex. Comparing the functional connectivity patterns during different stimulus types, we made several nontrivial observations: (1) while the frequencies of different functional motifs were preserved across stimuli, the identities of the neurons within those motifs changed; (2) the degree to which functional modules are contained within a single brain region increases with stimulus complexity. Altogether, our work reveals unexpected stimulus-dependence to the way groups of neurons interact to process incoming sensory information.

Tonically active GABAergic neurons in the dorsal periaqueductal gray control instinctive escape in mice

Escape behavior is a set of locomotor actions that move an animal away from threat. While these actions can be stereotyped, it is advantageous for survival that they are flexible.1,2,3 For example, escape probability depends on predation risk and competing motivations,4,5,6,7,8,9,10,11 and flight to safety requires continuous adjustments of trajectory and must terminate at the appropriate place and time.12,13,14,15,16 This degree of flexibility suggests that modulatory components, like inhibitory networks, act on the neural circuits controlling instinctive escape.17,18,19,20,21,22 In mice, the decision to escape from imminent threats is implemented by a feedforward circuit in the midbrain, where excitatory vesicular glutamate transporter 2-positive (VGluT2+) neurons in the dorsal periaqueductal gray (dPAG) compute escape initiation and escape vigor.23,24,25 Here we tested the hypothesis that local GABAergic neurons within the dPAG control escape behavior by setting the excitability of the dPAG escape network. Using in vitro patch-clamp and in vivo neural activity recordings, we found that vesicular GABA transporter-positive (VGAT+) dPAG neurons fire action potentials tonically in the absence of synaptic inputs and are a major source of inhibition to VGluT2+ dPAG neurons. Activity in VGAT+ dPAG cells transiently decreases at escape onset and increases during escape, peaking at escape termination. Optogenetically increasing or decreasing VGAT+ dPAG activity changes the probability of escape when the stimulation is delivered at threat onset and the duration of escape when delivered after escape initiation. We conclude that the activity of tonically firing VGAT+ dPAG neurons sets a threshold for escape initiation and controls the execution of the flight action.

Neural population dynamics in songbird RA and HVC during learned motor-vocal behavior

Complex, learned motor behaviors involve the coordination of large-scale neural activity across multiple brain regions, but our understanding of the population-level dynamics within different regions tied to the same behavior remains limited. Here, we investigate the neural population dynamics underlying learned vocal production in awake-singing songbirds. We use Neuropixels probes to record the simultaneous extracellular activity of populations of neurons in two regions of the vocal motor pathway. In line with observations made in non-human primates during limb-based motor tasks, we show that the population-level activity in both the premotor nucleus HVC and the motor nucleus RA is organized on low-dimensional neural manifolds upon which coordinated neural activity is well described by temporally structured trajectories during singing behavior. Both the HVC and RA latent trajectories provide relevant information to predict vocal sequence transitions between song syllables. However, the dynamics of these latent trajectories differ between regions. Our state-space models suggest a unique and continuous-over-time correspondence between the latent space of RA and vocal output, whereas the corresponding relationship for HVC exhibits a higher degree of neural variability. We then demonstrate that comparable high-fidelity reconstruction of continuous vocal outputs can be achieved from HVC and RA neural latents and spiking activity. Unlike those that use spiking activity, however, decoding models using neural latents generalize to novel sub-populations in each region, consistent with the existence of preserved manifolds that confine vocal-motor activity in HVC and RA.

Open Data In Neurophysiology: Advancements, Solutions & Challenges

Across the life sciences, an ongoing effort over the last 50 years has made data and methods more reproducible and transparent. This openness has led to transformative insights and vastly accelerated scientific progress1,2. For example, structural biology3 and genomics4,5 have undertaken systematic collection and publication of protein sequences and structures over the past half-century, and these data have led to scientific breakthroughs that were unthinkable when data collection first began (e.g.6). We believe that neuroscience is poised to follow the same path, and that principles of open data and open science will transform our understanding of the nervous system in ways that are impossible to predict at the moment. To this end, new social structures along with active and open scientific communities are essential7 to facilitate and expand the still limited adoption of open science practices in our field8. Unified by shared values of openness, we set out to organize a symposium for Open Data in Neuroscience (ODIN) to strengthen our community and facilitate transformative neuroscience research at large. In this report, we share what we learned during this first ODIN event. We also lay out plans for how to grow this movement, document emerging conversations, and propose a path toward a better and more transparent science of tomorrow.

Low-intensity pulsed ultrasound stimulation (LIPUS) modulates microglial activation following intracortical microelectrode implantation

Microglia are important players in surveillance and repair of the brain. Implanting an electrode into the cortex activates microglia, produces an inflammatory cascade, triggers the foreign body response, and opens the blood-brain barrier. These changes can impede intracortical brain-computer interfaces performance. Using two-photon imaging of implanted microelectrodes, we test the hypothesis that low-intensity pulsed ultrasound stimulation can reduce microglia-mediated neuroinflammation following the implantation of microelectrodes. In the first week of treatment, we found that low-intensity pulsed ultrasound stimulation increased microglia migration speed by 128%, enhanced microglia expansion area by 109%, and a reduction in microglial activation by 17%, indicating improved tissue healing and surveillance. Microglial coverage of the microelectrode was reduced by 50% and astrocytic scarring by 36% resulting in an increase in recording performance at chronic time. The data indicate that low-intensity pulsed ultrasound stimulation helps reduce the foreign body response around chronic intracortical microelectrodes.

A Modular 512-Channel Neural Signal Acquisition ASIC for High-Density 4096 Channel Electrophysiology

The complexity of information processing in the brain requires the development of technologies that can provide spatial and temporal resolution by means of dense electrode arrays paired with high-channel-count signal acquisition electronics. In this work, we present an ultra-low noise modular 512-channel neural recording circuit that is scalable to up to 4096 simultaneously recording channels. The neural readout application-specific integrated circuit (ASIC) uses a dense 8.2 mm × 6.8 mm 2D layout to enable high-channel count, creating an ultra-light 350 mg flexible module. The module can be deployed on headstages for small animals like rodents and songbirds, and it can be integrated with a variety of electrode arrays. The chip was fabricated in a TSMC 0.18 µm 1.8 V CMOS technology and dissipates a total of 125 mW. Each DC-coupled channel features a gain and bandwidth programmable analog front-end along with 14 b analog-to-digital conversion at speeds up to 30 kS/s. Additionally, each front-end includes programmable electrode plating and electrode impedance measurement capability. We present both standalone and in vivo measurements results, demonstrating the readout of spikes and field potentials that are modulated by a sensory input.

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.

Brain implantation of tissue-level-soft bioelectronics via embryonic development

The design of bioelectronics capable of stably tracking brain-wide, single-cell, and millisecond-resolved neural activities in the developing brain is critical to the study of neuroscience and neurodevelopmental disorders. During development, the three-dimensional (3D) structure of the vertebrate brain arises from a 2D neural plate 1,2 . These large morphological changes previously posed a challenge for implantable bioelectronics to track neural activity throughout brain development 3-9 . Here, we present a tissue-level-soft, sub-micrometer-thick, stretchable mesh microelectrode array capable of integrating into the embryonic neural plate of vertebrates by leveraging the 2D-to-3D reconfiguration process of the tissue itself. Driven by the expansion and folding processes of organogenesis, the stretchable mesh electrode array deforms, stretches, and distributes throughout the entire brain, fully integrating into the 3D tissue structure. Immunostaining, gene expression analysis, and behavioral testing show no discernible impact on brain development or function. The embedded electrode array enables long-term, stable, brain-wide, single-unit-single-spike-resolved electrical mapping throughout brain development, illustrating how neural electrical activities and population dynamics emerge and evolve during brain development.

Auditory oddball responses in the human subthalamic nucleus and substantia nigra pars reticulata

The auditory oddball is a mainstay in research on attention, novelty, and sensory prediction. How this task engages subcortical structures like the subthalamic nucleus and substantia nigra pars reticulata is unclear. We administered an auditory OB task while recording single unit activity (35 units) and local field potentials (57 recordings) from the subthalamic nucleus and substantia nigra pars reticulata of 30 patients with Parkinson's disease undergoing deep brain stimulation surgery. We found tone modulated and oddball modulated units in both regions. Population activity differentiated oddball from standard trials from 200 ms to 1000 ms after the tone in both regions. In the substantia nigra, beta band activity in the local field potential was decreased following oddball tones. The oddball related activity we observe may underlie attention, sensory prediction, or surprise-induced motor suppression.

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.

Facilitating the Sharing of Electrophysiology Data Analysis Results Through In-Depth Provenance Capture

Scientific research demands reproducibility and transparency, particularly in data-intensive fields like electrophysiology. Electrophysiology data are typically analyzed using scripts that generate output files, including figures. Handling these results poses several challenges due to the complexity and iterative nature of the analysis process. These stem from the difficulty to discern the analysis steps, parameters, and data flow from the results, making knowledge transfer and findability challenging in collaborative settings. Provenance information tracks data lineage and processes applied to it, and provenance capture during the execution of an analysis script can address those challenges. We present Alpaca (Automated Lightweight Provenance Capture), a tool that captures fine-grained provenance information with minimal user intervention when running data analysis pipelines implemented in Python scripts. Alpaca records inputs, outputs, and function parameters and structures information according to the W3C PROV standard. We demonstrate the tool using a realistic use case involving multichannel local field potential recordings of a neurophysiological experiment, highlighting how the tool makes result details known in a standardized manner in order to address the challenges of the analysis process. Ultimately, using Alpaca will help to represent results according to the FAIR principles, which will improve research reproducibility and facilitate sharing the results of data analyses.

Ultraflexible electrodes for recording neural activity in the mouse spinal cord during motor behavior

Implantable electrode arrays are powerful tools for directly interrogating neural circuitry in the brain, but implementing this technology in the spinal cord in behaving animals has been challenging due to the spinal cord's significant motion with respect to the vertebral column during behavior. Consequently, the individual and ensemble activity of spinal neurons processing motor commands remains poorly understood. Here, we demonstrate that custom ultraflexible 1-μm-thick polyimide nanoelectronic threads can conduct laminar recordings of many neuronal units within the lumbar spinal cord of unrestrained, freely moving mice. The extracellular action potentials have high signal-to-noise ratio, exhibit well-isolated feature clusters, and reveal diverse patterns of activity during locomotion. Furthermore, chronic recordings demonstrate the stable tracking of single units and their functional tuning over multiple days. This technology provides a path for elucidating how spinal circuits compute motor actions.

Spike sorting with Kilosort4

Spike sorting is the computational process of extracting the firing times of single neurons from recordings of local electrical fields. This is an important but hard problem in neuroscience, made complicated by the nonstationarity of the recordings and the dense overlap in electrical fields between nearby neurons. To address the spike-sorting problem, we have been openly developing the Kilosort framework. Here we describe the various algorithmic steps introduced in different versions of Kilosort. We also report the development of Kilosort4, a version with substantially improved performance due to clustering algorithms inspired by graph-based approaches. To test the performance of Kilosort, we developed a realistic simulation framework that uses densely sampled electrical fields from real experiments to generate nonstationary spike waveforms and realistic noise. We found that nearly all versions of Kilosort outperformed other algorithms on a variety of simulated conditions and that Kilosort4 performed best in all cases, correctly identifying even neurons with low amplitudes and small spatial extents in high drift conditions.

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/.

Theta Phase Entrainment of Single-Cell Spiking in Rat Somatosensory Barrel Cortex and Secondary Visual Cortex Is Enhanced during Multisensory Discrimination Behavior

Phase entrainment of cells by theta oscillations is thought to globally coordinate the activity of cell assemblies across different structures, such as the hippocampus and neocortex. This coordination is likely required for optimal processing of sensory input during recognition and decision-making processes. In quadruple-area ensemble recordings from male rats engaged in a multisensory discrimination task, we investigated phase entrainment of cells by theta oscillations in areas along the corticohippocampal hierarchy: somatosensory barrel cortex (S1BF), secondary visual cortex (V2L), perirhinal cortex (PER), and dorsal hippocampus (dHC). Rats discriminated between two 3D objects presented in tactile-only, visual-only, or both tactile and visual modalities. During task engagement, S1BF, V2L, PER, and dHC LFP signals showed coherent theta-band activity. We found phase entrainment of single-cell spiking activity to locally recorded as well as hippocampal theta activity in S1BF, V2L, PER, and dHC. While phase entrainment of hippocampal spikes to local theta oscillations occurred during sustained epochs of task trials and was nonselective for behavior and modality, somatosensory and visual cortical cells were only phase entrained during stimulus presentation, mainly in their preferred modality (S1BF, tactile; V2L, visual), with subsets of cells selectively phase-entrained during cross-modal stimulus presentation (S1BF: visual; V2L: tactile). This effect could not be explained by modulations of firing rate or theta amplitude. Thus, hippocampal cells are phase entrained during prolonged epochs, while sensory and perirhinal neurons are selectively entrained during sensory stimulus presentation, providing a brief time window for coordination of activity.

Cyborg islets: implanted flexible electronics reveal principles of human islet electrical maturation

Flexible electronics implanted during tissue formation enable chronic studies of tissue-wide electrophysiology. Here, we integrate tissue-like stretchable electronics during organogenesis of human stem cell-derived pancreatic islets, stably tracing single-cell extracellular spike bursting dynamics over months of functional maturation. Adapting spike sorting methods from neural studies reveals maturation-dependent electrical patterns of α and β-like (SC-α and β) cells, and their stimulus-coupled dynamics. We identified two major electrical states for both SC-α and β cells, distinguished by their glucose threshold for action potential firing. We find that improved hormone stimulation capacity during extended culture reflects increasing numbers of SC-α/β cells in low basal firing states, linked to energy and hormone metabolism gene upregulation. Continuous recording during further maturation by entrainment to daily feeding cycles reveals that circadian islet-level hormone secretion rhythms reflect sustained and coordinate oscillation of cell-level SC-α and β electrical activities. We find that this correlates with cell-cell communication and exocytic network induction, indicating a role for circadian rhythms in coordinating system-level stimulus-coupled responses. Cyborg islets thus reveal principles of electrical maturation that will be useful to build fully functional in vitro islets for research and therapeutic applications.

Parahippocampal neurons encode task-relevant information for goal-directed navigation

A behavioral strategy crucial to survival is directed navigation to a goal, such as a food or home location. One potential neural substrate for supporting goal-directed navigation is the parahippocampus, which contains neurons that represent an animal's position, orientation, and movement through the world, and that change their firing activity to encode behaviorally relevant variables such as reward. However, little prior work on the parahippocampus has considered how neurons encode variables during goal-directed navigation in environments that dynamically change. Here, we recorded single units from rat parahippocampal cortex while subjects performed a goal-directed task. The maze dynamically changed goal-locations via a visual cue on a trial-to-trial basis, requiring subjects to use cue-location associations to receive reward. We observed a mismatch-like signal, with elevated neural activity on incorrect trials, leading to rate-remapping. The strength of this remapping correlated with task performance. Recordings during open-field foraging allowed us to functionally define navigational coding for a subset of the neurons recorded in the maze. This approach revealed that head-direction coding units remapped more than other functional-defined units. Taken together, this work thus raises the possibility that during goal-directed navigation, parahippocampal neurons encode error information reflective of an animal's behavioral performance.

DeePhys: A machine learning-assisted platform for electrophysiological phenotyping of human neuronal networks

Reproducible functional assays to study in vitro neuronal networks represent an important cornerstone in the quest to develop physiologically relevant cellular models of human diseases. Here, we introduce DeePhys, a MATLAB-based analysis tool for data-driven functional phenotyping of in vitro neuronal cultures recorded by high-density microelectrode arrays. DeePhys is a modular workflow that offers a range of techniques to extract features from spike-sorted data, allowing for the examination of functional phenotypes both at the individual cell and network levels, as well as across development. In addition, DeePhys incorporates the capability to integrate novel features and to use machine-learning-assisted approaches, which facilitates a comprehensive evaluation of pharmacological interventions. To illustrate its practical application, we apply DeePhys to human induced pluripotent stem cell-derived dopaminergic neurons obtained from both patients and healthy individuals and showcase how DeePhys enables phenotypic screenings.

A model of human neural networks reveals NPTX2 pathology in ALS and FTLD

Human cellular models of neurodegeneration require reproducibility and longevity, which is necessary for simulating age-dependent diseases. Such systems are particularly needed for TDP-43 proteinopathies1, which involve human-specific mechanisms2-5 that cannot be directly studied in animal models. Here, to explore the emergence and consequences of TDP-43 pathologies, we generated induced pluripotent stem cell-derived, colony morphology neural stem cells (iCoMoNSCs) via manual selection of neural precursors6. Single-cell transcriptomics and comparison to independent neural stem cells7 showed that iCoMoNSCs are uniquely homogenous and self-renewing. Differentiated iCoMoNSCs formed a self-organized multicellular system consisting of synaptically connected and electrophysiologically active neurons, which matured into long-lived functional networks (which we designate iNets). Neuronal and glial maturation in iNets was similar to that of cortical organoids8. Overexpression of wild-type TDP-43 in a minority of neurons within iNets led to progressive fragmentation and aggregation of the protein, resulting in a partial loss of function and neurotoxicity. Single-cell transcriptomics revealed a novel set of misregulated RNA targets in TDP-43-overexpressing neurons and in patients with TDP-43 proteinopathies exhibiting a loss of nuclear TDP-43. The strongest misregulated target encoded the synaptic protein NPTX2, the levels of which are controlled by TDP-43 binding on its 3' untranslated region. When NPTX2 was overexpressed in iNets, it exhibited neurotoxicity, whereas correcting NPTX2 misregulation partially rescued neurons from TDP-43-induced neurodegeneration. Notably, NPTX2 was consistently misaccumulated in neurons from patients with amyotrophic lateral sclerosis and frontotemporal lobar degeneration with TDP-43 pathology. Our work directly links TDP-43 misregulation and NPTX2 accumulation, thereby revealing a TDP-43-dependent pathway of neurotoxicity.

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.

Sleep restores an optimal computational regime in cortical networks

Sleep is assumed to subserve homeostatic processes in the brain; however, the set point around which sleep tunes circuit computations is unknown. Slow-wave activity (SWA) is commonly used to reflect the homeostatic aspect of sleep; although it can indicate sleep pressure, it does not explain why animals need sleep. This study aimed to assess whether criticality may be the computational set point of sleep. By recording cortical neuron activity continuously for 10-14 d in freely behaving rats, we show that normal waking experience progressively disrupts criticality and that sleep functions to restore critical dynamics. Criticality is perturbed in a context-dependent manner, and waking experience is causal in driving these effects. The degree of deviation from criticality predicts future sleep/wake behavior more accurately than SWA, behavioral history or other neural measures. Our results demonstrate that perturbation and recovery of criticality is a network homeostatic mechanism consistent with the core, restorative function of sleep.

Regional specialization manifests in the reliability of neural population codes

The brain has the remarkable ability to learn and guide the performance of complex tasks. Decades of lesion studies suggest that different brain regions perform specialized functions in support of complex behaviors1-3. Yet recent large-scale studies of neural activity reveal similar patterns of activity and encoding distributed widely throughout the brain4-6. How these distributed patterns of activity and encoding are compatible with regional specialization of brain function remains unclear. Two frontal brain regions, the dorsal medial prefrontal cortex (dmPFC) and orbitofrontal cortex (OFC), are a paradigm of this conundrum. In the setting complex behaviors, the dmPFC is necessary for choosing optimal actions2,7,8, whereas the OFC is necessary for waiting for3,9 and learning from2,7,9-12 the outcomes of those actions. Yet both dmPFC and OFC encode both choice- and outcome-related quantities13-20. Here we show that while ensembles of neurons in the dmPFC and OFC of rats encode similar elements of a cognitive task with similar patterns of activity, the two regions differ in when that coding is consistent across trials ("reliable"). In line with the known critical functions of each region, dmPFC activity is more reliable when animals are making choices and less reliable preceding outcomes, whereas OFC activity shows the opposite pattern. Our findings identify the dynamic reliability of neural population codes as a mechanism whereby different brain regions may support distinct cognitive functions despite exhibiting similar patterns of activity and encoding similar quantities.

A dynamic subset of network interactions underlies tuning to natural movements in marmoset sensorimotor cortex

Mechanisms of computation in sensorimotor cortex must be flexible and robust to support skilled motor behavior. Patterns of neuronal coactivity emerge as a result of computational processes. Pairwise spike-time statistical relationships, across the population, can be summarized as a functional network (FN) which retains single-unit properties. We record populations of single-unit neural activity in forelimb sensorimotor cortex during prey-capture and spontaneous behavior and use an encoding model incorporating kinematic trajectories and network features to predict single-unit activity during forelimb movements. The contribution of network features depends on structured connectivity within strongly connected functional groups. We identify a context-specific functional group that is highly tuned to kinematics and reorganizes its connectivity between spontaneous and prey-capture movements. In the remaining context-invariant group, interactions are comparatively stable across behaviors and units are less tuned to kinematics. This suggests different roles in producing natural forelimb movements and contextualizes single-unit tuning properties within population dynamics.

Assessing cross-contamination in spike-sorted electrophysiology data

Recent advances in extracellular electrophysiology now facilitate the recording of spikes from hundreds or thousands of neurons simultaneously. This has necessitated both the development of new computational methods for spike sorting and better methods to determine spike sorting accuracy. One longstanding method of assessing the false discovery rate (FDR) of spike sorting - the rate at which spikes are misassigned to the wrong cluster - has been the rate of inter-spike-interval (ISI) violations. Despite their near ubiquitous usage in spike sorting, our understanding of how exactly ISI violations relate to FDR, as well as best practices for using ISI violations as a quality metric, remain limited. Here, we describe an analytical solution that can be used to predict FDR from ISI violation rate. We test this model in silico through Monte Carlo simulation, and apply it to publicly available spike-sorted electrophysiology datasets. We find that the relationship between ISI violation rate and FDR is highly nonlinear, with additional dependencies on firing rate, the correlation in activity between neurons, and contaminant neuron count. Predicted median FDRs in public datasets were found to range from 3.1% to 50.0%. We find that stochasticity in the occurrence of ISI violations as well as uncertainty in cluster-specific parameters make it difficult to predict FDR for single clusters with high confidence, but that FDR can be estimated accurately across a population of clusters. Our findings will help the growing community of researchers using extracellular electrophysiology assess spike sorting accuracy in a principled manner.

Adaptive frequency-domain filtering for neural signal preprocessing

Electrical interference from various sources is a common issue for experimental extracellular electrophysiology recordings collected using multi-electrode array neural recording systems. This interference deteriorates the signal-to-noise ratio (SNR) of the raw electrophysiology signals and hampers the accuracy of data post-processing using techniques such as spike-sorting. Traditional signal processing methods to digitally remove electrical interference during post-processing include bandpass filtering to limit the signal to the relevant spectral range of the biological data, e.g., the spikes band (300 Hz - 7 kHz), targeted notch filtering to remove power line interference from standard alternating current mains electricity and common reference removal to minimize noise common to all electrodes. These methods require a priori knowledge of the frequency of the interfering signal source to address the unique electromagnetic interference environment of each experimental setup. We discuss an adaptive method for automatically removing narrow-band electrical interference through a spectral peak detection and removal (SPDR) step that can be applied during post-processing of the recorded data, based on the intuition that tall, narrowband signals localized in the signal spectrum correspond to interference, rather than the activity of neurons. A spectral peak prominence (SPP) threshold is used to detect these peaks in the frequency domain, which will then be removed via notch filtering. We applied this method to simulated waveforms and also experimental electrophysiology data collected from cerebral organoids to demonstrate its effectiveness for removing unwanted interference without significantly distorting the neural signals. We discuss that proper selection of the SPP threshold is required to avoid over-filtering, which can result in distortion of the electrophysiology data. We also compare the firing-rate activity in the filtered electrophysiology with fluorescence calcium imaging, a secondary cellular activity marker, to quantify signal distortion and provide bounds on SNR-based optimization of the SPP threshold. The adaptive filtering technique demonstrated in this paper is a powerful method that can automatically detect and remove interband interference in recorded neural signals, potentially enabling data collection in more naturalistic settings where external interference signals are difficult to eliminate.

In vitro modeling of the human dopaminergic system using spatially arranged ventral midbrain-striatum-cortex assembloids

Ventral midbrain dopaminergic neurons project to the striatum as well as the cortex and are involved in movement control and reward-related cognition. In Parkinson's disease, nigrostriatal midbrain dopaminergic neurons degenerate and cause typical Parkinson's disease motor-related impairments, while the dysfunction of mesocorticolimbic midbrain dopaminergic neurons is implicated in addiction and neuropsychiatric disorders. Study of the development and selective neurodegeneration of the human dopaminergic system, however, has been limited due to the lack of an appropriate model and access to human material. Here, we have developed a human in vitro model that recapitulates key aspects of dopaminergic innervation of the striatum and cortex. These spatially arranged ventral midbrain-striatum-cortical organoids (MISCOs) can be used to study dopaminergic neuron maturation, innervation and function with implications for cell therapy and addiction research. We detail protocols for growing ventral midbrain, striatal and cortical organoids and describe how they fuse in a linear manner when placed in custom embedding molds. We report the formation of functional long-range dopaminergic connections to striatal and cortical tissues in MISCOs, and show that injected, ventral midbrain-patterned progenitors can mature and innervate the tissue. Using these assembloids, we examine dopaminergic circuit perturbations and show that chronic cocaine treatment causes long-lasting morphological, functional and transcriptional changes that persist upon drug withdrawal. Thus, our method opens new avenues to investigate human dopaminergic cell transplantation and circuitry reconstruction as well as the effect of drugs on the human dopaminergic system.

Using Human-Induced Pluripotent Stem Cell Derived Neurons on Microelectrode Arrays to Model Neurological Disease: A Review

In situ physiological signals of in vitro neural disease models are essential for studying pathogenesis and drug screening. Currently, an increasing number of in vitro neural disease models are established using human-induced pluripotent stem cell (hiPSC) derived neurons (hiPSC-DNs) to overcome interspecific gene expression differences. Microelectrode arrays (MEAs) can be readily interfaced with two-dimensional (2D), and more recently, three-dimensional (3D) neural stem cell-derived in vitro models of the human brain to monitor their physiological activity in real time. Therefore, MEAs are emerging and useful tools to model neurological disorders and disease in vitro using human iPSCs. This is enabling a real-time window into neuronal signaling at the network scale from patient derived. This paper provides a comprehensive review of MEA's role in analyzing neural disease models established by hiPSC-DNs. It covers the significance of MEA fabrication, surface structure and modification schemes for hiPSC-DNs culturing and signal detection. Additionally, this review discusses advances in the development and use of MEA technology to study in vitro neural disease models, including epilepsy, autism spectrum developmental disorder (ASD), and others established using hiPSC-DNs. The paper also highlights the application of MEAs combined with hiPSC-DNs in detecting in vitro neurotoxic substances. Finally, the future development and outlook of multifunctional and integrated devices for in vitro medical diagnostics and treatment are discussed.

DREDge: robust motion correction for high-density extracellular recordings across species

High-density microelectrode arrays (MEAs) have opened new possibilities for systems neuroscience in human and non-human animals, but brain tissue motion relative to the array poses a challenge for downstream analyses, particularly in human recordings. We introduce DREDge (Decentralized Registration of Electrophysiology Data), a robust algorithm which is well suited for the registration of noisy, nonstationary extracellular electrophysiology recordings. In addition to estimating motion from spikes in the action potential (AP) frequency band, DREDge enables automated tracking of motion at high temporal resolution in the local field potential (LFP) frequency band. In human intraoperative recordings, which often feature fast (period <1s) motion, DREDge correction in the LFP band enabled reliable recovery of evoked potentials, and significantly reduced single-unit spike shape variability and spike sorting error. Applying DREDge to recordings made during deep probe insertions in nonhuman primates demonstrated the possibility of tracking probe motion of centimeters across several brain regions while simultaneously mapping single unit electrophysiological features. DREDge reliably delivered improved motion correction in acute mouse recordings, especially in those made with an recent ultra-high density probe. We also implemented a procedure for applying DREDge to recordings made across tens of days in chronic implantations in mice, reliably yielding stable motion tracking despite changes in neural activity across experimental sessions. Together, these advances enable automated, scalable registration of electrophysiological data across multiple species, probe types, and drift cases, providing a stable foundation for downstream scientific analyses of these rich datasets.

Modified Neuropixels probes for recording human neurophysiology in the operating room

Neuropixels are silicon-based electrophysiology-recording probes with high channel count and recording-site density. These probes offer a turnkey platform for measuring neural activity with single-cell resolution and at a scale that is beyond the capabilities of current clinically approved devices. Our team demonstrated the first-in-human use of these probes during resection surgery for epilepsy or tumors and deep brain stimulation electrode placement in patients with Parkinson's disease. Here, we provide a better understanding of the capabilities and challenges of using Neuropixels as a research tool to study human neurophysiology, with the hope that this information may inform future efforts toward regulatory approval of Neuropixels probes as research devices. In perioperative procedures, the major concerns are the initial sterility of the device, maintaining a sterile field during surgery, having multiple referencing and grounding schemes available to de-noise recordings (if necessary), protecting the silicon probe from accidental contact before insertion and obtaining high-quality action potential and local field potential recordings. The research team ensures that the device is fully operational while coordinating with the surgical team to remove sources of electrical noise that could otherwise substantially affect the signals recorded by the sensitive hardware. Prior preparation using the equipment and training in human clinical research and working in operating rooms maximize effective communication within and between the teams, ensuring high recording quality and minimizing the time added to the surgery. The perioperative procedure requires ~4 h, and the entire protocol requires multiple weeks.

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.

NeuSort: an automatic adaptive spike sorting approach with neuromorphic models

Objective.Spike sorting, a critical step in neural data processing, aims to classify spiking events from single electrode recordings based on different waveforms. This study aims to develop a novel online spike sorter, NeuSort, using neuromorphic models, with the ability to adaptively adjust to changes in neural signals, including waveform deformations and the appearance of new neurons.Approach.NeuSort leverages a neuromorphic model to emulate template-matching processes. This model incorporates plasticity learning mechanisms inspired by biological neural systems, facilitating real-time adjustments to online parameters.Results.Experimental findings demonstrate NeuSort's ability to track neuron activities amidst waveform deformations and identify new neurons in real-time. NeuSort excels in handling non-stationary neural signals, significantly enhancing its applicability for long-term spike sorting tasks. Moreover, its implementation on neuromorphic chips guarantees ultra-low energy consumption during computation.Significance.NeuSort caters to the demand for real-time spike sorting in brain-machine interfaces through a neuromorphic approach. Its unsupervised, automated spike sorting process makes it a plug-and-play solution for online spike sorting.