EEG/ECoG/EMG publications

The mediodorsal thalamic nucleus (MD) forms neural circuits with various brain regions, including the prefrontal cortex (PFC), the reticular thalamic nucleus (TRN), and the midbrain reticular nucleus (MRN). However, the specific roles and underlying mechanisms in auditory information processing remain unclear. Here, we perform multi-channel electrophysiological recordings in awake mice to investigate the response patterns of the MD to auditory stimuli, as well as the regulatory effects of PFC, MRN, and TRN inputs. We identify two distinct types of sound-evoked responses. The Phasic-response features a transient burst firing to the stimulus with short latency, rapidly adapting to baseline and corresponding to the onset fluctuation of the local field potential. The Sustained-response is marked by prolonged firing with longer latency and is accompanied by persistent enhancement of oscillatory power following stimulus offset. The response patterns of MD neurons remain consistent across different types of auditory stimuli. Optogenetic inactivation of the MRN suppresses both response types in the MD. The Sustained-response is attenuated by PFC inactivation but enhanced by TRN inactivation, while the Phasic-response remains unaffected by inactivation of either the PFC or TRN. Our findings expand the understanding of the MD's role in sound information integration and auditory cognitive regulation.

Objective Cortical delta band (1–4 Hz) activity is considered a biomarker for states of altered consciousness, with increased delta power observed during anesthesia, sleep, coma, and delirium. The current study sought to characterize delta power following electrode implantation with respect to patient demographics and clinical characteristics as well as type and duration of surgery. Methods Participants were 25 adult neurosurgical patients implanted with intracranial electrodes for clinical monitoring of their epilepsy. Resting state cortical activity was recorded at multiple occasions over the course of the monitoring period. The initial time point was defined as the first recording within 72 h following surgery. Analyses of cortical activity were conducted using a linear mixed effects modeling approach to account for within‐participant correlations and between‐participant heterogeneity. Results Throughout the monitoring period, delta power decreased in frontal, occipital, parietal, and temporal regions, indicating a global phenomenon. By contrast, beta (14–30 Hz) power remained stable. Delta power was higher following surgical cases that required craniotomy compared to stereoelectroencephalography cases. Surgery duration and anesthesia emergence duration were associated with higher delta power. Recordings from depth electrodes showed higher delta power compared to subdural electrodes. No significant effects of patients' age, sex, white blood cell count, antiseizure medication, and opioid medication dosage on postoperative delta power were found. Significance The results are consistent with a postoperative elevation in delta power that resolves over the course of the monitoring period and indicate an association between increased delta power and craniotomy surgery, as well as longer surgery and emergence durations. The current work provides a comprehensive analysis of surgical, clinical, and physiological factors, suggests risk factors, and lays fundamental groundwork for future studies.

Restoring upper limb function is critical in individuals with central paralysis, and hand control is a priority in patients with neurological impairments. Functional electrical stimulation with implantable electrodes targeting the peripheral nervous system has the potential to selectively recruit hand muscles and generate multiple functional hand movements. However, the implantation of electrodes in the forearm or elbow areas requires multiple incisions for surgery, and elbow joint movements cannot be performed. In this study, we designed and implanted two epineural cuffs on the median and radial nerves in the distal brachial plexus of a single Japanese macaque ( Macaca fuscata ) monkey. The cuffs were successfully placed via an axillary approach using a single incision. Electrical stimuli were applied to innervate the contraction patterns of the hand, forearm, and triceps muscles relevant to the median and radial nerves. The evoked potentials of the target muscles electrically stimulated the distal brachial plexus to reliably and selectively innervate the upper limb muscles at the functional group level. Our results demonstrated that the distal brachial plexus can be a useful stimulation site for upper limb muscle contraction and that the axillary approach enables electrode placement to peripheral nerves required for upper limb control.

Closed-loop neural control is a powerful tool for both the scientific exploration of neural function and for mitigating deficiencies found in open-loop deep brain stimulation (DBS). Here, we present a protocol for artificial intelligence-guided neural control in rats using deep reinforcement learning (RL) and infrared neural stimulation (INS). We describe steps for integrating RL closed-loop control into neuroscience and neuromodulation studies. We then detail procedures for using and deploying infrared INS in chronic DBS applications. For complete details on the use and execution of this protocol, please refer to Coventry et al.1 and Coventry and Bartlett.2.

In this paper we aim to examine the correlation between diffusion tensor imaging parameters of anatomical connectivity and characteristics of signal transmission obtained from patient-specific transfer function models. Here, we focused on elucidating the correlation between structural and functional neural connectivity within a cohort of pediatric patients diagnosed with dystonia. Diffusion tractography images were obtained from 12 patients with dystonia prior to the deep brain stimulation surgery. For each patient, we processed the imaging data to estimate anatomical measures including fractional anisotropy, axial diffusivity, number of fibre tracts per unit area, and fibre tract length. After the implantation of temporary depth leads for each patient as part of their treatment plan, intracranial signals were recorded. Transfer function models of local field potential recordings and the corresponding measures of functional connectivity were computed for each patient. Linear mixed effect analysis was then employed to determine the relationship between transfer function measures and diffusion tractography parameters. Our results illustrate a positive correlation between fractional anisotropy, AD, and intrinsic neural transmission measures, representing amplification and spread of intrinsic neural oscillations, obtained from the transfer functions models. However, no significant correlation was found between the functional connectivity and number of fibre tracts or fibre lengths. Our findings suggest that white matter integrity, as measured by fractional anisotropy and AD, can potentially reflect the amplification and spread of intrinsic brain signals throughout the network. This study underscores the significant relationship between structural and functional connectivity, offering valuable insights into propagation of neural activity in the brain network and potential implications for optimizing non-invasive treatments and planning for neurological disorders.

Micro-electrocorticography (µECoG) electrodes have emerged to balance the trade-off between invasiveness and signal quality in brain recordings. However, its large-scale applicability is still hindered by a lack of comparative studies assessing the relationship between ECoG and traditional recording methods such as penetrating electrodes. This study aimed to compare somatosensory evoked potentials (SEPs) through the lenses of a µECoG and an intracortical microelectrode array (MEA). The electrodes were implanted in the pig’s primary somatosensory cortex, while SEPs were generated by applying electrical stimulation to the ulnar nerve. The SEP amplitude, signal-to-noise ratio (SNR), power spectral density (PSD), and correlation structure were analysed. Overall, SEPs resulting from MEA recordings had higher amplitudes and contained significantly more spectral power, especially at higher frequencies. However, the SNRs were similar between the interfaces. These results demonstrate the feasibility of using µECoG to decode SEPs with wide-range applications in physiology monitoring and brain–computer interfaces.

Here, we characterized the dynamics of sleep spindles, focusing on their damping, which we estimated using a metric called oscillatory-Quality (o-Quality), derived by fitting an autoregressive model to electrophysiological signals, recorded from the cortex in mice. The o-Quality of sleep spindles correlates weakly with their amplitude, shows marked laminar differences and regional topography across cortical regions, reflects the level of synchrony within and between cortical networks, is strongly modulated by sleep-wake history, reflects the degree of sensory disconnection, and correlates with the strength of coupling between spindles and slow waves. As most spindle events are highly localized and not detectable with conventional low-density recording approaches, o-Quality thus emerges as a valuable metric that allows us to infer the spread and dynamics of spindle activity across the brain and directly links their spatiotemporal dynamics with local and global regulation of brain states, sleep regulation, and function.

In autism spectrum disorder (ASD), atypical sensory experiences are often associated with irregularities in predictive coding, which proposes that the brain creates hierarchical sensory models via a bidirectional process of predictions and prediction errors. However, it remains unclear how these irregularities manifest across different functional hierarchies in the brain. To address this, we study a marmoset model of ASD induced by valproic acid (VPA) treatment. We record high-density electrocorticography (ECoG) during an auditory task with two layers of temporal control, and applied a quantitative model to quantify the integrity of predictive coding across two distinct hierarchies. Our results demonstrate a persistent pattern of sensory hypersensitivity and unstable predictions across two brain hierarchies in VPA-treated animals, and reveal the associated spatio-spectro-temporal neural signatures. Despite the regular occurrence of imprecise predictions in VPA-treated animals, we observe diverse configurations of underestimation or overestimation of sensory regularities within the hierarchies. Our results demonstrate the coexistence of the two primary Bayesian accounts of ASD: overly-precise sensory observations and weak prior beliefs, and offer a potential multi-layered biomarker for ASD, which could enhance our understanding of its diverse symptoms.

Converging electrophysiological, molecular and ultrastructural evidence supports the hypothesis that sleep promotes a net decrease in excitatory synaptic strength, counteracting the net synaptic potentiation caused by ongoing learning during waking. However, several outstanding questions about sleep‐dependent synaptic weakening remain. Here, we address some of these questions by using two established molecular markers of synaptic strength, the levels of the AMPA (alpha‐amino‐3‐hydroxy‐5‐methyl‐4‐isoxazole propionic acid) receptors containing the GluA1 subunit and the phosphorylation of GluA1 at serine 845 (p‐GluA1(845)). We previously found that, in the rat cortex and hippocampus, these markers are lower after 6–8 h of sleep than after the same time spent awake. Here, we measure GluA1 and p‐GluA1(845) levels in synaptosomes of mouse cortex after 5 h of either sleep, sleep deprivation, recovery sleep after sleep deprivation or selective REM sleep deprivation (32 C57BL/B6 adult mice, 16 females). We find that relative to after sleep deprivation, these synaptic markers are lower after sleep independent of whether the mice were allowed to enter REM sleep. Moreover, 5 h of recovery sleep following acute sleep deprivation is enough to renormalize their expression. Thus, the renormalization of GluA1 and p‐GluA1(845) expression crucially relies on NREM sleep and can occur in a few hours of sleep after acute sleep deprivation.

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.

The cerebral cortex is widely considered part of the neural substrate of consciousness, but direct causal evidence is missing. Here, we tested in mice whether optogenetic activation of cortical neurons in posterior parietal cortex (PtA) or medial prefrontal cortex (mPFC) is sufficient for arousal from three behavioral states characterized by progressively deeper unresponsiveness: sleep, a coma-like state induced by muscimol injection in the midbrain, and deep sevoflurane-dexmedetomidine anesthesia. We find that cortical stimulation always awakens the mice from both NREM sleep and REM sleep, with PtA requiring weaker/shorter light pulses than mPFC. Moreover, in most cases light pulses produce both cortical activation (decrease in low frequencies) and behavioral arousal (recovery of the righting reflex) from brainstem coma, as well as cortical activation from anesthesia. These findings provide evidence that direct activation of cortical neurons is sufficient for behavioral and/or cortical arousal from sleep, brainstem coma, and anesthesia.

Introduction: Evoked compound action potentials (ECAPs) measured using epidural spinal recordings (ESRs) during epidural spinal cord stimulation (SCS) can help elucidate fundamental mechanisms for the treatment of pain, as well as inform closed-loop control of SCS. Previous studies have used ECAPs to characterize the neural response to various neuromodulation therapies and have demonstrated that ECAPs are highly prone to multiple sources of artifact, including post-stimulus pulse capacitive artifact, electromyography (EMG) bleed-through, and motion artifact resulting from disturbance of the electrode/tissue interface during normal behavior. However, a thorough characterization has yet to be performed for how these sources of artifact may contaminate recordings within the temporal window commonly used to determine activation of A-beta fibers in a large animal model. Methods: We characterized the sources of artifacts that can contaminate the recording of ECAPs in an epidural SCS swine model using the Abbott Octrode™ lead. Muscle paralytics were administered to block muscle activation preventing EMG from contaminating the recorded ECAPs. Concurrent EMG recordings of the longissimus, a long muscle of the back, were used to confirm a 2-4 millisecond (ms) latency source of EMG bleed-through that frequently contaminated the A-beta temporal window. Additionally, we obtained recordings approximately 5-10 minutes post-mortem after clear evoked A-beta and associated EMG responses ceased to characterize the representation of stimulation artifact across the array. Results: Spinal ECAP recordings can be contaminated by capacitive artifact, short latency EMG from nearby long muscles of the back, and motion artifact from multiple sources. In many cases, the capacitive artifact can appear nearly identical in duration and waveshape to evoked A-beta responses. These sources of EMG can have phase shifts across the electrode array, very similar to the phase shift anticipated by propagation of an evoked A-beta fiber response across the array. This short latency EMG is often evident at currents similar to those needed to activate A-beta fibers associated with the treatment of pain. Changes in cerebrospinal fluid between the cord and dura, and motion induced during breathing created a cyclic oscillation in all evoked components of the recorded ECAP signal. Conclusion: Careful controls must be implemented to accurately separate neural signal from the sources of artifact in spinal cord ECAPs. To address this, we suggest experimental procedures and associated reporting requirements necessary to disambiguate the underlying neural response from these confounds. These data are important to better understand the conceptual framework for recorded ESRs, with components such as ECAPs, EMG responses and artifacts, and have important implications for closed-loop control algorithms to account for transient motion such as postural changes and cough.

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.

Despite the widespread application of general anesthesia since the 19th century, a fundamental understanding of its operating principle and a reliable metric of consciousness remains elusive. We here analyze Loss of Responsiveness (LoR) induced by propofol anesthesia in humans and macaque monkeys. Through application of “highly comparable time series analysis” (hctsa) over 7000 features, we aim to identify LoR metrics common to both humans and macaque monkeys. In Stage 1 analysis, we trained a nearest median classifier using 200 ms data epochs using one macaque individual, and evaluated its classification accuracy on human intracranial data, recorded from one participant. This analysis suggested the existence of univariate features capable of discriminating between awake and propofol-induced unresponsive states in both species. Stage 2 extended the analysis to additional macaque and human recordings to evaluate cross-session and cross-species generalization. We identified many hctsa features that commonly distinguished awake from unresponsive states across sessions and species, in the frontal, but not in parietal and temporal, recordings. In exploratory analyses, we aimed to distinguish loss of responsiveness from consciousness, by additionally analyzing human data under sedated yet conscious state in a within-subject manner. This analysis revealed that frontal, temporal and parietal recordings contain similarly sensitive features and channels as signatures of consciousness. These results support the existence of ECoG time-series features associated with propofol-induced loss of responsiveness across humans and macaques, while also highlighting the challenge to determine whether such features specifically track conscious experience rather than responsiveness.

Effective neural stimulation requires adequate parametrization. Gaussian-process (GP)-based Bayesian optimization (BO) offers a framework to discover optimal stimulation parameters in real time. Here, we first provide a general protocol to deploy this framework in neurostimulation interventions and follow by exemplifying its use in detail. Specifically, we describe the steps to implant rats with multi-channel electrode arrays in the hindlimb motor cortex. We then detail how to utilize the GP-BO algorithm to maximize evoked target movements, measured as electromyographic responses. For complete details on the use and execution of this protocol, please refer to Bonizzato and colleagues (2023).1.

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.

Introduction Cochlear implants (CIs) are the treatment of choice for severe to profound hearing loss. Variability in CI outcomes remains despite advances in technology and is attributed in part to differences in cortical processing. Studying these differences in CI users is technically challenging. Spectrally degraded stimuli presented to normal-hearing individuals approximate input to the central auditory system in CI users. This study used intracranial electroencephalography (iEEG) to investigate cortical processing of spectrally degraded speech. Methods Participants were adult neurosurgical epilepsy patients. Stimuli were utterances /aba/ and /ada/, spectrally degraded using a noise vocoder (1–4 bands) or presented without vocoding. The stimuli were presented in a two-alternative forced choice task. Cortical activity was recorded using depth and subdural iEEG electrodes. Electrode coverage included auditory core in posteromedial Heschl’s gyrus (HGPM), superior temporal gyrus (STG), ventral and dorsal auditory-related areas, and prefrontal and sensorimotor cortex. Analysis focused on high gamma (70–150 Hz) power augmentation and alpha (8–14 Hz) suppression. Results Chance task performance occurred with 1–2 spectral bands and was near-ceiling for clear stimuli. Performance was variable with 3–4 bands, permitting identification of good and poor performers. There was no relationship between task performance and participants demographic, audiometric, neuropsychological, or clinical profiles. Several response patterns were identified based on magnitude and differences between stimulus conditions. HGPM responded strongly to all stimuli. A preference for clear speech emerged within non-core auditory cortex. Good performers typically had strong responses to all stimuli along the dorsal stream, including posterior STG, supramarginal, and precentral gyrus; a minority of sites in STG and supramarginal gyrus had a preference for vocoded stimuli. In poor performers, responses were typically restricted to clear speech. Alpha suppression was more pronounced in good performers. In contrast, poor performers exhibited a greater involvement of posterior middle temporal gyrus when listening to clear speech. Discussion Responses to noise-vocoded speech provide insights into potential factors underlying CI outcome variability. The results emphasize differences in the balance of neural processing along the dorsal and ventral stream between good and poor performers, identify specific cortical regions that may have diagnostic and prognostic utility, and suggest potential targets for neuromodulation-based CI rehabilitation strategies.

Electronic retinal prostheses for stimulating retinal neurons are promising for vision restoration. However, the rigid electrodes of conventional retinal implants can inflict damage on the soft retina tissue. They also have limited selectivity due to their poor proximity to target cells in the degenerative retina. Here we present a soft artificial retina (thickness, 10 μm) where flexible ultrathin photosensitive transistors are integrated with three-dimensional stimulation electrodes of eutectic gallium–indium alloy. Platinum nanoclusters locally coated only on the tip of these three-dimensional liquid-metal electrodes show advantages in reducing the impedance of the stimulation electrodes. These microelectrodes can enhance the proximity to the target retinal ganglion cells and provide effective charge injections (72.84 mC cm −2 ) to elicit neural responses in the retina. Their low Young’s modulus (234 kPa), owing to their liquid form, can minimize damage to the retina. Furthermore, we used an unsupervised machine learning approach to effectively identify the evoked spikes to grade neural activities within the retinal ganglion cells. Results from in vivo experiments on a retinal degeneration mouse model reveal that the spatiotemporal distribution of neural responses on their retina can be mapped under selective localized illumination areas of light, suggesting the restoration of their vision.

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.

Understanding central auditory processing critically depends on defining underlying auditory cortical networks and their relationship to the rest of the brain. We addressed these questions using resting state functional connectivity derived from human intracranial electroencephalography. Mapping recording sites into a low-dimensional space where proximity represents functional similarity revealed a hierarchical organization. At a fine scale, a group of auditory cortical regions excluded several higher-order auditory areas and segregated maximally from the prefrontal cortex. On mesoscale, the proximity of limbic structures to the auditory cortex suggested a limbic stream that parallels the classically described ventral and dorsal auditory processing streams. Identities of global hubs in anterior temporal and cingulate cortex depended on frequency band, consistent with diverse roles in semantic and cognitive processing. On a macroscale, observed hemispheric asymmetries were not specific for speech and language networks. This approach can be applied to multivariate brain data with respect to development, behavior, and disorders.