High-density extracellular probes reveal dendritic backpropagation and facilitate neuron classification

Real-Time Precise Targeting of the Subthalamic Nucleus via Transfer Learning in a Rat Model of Parkinson's Disease Based on Microelectrode Arrays

In neurodegenerative disorders, neuronal firing patterns and oscillatory activity are remarkably altered in specific brain regions, which can serve as valuable biomarkers for the identification of deep brain regions. The subthalamic nucleus (STN) has been the primary target for DBS in patients with Parkinson's disease (PD). In this study, changes in the spike firing patterns and spectral power of local field potentials (LFPs) in the pre-STN (zona incerta, ZI) and post-STN (cerebral peduncle, cp) regions were investigated in PD rats, providing crucial evidence for the functional localization of the STN. Sixteen-channel microelectrode arrays (MEAs) with sites distributed at different depths and widths were utilized to record neuronal activities. The spikes in the STN exhibited higher firing rates than those in the ZI and cp. Furthermore, the LFP power in the delta band in the STN was the greatest, followed by that in the ZI, and was greater than that in the cp. Additionally, increased LFP power was observed in the beta bands in the STN. To identify the best performing classification model, we applied various convolutional neural networks (CNNs) based on transfer learning to analyze the recorded raw data, which were processed using the Gram matrix of the spikes and the fast Fourier transform of the LFPs. The best transfer learning model achieved an accuracy of 95.16%. After fusing the spike and LFP classification results, the time precision for processing the raw data reached 500 ms. The pretrained model, utilizing raw data, demonstrated the feasibility of employing transfer learning for training models on neural activity. This approach highlights the potential for functional localization within deep brain regions.

Development of wafer-scale multifunctional nanophotonic neural probes for brain activity mapping

Optical techniques, such as optogenetic stimulation and functional fluorescence imaging, have been revolutionary for neuroscience by enabling neural circuit analysis with cell-type specificity. To probe deep brain regions, implantable light sources are crucial. Silicon photonics, commonly used for data communications, shows great promise in creating implantable devices with complex optical systems in a compact form factor compatible with high volume manufacturing practices. This article reviews recent developments of wafer-scale multifunctional nanophotonic neural probes. The probes can be realized on 200 or 300 mm wafers in commercial foundries and integrate light emitters for photostimulation, microelectrodes for electrophysiological recording, and microfluidic channels for chemical delivery and sampling. By integrating active optical devices to the probes, denser emitter arrays, enhanced on-chip biosensing, and increased ease of use may be realized. Silicon photonics technology makes possible highly versatile implantable neural probes that can transform neuroscience experiments.

Nerve injury disrupts temporal processing in the spinal cord dorsal horn through alterations in PV+ interneurons

How mechanical allodynia following nerve injury is encoded in patterns of neural activity in the spinal cord dorsal horn (DH) remains incompletely understood. We address this in mice using the spared nerve injury model of neuropathic pain and in vivo electrophysiological recordings. Surprisingly, despite dramatic behavioral over-reactivity to mechanical stimuli following nerve injury, an overall increase in sensitivity or reactivity of DH neurons is not observed. We do, however, observe a marked decrease in correlated neural firing patterns, including the synchrony of mechanical stimulus-evoked firing, across the DH. Alterations in DH temporal firing patterns are recapitulated by silencing DH parvalbumin+ (PV+) interneurons, previously implicated in mechanical allodynia, as are allodynic pain-like behaviors. These findings reveal decorrelated DH network activity, driven by alterations in PV+ interneurons, as a prominent feature of neuropathic pain and suggest restoration of proper temporal activity as a potential therapeutic strategy to treat chronic neuropathic pain.

Discovering optimal features for neuron-type identification from extracellular recordings

Advancements in multichannel recordings of single-unit activity (SUA) in vivo present an opportunity to discover novel features of spatially-varying extracellularly-recorded action potentials (EAPs) that are useful for identifying neuron-types. Traditional approaches to classifying neuron-types often rely on computing EAP waveform features based on conventions of single-channel recordings and thus inherit their limitations. However, spatiotemporal EAP waveforms are the product of signals from underlying current sources being mixed within the extracellular space. We introduce a machine learning approach to demix the underlying sources of spatiotemporal EAP waveforms. Using biophysically realistic computational models, we simulate EAP waveforms and characterize them by the relative prevalence of these sources, which we use as features for identifying the neuron-types corresponding to recorded single units. These EAP sources have distinct spatial and multi-resolution temporal patterns that are robust to various sampling biases. EAP sources also are shared across many neuron-types, are predictive of gross morphological features, and expose underlying morphological domains. We then organize known neuron-types into a hierarchy of latent morpho-electrophysiological types based on differences in the source prevalences, which provides a multi-level classification scheme. We validate the robustness, accuracy, and interpretations of our demixing approach by analyzing simulated EAPs from morphologically detailed models with classification and clustering methods. This simulation-based approach provides a machine learning strategy for neuron-type identification.

High-Resolution Laminar Identification in Macaque Primary Visual Cortex Using Neuropixels Probes

Laminar electrode arrays allow simultaneous recording of activity of many cortical neurons and assignment to correct layers using current source density (CSD) analyses. Electrode arrays with 100-micron contact spacing can estimate borders between layer 4 versus superficial or deep layers, but in macaque primary visual cortex (V1) there are far more layers, such as 4A which is only 50-100 microns thick. Neuropixels electrode arrays have 20-micron spacing, and thus could potentially discern thinner layers and more precisely identify laminar borders. Here we show that CSD signals lack the spatial resolution required to take advantage of high density Neuropixels arrays and describe the development of approaches based on higher resolution electrical signals and analyses, including spike waveforms and spatial spread, unit density, high-frequency action potential (AP) power spectrum, temporal power change, and coherence spectrum, that afford far higher resolution of laminar distinctions, including the ability to precisely detect the borders of even the thinnest layers of V1.

Electrode Development for Epilepsy Diagnosis and Treatment

Recording neural activity has been a critical aspect in the diagnosis and treatment of patients with epilepsy. For those with intractable epilepsy, intracranial neural monitoring has been of substantial importance. Clinically, however, methods for recording neural information have remained essentially unchanged for decades. Over the last decade or so, rapid advances in electrode technology have begun to change this landscape. New systems allow for the observation of neural activity with high spatial resolution and, in some cases, at the level of the activity of individual neurons. These new tools have contributed greatly to our understanding of brain function and dysfunction. Here, the authors review the primary technologies currently in use in humans. The authors discuss other possible systems, some of the challenges which come along with these devices, and how they will become incorporated into the clinical workflow. Ultimately, the expectation is that these new, high-density, high-spatial-resolution recording systems will become a valuable part of the clinical arsenal used in the diagnosis and surgical management of epilepsy.

Sensory cortical ensembles exhibit differential coupling to ripples in distinct hippocampal subregions

Cortical neurons activated during recent experiences often reactivate with dorsal hippocampal CA1 ripples during subsequent rest. Less is known about cortical interactions with intermediate hippocampal CA1, whose connectivity, functions, and ripple events differ from dorsal CA1. We identified three clusters of putative excitatory neurons in mouse visual cortex that are preferentially excited together with either dorsal or intermediate CA1 ripples or suppressed before both ripples. Neurons in each cluster were evenly distributed across primary and higher visual cortices and co-active even in the absence of ripples. These ensembles exhibited similar visual responses but different coupling to thalamus and pupil-indexed arousal. We observed a consistent activity sequence preceding and predicting ripples: (1) suppression of ripple-suppressed cortical neurons, (2) thalamic silence, and (3) activation of intermediate CA1-ripple-activated cortical neurons. We propose that coordinated dynamics of these ensembles relay visual experiences to distinct hippocampal subregions for incorporation into different cognitive maps.

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.

Positive and biphasic extracellular waveforms correspond to return currents and axonal spikes

Multiple biophysical mechanisms may generate non-negative extracellular waveforms during action potentials, but the origin and prevalence of positive spikes and biphasic spikes in the intact brain are unknown. Using extracellular recordings from densely-connected cortical networks in freely-moving mice, we find that a tenth of the waveforms are non-negative. Positive phases of non-negative spikes occur in synchrony or just before wider same-unit negative spikes. Narrow positive spikes occur in isolation in the white matter. Isolated biphasic spikes are narrower than negative spikes, occurring right after spikes of verified inhibitory units. In CA1, units with dominant non-negative spikes exhibit place fields, phase precession, and phase-locking to ripples. Thus, near-somatic narrow positive extracellular potentials correspond to return currents, and isolated non-negative spikes correspond to axonal potentials. Identifying non-negative extracellular waveforms that correspond to non-somatic compartments during spikes can enhance the understanding of physiological and pathological neural mechanisms in intact animals.

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.

Characterization of extracellular spike waveforms recorded in wallaby primary visual cortex

Extracellular recordings were made from 642 units in the primary visual cortex (V1) of a highly visual marsupial, the Tammar wallaby. The receptive field (RF) characteristics of the cells were objectively estimated using the non-linear input model (NIM), and these were correlated with spike shapes. We found that wallaby cortical units had 68% regular spiking (RS), 12% fast spiking (FS), 4% triphasic spiking (TS), 5% compound spiking (CS) and 11% positive spiking (PS). RS waveforms are most often associated with recordings from pyramidal or spiny stellate cell bodies, suggesting that recordings from these cell types dominate in the wallaby cortex. In wallaby, 70-80% of FS and RS cells had orientation selective RFs and had evenly distributed linear and nonlinear RFs. We found that 47% of wallaby PS units were non-orientation selective and they were dominated by linear RFs. Previous studies suggest that the PS units represent recordings from the axon terminals of non-orientation selective cells originating in the lateral geniculate nucleus (LGN). If this is also true in wallaby, as strongly suggested by their low response latencies and bursty spiking properties, the results suggest that significantly more neurons in wallaby LGN are already orientation selective. In wallaby, less than 10% of recorded spikes had triphasic (TS) or sluggish compound spiking (CS) waveforms. These units had a mixture of orientation selective and non-oriented properties, and their cellular origins remain difficult to classify.

Retinal input integration in excitatory and inhibitory neurons in the mouse superior colliculus in vivo

The superior colliculus (SC) is a midbrain structure that receives inputs from retinal ganglion cells (RGCs). The SC contains one of the highest densities of inhibitory neurons in the brain but whether excitatory and inhibitory SC neurons differentially integrate retinal activity in vivo is still largely unknown. We recently established a recording approach to measure the activity of RGCs simultaneously with their postsynaptic SC targets in vivo, to study how SC neurons integrate RGC activity. Here, we employ this method to investigate the functional properties that govern retinocollicular signaling in a cell type-specific manner by identifying GABAergic SC neurons using optotagging in VGAT-ChR2 mice. Our results demonstrate that both excitatory and inhibitory SC neurons receive comparably strong RGC inputs and similar wiring rules apply for RGCs innervation of both SC cell types, unlike the cell type-specific connectivity in the thalamocortical system. Moreover, retinal activity contributed more to the spiking activity of postsynaptic excitatory compared to inhibitory SC neurons. This study deepens our understanding of cell type-specific retinocollicular functional connectivity and emphasizes that the two major brain areas for visual processing, the visual cortex and the SC, differently integrate sensory afferent inputs.

Opto-juxtacellular interrogation of neural circuits in freely moving mice

Neural circuits are assembled from an enormous variety of neuronal cell types. Although significant advances have been made in classifying neurons on the basis of morphological, molecular and electrophysiological properties, understanding how this diversity contributes to brain function during behavior has remained a major experimental challenge. Here, we present an extension to our previous protocol, in which we describe the technical procedures for performing juxtacellular opto-tagging of single neurons in freely moving mice by using Channelrhodopsin-2-expressing viral vectors. This method allows one to selectively target molecularly defined cell classes for in vivo single-cell recordings. The targeted cells can be labeled via juxtacellular procedures and further characterized via post-hoc morphological and molecular analysis. In its current form, the protocol allows multiple recording and labeling attempts to be performed within individual animals, by means of a mechanical pipette micropositioning system. We provide proof-of-principle validation of this technique by recording from Calbindin-positive pyramidal neurons in the mouse hippocampus during spatial exploration; however, this approach can easily be extended to other behaviors and cortical or subcortical areas. The procedures described here, from the viral injection to the histological processing of brain sections, can be completed in ~4-5 weeks.This protocol is an extension to: Nat. Protoc. 9, 2369-2381 (2014): https://doi.org/10.1038/nprot.2014.161.

Cortico-thalamo-cortical interactions modulate electrically evoked EEG responses in mice

Perturbational complexity analysis predicts the presence of consciousness in volunteers and patients by stimulating the brain with brief pulses, recording EEG responses, and computing their spatiotemporal complexity. We examined the underlying neural circuits in mice by directly stimulating cortex while recording with EEG and Neuropixels probes during wakefulness and isoflurane anesthesia. When mice are awake, stimulation of deep cortical layers reliably evokes locally a brief pulse of excitation, followed by a biphasic sequence of 120 ms profound off period and a rebound excitation. A similar pattern, partially attributed to burst spiking, is seen in thalamic nuclei and is associated with a pronounced late component in the evoked EEG. We infer that cortico-thalamo-cortical interactions drive the long-lasting evoked EEG signals elicited by deep cortical stimulation during the awake state. The cortical and thalamic off period and rebound excitation, and the late component in the EEG, are reduced during running and absent during anesthesia.

WaveMAP for identifying putative cell types from in vivo electrophysiology

Action potential spike widths are used to classify cell types as either excitatory or inhibitory; however, this approach obscures other differences in waveform shape useful for identifying more fine-grained cell types. Here, we present a protocol for using WaveMAP to generate nuanced average waveform clusters more closely linked to underlying cell types. We describe steps for installing WaveMAP, preprocessing data, and clustering waveform into putative cell types. We also detail cluster evaluation for functional differences and interpretation of WaveMAP output. For complete details on the use and execution of this protocol, please refer to Lee et al. (2021).¹.

Large-scale high-density brain-wide neural recording in nonhuman primates

High-density, integrated silicon electrodes have begun to transform systems neuroscience, by enabling large-scale neural population recordings with single cell resolution. Existing technologies, however, have provided limited functionality in nonhuman primate species such as macaques, which offer close models of human cognition and behavior. Here, we report the design, fabrication, and performance of Neuropixels 1.0-NHP, a high channel count linear electrode array designed to enable large-scale simultaneous recording in superficial and deep structures within the macaque or other large animal brain. These devices were fabricated in two versions: 4416 electrodes along a 45 mm shank, and 2496 along a 25 mm shank. For both versions, users can programmatically select 384 channels, enabling simultaneous multi-area recording with a single probe. We demonstrate recording from over 3000 single neurons within a session, and simultaneous recordings from over 1000 neurons using multiple probes. This technology represents a significant increase in recording access and scalability relative to existing technologies, and enables new classes of experiments involving fine-grained electrophysiological characterization of brain areas, functional connectivity between cells, and simultaneous brain-wide recording at scale.

Multimodal charting of molecular and functional cell states via in situ electro-sequencing

Paired mapping of single-cell gene expression and electrophysiology is essential to understand gene-to-function relationships in electrogenic tissues. Here, we developed in situ electro-sequencing (electro-seq) that combines flexible bioelectronics with in situ RNA sequencing to stably map millisecond-timescale electrical activity and profile single-cell gene expression from the same cells across intact biological networks, including cardiac and neural patches. When applied to human-induced pluripotent stem-cell-derived cardiomyocyte patches, in situ electro-seq enabled multimodal in situ analysis of cardiomyocyte electrophysiology and gene expression at the cellular level, jointly defining cell states and developmental trajectories. Using machine-learning-based cross-modal analysis, in situ electro-seq identified gene-to-electrophysiology relationships throughout cardiomyocyte development and accurately reconstructed the evolution of gene expression profiles based on long-term stable electrical measurements. In situ electro-seq could be applicable to create spatiotemporal multimodal maps in electrogenic tissues, potentiating the discovery of cell types and gene programs responsible for electrophysiological function and dysfunction.

Transcriptomic cell type structures in vivo neuronal activity across multiple timescales

Cell type is hypothesized to be a key determinant of a neuron's role within a circuit. Here, we examine whether a neuron's transcriptomic type influences the timing of its activity. We develop a deep-learning architecture that learns features of interevent intervals across timescales (ms to >30 min). We show that transcriptomic cell-class information is embedded in the timing of single neuron activity in the intact brain of behaving animals (calcium imaging and extracellular electrophysiology) as well as in a bio-realistic model of the visual cortex. Further, a subset of excitatory cell types are distinguishable but can be classified with higher accuracy when considering cortical layer and projection class. Finally, we show that computational fingerprints of cell types may be universalizable across structured stimuli and naturalistic movies. Our results indicate that transcriptomic class and type may be imprinted in the timing of single neuron activity across diverse stimuli.

Associations between in vitro, in vivo and in silico cell classes in mouse primary visual cortex

The brain consists of many cell classes yet in vivo electrophysiology recordings are typically unable to identify and monitor their activity in the behaving animal. Here, we employed a systematic approach to link cellular, multi-modal in vitro properties from experiments with in vivo recorded units via computational modeling and optotagging experiments. We found two one-channel and six multi-channel clusters in mouse visual cortex with distinct in vivo properties in terms of activity, cortical depth, and behavior. We used biophysical models to map the two one- and the six multi-channel clusters to specific in vitro classes with unique morphology, excitability and conductance properties that explain their distinct extracellular signatures and functional characteristics. These concepts were tested in ground-truth optotagging experiments with two inhibitory classes unveiling distinct in vivo properties. This multi-modal approach presents a powerful way to separate in vivo clusters and infer their cellular properties from first principles.

Associations between in vitro , in vivo and in silico cell classes in mouse primary visual cortex

The brain consists of many cell classes yet in vivo electrophysiology recordings are typically unable to identify and monitor their activity in the behaving animal. Here, we employed a systematic approach to link cellular, multi-modal in vitro properties from experiments with in vivo recorded units via computational modeling and optotagging experiments. We found two one-channel and six multi-channel clusters in mouse visual cortex with distinct in vivo properties in terms of activity, cortical depth, and behavior. We used biophysical models to map the two one- and the six multi-channel clusters to specific in vitro classes with unique morphology, excitability and conductance properties that explain their distinct extracellular signatures and functional characteristics. These concepts were tested in ground-truth optotagging experiments with two inhibitory classes unveiling distinct in vivo properties. This multi-modal approach presents a powerful way to separate in vivo clusters and infer their cellular properties from first principles.

Nerve injury disrupts temporal processing in the spinal cord dorsal horn through alterations in PV + interneurons

How mechanical allodynia following nerve injury is encoded in patterns of neural activity in the spinal cord dorsal horn (DH) is not known. We addressed this using the spared nerve injury model of neuropathic pain and in vivo electrophysiological recordings. Surprisingly, despite dramatic behavioral over-reactivity to mechanical stimuli following nerve injury, an overall increase in sensitivity or reactivity of DH neurons was not observed. We did, however, observe a marked decrease in correlated neural firing patterns, including the synchrony of mechanical stimulus-evoked firing, across the DH. Alterations in DH temporal firing patterns were recapitulated by silencing DH parvalbumin + (PV + ) inhibitory interneurons, previously implicated in mechanical allodynia, as were allodynic pain-like behaviors in mice. These findings reveal decorrelated DH network activity, driven by alterations in PV + interneurons, as a prominent feature of neuropathic pain, and suggest that restoration of proper temporal activity is a potential treatment for chronic neuropathic pain.

Tracking neural activity from the same cells during the entire adult life of mice

Stably recording the electrical activity of the same neurons over the adult life of an animal is important to neuroscience research and biomedical applications. Current implantable devices cannot provide stable recording on this timescale. Here, we introduce a method to precisely implant electronics with an open, unfolded mesh structure across multiple brain regions in the mouse. The open mesh structure forms a stable interwoven structure with the neural network, preventing probe drifting and showing no immune response and neuron loss during the year-long implantation. Rigorous statistical analysis, visual stimulus-dependent measurement and unbiased, machine-learning-based analysis demonstrated that single-unit action potentials have been recorded from the same neurons of behaving mice in a very long-term stable manner. Leveraging this stable structure, we demonstrated that the same neurons can be recorded over the entire adult life of the mouse, revealing the aging-associated evolution of single-neuron activities.

Nanomaterial-based microelectrode arrays for in vitro bidirectional brain-computer interfaces: a review

A bidirectional in vitro brain-computer interface (BCI) directly connects isolated brain cells with the surrounding environment, reads neural signals and inputs modulatory instructions. As a noninvasive BCI, it has clear advantages in understanding and exploiting advanced brain function due to the simplified structure and high controllability of ex vivo neural networks. However, the core of ex vivo BCIs, microelectrode arrays (MEAs), urgently need improvements in the strength of signal detection, precision of neural modulation and biocompatibility. Notably, nanomaterial-based MEAs cater to all the requirements by converging the multilevel neural signals and simultaneously applying stimuli at an excellent spatiotemporal resolution, as well as supporting long-term cultivation of neurons. This is enabled by the advantageous electrochemical characteristics of nanomaterials, such as their active atomic reactivity and outstanding charge conduction efficiency, improving the performance of MEAs. Here, we review the fabrication of nanomaterial-based MEAs applied to bidirectional in vitro BCIs from an interdisciplinary perspective. We also consider the decoding and coding of neural activity through the interface and highlight the various usages of MEAs coupled with the dissociated neural cultures to benefit future developments of BCIs.

Strong tuning for stereoscopic depth indicates orientation-specific recurrent circuitry in tree shrew V1

Neurons in the primary visual cortex (V1) are tuned to specific disparities between the two retinal images, which form the neural substrate for stereoscopic vision. We show that V1 neurons in tree shrews, but not in mice, display highly selective responses to narrow ranges of disparity in random-dot stereograms. Surprisingly, V1 neurons in both species show similarly strong tuning to gratings of varying interocular phase differences. This stimulus-dependent dissociation of disparity tuning can be explained by a network model that combines both feedforward and recurrent connections. The features of the model connections are supported by cortical organizations specific to each species. We validate this model by identifying putative inhibitory neurons and confirming their predicted disparity tuning in both species. Together, our studies establish a foundation for using tree shrews in studying binocular vision and raise an exciting possibility of how cortical columns could be uniquely important in computing stereoscopic depth.

Functional imaging of brain organoids using high-density microelectrode arrays

ABSTRACT

Studies have provided evidence that human cerebral organoids (hCOs) recapitulate fundamental milestones of early brain development, but many important questions regarding their functionality and electrophysiological properties persist. High-density microelectrode arrays (HD-MEAs) represent an attractive analysis platform to perform functional studies of neuronal networks at the cellular and network scale. Here, we use HD-MEAs to derive large-scale electrophysiological recordings from sliced hCOs. We record the activity of hCO slices over several weeks and probe observed neuronal dynamics pharmacologically. Moreover, we present results on how the obtained recordings can be spike-sorted and subsequently studied across scales. For example, we show how to track single neurons across several days on the HD-MEA and how to infer axonal action potential velocities. We also infer putative functional connectivity from hCO recordings. The introduced methodology will contribute to a better understanding of developing neuronal networks in brain organoids and provide new means for their functional characterization.

IMPACT STATEMENT

Human cerebral organoids (hCOs) represent an attractive in vitro model system to study key physiological mechanisms underlying early neuronal network formation in tissue with healthy or disease-related genetic backgrounds. Despite remarkable advances in the generation of brain organoids, knowledge on the functionality of their neuronal circuits is still scarce. Here, we used complementary metal-oxide-semiconductor (CMOS)-based high-density microelectrode arrays (HD-MEAs) to perform large-scale recordings from sliced hCOs over several weeks and quantified their activity across scales. Using single-cell and network metrics, we were able to probe aspects of hCO neurophysiology that are more difficult to obtain with other techniques, such as patch clamping (lower yield) and calcium imaging (lower temporal resolution). These metrics included, for example, extracellular action potential (AP) waveform features and axonal AP velocity at the cellular level, as well as functional connectivity at the network level. Analysis was enabled by the large sensing area and the high spatiotemporal resolution provided by HD-MEAs, which allowed recordings from hundreds of neurons and spike sorting of their activity. Our results demonstrate that HD-MEAs provide a multi-purpose platform for the functional characterization of hCOs, which will be key in improving our understanding of this model system and assessing its relevance for translational research.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1557/s43577-022-00282-w.

Multi-regional module-based signal transmission in mouse visual cortex

The visual cortex is hierarchically organized, yet the presence of extensive recurrent and parallel pathways make it challenging to decipher how signals flow between neuronal populations. Here, we tracked the flow of spiking activity recorded from six interconnected levels of the mouse visual hierarchy. By analyzing leading and lagging spike-timing relationships among pairs of simultaneously recorded neurons, we created a cellular-scale directed network graph. Using a module-detection algorithm to cluster neurons based on shared connectivity patterns, we uncovered two multi-regional communication modules distributed across the hierarchy. The direction of signal flow both between and within these modules, differences in layer and area distributions, and distinct temporal dynamics suggest that one module transmits feedforward sensory signals, whereas the other integrates inputs for recurrent processing. These results suggest that multi-regional functional modules may be a fundamental feature of organization beyond cortical areas that supports signal propagation across hierarchical recurrent networks.

Large-scale neural recordings with single neuron resolution using Neuropixels probes in human cortex

Recent advances in multi-electrode array technology have made it possible to monitor large neuronal ensembles at cellular resolution in animal models. In humans, however, current approaches restrict recordings to a few neurons per penetrating electrode or combine the signals of thousands of neurons in local field potential (LFP) recordings. Here we describe a new probe variant and set of techniques that enable simultaneous recording from over 200 well-isolated cortical single units in human participants during intraoperative neurosurgical procedures using silicon Neuropixels probes. We characterized a diversity of extracellular waveforms with eight separable single-unit classes, with differing firing rates, locations along the length of the electrode array, waveform spatial spread and modulation by LFP events such as inter-ictal discharges and burst suppression. Although some challenges remain in creating a turnkey recording system, high-density silicon arrays provide a path for studying human-specific cognitive processes and their dysfunction at unprecedented spatiotemporal resolution.

A Multimodal Multi-Shank Fluorescence Neural Probe for Cell-Type-Specific Electrophysiology in Multiple Regions across a Neural Circuit

Cell-type-specific, activity-dependent electrophysiology can allow in-depth analysis of functional connectivity inside complex neural circuits composed of various cell types. To date, optics-based fluorescence recording devices enable monitoring cell-type-specific activities. However, the monitoring is typically limited to a single brain region, and the temporal resolution is significantly low. Herein, a multimodal multi-shank fluorescence neural probe that allows cell-type-specific electrophysiology from multiple deep-brain regions at a high spatiotemporal resolution is presented. A photodiode and an electrode-array pair are monolithically integrated on each tip of a minimal-form-factor silicon device. Both fluorescence and electrical signals are successfully measured simultaneously in GCaMP6f expressing mice, and the cell type from sorted neural spikes is identified. The probe's capability of combined electro-optical recordings for cell-type-specific electrophysiology at multiple brain regions within a neural circuit is demonstrated. The new experimental paradigm to enable the precise investigation of functional connectivity inside and across complex neural circuits composed of various cell types is expected.

CellExplorer: A framework for visualizing and characterizing single neurons

The large diversity of neuron types provides the means by which cortical circuits perform complex operations. Neuron can be described by biophysical and molecular characteristics, afferent inputs, and neuron targets. To quantify, visualize, and standardize those features, we developed the open-source, MATLAB-based framework CellExplorer. It consists of three components: a processing module, a flexible data structure, and a powerful graphical interface. The processing module calculates standardized physiological metrics, performs neuron-type classification, finds putative monosynaptic connections, and saves them to a standardized, yet flexible, machine-readable format. The graphical interface makes it possible to explore the computed features at the speed of a mouse click. The framework allows users to process, curate, and relate their data to a growing public collection of neurons. CellExplorer can link genetically identified cell types to physiological properties of neurons collected across laboratories and potentially lead to interlaboratory standards of single-cell metrics.

Non-linear dimensionality reduction on extracellular waveforms reveals cell type diversity in premotor cortex

Cortical circuits are thought to contain a large number of cell types that coordinate to produce behavior. Current in vivo methods rely on clustering of specified features of extracellular waveforms to identify putative cell types, but these capture only a small amount of variation. Here, we develop a new method (WaveMAP) that combines non-linear dimensionality reduction with graph clustering to identify putative cell types. We apply WaveMAP to extracellular waveforms recorded from dorsal premotor cortex of macaque monkeys performing a decision-making task. Using WaveMAP, we robustly establish eight waveform clusters and show that these clusters recapitulate previously identified narrow- and broad-spiking types while revealing previously unknown diversity within these subtypes. The eight clusters exhibited distinct laminar distributions, characteristic firing rate patterns, and decision-related dynamics. Such insights were weaker when using feature-based approaches. WaveMAP therefore provides a more nuanced understanding of the dynamics of cell types in cortical circuits.

Dentate spikes and external control of hippocampal function

Mouse hippocampus CA1 place-cell discharge typically encodes current location, but during slow gamma dominance (SGdom), when SG oscillations (30-50 Hz) dominate mid-frequency gamma oscillations (70-90 Hz) in CA1 local field potentials, CA1 discharge switches to represent distant recollected locations. We report that dentate spike type 2 (DSM) events initiated by medial entorhinal cortex II (MECII)→ dentate gyrus (DG) inputs promote SGdom and change excitation-inhibition coordinated discharge in DG, CA3, and CA1, whereas type 1 (DSL) events initiated by lateral entorhinal cortex II (LECII)→DG inputs do not. Just before SGdom, LECII-originating SG oscillations in DG and CA3-originating SG oscillations in CA1 phase and frequency synchronize at the DSM peak when discharge within DG and CA3 increases to promote excitation-inhibition cofiring within and across the DG→CA3→CA1 pathway. This optimizes discharge for the 5-10 ms DG-to-CA1 neuro-transmission that SGdom initiates. DSM properties identify extrahippocampal control of SGdom and a cortico-hippocampal mechanism that switches between memory-related modes of information processing.

Reconciling functional differences in populations of neurons recorded with two-photon imaging and electrophysiology

Extracellular electrophysiology and two-photon calcium imaging are widely used methods for measuring physiological activity with single-cell resolution across large populations of cortical neurons. While each of these two modalities has distinct advantages and disadvantages, neither provides complete, unbiased information about the underlying neural population. Here, we compare evoked responses in visual cortex recorded in awake mice under highly standardized conditions using either imaging of genetically expressed GCaMP6f or electrophysiology with silicon probes. Across all stimulus conditions tested, we observe a larger fraction of responsive neurons in electrophysiology and higher stimulus selectivity in calcium imaging, which was partially reconciled by applying a spikes-to-calcium forward model to the electrophysiology data. However, the forward model could only reconcile differences in responsiveness when restricted to neurons with low contamination and an event rate above a minimum threshold. This work established how the biases of these two modalities impact functional metrics that are fundamental for characterizing sensory-evoked responses.

Dataset of cortical activity recorded with high spatial resolution from anesthetized rats

Publicly available neural recordings obtained with high spatial resolution are scarce. Here, we present an electrophysiological dataset recorded from the neocortex of twenty rats anesthetized with ketamine/xylazine. The wideband, spontaneous recordings were acquired with a single-shank silicon-based probe having 128 densely-packed recording sites arranged in a 32 × 4 array. The dataset contains the activity of a total of 7126 sorted single units extracted from all layers of the cortex. Here, we share raw neural recordings, as well as spike times, extracellular spike waveforms and several properties of units packaged in a standardized electrophysiological data format. For technical validation of our dataset, we provide the distributions of derived single unit properties along with various spike sorting quality metrics. This large collection of in vivo data enables the investigation of the high-resolution electrical footprint of cortical neurons which in turn may aid their electrophysiology-based classification. Furthermore, the dataset might be used to study the laminar-specific neuronal activity during slow oscillation, a brain rhythm strongly involved in neural mechanisms underlying memory consolidation and sleep.

Microscale Physiological Events on the Human Cortical Surface

Despite ongoing advances in our understanding of local single-cellular and network-level activity of neuronal populations in the human brain, extraordinarily little is known about their "intermediate" microscale local circuit dynamics. Here, we utilized ultra-high-density microelectrode arrays and a rare opportunity to perform intracranial recordings across multiple cortical areas in human participants to discover three distinct classes of cortical activity that are not locked to ongoing natural brain rhythmic activity. The first included fast waveforms similar to extracellular single-unit activity. The other two types were discrete events with slower waveform dynamics and were found preferentially in upper cortical layers. These second and third types were also observed in rodents, nonhuman primates, and semi-chronic recordings from humans via laminar and Utah array microelectrodes. The rates of all three events were selectively modulated by auditory and electrical stimuli, pharmacological manipulation, and cold saline application and had small causal co-occurrences. These results suggest that the proper combination of high-resolution microelectrodes and analytic techniques can capture neuronal dynamics that lay between somatic action potentials and aggregate population activity. Understanding intermediate microscale dynamics in relation to single-cell and network dynamics may reveal important details about activity in the full cortical circuit.

Dendritic calcium signals in rhesus macaque motor cortex drive an optical brain-computer interface

Calcium imaging is a powerful tool for recording from large populations of neurons in vivo. Imaging in rhesus macaque motor cortex can enable the discovery of fundamental principles of motor cortical function and can inform the design of next generation brain-computer interfaces (BCIs). Surface two-photon imaging, however, cannot presently access somatic calcium signals of neurons from all layers of macaque motor cortex due to photon scattering. Here, we demonstrate an implant and imaging system capable of chronic, motion-stabilized two-photon imaging of neuronal calcium signals from macaques engaged in a motor task. By imaging apical dendrites, we achieved optical access to large populations of deep and superficial cortical neurons across dorsal premotor (PMd) and gyral primary motor (M1) cortices. Dendritic signals from individual neurons displayed tuning for different directions of arm movement. Combining several technical advances, we developed an optical BCI (oBCI) driven by these dendritic signalswhich successfully decoded movement direction online. By fusing two-photon functional imaging with CLARITY volumetric imaging, we verified that many imaged dendrites which contributed to oBCI decoding originated from layer 5 output neurons, including a putative Betz cell. This approach establishes new opportunities for studying motor control and designing BCIs via two photon imaging.

Cellular Classes in the Human Brain Revealed In Vivo by Heartbeat-Related Modulation of the Extracellular Action Potential Waveform

Determining cell types is critical for understanding neural circuits but remains elusive in the living human brain. Current approaches discriminate units into putative cell classes using features of the extracellular action potential (EAP); in absence of ground truth data, this remains a problematic procedure. We find that EAPs in deep structures of the brain exhibit robust and systematic variability during the cardiac cycle. These cardiac-related features refine neural classification. We use these features to link bio-realistic models generated from in vitro human whole-cell recordings of morphologically classified neurons to in vivo recordings. We differentiate aspiny inhibitory and spiny excitatory human hippocampal neurons and, in a second stage, demonstrate that cardiac-motion features reveal two types of spiny neurons with distinct intrinsic electrophysiological properties and phase-locking characteristics to endogenous oscillations. This multi-modal approach markedly improves cell classification in humans, offers interpretable cell classes, and is applicable to other brain areas and species.

Growing evidence for separate neural mechanisms for attention and consciousness

Our conscious experience of the world seems to go in lockstep with our attentional focus: We tend to see, hear, taste, and feel what we attend to, and vice versa. This tight coupling between attention and consciousness has given rise to the idea that these two phenomena are indivisible. In the late 1950s, the honoree of this special issue, Charles Eriksen, was among a small group of early pioneers that sought to investigate whether a transient increase in overall level of attention (alertness) in response to a noxious stimulus can be decoupled from conscious perception using experimental techniques. Recent years saw a similar debate regarding whether attention and consciousness are two dissociable processes. Initial evidence that attention and consciousness are two separate processes primarily rested on behavioral data. However, the past couple of years witnessed an explosion of studies aimed at testing this conjecture using neuroscientific techniques. Here we provide an overview of these and related empirical studies on the distinction between the neuronal correlates of attention and consciousness, and detail how advancements in theory and technology can bring about a more detailed understanding of the two. We argue that the most promising approach will combine ever-evolving neurophysiological and interventionist tools with quantitative, empirically testable theories of consciousness that are grounded in a mathematically formalized understanding of phenomenology.

Survey of spiking in the mouse visual system reveals functional hierarchy

The anatomy of the mammalian visual system, from the retina to the neocortex, is organized hierarchically¹. However, direct observation of cellular-level functional interactions across this hierarchy is lacking due to the challenge of simultaneously recording activity across numerous regions. Here we describe a large, open dataset-part of the Allen Brain Observatory²-that surveys spiking from tens of thousands of units in six cortical and two thalamic regions in the brains of mice responding to a battery of visual stimuli. Using cross-correlation analysis, we reveal that the organization of inter-area functional connectivity during visual stimulation mirrors the anatomical hierarchy from the Allen Mouse Brain Connectivity Atlas³. We find that four classical hierarchical measures-response latency, receptive-field size, phase-locking to drifting gratings and response decay timescale-are all correlated with the hierarchy. Moreover, recordings obtained during a visual task reveal that the correlation between neural activity and behavioural choice also increases along the hierarchy. Our study provides a foundation for understanding coding and signal propagation across hierarchically organized cortical and thalamic visual areas.

Synchrony of Spontaneous Burst Firing between Retinal Ganglion Cells Across Species

Neurons communicate with other neurons in response to environmental changes. Their goal is to transmit information to their targets reliably. A burst, which consists of multiple spikes within a short time interval, plays an essential role in enhancing the reliability of information transmission through synapses. In the visual system, retinal ganglion cells (RGCs), the output neurons of the retina, show bursting activity and transmit retinal information to the lateral geniculate neuron of the thalamus. In this study, to extend our interest to the population level, the burstings of multiple RGCs were simultaneously recorded using a multi-channel recording system. As the first step in network analysis, we focused on investigating the pairwise burst correlation between two RGCs. Furthermore, to assess if the population bursting is preserved across species, we compared the synchronized bursting of RGCs between marmoset monkey (callithrix jacchus), one species of the new world monkeys and mouse (C57BL/6J strain). First, monkey RGCs showed a larger number of spikes within a burst, while the inter-spike interval, burst duration, and inter-burst interval were smaller compared with mouse RGCs. Monkey RGCs showed a strong burst synchronization between RGCs, whereas mouse RGCs showed no correlated burst firing. Monkey RGC pairs showed significantly higher burst synchrony and mutual information than mouse RGC pairs did. Comprehensively, through this study, we emphasize that two species have a different bursting activity of RGCs and different burst synchronization suggesting two species have distinctive retinal processing.

Human Cerebrospinal Fluid Induces Neuronal Excitability Changes in Resected Human Neocortical and Hippocampal Brain Slices

Human cerebrospinal fluid (hCSF) has proven advantageous over conventional medium for culturing both rodent and human brain tissue. In addition, increased activity and synchrony, closer to the dynamic states exclusively recorded in vivo, were reported in rodent slices and cell cultures switching from artificial cerebrospinal fluid (aCSF) to hCSF. This indicates that hCSF possesses properties that are not matched by the aCSF, which is generally used for most electrophysiological recordings. To evaluate the possible significance of using hCSF as an electrophysiological recording medium, also for human brain tissue, we compared the network and single-cell firing properties of human brain slice cultures during perfusion with hCSF and aCSF. For measuring the overall activity from a majority of neurons within neocortical and hippocampal human slices, we used a microelectrode array (MEA) recording technique with 252 electrodes covering an area of 3.2 × 3.2 mm2. A second CMOS-based MEA with 4225 sensors on a 2 × 2 mm2 area was used for detailed mapping of action potential waveforms and cell identification. We found that hCSF increased the number of active electrodes and neurons and the firing rate of the neurons in the slices and induced an increase in the numbers of single channel and population bursts. Interestingly, not only an increase in the overall activity in the slices was observed, but a reconfiguration of the network could also be detected with specific activation and inactivation of subpopulations of neuronal ensembles. In conclusion, hCSF is an important component to consider for future human brain slice studies, especially for experiments designed to mimic parts of physiology and disease observed in vivo.