Electrophysiology publications

Sleep homeostasis manifests as a relative constancy of its daily amount and intensity. Theoretical descriptions define ‘Process S’, a variable with dynamics dependent on global sleep-wake history, and reflected in electroencephalogram (EEG) slow wave activity (SWA, 0.5–4 Hz) during sleep. The notion of sleep as a local, activity-dependent process suggests that activity history must be integrated to determine the dynamics of global Process S. Here, we developed novel mathematical models of Process S based on cortical activity recorded in freely behaving mice, describing local Process S as a function of the deviation of neuronal firing rates from a locally defined set-point, independent of global sleep-wake state. Averaging locally derived Processes S and their rate parameters yielded values resembling those obtained from EEG SWA and global vigilance states. We conclude that local Process S dynamics reflects neuronal activity integrated over time, and global Process S reflects local processes integrated over space.

The cerebellum is involved in the control of voluntary and autonomic rhythmic behaviors, yet it is unclear to what extent it coordinates these in concert. We studied Purkinje cell activity during unperturbed and perturbed respiration in lobules simplex, crus 1, and crus 2. During unperturbed (eupneic) respiration, complex spike and simple spike activity encode the phase of ongoing sensorimotor processing. In contrast, when the respiratory cycle is perturbed by whisker stimulation, mice concomitantly protract their whiskers and advance their inspiration in a phase-dependent manner, preceded by increased simple spike activity. This phase advancement of respiration in response to whisker stimulation can be mimicked by optogenetic stimulation of Purkinje cells and prevented by cell-specific genetic modification of their AMPA receptors, hampering increased simple spike firing. Thus, the impact of Purkinje cell activity on respiratory control is context and phase dependent, highlighting a coordinating role for the cerebellar hemispheres in aligning autonomic and sensorimotor behaviors.

Knowing the depth and laminar location of the microelectrodes used to record neuronal activity from the cerebral cortex is crucial to properly interpret the recorded patterns of neuronal responses. Here, we present an innovative approach that allows inferring such properties with high accuracy and in an automated way (i.e., without the need of visual inspection and manual annotation) from the evoked response potentials elicited by sensory (e.g., visual) stimuli.

The peripheral nervous system (PNS) is an attractive target for modulation of afferent input (e.g., nociceptive input signaling tissue damage) to the central nervous system. To advance mechanistic understanding of PNS neural encoding and modulation requires single-unit recordings from individual peripheral neurons or axons. This is challenged by multiple connective tissue layers surrounding peripheral nerve fibers that prevent electrical recordings by existing electrodes or electrode arrays. In this study, we developed a novel microelectrode array (MEA) via silicon-based microfabrication that consists of 5 parallel hydrophilic gold electrodes surrounded by silanized hydrophobic surfaces. This novel hydrophilic/hydrophobic surface pattern guides the peripheral nerve filaments to self-align towards the hydrophilic electrodes, which dramatically reduces the technical challenges in conducting single-unit recordings. We validated our MEA by recording simultaneous single-unit action potentials from individual axons in mouse sciatic nerves, including both myelinated A-fibers and unmyelinated C-fibers. We confirmed that our recordings were single units from individual axons by increasing nerve trunk electrical stimulus intensity, which did not alter the spike shape or amplitude. By reducing the technical challenges, our novel MEA will likely allow peripheral single-unit recordings to be adopted by a larger research community and thus expedite our mechanistic understanding of peripheral neural encoding and modulation.

Basal forebrain modulation of central circuits is associated with active sensation, attention, and learning. While cholinergic modulations have been studied extensively the effect of non-cholinergic basal forebrain subpopulations on sensory processing remains largely unclear. Here, we directly compare optogenetic manipulation effects of two major basal forebrain subpopulations on principal neuron activity in an early sensory processing area, i.e. mitral/tufted cells (MTCs) in the olfactory bulb. In contrast to cholinergic projections, which consistently increased MTC firing, activation of GABAergic fibers from basal forebrain to the olfactory bulb leads to differential modulation effects: while spontaneous MTC activity is mainly inhibited, odor-evoked firing is predominantly enhanced. Moreover, sniff-triggered averages revealed an enhancement of maximal sniff evoked firing amplitude and an inhibition of firing rates outside the maximal sniff phase. These findings demonstrate that GABAergic neuromodulation affects MTC firing in a bimodal, sensory-input dependent way, suggesting that GABAergic basal forebrain modulation could be an important factor in attention mediated filtering of sensory information to the brain.

Proprioceptive feedback is a critical component of voluntary movement planning and execution. Neuroprosthetic technologies aiming at restoring movement must interact with it to restore accurate motor control. Optimization and design of such technologies depends on the availability of quantitative insights into the neural dynamics of proprioceptive afferents during functional movements. However, recording proprioceptive neural activity during unconstrained movements in clinically relevant animal models presents formidable challenges. In this work, we developed a computational framework to estimate the spatiotemporal patterns of proprioceptive inputs to the cervical spinal cord during three-dimensional arm movements in monkeys. We extended a biomechanical model of the monkey arm with ex-vivo measurements, and combined it with models of mammalian group-Ia, Ib and II afferent fibers. We then used experimental recordings of arm kinematics and muscle activity of two monkeys performing a reaching and grasping task to estimate muscle stretches and forces with computational biomechanics. Finally, we projected the simulated proprioceptive firing rates onto the cervical spinal roots, thus obtaining spatiotemporal maps of spinal proprioceptive inputs during voluntary movements. Estimated maps show complex and markedly distinct patterns of neural activity for each of the fiber populations spanning the spinal cord rostro-caudally. Our results indicate that reproducing the proprioceptive information flow to the cervical spinal cord requires complex spatio-temporal modulation of each spinal root. Our model can support the design of neuroprosthetic technologies as well as in-silico investigations of the primate sensorimotor system.

The ventral tegmental area (VTA) has dopamine, GABA, and glutamate neurons, which have been implicated in reward and aversion. Here, we determined whether VTA-glutamate or -GABA neurons play a role in innate defensive behavior. By VTA cell-type-specific genetic ablation, we found that ablation of glutamate, but not GABA, neurons abolishes escape behavior in response to threatening stimuli. We found that escape behavior is also decreased by chemogenetic inhibition of VTA-glutamate neurons and detected increases in activity in VTA-glutamate neurons in response to the threatening stimuli. By ultrastructural and electrophysiological analysis, we established that VTA-glutamate neurons receive a major monosynaptic glutamatergic input from the lateral hypothalamic area (LHA) and found that photoinhibition of this input decreases escape responses to threatening stimuli. These findings indicate that VTA-glutamate neurons are activated by and required for innate defensive responses and that information on threatening stimuli to VTA-glutamate neurons is relayed by LHA-glutamate neurons.

Odors are well known to elicit strong emotional and behavioral responses that become strengthened throughout learning, yet the specific cellular systems involved in odor learning and the direct influence of these on behavior are unclear. Here, we investigate the representation of odor-reward associations within two areas recipient of dense olfactory input, the posterior piriform cortex (pPCX) and the olfactory tubercle (OT), using electrophysiological recordings from mice engaged in reward-based learning. Neurons in both regions represent conditioned odors and do so with similar information content, yet the proportion of neurons recruited by conditioned rewarded odors and the magnitudes and durations of their responses are greater in the OT. Using fiber photometry, we find that OT D1-type dopamine-receptor-expressing neurons flexibly represent odors based on reward associations, and using optogenetics, we show that these neurons influence behavioral engagement. These findings contribute to a model whereby OT D1 neurons support odor-guided motivated behaviors.

Traumatic stress can lead to heightened aggression which may be a symptom of psychiatric diseases such as PTSD and intermittent explosive disorder. The medial amygdala (MeA) is an evolutionarily conserved subnucleus of the amygdala that regulates attack behavior and behavioral responses to stressors. The precise contribution of the MeA in traumatic stress-induced aggression, however, requires further elucidation. In this study, we used foot shock to induce traumatic stress in mice and examine the mechanisms of prolonged aggression increase associated with it. Foot shock causes a prolonged increase in aggression that lasts at least one week. In vivo electrophysiological recordings revealed that foot shock induces potentiation of synapses formed between the MeA and the ventromedial hypothalamus (VmH) and bed nucleus of the stria terminalis (BNST). This synaptic potentiation lasts at least one week. Induction of synaptic depotentiation with low-frequency photostimulation (LFPS) immediately after foot shock suppresses the prolonged aggression increase without affecting non-aggressive social behavior, anxiety-like and depression-like behaviors, or fear learning. These results show that potentiation of the MeA-VmH and MeA-BNST circuits is essential for traumatic stress to cause a prolonged increase in aggression. These circuits may be potential targets for the development of therapeutic strategies to treat the aggression symptom associated with psychiatric diseases.

In vivo loose patch and breakthrough whole-cell recordings are useful tools for investigating the intrinsic and synaptic properties of neurons. However, the correlation among pipette resistance, seal condition, and recording time is not thoroughly clear. Presently, we investigated the recording time of different pipette resistances and seal conditions in loose patch and breakthrough whole-cell recordings. The recording time did not change with pipette resistance for loose patch recording (Rp-loose) and first increased and then decreased as seal resistance for loose patch recording (Rs-loose) increased. For a high probability of a recording time ≥30 min, the low and high cutoff values of Rs-loose were 21.5 and 36 MΩ, respectively. For neurons with Rs-loose values of 21.5-36 MΩ, the action potential (AP) amplitudes changed slightly 30 min after the seal. The recording time increased as seal resistance for whole-cell recording (Rs-tight) increased and the zero-current membrane potential for breakthrough whole-cell recording (MPzero-current) decreased. For a high probability of a recording time ≥30 min, the cutoff values of Rs-tight and MPzero-current were 2.35 GΩ and -53.5 mV, respectively. The area under the curve (AUC) of the MPzero-current receiver operating characteristic (ROC) curve was larger than that of the Rs-tight ROC curve. For neurons with MPzero-current values ≤ -53.5 mV, the inhibitory or excitatory postsynaptic current amplitudes did not show significant changes 30 min after the seal. In neurons with Rs-tight values ≥2.35 GΩ, the recording time gradually increased and then decreased as the pipette resistance for whole-cell recording (Rp-tight) increased. For the high probability of a recording time ≥30 min, the low and high cutoff values of Rp-tight were 6.15 and 6.45 MΩ, respectively. Together, we concluded that the optimal Rs-loose range is 21.5-36 MΩ, the optimal Rp-tight range is 6.15-6.45 MΩ, and the optimal Rs-tight and MPzero-current values are ≥2.35 GΩ and ≤ -53.5 mV, respectively. Compared with Rs-tight, the MPzero-current value can more accurately discriminate recording times ≥30 min and <30 min.

Our understanding of nonlinear stimulus transformations by neural circuits is hindered by the lack of comprehensive yet interpretable computational modeling frameworks. Here, we propose a data-driven approach based on deep neural networks to directly model arbitrarily nonlinear stimulus-response mappings. Reformulating the exact function of a trained neural network as a collection of stimulus-dependent linear functions enables a locally linear receptive field interpretation of the neural network. Predicting the neural responses recorded invasively from the auditory cortex of neurosurgical patients as they listened to speech, this approach significantly improves the prediction accuracy of auditory cortical responses, particularly in nonprimary areas. Moreover, interpreting the functions learned by neural networks uncovered three distinct types of nonlinear transformations of speech that varied considerably from primary to nonprimary auditory regions. The ability of this framework to capture arbitrary stimulus-response mappings while maintaining model interpretability leads to a better understanding of cortical processing of sensory signals.

The mechanical mismatch between implantable interfaces and neural tissues may be reduced by employing soft polymeric materials. Here, we report on a simple strategy to prepare and pattern a soft electrode coating of neural interfacing devices based on a screen-printable conducting hydrogel. The coating formulation, based on polyacrylamide and poly(3,4-ethylenedioxythiophene) polystyrene sulfonate, is suitable to additive manufacturing and exhibits excellent adhesion to polydimethylsiloxane, an elastomer commonly used as a substrate in soft neural interfaces. The soft conductive coating displays a tunable elastic modulus in the 10–100 kPa range and electrochemical properties on a par with stiff conductive inks while supporting good neural cell attachment and proliferation in vitro. Next, the soft printable hydrogel is integrated within a 4 × 4 microelectrode array for electrocorticography with 250 μm-diameter contacts. Acute recording of cortical local field potentials and electrochemical characterization preimplantation and postimplantation highlight the stability of the soft organic conductor. The overall robustness of the soft coating and its patterning method provide a promising route for a range of implantable bioelectronic applications.

Repeated stimulus presentations commonly produce decreased neural responses—a phenomenon known as repetition suppression (RS) or adaptation—in ventral temporal cortex (VTC) of humans and nonhuman primates. However, the temporal features of RS in human VTC are not well understood. To fill this gap in knowledge, we utilized the precise spatial localization and high temporal resolution of electrocorticography (ECoG) from nine human subjects implanted with intracranial electrodes in the VTC. The subjects viewed nonrepeated and repeated images of faces with long-lagged intervals and many intervening stimuli between repeats. We report three main findings: 1) robust RS occurs in VTC for activity in high-frequency broadband (HFB), but not lower-frequency bands; 2) RS of the HFB signal is associated with lower peak magnitude (PM), lower total responses, and earlier peak responses; and 3) RS effects occur early within initial stages of stimulus processing and persist for the entire stimulus duration. We discuss these findings in the context of early and late components of visual perception, as well as theoretical models of repetition suppression.

Understanding the mechanisms by which neuronal population activity gives rise to perception and behavior remains a central question in systems neuroscience. Such understanding is complicated by the fact that natural stimuli often have complex structure. Here we investigated how heterogeneities within a sensory neuron population influence the coding of a noisy stimulus waveform (i.e., the noise) and its behaviorally relevant envelope signal (i.e., the signal). We found that On- and Off-type neurons displayed more heterogeneities in their responses to the noise than in their responses to the signal. These differences in heterogeneities had important consequences when quantifying response similarity between pairs of neurons. Indeed, the larger response heterogeneity displayed by On- and Off-type neurons made their pairwise responses to the noise on average more independent than when instead considering pairs of On-type or Off-type neurons. Such relative independence allowed for better averaging out of the noise response when pooling neural activities in a mixed-type (i.e., On- and Off-type) than for same-type (i.e., only On-type or only Off-type), thereby leading to greater information transmission about the signal. Our results thus reveal a function for the combined activities of On- and Off-type neurons towards improving information transmission of envelope stimuli at the population level. Our results will likely generalize because natural stimuli across modalities are characterized by a stimulus waveform whose envelope varies independently as well as because On- and Off-type neurons are observed across systems and species.

Dopamine guides behavior and learning through pleasure, according to classic understanding. Dopaminergic neurons are traditionally thought to signal positive or negative prediction errors (PEs) when reward expectations are, respectively, exceeded or not matched. These signed PEs are quite different from the unsigned PEs, which report surprise during sensory processing. But mounting theoretical accounts from the predictive processing framework postulate that dopamine, as a neuromodulator, could potentially regulate the postsynaptic gain of sensory neurons, thereby scaling unsigned PEs according to their expected precision or confidence. Despite ample modeling work, the physiological effects of dopamine on the processing of surprising sensory information are yet to be addressed experimentally. In this study, we tested how dopamine modulates midbrain processing of unexpected tones. We recorded extracellular responses from the rat inferior colliculus to oddball and cascade sequences, before, during, and after the microiontophoretic application of dopamine or eticlopride (a D2-like receptor antagonist). Results demonstrate that dopamine reduces the net neuronal responsiveness exclusively to unexpected sensory input without significantly altering the processing of expected input. We conclude that dopaminergic projections from the thalamic subparafascicular nucleus to the inferior colliculus could encode the expected precision of unsigned PEs, attenuating via D2-like receptors the postsynaptic gain of sensory inputs forwarded by the auditory midbrain neurons. This direct dopaminergic modulation of sensory PE signaling has profound implications for both the predictive coding framework and the understanding of dopamine function.

Recognizing reward-related stimuli is crucial for survival. Neuronal projections from the basolateral amygdala (BLA) to the nucleus accumbens (NAc) play an important role in processing reward-related cues. Previous studies revealed synchronization between distant brain regions in reward-sensitive neurocircuits; however, whether the NAc synchronizes with the BLA is unknown. Here, we recorded local field potentials simultaneously from the BLA and NAc of rats during social preference tests and an appetitive conditioning test in which explicit stimuli were associated with food. BLA-NAc coherence in the theta band (5-8 Hz) increased in response to food-associated cues. Meanwhile, the modulatory strength of theta-high gamma (50-110 Hz) phase-amplitude cross-frequency coupling (PAC) in the NAc decreased. Importantly, both of these neuromodulations disappeared upon extinction. In contrast, both theta and gamma power oscillations in each region increased in the presence of social conspecifics or contexts associated with conspecifics, but coherence did not change. To potentially disrupt behavior and associated neural activity, a subgroup of rats was exposed prenatally to valproic acid (VPA), which has been shown to disrupt transcriptome and excitatory/inhibitory balance in the amygdala. VPA-exposed rats demonstrated impulsive-like behavior, but VPA did not affect BLA-NAc coherence. These findings reveal changes in BLA-NAc coherence in response to select reward-related stimuli (i.e., food-predictive cues); the differences between the tasks used here could shed light onto the functional nature of BLA-NAc coherence and are discussed.

Listeners differ in their ability to attend to a speech stream in the presence of a competing sound. Differences in speech intelligibility in noise cannot be fully explained by the hearing ability which suggests the involvement of additional cognitive factors. A better understanding of the temporal fluctuations in the ability to pay selective auditory attention to a desired speech stream may help in explaining these variabilities. In order to better understand the temporal dynamics of selective auditory attention, we developed an online auditory attention decoding (AAD) processing pipeline based on speech envelope tracking in the electroencephalogram (EEG). Participants had to attend to one audiobook story while a second one had to be ignored. Online AAD was applied to track the attention toward the target speech signal. Individual temporal attention profiles were computed by combining an established AAD method with an adaptive staircase procedure. The individual decoding performance over time was analyzed and linked to behavioral performance as well as subjective ratings of listening effort, motivation, and fatigue. The grand average attended speaker decoding profile derived in the online experiment indicated performance above chance level. Parameters describing the individual AAD performance in each testing block indicated significant differences in decoding performance over time to be closely related to the behavioral performance in the selective listening task. Further, an exploratory analysis indicated that subjects with poor decoding performance reported higher listening effort and fatigue compared to good performers. Taken together our results show that online EEG based AAD in a complex listening situation is feasible. Adaptive attended speaker decoding profiles over time could be used as an objective measure of behavioral performance and listening effort. The developed online processing pipeline could also serve as a basis for future EEG based near real-time auditory neurofeedback systems.

Cortical injuries (e.g. – strokes or traumatic brain injuries) can create a host of secondary events that further impair the brain’s sensory, motor, or cognitive capabilities. Here, we attempted to isolate the acute effects of the primary injury – the loss of cortical activity – on rodent motor cortex (caudal forelimb area, CFA) without the secondary effects that arise from damage to cortical tissue. We then observed the effects of this loss of activity on the rodent premotor cortex (rostral forelimb area, RFA). In anesthetized rats, CFA was temporarily inactivated with the GABA-A agonist muscimol, disrupting motor network function while leaving neural connectivity intact. Using intracortical microstimulation (ICMS) techniques, we found that CFA inactivation completely abolished ICMS-evoked forelimb movement from RFA yet spared some CFA evoked-movement. Neural recordings confirmed that neural suppression by muscimol was isolated to CFA and did not spread into RFA. We next observed how CFA inactivation suppressed RFA influence on forelimb muscles by obtaining intramuscular electromyographical (EMG) recordings from forelimb muscles during ICMS. EMG recordings showed that despite the presence of evoked movement in CFA, but not RFA, muscle activation in both areas were similarly reduced. These results suggest that the primary reason for the loss of ICMS-evoked movement in RFA is not reduced forelimb muscle activity, but rather a loss of the specific activity between RFA and CFA. Therefore, within the intact motor network of the rat, RFA’s influence on forelimb movement is mediated by CFA, similar to the premotor and motor organization observed in non-human primates.

Brain injury affecting the isocortical frontal cortex is a common pathological occurrence. Many patients report severe deficits to functions of daily living. However, there is a variable degree of motor recovery that occurs with some individuals recovering astounding degrees of motor recovery while others have not. This variability has led researchers into investigating the possible mechanisms for this variability. Recently, several non-human primate studies have shed light on the possibility of spared, ipsilesional motor area taken over the lost function to the damaged cortex. Unfortunately, these studies have focused on long-term adaption ranging from 5months to one year post injury. In this present study, we are the first use rigorous stereological quantification to show that significant neuroplastic changes in the form of changes to neuroanatomical connections between distant cortical area occurs at a very early time point of 4 weeks post injury. Much like the Dancause study in 2005, we found that ishemic damage to the distal forelimb area (DFL) of the primary motor cortex (M1) induced plastic changes between the DFL of the ventral premotor cortex (PMv) and area 1/2 of the somatosensory cortex. Indeed, we found a nearly 2 fold increase in the number of boutons between PMV and area 1/2. Additionally, labeled fibers from PMv change direction from their normal termination within M1 and traveled in a ventral posterior direction toward the somatosensory cortex. Also of interest, several labeled fibers actually traveled through the glial scar of M1 toward the somatosensory cortex. These data demonstrate that a massive neuroplastic response has occurred following an ischemic insult to the DFL of M1. These data may suggest that the brain may be undergoing an attempt to re-establish a degree of motor and or sensory control to compensate for the lost function due to the injury.

Key points Epidural electrical stimulation (ES) of the spinal cord restores/improves locomotion in patients. ES‐evoked locomotor movements differ to some extent from the normal ones. Operation of the locomotor network during ES is unknown. We compared the activity of individual spinal neurons during locomotion initiated by signals from the brainstem and by ES. We demonstrated that the spinal network generating locomotion under each of the two conditions is formed by the same neurons. A part of this network operates similarly under the two conditions, suggesting that it is essential for generation of locomotion under both conditions. Another part of this network operates differently under the two conditions, suggesting that it is responsible for differences in the movement kinematics observed under the two conditions. Abstract Locomotion is a vital motor function for both animals and humans. Epidural electrical stimulation (ES) of the spinal cord is used to restore/improve locomotor movements in patients. However, operation of locomotor networks during ES has never been studied. Here we compared the activity of individual spinal neurons recorded in decerebrate cats of either sex during locomotion initiated by supraspinal commands (caused by stimulation of the mesencephalic locomotor region, MLR) and by ES. We found that under both conditions, the same neurons had modulation of their activity related to the locomotor rhythm, suggesting that the network generating locomotion under the two conditions is formed by the same neurons. About 40% of these neurons had stable modulation (i.e. small dispersion of their activity phase in sequential cycles), as well as a similar phase and shape of activity burst in MLR‐ and ES‐evoked locomotor cycles. We suggest that these neurons form a part of the locomotor network that operates similarly under the two conditions, and are critical for generation of locomotion. About 23% of the modulated neurons had stable modulation only during MLR‐evoked locomotion. We suggest that these neurons are responsible for some differences in kinematics of MLR‐ and ES‐evoked locomotor movements. Finally, 25% of the modulated neurons had unstable modulation during both MLR‐ and ES‐evoked locomotion. One can assume that these neurons contribute to maintenance of the excitability level of locomotor networks necessary for generation of stepping, or belong to postural networks, activated simultaneously with locomotor networks by both MLR stimulation and ES.