Extracellularly Recorded Somatic and Neuritic Signal Shapes and Classification Algorithms for High-Density Microelectrode Array Electrophysiology

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.

Portrait of intense communications within microfluidic neural networks

In vitro model networks could provide cellular models of physiological relevance to reproduce and investigate the basic function of neural circuits on a chip in the laboratory. Several tools and methods have been developed since the past decade to build neural networks on a chip; among them, microfluidic circuits appear to be a highly promising approach. One of the numerous advantages of this approach is that it preserves stable somatic and axonal compartments over time due to physical barriers that prevent the soma from exploring undesired areas and guide neurites along defined pathways. As a result, neuron compartments can be identified and isolated, and their interconnectivity can be modulated to build a topological neural network (NN). Here, we have assessed the extent to which the confinement imposed by the microfluidic environment can impact cell development and shape NN activity. Toward that aim, microelectrode arrays have enabled the monitoring of the short- and mid-term evolution of neuron activation over the culture period at specific locations in organized (microfluidic) and random (control) networks. In particular, we have assessed the spike and burst rate, as well as the correlations between the extracted spike trains over the first stages of maturation. This study enabled us to observe intense neurite communications that would have been weaker and more delayed within random networks; the spiking rate, burst and correlations being reinforced over time in terms of number and amplitude, exceeding the electrophysiological features of standard cultures. Beyond the enhanced detection efficiency that was expected from the microfluidic channels, the confinement of cells seems to reinforce neural communications and cell development throughout the network.

Visuomotor learning promotes visually evoked activity in the medial prefrontal cortex

The medial prefrontal cortex (mPFC) is necessary for executing many learned associations between stimuli and movement. It is unclear, however, how activity in the mPFC evolves across learning, and how this activity correlates with sensory stimuli and the learned movements they evoke. To address these questions, we record cortical activity with widefield calcium imaging while mice learned to associate a visual stimulus with a forelimb movement. After learning, the mPFC shows stimulus-evoked activity both during task performance and during passive viewing, when the stimulus evokes no action. This stimulus-evoked activity closely tracks behavioral performance across training, with both exhibiting a marked increase between days when mice first learn the task, followed by a steady increase with further training. Electrophysiological recordings localized this activity to the secondary motor and anterior cingulate cortex. We conclude that learning a visuomotor task promotes a route for visual information to reach the prefrontal cortex.

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.

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.

Striatal activity topographically reflects cortical activity

The cortex projects to the dorsal striatum topographically^1,2 to regulate behaviour^3-5, but spiking activity in the two structures has previously been reported to have markedly different relations to sensorimotor events^6-9. Here we show that the relationship between activity in the cortex and striatum is spatiotemporally precise, topographic, causal and invariant to behaviour. We simultaneously recorded activity across large regions of the cortex and across the width of the dorsal striatum in mice that performed a visually guided task. Striatal activity followed a mediolateral gradient in which behavioural correlates progressed from visual cue to response movement to reward licking. The summed activity in each part of the striatum closely and specifically mirrored activity in topographically associated cortical regions, regardless of task engagement. This relationship held for medium spiny neurons and fast-spiking interneurons, whereas the activity of tonically active neurons differed from cortical activity with stereotypical responses to sensory or reward events. Inactivation of the visual cortex abolished striatal responses to visual stimuli, supporting a causal role of cortical inputs in driving the striatum. Striatal visual responses were larger in trained mice than untrained mice, with no corresponding change in overall activity in the visual cortex. Striatal activity therefore reflects a consistent, causal and scalable topographical mapping of cortical activity.

The Claustrum is Involved in Cognitive Processes Related to the Classical Conditioning of Eyelid Responses in Behaving Rabbits

It is assumed that the claustrum (CL) is involved in sensorimotor integration and cognitive processes. We recorded the firing activity of identified CL neurons during classical eyeblink conditioning in rabbits, using a delay paradigm in which a tone was presented as conditioned stimulus (CS), followed by a corneal air puff as unconditioned stimulus (US). Neurons were identified by their activation from motor (MC), cingulate (CC), and medial prefrontal (mPFC) cortices. CL neurons were rarely activated by single stimuli of any modality. In contrast, their firing was significantly modulated during the first sessions of paired CS/US presentations, but not in well-trained animals. Neuron firing rates did not correlate with the kinematics of conditioned responses (CRs). CL local field potentials (LFPs) changed their spectral power across learning and presented well-differentiated CL–mPFC/CL–MC network dynamics, as shown by crossfrequency spectral measurements. CL electrical stimulation did not evoke eyelid responses, even in trained animals. Silencing of synaptic transmission of CL neurons by the vINSIST method delayed the acquisition of CRs but did not affect their presentation rate. The CL plays an important role in the acquisition of associative learning, mostly in relation to the novelty of CS/US association, but not in the expression of CRs.

Large-Scale Mapping of Axonal Arbors Using High-Density Microelectrode Arrays

Understanding the role of axons in neuronal information processing is a fundamental task in neuroscience. Over the last years, sophisticated patch-clamp investigations have provided unexpected and exciting data on axonal phenomena and functioning, but there is still a need for methods to investigate full axonal arbors at sufficient throughput. Here, we present a new method for the simultaneous mapping of the axonal arbors of a large number of individual neurons, which relies on their extracellular signals that have been recorded with high-density microelectrode arrays (HD-MEAs). The segmentation of axons was performed based on the local correlation of extracellular signals. Comparison of the results with both, ground truth and receiver operator characteristics, shows that the new segmentation method outperforms previously used methods. Using a standard HD-MEA, we mapped the axonal arbors of 68 neurons in <6 h. The fully automated method can be extended to new generations of HD-MEAs with larger data output and is estimated to provide data of axonal arbors of thousands of neurons within recording sessions of a few hours.