Publications using the RZ2

Theta oscillations in the primary visual cortex (VC) have been observed during running tasks, but the mechanism behind their generation is not well understood. Some studies have suggested that theta in the VC is locally generated, while others have proposed that it is volume conducted from the hippocampus. The present study aimed to investigate the relationship between hippocampal and VC LFP dynamics. Analysis of power spectral density revealed that LFP in the VC was similar to that in the hippocampus, but with lower overall magnitude. As running velocity increased, both the power and frequency of theta and its harmonics increased in the VC, similarly to what is observed in the hippocampus. Current source density analysis triggered to theta did not identify distinct current sources and sinks in the VC, supporting the idea that theta in the VC is conducted from the adjacent hippocampus. Phase coupling between theta, its harmonics, and gamma is a notable feature in the hippocampus, particularly in the lacunosum moleculare. While some evidence of coupling between theta and its harmonics in the VC was found, bicoherence estimates did not reveal significant phase coupling between theta and gamma. Similar results were seen in the cross-region bicoherence analysis, where theta showed strong coupling with its harmonics with increasing velocity. Thus, theta oscillations observed in the VC during running tasks are likely due to volume conduction from the hippocampus.

The perception of pitch is a fundamental percept, which is mediated by the auditory system, requiring the abstraction of stimulus properties related to the spectro-temporal structure of sound. Despite its importance, there is still debate as to the precise areas responsible for its encoding, which may be due to species differences or differences in the recording measures and choices of stimuli used in previous studies. Moreover, it was unknown whether the human brain contains pitch neurons and how distributed such neurons might be. Here, we present the first study to measure multiunit neural activity in response to pitch stimuli in the auditory cortex of intracranially implanted humans. The stimulus sets were regular-interval noise with a pitch strength that is related to the temporal regularity and a pitch value determined by the repetition rate and harmonic complexes. Specifically, we demonstrate reliable responses to these different pitch-inducing paradigms that are distributed throughout Heschl’s gyrus, rather than being localized to a particular region, and this finding was evident regardless of the stimulus presented. These data provide a bridge across animal and human studies and aid our understanding of the processing of a critical percept associated with acoustic stimuli.

The subthalamic nucleus (STN) is hypothesized to play a central role in neural processes that regulate self-control. Still uncertain, however, is how that brain structure participates in the dynamically evolving estimation of value that underlies the ability to delay gratification and wait patiently for a gain. To address that gap in knowledge, we studied the spiking activity of neurons in the STN of monkeys during a task in which animals were required to remain motionless for varying periods of time in order to obtain food reward. At the single-neuron and population levels, we found a cost–benefit integration between the desirability of the expected reward and the imposed delay to reward delivery, with STN signals that dynamically combined both attributes of the reward to form a single integrated estimate of value. This neural encoding of subjective value evolved dynamically across the waiting period that intervened after instruction cue. Moreover, this encoding was distributed inhomogeneously along the antero-posterior axis of the STN such that the most dorso-posterior-placed neurons represented the temporal discounted value most strongly. These findings highlight the selective involvement of the dorso-posterior STN in the representation of temporally discounted rewards. The combination of rewards and time delays into an integrated representation is essential for self-control, the promotion of goal pursuit, and the willingness to bear the costs of time delays.

Objective. Neuropixels (NP) probes are a significant advance in electrophysiological recording technology that enable monitoring of hundreds of neurons in the brain simultaneously at different depths. Application of this technology has been predominately in rodents, however widespread use in non-human primates (NHPs) such as rhesus macaques has been limited. In this study we sought to overcome two overarching challenges that impede acute NP implantation in NHPs: (1) traditional microdrive systems that mount to cephalic chambers are commonly used to access cortical areas for microelectrode recordings but are not designed to accommodate NP probes, and (2) NHPs have thick dura mater and tissue growth within the cephalic chambers which poses a challenge for insertion of the extremely fragile NP probe. Approach. In this study we present a novel NP guide tube system that can be adapted to commercial microdrive systems and demonstrate an implant method using the NP guide tube system. This system was developed using a combination of CAD design, 3D printing, and small part machining. Software programs, 3D Slicer and SolidWorks were used to target cortical areas, approximate recording depths and locations, and for in-silico implant testing. Main results. We performed in vivo testing to validate our methodology, successfully implanting, explanting, and reimplanting NP probes. We collected stable neurophysiological recordings in the premotor cortex of a rhesus macaque at rest and during performance of a reaching task. Significance. In this study we demonstrate a robust Neuropixels implant system that allows multiple penetrations with the same NP probe and share design files that will facilitate the adoption of this powerful recording technology for NHP studies.

Biocompatible and plastic neural interface devices allow for minimally invasive recording of brain activity. Increasing electrode density in such devices is essential for high-resolution neural recordings. Superimposing conductive leads in devices can help multiply the number of recording sites while keeping probes width small and suitable for implantation. However, because of leads' vertical proximity, this can create capacitive coupling (CC) between overlapping channels, which leads to crosstalk. Here, we present a thorough investigation of CC phenomenon in multi-gold layer thin-film multi-electrode arrays with a parylene C (PaC) insulation layer between superimposed leads. We also propose a guideline on the design, fabrication, and characterization of such type of neural interface devices for high spatial resolution recording. Our results demonstrate that the capacitance created through CC between superimposed tracks decreases non-linearly and then linearly with the increase of insulation thickness. We identify an optimal PaC insulation thickness that leads to a drastic reduction of CC between superimposed gold channels while not significantly increasing the overall device thickness. Finally, we show that double gold layer electrocorticography probes with the optimal insulation thickness exhibit similar performances in vivo when compared to single-layer devices. This confirms that these probes are adequate for high-quality neural recordings.

We report on the power spectral analysis of multiunit activity (MUA) of lumbar spinal interneurons during a locomotor task. In line with prior studies, we found evidence of greater power at the frequency of locomotion in high lumbar segments and during the flexion phase. Our results also confirm prior observations from our laboratory that the rhythmically active MUA behaves as a longitudinal standing wave of neural activation that is flexor dominant.

Application of closed-loop approaches in systems neuroscience and brain-computer interfaces holds great promise for revolutionizing our understanding of the brain and for developing novel neuromodulation strategies to restore lost function. The anterior forebrain mesocircuit (AFM) of the mammalian brain is hypothesized to underlie arousal regulation of the cortex and striatum, and support cognitive functions during wakefulness. Dysfunction of arousal regulation is hypothesized to contribute to cognitive dysfunctions in various neurological disorders, and most prominently in patients following traumatic brain injury (TBI). Several clinical studies have explored the use of daily central thalamic deep brain stimulation (CT-DBS) within the AFM to restore consciousness and executive attention in TBI patients. In this study, we explored the use of closed-loop CT-DBS in order to episodically regulate arousal of the AFM of a healthy non-human primate (NHP) with the goal of restoring behavioral performance. We used pupillometry and near real-time analysis of ECoG signals to episodically initiate closed-loop CT-DBS and here we report on our ability to enhance arousal and restore the animal's performance. The initial computer based approach was then experimentally validated using a customized clinical-grade DBS device, the DyNeuMo-X, a bi-directional research platform used for rapidly testing closed-loop DBS. The successful implementation of the DyNeuMo-X in a healthy NHP supports ongoing clinical trials employing the internal DyNeuMo system (NCT05437393, NCT05197816) and our goal of developing and accelerating the deployment of novel neuromodulation approaches to treat cognitive dysfunction in patients with structural brain injuries and other etiologies.

Intracortical microstimulation (ICMS) is commonly used in many experimental and clinical paradigms; however, its effects on the activation of neurons are still not completely understood. To document the responses of cortical neurons in awake nonhuman primates to stimulation, we recorded single-unit activity while delivering single-pulse stimulation via Utah arrays implanted in primary motor cortex (M1) of three macaque monkeys. Stimuli between 5 and 50 μA delivered to single channels reliably evoked spikes in neurons recorded throughout the array with delays of up to 12 ms. ICMS pulses also induced a period of inhibition lasting up to 150 ms that typically followed the initial excitatory response. Higher current amplitudes led to a greater probability of evoking a spike and extended the duration of inhibition. The likelihood of evoking a spike in a neuron was dependent on the spontaneous firing rate as well as the delay between its most recent spike time and stimulus onset. Tonic repetitive stimulation between 2 and 20 Hz often modulated both the probability of evoking spikes and the duration of inhibition; high-frequency stimulation was more likely to change both responses. On a trial-by-trial basis, whether a stimulus evoked a spike did not affect the subsequent inhibitory response; however, their changes over time were often positively or negatively correlated. Our results document the complex dynamics of cortical neural responses to electrical stimulation that need to be considered when using ICMS for scientific and clinical applications.

In bistable perception, observers experience alternations between two interpretations of an unchanging stimulus. Neurophysiological studies of bistable perception typically partition neural measurements into stimulus-based epochs and assess neuronal differences between epochs based on subjects' perceptual reports. Computational studies replicate statistical properties of percept durations with modeling principles like competitive attractors or Bayesian inference. However, bridging neuro-behavioral findings with modeling theory requires the analysis of single-trial dynamic data. Here, we propose an algorithm for extracting nonstationary timeseries features from single-trial electrocorticography (ECoG) data. We applied the proposed algorithm to 5-min ECoG recordings from human primary auditory cortex obtained during perceptual alternations in an auditory triplet streaming task (six subjects: four male, two female). We report two ensembles of emergent neuronal features in all trial blocks. One ensemble consists of periodic functions that encode a stereotypical response to the stimulus. The other comprises more transient features and encodes dynamics associated with bistable perception at multiple time scales: minutes (within-trial alternations), seconds (duration of individual percepts), and milliseconds (switches between percepts). Within the second ensemble, we identified a slowly drifting rhythm that correlates with the perceptual states and several oscillators with phase shifts near perceptual switches. Projections of single-trial ECoG data onto these features establish low-dimensional attractor-like geometric structures invariant across subjects and stimulus types. These findings provide supporting neural evidence for computational models with oscillatory-driven attractor-based principles. The feature extraction techniques described here generalize across recording modality and are appropriate when hypothesized low-dimensional dynamics characterize an underlying neural system. SIGNIFICANCE STATEMENT Irrespective of the sensory modality, neurophysiological studies of multistable perception have typically investigated events time-locked to the perceptual switching rather than the time course of the perceptual states per se. Here, we propose an algorithm that extracts neuronal features of bistable auditory perception from largescale single-trial data while remaining agnostic to the subject's perceptual reports. The algorithm captures the dynamics of perception at multiple timescales, minutes (within-trial alternations), seconds (durations of individual percepts), and milliseconds (timing of switches), and distinguishes attributes of neural encoding of the stimulus from those encoding the perceptual states. Finally, our analysis identifies a set of latent variables that exhibit alternating dynamics along a low-dimensional manifold, similar to trajectories in attractor-based models for perceptual bistability.

The primate brain is equipped to learn and remember newly encountered visual stimuli such as faces and objects. In the macaque inferior temporal (IT) cortex, neurons mark the familiarity of a visual stimulus through response modification, often involving a decrease in spiking rate. Here, we investigate the emergence of this neural plasticity by longitudinally tracking IT neurons during several weeks of familiarization with face images. We found that most neurons in the anterior medial (AM) face patch exhibited a gradual decline in their late-phase visual responses to multiple stimuli. Individual neurons varied from days to weeks in their rates of plasticity, with time constants determined by the number of days of exposure rather than the cumulative number of presentations. We postulate that the sequential recruitment of neurons with experience-modified responses may provide an internal and graded measure of familiarity strength, which is a key mnemonic component of visual recognition.

Sleep is ubiquitous and essential, but its mechanisms remain unclear. Studies in animals and humans have provided insights of sleep at vastly different spatiotemporal scales. However, challenges remain to integrate local and global information of sleep. Therefore, we developed sleep fMRI based on simultaneous electrophysiology at 9.4 T in male mice. Optimized un-anesthetized mouse fMRI setup allowed manifestation of NREM and REM sleep, and a large sleep fMRI dataset was collected and openly accessible. State dependent global patterns were revealed, and state transitions were found to be global, asymmetrical and sequential, which can be predicted up to 17.8 s using LSTM models. Importantly, sleep fMRI with hippocampal recording revealed potentiated sharp-wave ripple triggered global patterns during NREM than awake state, potentially attributable to co-occurrence of spindle events. To conclude, we established mouse sleep fMRI with simultaneous electrophysiology, and demonstrated its capability by revealing global dynamics of state transitions and neural events.

Fear extinction is a form of inhibitory learning that suppresses the expression of aversive memories and plays a key role in the recovery of anxiety and trauma-related disorders. Here, using male mice, we identify a cerebello-thalamo-cortical pathway regulating fear extinction. The cerebellar fastigial nucleus (FN) projects to the lateral subregion of the mediodorsal thalamic nucleus (MD), which is reciprocally connected with the dorsomedial prefrontal cortex (dmPFC). The inhibition of FN inputs to MD in male mice impairs fear extinction in animals with high fear responses and increases the bursting of MD neurons, a firing pattern known to prevent extinction learning. Indeed, this MD bursting is followed by high levels of the dmPFC 4 Hz oscillations causally associated with fear responses during fear extinction, and the inhibition of FN-MD neurons increases the coherence of MD bursts and oscillations with dmPFC 4 Hz oscillations. Overall, these findings reveal a regulation of fear-related thalamo-cortical dynamics by the cerebellum and its contribution to fear extinction.

In humans, the level of consciousness is assessed by quantifying the spatiotemporal complexity of cortical responses using Perturbational Complexity Index (PCI) and related PCIst (st, state transitions). Here we validate PCIst in freely moving rats and mice by showing that it is lower in NREM sleep and slow wave anesthesia than in wake or REM sleep, as in humans. We then show that (1) low PCIst is associated with the occurrence of an OFF period of neuronal silence; (2) stimulation of deep, but not superficial, cortical layers leads to reliable PCIst changes across sleep/wake and anesthesia; (3) consistent PCIst changes are independent of which single area is being stimulated or recorded, except for recordings in mouse prefrontal cortex. These experiments show that PCIst can reliably measure vigilance states in unresponsive animals and support the hypothesis that it is low when an OFF period disrupts causal interactions in cortical networks.

Background During non-rapid eye movement sleep (NREM), alternating periods of synchronised high (ON period) and low (OFF period) neuronal activity are associated with high amplitude delta band (0.5–4 Hz) oscillations in neocortical electrophysiological signals termed slow waves. As this oscillation is dependent crucially on hyperpolarisation of cortical cells, there is an interest in understanding how neuronal silencing during OFF periods leads to the generation of slow waves and whether this relationship changes between cortical layers. A formal, widely adopted definition of OFF periods is absent, complicating their detection. Here, we grouped segments of high frequency neural activity containing spikes, recorded as multiunit activity from the neocortex of freely behaving mice, on the basis of amplitude and asked whether the population of low amplitude (LA) segments displayed the expected characteristics of OFF periods. Results Average LA segment length was comparable to previous reports for OFF periods but varied considerably, from as short as 8 ms to > 1 s. LA segments were longer and occurred more frequently in NREM but shorter LA segments also occurred in half of rapid eye movement sleep (REM) epochs and occasionally during wakefulness. LA segments in all states were associated with a local field potential (LFP) slow wave that increased in amplitude with LA segment duration. We found that LA segments > 50 ms displayed a homeostatic rebound in incidence following sleep deprivation whereas short LA segments (< 50 ms) did not. The temporal organisation of LA segments was more coherent between channels located at a similar cortical depth. Conclusion We corroborate previous studies showing neural activity signals contain uniquely identifiable periods of low amplitude with distinct characteristics from the surrounding signal known as OFF periods and attribute the new characteristics of vigilance-state-dependent duration and duration-dependent homeostatic response to this phenomenon. This suggests that ON/OFF periods are currently underdefined and that their appearance is less binary than previously considered, instead representing a continuum.

The common marmoset has emerged as a key model in neuroscience. Marmosets are small in size, show great potential for genetic modification and exhibit complex behaviors. Thus, it is necessary to develop technology that enables monitoring and manipulation of the underlying neural circuits. Here, we describe a novel approach to record and optogenetically manipulate neural activity in awake, behaving marmosets. Our design utilizes a light-weight, 3D printed titanium chamber that can house several high-density silicon probes for semi-chronic recordings, while enabling simultaneous optogenetic stimulation. We demonstrate the application of our method in male marmosets by recording multi- and single-unit data from areas V1 and V6 with 192 channels simultaneously, and show that optogenetic activation of excitatory neurons in area V6 can influence behavior in a detection task. This method may enable future studies to investigate the neural basis of perception and behavior in the marmoset.

Episodic Ataxia Type 2 (EA2) is a rare neurological disorder caused by a mutation in the CACNA1A gene, encoding the P/Q-type voltage-gated Ca 2+ channel important for neurotransmitter release. Patients with this channelopathy exhibit both cerebellar and cerebral pathologies, suggesting the condition affects both regions. The tottering ( tg / tg) mouse is the most commonly used EA2 model due to an orthologous mutation in the cacna1a gene. The tg / tg mouse has three prominent behavioral phenotypes: a dramatic episodic dystonia; absence seizures with generalized spike and wave discharges (GSWDs); and mild ataxia. We previously observed a novel brain state, transient low-frequency oscillations (LFOs) in the cerebellum and cerebral cortex under anesthesia. In this study, we examine the relationships among the dystonic attack, GSWDs, and LFOs in the cerebral cortex. Previous studies characterized LFOs in the motor cortex of anesthetized tg / tg mice using flavoprotein autofluorescence imaging testing the hypothesis that LFOs provide a mechanism for the paroxysmal dystonia. We sought to obtain a more direct understanding of motor cortex (M1) activity during the dystonic episodes. Using two-photon Ca 2+ imaging to investigate neuronal activity in M1 before, during, and after the dystonic attack, we show that there is not a significant change in the activity of M1 neurons from baseline through the attack. We also conducted simultaneous, multi-electrode recordings to further understand how M1 cellular activity and local field potentials change throughout the progression of the dystonic attack. Neither putative pyramidal nor inhibitory interneuron firing rate changed during the dystonic attack. However, we did observe a near complete loss of GSWDs during the dystonic attack in M1. Finally, using spike triggered averaging to align simultaneously recorded limb kinematics to the peak Ca 2+ response, and vice versa , revealed a reduction in the spike triggered average during the dystonic episodes. Both the loss of GSWDs and the reduction in the coupling suggest that, during the dystonic attack, M1 is effectively decoupled from other structures. Overall, these results indicate that the attack is not initiated or controlled in M1, but elsewhere in the motor circuitry. The findings also highlight that LFOs, GSWDs, and dystonic attacks represent three brain states in tg / tg mice.

The motor cortex controls skilled arm movement by recruiting a variety of targets in the nervous system, and it is important to understand the emergent activity in these regions as refinement of a motor skill occurs. One fundamental projection of the motor cortex (M1) is to the cerebellum. However, the emergent activity in the motor cortex and the cerebellum that appears as a dexterous motor skill is consolidated is incompletely understood. Here, we report on low-frequency oscillatory (LFO) activity that emerges in cortico-cerebellar networks with learning the reach-to-grasp motor skill. We chronically recorded the motor and the cerebellar cortices in rats, which revealed the emergence of coordinated movement-related activity in the local-field potentials as the reaching skill consolidated. Interestingly, we found this emergent activity only in the rats that gained expertise in the task. We found that the local and cross-area spiking activity was coordinated with LFOs in proficient rats. Finally, we also found that these neural dynamics were more prominently expressed during accurate behavior in the M1. This work furthers our understanding on emergent dynamics in the cortico-cerebellar loop that underlie learning and execution of precise skilled movement.

Systems consolidation—a process for long-term memory stabilization—has been hypothesized to occur in two stages 1–4 . Whereas new memories require the hippocampus 5–9 , they become integrated into cortical networks over time 10–12 , making them independent of the hippocampus. How hippocampal–cortical dialogue precisely evolves during this and how cortical representations change in concert is unknown. Here, we use a skill learning task 13,14 to monitor the dynamics of cross-area coupling during non-rapid eye movement sleep along with changes in primary motor cortex (M1) representational stability. Our results indicate that precise cross-area coupling between hippocampus, prefrontal cortex and M1 can demarcate two distinct stages of processing. We specifically find that each animal demonstrates a sharp increase in prefrontal cortex and M1 sleep slow oscillation coupling with stabilization of performance. This sharp increase then predicts a drop in hippocampal sharp-wave ripple (SWR)–M1 slow oscillation coupling—suggesting feedback to inform hippocampal disengagement and transition to a second stage. Notably, the first stage shows significant increases in hippocampal SWR–M1 slow oscillation coupling in the post-training sleep and is closely associated with rapid learning and variability of the M1 low-dimensional manifold. Strikingly, even after consolidation, inducing new manifold exploration by changing task parameters re-engages hippocampal–M1 coupling. We thus find evidence for dynamic hippocampal–cortical dialogue associated with manifold exploration during learning and adaptation.

From myriads of ongoing stimuli, the brain creates a fused percept of the environment. This process, which culminates in perceptual binding, is presumed to occur through the operations of multisensory neurons that occur throughout the brain. However, because different brain areas receive different inputs and have different cytoarchitechtonics, it would be expected that local multisensory features would also vary across regions. The present study investigated that hypothesis using multiple single-unit recordings from anesthetized cats in response to controlled, electronically-generated separate and combined auditory, visual, and somatosensory stimulation. These results were used to compare the multisensory features of neurons in cat primary auditory cortex (A1) with those identified in the nearby higher-order auditory region, the Dorsal Zone (DZ). Both regions exhibited the same forms of multisensory neurons, albeit in different proportions. Multisensory neurons exhibiting excitatory or inhibitory properties occurred in similar proportions in both areas. Also, multisensory neurons in both areas expressed similar levels of multisensory integration. Because responses to auditory cues alone were so similar to those that included non-auditory stimuli, it is proposed that this effect represents a mechanism by which multisensory neurons subserve the process of perceptual binding.

Although sleep and anesthesia are predicted to share common neural signatures of transitions into and out of unconsciousness, supportive evidence has been elusive. We identified these signatures using intracranial electroencephalography in neurosurgical patients. We applied diffusion map embedding to map cortical location into a space where proximity indicates functional similarity using a normalized connectivity (‘diffusion’) matrix, itself a rich source of information about network properties. During reduced consciousness, diffusion matrices exhibited decreased effective dimensionality, reflecting reduced network entropy. Furthermore, functional brain regions exhibited tighter clustering in embedding space with greater distances between regions, corresponding to decreased differentiation and functional integration. These changes were not region-specific, suggesting global network reorganization. These results strongly suggest common neural substrates for loss and recovery of consciousness during anesthesia and sleep, providing a systems-level mechanistic understanding within an intuitive geometric context and laying the foundation for evaluation of cortical state transitions in clinical settings.