EEG/ECoG/EMG publications

Brain function relies on coordinated activity across spatial and temporal scales. Single neurons integrate local and long-range connectivity and reflect activity across brain-wide networks. Understanding integrated brain function requires tools capable of recording from anatomically connected populations in distributed brain areas to bridge local and global dynamics. Here, we present high-density, micro-electrocorticography arrays that facilitate multi-scale studies of brain activity. The hybrid arrays integrate the desirable features of silicone elastomers and polyimide films: silicone provides optical transparency and permits repeated penetration with intracortical arrays, while polyimide enables fine photolithographic feature definition. This combination facilitates high-throughput functional mapping to identify targets and insertion of intracortical arrays for dense local sampling. We demonstrate functional mapping in rats, cats, and marmosets, showing how functional maps guide multi-area laminar recordings. Finally, we demonstrate local and feedforward optogenetic stimulation to investigate cortico-cortical interactions. These capabilities establish the hybrid μECoG as a compelling tool for systems neuroscience.

Purpose: In anesthesia research and clinical practice, the level of consciousness is often assessed by stimulating subjects and observing their responses. It is possible that such stimulation alters the very state of consciousness being assessed. Whether this phenomenon is accompanied by specific cortical activities remains unknown. Here, we tested the effects of behavioral assessments on neuronal activity in the rat primary auditory cortex (A1) and perirhinal cortex (PrC). Methods: Neuronal spiking activity was recorded from microwire arrays targeting A1 and/or PrC in rats (n = 9), whilst presenting 40 Hz click-trains and gradually increasing infusion rates of propofol, remifentanil, ketamine or dexmedetomidine. We assessed righting reflex, pain withdrawal, whisking reflex, and corneal reflex every 12 minutes (every 360 auditory trials). We searched for significant changes in baseline firing rate and stimulus-evoked activities (where present) due to the behavioral assessments, by comparing 4 min intervals before and after the behavioral tests. Results: Out of 1004 units in 20 sessions, we found significantly elevated activity (p < 0.01) following behavioral tests in 32% and 20% of A1 auditory-responsive units and 8.6% and 5.3% of auditory-unresponsive units under propofol and dexmedetomidine, respectively, but not under remifentanil or ketamine. A higher proportion of units modulated their firing following behavioral tests in A1 compared to PrC, with median gains ranging from +0.23 to +1.06, depending on the specific response feature. Conclusion: Behavioral testing increases ongoing firing rate and responses to sounds in some cortical neurons. This increase is more prevalent in A1 than in PrC, mirroring clinical observations that patients exhibit basic responses under sedation, while complex processing, such as memory encoding, is impaired. The degree of neuronal modulation due to behavioral testing varies between anesthetic agents, mirroring their clinical phenotypes.

Pain management strategies have progressed beyond traditional pharmacologic and physical interventions, integrating advanced neuromodulation techniques such as deep brain stimulation, peripheral nerve stimulation, and high‐frequency spinal cord stimulation (hSCS). Despite its clinical efficacy, the supraspinal mechanisms underlying hSCS remain poorly understood. Prior work in sheep demonstrated that hSCS modulates gamma (γ) band (70–150 Hz) activity in the primary somatosensory and association cortices, implicating cortical involvement in pain modulation. Given the interaction between low‐ and high‐oscillations, we hypothesized that hSCS modulates γ activity in a region‐ and time‐dependent manner through specific coupling with theta rhythms (θ, 4–8 Hz). To test this, we quantified θ–γ phase‐amplitude coupling (PAC) during tibial nerve stimulation (TNS) before (TNS 0 ) and after (TNS 1 ) hSCS to assess cortical processing of TNS‐evoked activity, while recording cortical (association and somatosensory cortices) neural activity using 96‐channel subdural electrocorticography (ECoG). We then computed the corresponding modulation index (MI) to quantify PAC strength. While the preferred θ phase of γ activity remained consistent across conditions in both cortical areas, θ‐γ coupling was more robust and stable in the association than the somatosensory cortex. However, MI increased significantly post‐hSCS in both cortices, indicating enhanced θ‐γ coupling during TNS 1 . Temporally, both cortices showed distinct response patterns: the somatosensory cortex exhibited strongest coupling increase early post‐hSCS, followed by a progressive decline, whereas the association cortex showed a more delayed and temporally distributed enhancement. These findings indicate that hSCS modulates cortical cross‐frequency interactions during TNS‐evoked sensory processing, wherein somatosensory and association cortices show distinct temporal coupling dynamics that may contribute to supraspinal mechanisms of pain modulation.

The cortex generates diverse neural dynamics, ranging from broadband fluctuations to narrowband oscillations at specific frequencies. Here, we investigated whether broadband and oscillatory dynamics play different roles in the encoding and transmission of visual information. We used information-theoretical measures to dissociate neural signals sharing common information (i.e., redundancy) from signals encoding complementary information (i.e., synergy). We analyzed electrocorticography (ECoG) and local field potentials (LFP) in the visual cortex of human and non-human primates (macaque) to investigate the extent to which broadband signals (BB) and narrowband gamma (NBG) oscillations conveyed synergistic or redundant information about images. In both species, the information conveyed by BB signals was highly synergistic within and between visual areas. By contrast, the information carried by NBG was primarily redundant within and between the same visual areas. Finally, the information conveyed by BB signals emerged early after stimulus onset, while NBG sustained information at later time points. These results suggest a potential dual role of BB and NBG cortical dynamics in visual processing, with broadband dynamics supporting nonlinear pattern recognition and oscillations facilitating information maintenance across the cortex.

Speech production entails several processing steps that encode linguistic and articulatory structure, but whether these computations correspond to spatiotemporally discrete patterns of neural activity is unclear. To address this issue, we use electrocorticography to directly measure the brains of neurosurgical participants performing an interactive speech paradigm. We observe a broad range of cortical modulation profiles, and subsequent clustering analyses establish that responses comprised distinct classes associated with sensory perception, planning, motor execution, and task-related suppression. These activity classes are also localized to separate neural substrates, indicating their status as specialized networks. We then parse dynamics in the planning and motor networks using unsupervised dimensionality reduction, which reveals subnetworks that are sequentially active throughout preparation and articulation. These results therefore support and extend a localizationist model of speech production where cortical activity "flows" within and across discrete pathways during language use.

The non-rapid-eye movement (NREM) sleep power spectrum is composed of rhythmic and arrhythmic components, respectively associated with brain rhythms and scale-free dynamics. Both components are hypothesized to represent distinct processes underlying sleep-dependent memory consolidation as well as other brain complex processes. Recent advancements in spectral parametrization techniques and the use of multifractal models have enabled new insights into these components and their connection to brain networks. Aging impacts NREM sleep oscillations, but no animal studies considered the impact of age on rhythmic and arrhythmic components. In this study, we assessed the effects of age on the power spectrum and its two components using local field potential recordings in mice. We recorded across the cerebral cortex and within the hippocampus—a central brain hub involved in NREM sleep-dependent cognitive processing. Ten younger (7.6 months) and eleven older (15.7 months) C57BL/6 J male and female mice were continuously recorded over a 24-hour period. We extracted the NREM sleep standard and rhythmic spectra, controlled for scale-free activity, and estimated multifractal arrhythmic properties during NREM sleep using the Wavelet Leader and Bootstrap based MultiFractal analysis toolbox. Older mice showed specific alterations in both components compared to younger animals: reduced rhythmic gamma power in the anterior cortex, greater regional differentiation of scaling exponents, and higher multifractal dispersion of the arrhythmic component in the hippocampus. These findings demonstrate that aging alters rhythmic and arrhythmic properties of NREM sleep in distinct brain regions, suggesting that both contribute to age-related changes in sleep-dependent cognition.

The locus coeruleus–norepinephrine (LC-NE) system regulates arousal and awakening; however, it remains unclear whether the LC does this in a global or circuit-specific manner. We hypothesized that sensory-evoked awakenings are predominantly regulated by specific LC-NE efferent pathways. Anatomical, physiological, and functional modularities of LC-NE pathways involving the mouse basal forebrain (BF) and pontine reticular nucleus (PRN) were tested. We found partial anatomical segregation between the LC → PRN and LC → BF circuits. Extracellular NE dynamics in BF and PRN exhibited distinct sound-evoked activation during sleep, including a fast sound-evoked NE peak specific to PRN. Causal optogenetic interrogation of LC efferent pathways, by retro-channelrhodopsin (ChR2) activation or Platynereis dumerilii ciliary opsin (PdCO) silencing of synapses in target regions, revealed a role for early LC → PRN activity in driving arousal and sound-evoked awakenings. Together, our results uncover a role for early LC-NE PRN activity in connecting sensory and arousal pathways and establish LC heterogeneity in regulating arousal.

Mammals can execute intended limb movements despite the fact that spinal reflexes involuntarily modulate muscle activity. To generate appropriate muscle activity, the cortical descending motor output must coordinate with spinal reflexes, yet the underlying neural mechanism remains unclear. We simultaneously recorded activities in motor-related cortical areas, afferent neurons, and forelimb muscles of monkeys performing reaching movements. Motor-related cortical areas, predominantly primary motor cortex (M1), encode subsequent afferent activities attributed to forelimb movement. M1 also encodes a subcomponent of muscle activity evoked by these afferent activities, corresponding to spinal reflexes. Furthermore, selective disruption of the afferent pathway specifically reduced this subcomponent of muscle activity, suggesting that M1 output drives muscle activity not only through direct descending pathways but also through the “transafferent” pathway composed of descending plus subsequent spinal reflex pathways. Thus, M1 provides optimal motor output based on an internal forward model that prospectively computes future spinal reflexes.

Animacy perception, the ability to discern living from non-living entities, is crucial for survival and social interaction, as it includes recognizing abstract concepts such as movement, purpose, and intentions. This process involves interpreting cues that may suggest the intentions or actions of others. It engages the temporal cortex (TC), particularly the superior temporal sulcus (STS) and the adjacent region of the inferior temporal cortex (ITC), as well as the dorsomedial prefrontal cortex (dmPFC). However, it remains unclear how animacy is dynamically encoded over time in these brain areas and whether its processing is distributed or localized. In this study, we addressed these questions by employing a symbolic categorization task involving animate and inanimate objects using natural movie stimuli. Simultaneously, electrocorticography were conducted in both the TC and dmPFC. Time-frequency analysis revealed region-specific frequency representations throughout the observation of the movies. Spatial searchlight decoding analysis demonstrated that animacy processing is represented in a distributed manner. Regions encoding animacy information were found to be dispersed across the fundus and lip of the STS, as well as in the ITC. Next, we examined whether these dispersed regions form functional networks. Independent component analysis revealed that the spatial distribution of the component with the most significant animacy information corresponded with the dispersed regions identified by the spatial decoding analysis. Furthermore, Granger causality analysis indicated that these regions exhibit frequency-specific directional functional connectivity, with a general trend of causal influence from the ITC to STS across multiple frequency bands. Notably, a prominent feedback flow in the alpha band from the ITC to both the ventral bank and fundus of the STS was identified. These findings suggest a distributed and functionally interconnected neural substrate for animacy processing across the STS and ITC.

The claustrum has been linked to attention and sleep. We hypothesized that this reflects a shared function, determining responsiveness to stimuli, which spans the axis of engagement. To test this hypothesis, we recorded claustrum population dynamics from male mice during both sleep and an attentional task (‘ENGAGE’). Heightened activity in claustrum neurons projecting to the anterior cingulate cortex (ACCp) corresponded to reduced sensory responsiveness during sleep. Similarly, in the ENGAGE task, heightened ACCp activity correlated with disengagement and behavioral lapses, while low ACCp activity correlated with hyper-engagement and impulsive errors. Chemogenetic elevation of ACCp activity reduced both awakenings during sleep and impulsive errors in the ENGAGE task. Furthermore, mice employing an exploration strategy in the task showed a stronger correlation between ACCp activity and performance compared to mice employing an exploitation strategy which reduced task complexity. Our results implicate ACCp claustrum neurons in restricting engagement during sleep and goal-directed behavior.

For nearly a century, evidence has accumulated indicating that the lateral hypothalamus (LH) contains neurons essential to sustain wakefulness. While lesion or inactivation of LH neurons produces a profound increase in sleep, stimulation of inhibitory LH neurons promotes wakefulness. To date, the primary wake-promoting cells that have been identified in the LH are the hypocretin/orexin (Hcrt) neurons, yet these neurons have little impact on total sleep or wake duration across the 24-h period. Recently, we and others have identified other LH populations that increase wakefulness. In the present study, we conducted microendoscopic calcium imaging in the LH concomitant with EEG and locomotor activity (LMA) recordings and found that a subset of LH neurons that express Ca 2+ /calmodulin-dependent protein kinase IIα (CaMKIIα) are preferentially active during wakefulness. Chemogenetic activation of these neurons induced sustained wakefulness and greatly increased LMA even in the absence of Hcrt signaling. Few LH CaMKIIα-expressing neurons are hypocretinergic or histaminergic while a small but significant proportion are GABAergic. Ablation of LH inhibitory neurons followed by activation of the remaining LH CaMKIIα neurons induced similar levels of wakefulness but blunted the LMA increase. Ablated animals showed no significant changes in sleep architecture but both spontaneous LMA and high theta (8 to 10 Hz) power during wakefulness were reduced. Together, these findings indicate the existence of two subpopulations of LH CaMKIIα neurons: an inhibitory population that promotes locomotion without affecting sleep architecture and an excitatory population that promotes prolonged wakefulness even in the absence of Hcrt signaling.

Melody is a core component of music in which discrete pitches are serially arranged to convey emotion and meaning. Perception varies along several pitch-based dimensions: (i) the absolute pitch of notes, (ii) the difference in pitch between successive notes, and (iii) the statistical expectation of each note given prior context. How the brain represents these dimensions and whether their encoding is specialized for music remains unknown. We recorded high-density neurophysiological activity directly from the human auditory cortex while participants listened to Western musical phrases. Pitch, pitch-change, and expectation were selectively encoded at different cortical sites, indicating a spatial map for representing distinct melodic dimensions. The same participants listened to spoken English, and we compared responses to music and speech. Cortical sites selective for music encoded expectation, while sites that encoded pitch and pitch-change in music used the same neural code to represent equivalent properties of speech. Findings reveal how the perception of melody recruits both music-specific and general-purpose sound representations.

Direct neural recordings from human auditory cortex have demonstrated encoding for acoustic-phonetic features of consonants and vowels. Neural responses also encode distinct acoustic amplitude cues related to timing, such as those that occur at the onset of a sentence after a silent period or the onset of the vowel in each syllable. Here, we used a group reduced rank regression model to show that distributed cortical responses support a low-dimensional latent state representation of temporal context in speech. The timing cues each capture more unique variance than all other phonetic features and exhibit rotational or cyclical dynamics in latent space from activity that is widespread over the superior temporal gyrus. We propose that these spatially distributed timing signals could serve to provide temporal context for, and possibly bind across time, the concurrent processing of individual phonetic features, to compose higher-order phonological (e.g. word-level) representations.

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.

Significance Bioelectronic medicine relies on electrical stimulation for most applications in the peripheral nervous system. It faces persistent challenges in selectively activating bundled nerve fibers. Here, we investigated ion-concentration modulation with ion-selective membranes and whether this modality may enhance the functional selectivity of peripheral nerve stimulation. We designed a multimodal stimulator that could control Ca 2+ concentrations within a focused volume. Acutely implanting it on the sciatic nerve of a rat, we demonstrated that Ca 2+ depletion could increase the sensitivity of the nerve to electrical stimulation in vivo. We provided evidence that it selectively influenced individual fascicles of the nerve, allowing selective activation by electrical current. Improved functional selectivity may improve outcomes for important therapeutic modalities.

Circuits of excitatory and inhibitory neurons generate gamma-rhythmic activity (30–80 Hz). Gamma-cycles show spontaneous variability in amplitude and duration. To investigate the mechanisms underlying this variability, we recorded local-field-potentials (LFPs) and spikes from awake macaque V1. We developed a noise-robust method to detect gamma-cycle amplitudes and durations, which showed a weak but positive correlation. This correlation, and the joint amplitude-duration distribution, is well reproduced by a noise-driven damped harmonic oscillator. This model accurately fits LFP power-spectra, is equivalent to a linear, noise-driven E-I circuit, and recapitulates two additional features of gamma: (1) Amplitude-duration correlations decrease with oscillation strength; (2) amplitudes and durations exhibit strong and weak autocorrelations, respectively, depending on oscillation strength. Finally, longer gamma-cycles are associated with stronger spike-synchrony, but lower spike-rates in both (putative) excitatory and inhibitory neurons. In sum, V1 gamma-dynamics are well described by the simplest possible model of gamma: A damped harmonic oscillator driven by noise.

We report on the ensemble organization of interneuronal activity in the spinal cord during locomotor movements and show that lumbar intermediate zone interneurons organize in two groups related to the two major phases of walking: stance and swing. Ensemble organization is also shown to better correlate with muscular output than single-cell activity, although ensemble membership does not appear to be somatotopically organized within the spinal cord.

The human hand is unique in the animal kingdom for unparalleled dexterity, ranging from complex prehension to fine finger individuation. How does the brain represent such a diverse repertoire of movements? We evaluated mesoscale neural dynamics across the human "grasp network," using electrocorticography and dimensionality reduction methods, for a repertoire of hand movements. Strikingly, we found that the grasp network represented both finger and grasping movements alike. Specifically, the manifold characterizing the multi-areal neural covariance structure was preserved during all movements across this distributed network. In contrast, latent neural dynamics within this manifold were surprisingly specific to movement type. Aligning latent activity to kinematics further uncovered distinct submanifolds despite similarities in synergistic coupling of joints between movements. We thus find that despite preserved neural covariance at the distributed network level, mesoscale dynamics are compartmentalized into movement-specific submanifolds; this mesoscale organization may allow flexible switching between a repertoire of hand movements.

When a visual stimulus is repeated, average neuronal responses typically decrease, yet they might maintain or even increase their impact through increased synchronization. Previous work has found that many repetitions of a grating lead to increasing gamma-band synchronization. Here, we show in awake macaque area V1 that both repetition-related reductions in firing rate and increases in gamma are specific to the repeated stimulus. These effects show some persistence on the timescale of minutes. Gamma increases are specific to the presented stimulus location. Further, repetition effects on gamma and on firing rates generalize to images of natural objects. These findings support the notion that gamma-band synchronization subserves the adaptive processing of repeated stimulus encounters.

Engagement is a major determinant of performance. Hyper-engagement risks impulsivity and is fatiguing over time, while hypo-engagement could lead to missed opportunities. Even in sleep, when engagement levels are minimal, sensory responsiveness varies. Thus, maintaining an optimal engagement level with the environment is a fundamental cognitive ability. The claustrum, and in particular its reciprocal connectivity with executive regions in the frontal cortex, has been associated with salience, attention and sleep. These apparently disparate roles can be consolidated within the context of engagement. Here we describe the activity of claustro-frontal circuits in a task imposing a tradeoff between response inhibition and sensory acuity (‘ENGAGE’). Recording calcium fiber photometry during >80,000 trials, we characterize claustrum recruitment during salient behavioral events, and find that a moderate level of activity in claustro-cingulate projections defines optimal engagement. Low activity of this pathway is associated with impulsive actions, while high activity is associated with behavioral lapses. Chemogenetic activation of cingulate-projecting claustrum neurons suppressed impulsive behavior and reduced the engagement of mice in the task. This relationship became even clearer upon addressing individual variability in the strategy mice employed during the ENGAGE task. Furthermore, this association of claustrum activity and engagement extends into sleep. Using simultaneous EEG and photometry recordings in the claustrum, we find that cingulate projecting claustrum neurons are most active during deep unresponsive slow-wave sleep, when mice are less prone to awakening by sensory stimuli.