Single Unit publications

To understand the underlying mechanisms of progressive neurophysiological phenomena, neural interfaces should interact bi-directionally with brain circuits over extended periods of time. However, such interfaces remain limited by the foreign body response that stems from the chemo-mechanical mismatch between the probes and the neural tissues. To address this challenge, we developed a multifunctional sensing and actuation platform consisting of multimaterial fibers intimately integrated within a soft hydrogel matrix mimicking the brain tissue. These hybrid devices possess adaptive bending stiffness determined by the hydration states of the hydrogel matrix. This enables their direct insertion into the deep brain regions, while minimizing tissue damage associated with the brain micromotion after implantation. The hydrogel hybrid devices permit electrophysiological, optogenetic, and behavioral studies of neural circuits with minimal foreign body responses and tracking of stable isolated single neuron potentials in freely moving mice over 6 months following implantation.

A dopaminergic neural circuit projecting from the VTA to the hindbrain regulates feeding by transducing meal satiation signals.

Background Cerebellar electrical stimulation has shown promise in improving motor recovery post-stroke in both rodent and human studies. Past studies have used motor evoked potentials (MEPs) to evaluate how cerebellar stimulation modulates ongoing activity in the cortex, but the underlying mechanisms are incompletely understood. Here we used invasive electrophysiological recordings from the intact and stroke-injured rodent primary motor cortex (M1) to assess how epidural cerebellar stimulation modulates neural dynamics at the level of single neurons as well as at the level of mesoscale dynamics. Methods We recorded single unit spiking and local field potentials (LFPs) in both the intact and acutely stroke-injured M1 contralateral to the stimulated cerebellum in adult Long-Evans rats under anesthesia. We analyzed changes in the firing rates of single units, the extent of synchronous spiking and power spectral density (PSD) changes in LFPs during and post-stimulation. Results Our results show that post-stimulation, the firing rates of a majority of M1 neurons changed significantly with respect to their baseline rates. These firing rate changes were diverse in character, as the firing rate of some neurons increased while others decreased. Additionally, these changes started to set in during stimulation. Furthermore, cross-correlation analysis showed a significant increase in coincident firing amongst neuronal pairs. Interestingly, this increase in synchrony was unrelated to the direction of firing rate change. We also found that neuronal ensembles derived through principal component analysis were more active post-stimulation. Lastly, these changes occurred without a significant change in the overall spectral power of LFPs post-stimulation. Conclusions Our results show that cerebellar stimulation caused significant, long-lasting changes in the activity patterns of M1 neurons by altering firing rates, boosting neural synchrony and increasing neuronal assemblies’ activation strength. Our study provides evidence that cerebellar stimulation can directly modulate cortical dynamics. Since these results are present in the perilesional cortex, our data might also help explain the facilitatory effects of cerebellar stimulation post-stroke.

Acoustic communication signals are inevitably challenged by ambient noise. In response to noise, many animals adjust their calls to maintain signal detectability. However, the mechanisms by which the auditory system adapts to the adjusted pulses are unclear. Our previous study revealed that the echolocating bat, Hipposideros pratti , increased its pulse intensity in the presence of background white noise. In vivo single-neuron recording demonstrated that the auditory midbrain neurons tuned to the second harmonic (H2 neurons) increased their minimal threshold (MT) to a similar degree as the increment of pulse intensity in the presence of the background noise. Furthermore, the H2 neurons exhibited consistent spike rates at their best amplitudes and sharper intensity tuning with background white noise compared with silent conditions. The previous data indicated that sound intensity analysis by auditory midbrain neurons was adapted to the increased pulse intensity in the same noise condition. This study further examined the echolocation pulse frequency and frequency analysis of auditory midbrain neurons with noise conditions. The data revealed that H. pratti did not shift the resting frequency in the presence of background noise. The auditory midbrain neuronal frequency analysis highly linked to processing the resting frequency with the presence of noise by presenting the constant best frequency (BF), frequency sensitivity, and frequency selectivity. Thus, our results suggested that auditory midbrain neuronal responses in background white noise are adapted to process echolocation pulses in the noise conditions.

Understanding how neural populations encode natural stimuli with complex spatiotemporal structure to give rise to perception remains a central problem in neuroscience. Here we investigated population coding of natural communication stimuli by hindbrain neurons within the electrosensory system of weakly electric fish Apteronotus leptorhynchus . Overall, we found that simultaneously recorded neural activities were correlated: signal but not noise correlations were variable depending on the stimulus waveform as well as the distance between neurons. Combining the neural activities using an equal-weight sum gave rise to discrimination performance between different stimulus waveforms that was limited by redundancy introduced by noise correlations. However, using an evolutionary algorithm to assign different weights to individual neurons before combining their activities (i.e., a weighted sum) gave rise to increased discrimination performance by revealing synergistic interactions between neural activities. Our results thus demonstrate that correlations between the neural activities of hindbrain electrosensory neurons can enhance information about the structure of natural communication stimuli that allow for reliable discrimination between different waveforms by downstream brain areas.

Rapid sensory adaptation is observed across all sensory systems, and strongly shapes sensory percepts in complex sensory environments. Yet despite its ubiquity and likely necessity for survival, the mechanistic basis is poorly understood. A wide range of primarily in vitro and anesthetized studies have demonstrated the emergence of adaptation at the level of primary sensory cortex, with only modest signatures in earlier stages of processing. The nature of rapid adaptation and how it shapes sensory representations during wakefulness, and thus the potential role in perceptual adaptation, is underexplored, as are the mechanisms that underlie this phenomenon. To address these knowledge gaps, we recorded spiking activity in primary somatosensory cortex (S1) and the upstream ventral posteromedial (VPm) thalamic nucleus in the vibrissa pathway of awake male and female mice, and quantified responses to whisker stimuli delivered in isolation and embedded in an adapting sensory background. We found that cortical sensory responses were indeed adapted by persistent sensory stimulation; putative excitatory neurons were profoundly adapted, and inhibitory neurons only modestly so. Further optogenetic manipulation experiments and network modeling suggest this largely reflects adaptive changes in synchronous thalamic firing combined with robust engagement of feedforward inhibition, with little contribution from synaptic depression. Taken together, these results suggest that cortical adaptation in the regime explored here results from changes in the timing of thalamic input, and the way in which this differentially impacts cortical excitation and feedforward inhibition, pointing to a prominent role of thalamic gating in rapid adaptation of primary sensory cortex. SIGNIFICANCE STATEMENT Rapid adaptation of sensory activity strongly shapes representations of sensory inputs across all sensory pathways over the timescale of seconds, and has profound effects on sensory perception. Despite its ubiquity and theoretical role in the efficient encoding of complex sensory environments, the mechanistic basis is poorly understood, particularly during wakefulness. In this study in the vibrissa pathway of awake mice, we show that cortical representations of sensory inputs are strongly shaped by rapid adaptation, and that this is mediated primarily by adaptive gating of the thalamic inputs to primary sensory cortex and the differential way in which these inputs engage cortical subpopulations of neurons.

Survival in a dynamic environment requires animals to plan future actions based on past sensory evidence, known as motor planning. However, the neuronal circuits underlying this crucial brain function remain elusive. Here, we employ projection-specific imaging and perturbation methods to investigate the direct pathway linking two key nodes in the motor planning network, the secondary motor cortex (M2) and the midbrain superior colliculus (SC), in mice performing a memory-dependent perceptual decision task. We find dynamic coding of choice information in SC-projecting M2 neurons during motor planning and execution, and disruption of this information by inhibiting M2 terminals in SC selectively impaired decision maintenance. Furthermore, we show that while both excitatory and inhibitory SC neurons receive synaptic inputs from M2, these SC subpopulations display differential temporal patterns in choice coding during behavior. Our results reveal the dynamic recruitment of the premotor-collicular pathway as a circuit mechanism for motor planning.

Transcranial focused ultrasound (tFUS) is a promising neuromodulation technique, but its mechanisms remain unclear. We hypothesize that if tFUS parameters exhibit distinct modulation effects in different neuron populations, then the mechanism can be understood through identifying unique features in these neuron populations. In this work, we investigate the effect of tFUS stimulation on different functional neuron types in in vivo anesthetized rodent brains. Single neuron recordings were separated into regular-spiking and fast-spiking units based on their extracellular spike shapes acquired through intracranial electrophysiological recordings, and further validated in transgenic optogenetic mice models of light-excitable excitatory and inhibitory neurons. We show that excitatory and inhibitory neurons are intrinsically different in response to ultrasound pulse repetition frequency (PRF). The results suggest that we can preferentially target specific neuron types noninvasively by tuning the tFUS PRF. Chemically deafened rats and genetically deafened mice were further tested for validating the directly local neural effects induced by tFUS without potential auditory confounds.

Primate social communication depends on the perceptual integration of visual and auditory cues, reflected in the multimodal mixing of sensory signals in certain cortical areas. The macaque cortical face patch network, identified through visual, face-selective responses measured with fMRI, is assumed to contribute to visual social interactions. However, whether face patch neurons are also influenced by acoustic information, such as the auditory component of a natural vocalization, remains unknown. Here, we recorded single-unit activity in the anterior fundus (AF) face patch, in the superior temporal sulcus, and anterior medial (AM) face patch, on the undersurface of the temporal lobe, in macaques presented with audiovisual, visual-only, and auditory-only renditions of natural movies of macaques vocalizing. The results revealed that 76% of neurons in face patch AF were significantly influenced by the auditory component of the movie, most often through enhancement of visual responses but sometimes in response to the auditory stimulus alone. By contrast, few neurons in face patch AM exhibited significant auditory responses or modulation. Control experiments in AF used an animated macaque avatar to demonstrate, first, that the structural elements of the face were often essential for audiovisual modulation and, second, that the temporal modulation of the acoustic stimulus was more important than its frequency spectrum. Together, these results identify a striking contrast between two face patches and specifically identify AF as playing a potential role in the integration of audiovisual cues during natural modes of social communication.

Working memory (WM), the ability to actively hold information in memory over a delay period of seconds, is a fundamental constituent of cognition. Delay-period activity in sensory cortices has been observed in WM tasks, but whether and when the activity plays a functional role for memory maintenance remains unclear. Here, we investigated the causal role of auditory cortex (AC) for memory maintenance in mice performing an auditory WM task. Electrophysiological recordings revealed that AC neurons were active not only during the presentation of the auditory stimulus but also early in the delay period. Furthermore, optogenetic suppression of neural activity in AC during the stimulus epoch and early delay period impaired WM performance, whereas suppression later in the delay period did not. Thus, AC is essential for information encoding and maintenance in auditory WM task, especially during the early delay period.

The LIM homeodomain transcription factor ISL1 is essential for the different aspects of neuronal development and maintenance. In order to study the role of ISL1 in the auditory system, we generated a transgenic mouse (Tg) expressing Isl1 under the Pax2 promoter control. We previously reported a progressive age-related decline in hearing and abnormalities in the inner ear, medial olivocochlear system, and auditory midbrain of these Tg mice. In this study, we investigated how Isl1 overexpression affects sound processing by the neurons of the inferior colliculus (IC). We recorded extracellular neuronal activity and analyzed the responses of IC neurons to broadband noise, clicks, pure tones, two-tone stimulation and frequency-modulated sounds. We found that Tg animals showed a higher inhibition as displayed by two-tone stimulation; they exhibited a wider dynamic range, lower spontaneous firing rate, longer first spike latency and, in the processing of frequency modulated sounds, showed a prevalence of high-frequency inhibition. Functional changes were accompanied by a decreased number of calretinin and parvalbumin positive neurons, and an increased expression of vesicular GABA/glycine transporter and calbindin in the IC of Tg mice, compared to wild type animals. The results further characterize abnormal sound processing in the IC of Tg mice and demonstrate that major changes occur on the side of inhibition.

The hippocampal local field potential (LFP) exhibits a strong correlation with behavior. During rest, the theta rhythm is not prominent, but during active behavior, there are strong rhythms in the theta, theta harmonics, and gamma ranges. With increasing running velocity, theta, theta harmonics and gamma increase in power and in cross-frequency coupling, suggesting that neural entrainment is a direct consequence of the total excitatory input. While it is common to study the parametric range between the LFP and its complementing power spectra between deep rest and epochs of high running velocity, it is also possible to explore how the spectra degrades as the energy is completely quenched from the system. Specifically, it is unknown whether the 1/f slope is preserved as synaptic activity becomes diminished, as low frequencies are generated by large pools of neurons while higher frequencies comprise the activity of more local neuronal populations. To test this hypothesis, we examined rat LFPs recorded from the hippocampus and entorhinal cortex during barbiturate overdose euthanasia. Within the hippocampus, the initial stage entailed a quasi-stationary LFP state with a power-law feature in the power spectral density. In the second stage, there was a successive erosion of power from high- to low-frequencies in the second stage that continued until the only dominant remaining power was <20 Hz. This stage was followed by a rapid collapse of power spectrum toward the absolute electrothermal noise background. As the collapse of activity occurred later in hippocampus compared with medial entorhinal cortex, it suggests that the ability of a neural network to maintain the 1/f slope with decreasing energy is a function of general connectivity. Broadly, these data support the energy cascade theory where there is a cascade of energy from large cortical populations into smaller loops, such as those that supports the higher frequency gamma rhythm. As energy is pulled from the system, neural entrainment at gamma frequency (and higher) decline first. The larger loops, comprising a larger population, are fault-tolerant to a point capable of maintaining their activity before a final collapse.

The marmoset has emerged as a promising primate model system, in particular for visual neuroscience. Many common experimental paradigms rely on head fixation and an extended period of eye fixation during the presentation of salient visual stimuli. Both of these behavioral requirements can be challenging for marmosets. Here, we present two methodological developments, each addressing one of these difficulties. First, we show that it is possible to use a standard eye-tracking system without head fixation to assess visual behavior in the marmoset. Eye-tracking quality from head-free animals is sufficient to obtain precise psychometric functions from a visual acuity task. Second, we introduce a novel method for efficient receptive field (RF) mapping that does not rely on moving stimuli but uses fast flashing annuli and wedges. We present data recorded during head-fixation in areas V1 and V6 and show that RF locations are readily obtained within a short period of recording time. Thus, the methodological advancements presented in this work will contribute to establish the marmoset as a valuable model in neuroscience.

In this work the impact of two widely used anesthetics on the electrical activity of auditory brainstem neurons was studied during postnatal development. Spontaneous electrical activity in neonate rats of either sex was analyzed through a ventral craniotomy in mechanically ventilated pups to carry out patch clamp and multi-electrode electrophysiology recordings in the medial region of the superior olivary complex (SOC) between birth (postnatal day 0, P0) and P12. Recordings were obtained in pups anesthetized with the injectable mix of ketamine/xylazine (K/X mix), with the volatile anesthetic isoflurane (ISO), or in pups anesthetized with K/X mix that were also exposed to ISO. The results of patch clamp recordings demonstrate for the first time that olivary and periolivary neurons in the medial region of the SOC fire bursts of action potentials. The results of multielectrode recordings suggest that the firing pattern of single units recorded in K/X mix is similar to that recorded in ISO anesthetized rat pups. Taken together, the results of this study provide a framework to use injectable and volatile anesthetics for future studies to obtain functional information on the activity of medial superior olivary neurons in vivo .

Brain asymmetry in the sensitivity to spectrotemporal modulation is an established functional feature that underlies the perception of speech and music. Recent studies in humans suggest separable neural mechanisms must underlie these cognitively complex functions. The left Auditory Cortex (ACx) is believed to specialize in processing fast temporal components of speech sounds, and the right ACx slower components. However, the circuit features and neural computations behind these lateralized spectrotemporal processes are poorly understood. To answer these mechanistic questions we use mice, an animal model that captures some relevant features of human communication systems. In this study, we screened for circuit features that could subserve temporal integration differences between the left and right ACx. We mapped excitatory input to principal neurons in all cortical layers and found significantly stronger recurrent connections in the superficial layers of the right ACx compared to the left. We hypothesized that the underlying recurrent neural dynamics would exhibit differential characteristic timescales corresponding to their hemispheric specialization. To investigate, we recorded spike trains from awake mice and estimated the network time constants using a quasi-Bayesian method to combine evidence from multiple weak signal-to-noise ratio neurons. We found longer temporal integration windows in the right ACx compared to the left as predicted by stronger recurrent excitation. Our study shows direct evidence that stronger recurrent synaptic connections lead to longer network time scales. These findings support speech processing theories that purport asymmetric integration time constants is a crucial feature of lateralization in auditory processing.

Development is a time of great opportunity. A heightened period of neural plasticity contributes to dramatic improvements in perceptual, motor, and cognitive skills. However, developmental plasticity poses a risk: greater malleability of neural circuits exposes them to environmental factors that may impede behavioral maturation. While these risks are well-established prior to sexual maturity (i.e., critical periods), the degree of neural vulnerability during adolescence remains uncertain. To address this question, we induced a transient period of hearing loss (HL) spanning adolescence in the gerbil, confirmed by assessment of circulating sex hormones, and asked whether behavioral and neural deficits are diminished. Wireless recordings were obtained from auditory cortex neurons during perceptual task performance, and within-session behavioral and neural sensitivity were compared. We found that a transient period of adolescent HL caused a significant perceptual deficit (i.e., amplitude modulation detection thresholds) that could be attributed to degraded auditory cortex processing, as confirmed with both single neuron and population-level analyses. In contrast, perceptual deficits did not occur when HL of the same duration was induced in adulthood. To determine whether degraded auditory cortex encoding was attributable to an intrinsic change, we obtained auditory cortex brain slices from adolescent HL animals, and recorded synaptic and discharge properties from auditory cortex pyramidal neurons. There was a clear and novel phenotype, distinct from critical period HL: excitatory postsynaptic potential amplitudes were elevated in adolescent HL animals, whereas inhibitory postsynaptic potentials were unchanged. This is in contrast to critical period deprivation, where there are large changes to synaptic inhibition. Taken together, these results show that diminished adolescent sensory experience can cause long-lasting behavioral deficits that originate, in part, from a dysfunctional cortical circuit. Abstract Figure Summary of experimental design and main findings.

Cortical spreading depolarization (CSD) is a propagating wave of tissue depolarization characterized by a large increase of extracellular potassium concentration and prolonged subsequent electrical silencing of neurons. Waves of CSD arise spontaneously in various acute neurological settings, including migraine aura and ischemic stroke. Recently, we have reported that pan-inhibition of adrenergic receptors (AdRs) facilitates the normalization of extracellular potassium after acute photothrombotic stroke in mice. Here, we have extended that mechanistic study to ask whether AdR antagonists also modify the dynamics of KCl-induced CSD and post-CSD recovery in vivo. Spontaneous neural activity and KCl-induced CSD were visualized by cortex-wide transcranial Ca 2+ imaging in G-CaMP7 transgenic mice. AdR antagonism decreased the recurrence of CSD waves and accelerated the post-CSD recovery of neural activity. Two-photon imaging revealed that astrocytes exhibited aberrant Ca 2+ signaling after passage of the CSD wave. This astrocytic Ca 2+ activity was diminished by the AdR antagonists. Furthermore, AdR pan-antagonism facilitated the normalization of the extracellular potassium level after CSD, which paralleled the recovery of neural activity. These observations add support to the proposal that neuroprotective effects of AdR pan-antagonism arise from accelerated normalization of extracellular K + levels in the setting of acute brain injury.

In sensory systems, representational features of increasing complexity emerge at successive stages of processing. In the mammalian auditory pathway, the clearest change from brainstem to cortex is defined by what is lost, not by what is gained, in that high-fidelity temporal coding becomes increasingly restricted to slower acoustic modulation rates.1,2 Here, we explore the idea that sluggish temporal processing is more than just an inability for fast processing, but instead reflects an emergent specialization for encoding sound features that unfold on very slow timescales.3,4 We performed simultaneous single unit ensemble recordings from three hierarchical stages of auditory processing in awake mice - the inferior colliculus (IC), medial geniculate body of the thalamus (MGB) and primary auditory cortex (A1). As expected, temporal coding of brief local intervals (0.001 - 0.1 s) separating consecutive noise bursts was robust in the IC and declined across MGB and A1. By contrast, slowly developing (∼1 s period) global rhythmic patterns of inter-burst interval sequences strongly modulated A1 spiking, were weakly captured by MGB neurons, and not at all by IC neurons. Shifts in stimulus regularity were not represented by changes in A1 spike rates, but rather in how the spikes were arranged in time. These findings show that low-level auditory neurons with fast timescales encode isolated sound features but not the longer gestalt, while the extended timescales in higher-level areas can facilitate sensitivity to slower contextual changes in the sensory environment.

Objective. The rapid acceleration of tools for recording neuronal populations and targeted optogenetic manipulation has enabled real-time, feedback control of neuronal circuits in the brain. Continuously-graded control of measured neuronal activity poses a wide range of technical challenges, which we address through a combination of optogenetic stimulation and a state-space optimal control framework implemented in the thalamocortical circuit of the awake mouse. Approach. Closed-loop optogenetic control of neurons was performed in real-time via stimulation of channelrhodopsin-2 expressed in the somatosensory thalamus of the head-fixed mouse. A state-space linear dynamical system model structure was used to approximate the light-to-spiking input-output relationship in both single-neuron as well as multi-neuron scenarios when recording from multielectrode arrays. These models were utilized to design state feedback controller gains by way of linear quadratic optimal control and were also used online for estimation of state feedback, where a parameter-adaptive Kalman filter provided robustness to model-mismatch. Main results. This model-based control scheme proved effective for feedback control of single-neuron firing rate in the thalamus of awake animals. Notably, the graded optical actuation utilized here did not synchronize simultaneously recorded neurons, but heterogeneity across the neuronal population resulted in a varied response to stimulation. Simulated multi-output feedback control provided better control of a heterogeneous population and demonstrated how the approach generalizes beyond single-neuron applications. Significance. To our knowledge, this work represents the first experimental application of state space model-based feedback control for optogenetic stimulation. In combination with linear quadratic optimal control, the approaches laid out and tested here should generalize to future problems involving the control of highly complex neural circuits. More generally, feedback control of neuronal circuits opens the door to adaptively interacting with the dynamics underlying sensory, motor, and cognitive signaling, enabling a deeper understanding of circuit function and ultimately the control of function in the face of injury or disease.

A spinal cord injury usually spares some components of the locomotor circuitry. Deep brain stimulation (DBS) of the midbrain locomotor region and epidural electrical stimulation of the lumbar spinal cord (EES) are being used to tap into this spared circuitry to enable locomotion in humans with spinal cord injury. While appealing, the potential synergy between DBS and EES remains unknown. Here, we report the synergistic facilitation of locomotion when DBS is combined with EES in a rat model of severe contusion spinal cord injury leading to leg paralysis. However, this synergy requires high amplitudes of DBS, which triggers forced locomotion associated with stress responses. To suppress these undesired responses, we link DBS to the intention to walk, decoded from cortical activity using a robust, rapidly calibrated unsupervised learning algorithm. This contingency amplifies the supraspinal descending command while empowering the rats into volitional walking. However, the resulting improvements may not outweigh the complex technological framework necessary to establish viable therapeutic conditions.