EEG/ECoG/EMG publications

A Rodent Sleep Spindle Detector (RSSD) application (app) was developed to assist researchers working with high volume studies examining the impact of sleep on neurological function. Our RSSD is a MATLAB-based software program with a user interface that automatically identifies sleep spindles from intracranial EEG (iEEG) recordings of rodents using two novel yet complementary algorithmic approaches, a primary and secondary one. To validate the program, 6,000 copies of real spindles of 5 different types, ranging from 11–17 Hz with a duration of at least 0.3 seconds, were randomly placed within a noisy simulated prefrontal cortex iEEG signal of 50,000 seconds in duration. When compared to the ground truth on a datapoint-by-datapoint basis (individual spindle detection), the program had an accuracy of 98.40 ± 5.62% (mean ± SD) with 95% C.I. [91.93, 100] and 96.90 ± 4.34% (mean ± SD) with 95% C.I. [91.91, 100] for the primary and secondary algorithmic approach, respectively. Evaluating total spindle count, the program had an accuracy of 93.68 ± 13.66% (mean ± SD) with 95% C.I. [81.71, 100], and of 99.85 ± 0.12% (mean ± SD) with 95% C.I. [99.71, 99.96] for the primary and secondary algorithmic approach, respectively. The robustness of the sleep spindle detection was further validated for a range of spindle's duration, amplitude, and frequency by embedding in the iEEG signal respective artificial spindles. Finally, the RSSD app further improves its performance by first processing available video recordings of rodents to identify periods of quiescence and then running the sleep spindle detection algorithms on the iEEG only for those periods.

Mice are increasingly used as models of human-acquired neurological or neurodevelopmental conditions, such as autism, schizophrenia, and Alzheimer’s disease. All these conditions involve central auditory processing disorders, which have been little investigated despite their potential for providing interesting insights into the mechanisms behind such disorders. Alterations of the auditory steady-state response to 40 Hz click trains are associated with an imbalance between neuronal excitation and inhibition, a mechanism thought to be common to many neurological disorders. Here, we demonstrate the value of presenting click trains at various rates to mice with chronically implanted pins above the inferior colliculus and the auditory cortex for obtaining easy, reliable, and long-lasting access to subcortical and cortical complex auditory processing in awake mice. Using this protocol on a mutant mouse model of autism with a defect of the Shank3 gene, we show that the neural response is impaired at high click rates (above 60 Hz) and that this impairment is visible subcortically—two results that cannot be obtained with classical protocols for cortical EEG recordings in response to stimulation at 40 Hz. These results demonstrate the value and necessity of a more complete investigation of central auditory processing disorders in mouse models of neurological or neurodevelopmental disorders.

In tonal languages, which are spoken by nearly one-third of the world’s population, speakers precisely control the tension of vocal folds in the larynx to modulate pitch in order to distinguish words with completely different meanings. The specific pitch trajectories for a given tonal language are called lexical tones. Here, we used high-density direct cortical recordings to determine the neural basis of lexical tone production in native Mandarin-speaking participants. We found that instead of a tone category-selective coding, local populations in the bilateral laryngeal motor cortex (LMC) encode articulatory kinematic information to generate the pitch dynamics of lexical tones. Using a computational model of tone production, we discovered two distinct patterns of population activity in LMC commanding pitch rising and lowering. Finally, we showed that direct electrocortical stimulation of different local populations in LMC evoked pitch rising and lowering during tone production, respectively. Together, these results reveal the neural basis of vocal pitch control of lexical tones in tonal languages.

The human auditory system extracts rich linguistic abstractions from speech signals. Traditional approaches to understanding this complex process have used linear feature-encoding models, with limited success. Artificial neural networks excel in speech recognition tasks and offer promising computational models of speech processing. We used speech representations in state-of-the-art deep neural network (DNN) models to investigate neural coding from the auditory nerve to the speech cortex. Representations in hierarchical layers of the DNN correlated well with the neural activity throughout the ascending auditory system. Unsupervised speech models performed at least as well as other purely supervised or fine-tuned models. Deeper DNN layers were better correlated with the neural activity in the higher-order auditory cortex, with computations aligned with phonemic and syllabic structures in speech. Accordingly, DNN models trained on either English or Mandarin predicted cortical responses in native speakers of each language. These results reveal convergence between DNN model representations and the biological auditory pathway, offering new approaches for modeling neural coding in the auditory cortex.

The kindling model has been used extensively by researchers to study the neurobiology of temporal lobe epilepsy (TLE) due to its capacity to induce intensification of seizures by the progressive recruitment of additional neuronal clusters into epileptogenic networks. We applied repetitive focal optogenetic activation of putative excitatory neurons in the dorsal CA1 area of the hippocampus of mice to investigate the role of inhibitory signaling during this process. This experimental protocol resulted in a kindling phenotype that was maintained for 2 weeks after the animals were fully kindled. As a result of the different phases of optogenetic kindling (OpK), key inhibitory signaling elements, such as KCC2 and NKCC1, exhibited distinct temporal and spatial dynamics of regulation. These alterations in protein expression were related to the distinct pattern of ictal activity propagation through the different hippocampal sublayers. Our results suggest the KCC2 disruption in the contralateral hippocampus of fully kindled animals progressively facilitated the creation of pathological pathways for seizure propagation through the hippocampal network. Upon completion of kindling, we observed animals that were restimulated after a rest period of 14-day showed, besides a persistent KCC2 downregulation, an NKCC1 upregulation in the bilateral dentate gyrus and hippocampus-wide loss of parvalbumin-positive interneurons. These alterations observed in the chronic phase of OpK suggest that the hippocampus of rekindled animals continued to undergo self-modifications during the rest period. The changes resulting from this period suggest the possibility of the development of a mirror focus on the hippocampus contralateral to the site of optical stimulations. Our results offer perspectives for preventing the recruitment and conversion of healthy neuronal networks into epileptogenic ones among patients with epilepsy.

Our brains have a propensity to integrate closely-timed auditory and visual stimuli into a unified percept; a phenomenon that is highly malleable based on prior sensory experiences, and is known to be altered in clinical populations. While the neural correlates of audiovisual temporal perception have been investigated using neuroimaging and electroencephalography techniques in humans, animal research will be required to uncover the underlying cellular and molecular mechanisms. Prior to conducting such mechanistic studies, it is important to first confirm the translational potential of any prospective animal model. Thus, in the present study, we conducted a series of experiments to determine if rats show the hallmarks of audiovisual temporal perception observed in neurotypical humans, and whether the rat behavioral paradigms could reveal when they experienced perceptual disruptions akin to those observed in neurodevelopmental disorders. After training rats to perform a temporal order judgment (TOJ) or synchrony judgment (SJ) task, we found that the rats’ perception was malleable based on their past and present sensory experiences. More specifically, passive exposure to asynchronous audiovisual stimulation in the minutes prior to behavioral testing caused the rats’ perception to predictably shift in the direction of the leading stimulus; findings which represent the first time that this form of audiovisual perceptual malleability has been reported in non-human subjects. Furthermore, rats performing the TOJ task also showed evidence of rapid recalibration, in which their audiovisual temporal perception on the current trial was predictably influenced by the timing lag between the auditory and visual stimuli in the preceding trial. Finally, by manipulating either experimental testing parameters or altering the rats’ neurochemistry with a systemic injection of MK-801, we showed that the TOJ and SJ tasks could identify when the rats had difficulty judging the timing of audiovisual stimuli. These findings confirm that the behavioral paradigms are indeed suitable for future testing of rats with perceptual disruptions in audiovisual processing. Overall, our collective results highlight that rats represent an excellent animal model to study the cellular and molecular mechanisms underlying the acuity and malleability of audiovisual temporal perception, as they showcase the perceptual hallmarks commonly observed in humans.

Rapid-eye-movement (REM) sleep is accompanied by intense cortical activity, underlying its wake-like electroencephalogram (EEG). The neural activity inducing REM sleep is thought to originate from subcortical circuits in brainstem and hypothalamus. However, whether cortical neurons can also trigger REM sleep has remained unknown. Here, we show in mice that the medial prefrontal cortex (mPFC) strongly promotes REM sleep. Bidirectional optogenetic manipulations demonstrate that excitatory mPFC neurons promote REM sleep through their projections to the lateral hypothalamus (LH) and regulate phasic events, reflected in accelerated EEG theta oscillations and increased eye-movement density during REM sleep. Calcium imaging reveals that the majority of LH-projecting mPFC neurons are maximally activated during REM sleep and a subpopulation is recruited during phasic theta accelerations. Our results delineate a cortico-hypothalamic circuit for the top-down control of REM sleep and identify a critical role of the mPFC in regulating phasic events during REM sleep.

The human brain extracts meaning using an extensive neural system for semantic knowledge. Whether broadly distributed systems depend on or can compensate after losing a highly interconnected hub is controversial. We report intracranial recordings from two patients during a speech prediction task, obtained minutes before and after neurosurgical treatment requiring disconnection of the left anterior temporal lobe (ATL), a candidate semantic knowledge hub. Informed by modern diaschisis and predictive coding frameworks, we tested hypotheses ranging from solely neural network disruption to complete compensation by the indirectly affected language-related and speech-processing sites. Immediately after ATL disconnection, we observed neurophysiological alterations in the recorded frontal and auditory sites, providing direct evidence for the importance of the ATL as a semantic hub. We also obtained evidence for rapid, albeit incomplete, attempts at neural network compensation, with neural impact largely in the forms stipulated by the predictive coding framework, in specificity, and the modern diaschisis framework, more generally. The overall results validate these frameworks and reveal an immediate impact and capability of the human brain to adjust after losing a brain hub.

The R47H variant of triggering receptor expressed on myeloid cells 2 (TREM2) increases the risk of Alzheimer's disease (AD). To investigate potential mechanisms, we analyzed knockin mice expressing human TREM2-R47H from one mutant mouse Trem2 allele. TREM2-R47H mice showed increased seizure activity in response to an acute excitotoxin challenge, compared to wildtype controls or knockin mice expressing the common variant of human TREM2. TREM2-R47H also increased spontaneous thalamocortical epileptiform activity in App knockin mice expressing amyloid precursor proteins bearing autosomal dominant AD mutations and a humanized amyloid-β sequence. In mice with or without such App modifications, TREM2-R47H increased the density of putative synapses in cortical regions without amyloid plaques. TREM2-R47H did not affect synaptic density in hippocampal regions with or without plaques. We conclude that TREM2-R47H increases AD-related network hyperexcitability and that it may do so, at least in part, by causing an imbalance in synaptic densities across brain regions.

Objective There is a complex interaction between nervous and cardiovascular systems, but sparse data exist on brain–heart electrophysiological responses to cardiac arrest resuscitation. Our aim was to investigate dynamic changes in autonomic and cortical function during hyperacute stage post‐resuscitation. Methods Ten rats were resuscitated from 7‐min cardiac arrest, as indicators of autonomic response, heart rate (HR), and its variability (HRV) were measured. HR was monitored through continuous electrocardiography, while HRV was assessed via spectral analysis, whereby the ratio of low−/high‐frequency (LF/HF) power indicates the balance between sympathetic/parasympathetic activities. Cortical response was evaluated by continuous electroencephalography and quantitative analysis. Parameters were quantified at 5‐min intervals over the first‐hour post‐resuscitation. Neurological outcome was assessed by Neurological Deficit Score (NDS, range 0–80, higher = better outcomes) at 4‐h post‐resuscitation. Results A significant increase in HR was noted over 15–30 min post‐resuscitation ( p < 0.01 vs.15‐min, respectively) and correlated with higher NDS ( r s = 0.56, p < 0.01). LF/HF ratio over 15–20 min was positively correlated with NDS ( r s = 0.75, p < 0.05). Gamma band power surged over 15–30 min post‐resuscitation ( p < 0.05 vs. 0–15 min, respectively), and gamma band fraction during this period was associated with NDS ( r s ≥0.70, p < 0.05, respectively). Significant correlations were identified between increased HR and gamma band power during 15–30 min ( r s ≥0.83, p < 0.01, respectively) and between gamma band fraction and LF/HF ratio over 15–20 min post‐resuscitation ( r s = 0.85, p < 0.01). Interpretations Hyperacute recovery of autonomic and cortical function is associated with favorable functional outcomes. While this observation needs further validation, it presents a translational opportunity for better autonomic and neurologic monitoring during early periods post‐resuscitation to develop novel interventions.

Epilepsy is characterized by spontaneous non-provoked seizures, yet the mechanisms that trigger a seizure and allow its evolution remain underexplored. To dissect out phases of ictogenesis, we evoked hypersynchronous activity with optogenetic stimulation. Focal optogenetic activation of putative excitatory neurons in the mouse hippocampal CA1 reliably evoked convulsive seizures in awake mice. A time-vs-time pulsogram plot characterized the evolution of the EEG pulse response from a light evoked response to induced seizure activity. Our results depict ictogenesis as a stepwise process comprised of three distinctive phases demarcated by two transition points. The induction phase undergoes the first transition to reverberant phase activity, followed by the second transition into the paroxysmal phase or a seizure. Non-seizure responses are confined to either induction or reverberant phases. The pulsogram was then constructed in seizures recorded from a murine model of temporal lobe epilepsy and it depicted a similar reverberance preceding spontaneous seizures. The discovery of these distinct phases of ictogenesis may offer means to abort a seizure before it develops.

The auditory system has exquisite temporal coding in the periphery which is transformed into a rate-based code in central auditory structures, like auditory cortex. However, the cortex is still able to synchronize, albeit at lower modulation rates, to acoustic fluctuations. The perceptual significance of this cortical synchronization is unknown. We estimated physiological synchronization limits of cortex (in humans with electroencephalography) and brainstem neurons (in chinchillas) to dynamic binaural cues using a novel system-identification technique, along with parallel perceptual measurements. We find that cortex can synchronize to dynamic binaural cues up to approximately 10 Hz, which aligns well with our measured limits of perceiving dynamic spatial information and utilizing dynamic binaural cues for spatial unmasking, i.e. measures of binaural sluggishness. We also find that the tracking limit for frequency modulation (FM) is similar to the limit for spatial tracking, demonstrating that this sluggish tracking is a more general perceptual limit that can be accounted for by cortical temporal integration limits.

Spoken language contains information at a broad range of timescales, from phonetic distinctions on the order of milliseconds to semantic contexts which shift over seconds to minutes. It is not well understood how the brain’s speech production systems combine features at these timescales into a coherent vocal output. We investigated the spatial and temporal representations in cerebral cortex of three phonological units with different durations: consonants, vowels, and syllables. Electrocorticography (ECoG) recordings were obtained from five participants while speaking single syllables. We developed a novel clustering and Kalman filter-based trend analysis procedure to sort electrodes into temporal response profiles. A linear discriminant classifier was used to determine how strongly each electrode’s response encoded phonological features. We found distinct time-courses of encoding phonological units depending on their duration: consonants were represented more during speech preparation, vowels were represented evenly throughout trials, and syllables during production. Locations of strongly speech-encoding electrodes (the top 30% of electrodes) likewise depended on phonological element duration, with consonant-encoding electrodes left-lateralized, vowel-encoding hemispherically balanced, and syllable-encoding right-lateralized. The lateralization of speech-encoding electrodes depended on onset time, with electrodes active before or after speech production favoring left hemisphere and those active during speech favoring the right. Single-electrode speech classification revealed cortical areas with preferential encoding of particular phonemic elements, including consonant encoding in the left precentral and postcentral gyri and syllable encoding in the right middle frontal gyrus. Our findings support neurolinguistic theories of left hemisphere specialization for processing short-timescale linguistic units and right hemisphere processing of longer-duration units.

The mechanisms underlying tinnitus perception are still under research. One of the proposed hypotheses involves an alteration in top-down processing of auditory activity. Low-frequency oscillations in the delta and theta bands have been recently described in brain and cochlear infrasonic signals during selective attention paradigms in normal hearing controls. Here, we propose that the top-down oscillatory activity observed in brain and cochlear signals during auditory and visual selective attention in normal subjects, is altered in tinnitus patients, reflecting an abnormal functioning of the corticofugal pathways that connect brain circuits with the cochlear receptor. Methods: To test this hypothesis, we used a behavioral task that alternates between auditory and visual top-down attention while we simultaneously measured electroencephalogram (EEG) and distortion-product otoacoustic emissions (DPOAE) signals in 14 tinnitus and 14 control subjects. Results: We found oscillatory activity in the delta and theta bands in cortical and cochlear channels in control and tinnitus patients. There were significant decreases in the DPOAE oscillatory amplitude during the visual attention period as compared to the auditory attention period in tinnitus and control groups, while frontal EEG oscillatory activity (Fz) was increased during the visual attention task only in the tinnitus group. The difference between auditory and visual DPOAE oscillatory power in the delta band was stronger in tinnitus individuals as compared to controls. Conclusions: These results confirm the presence of top-down infrasonic low-frequency cochlear oscillatory activity in the delta and theta bands in tinnitus patients. Our results show an increase in the EEG theta oscillatory power in frontal regions of tinnitus sufferers concomitant with a stronger corticofugal suppression of cochlear oscillations in the delta band during visual attention in tinnitus. These findings suggest an increased compensatory corticofugal inhibitory mechanism to the cochlear receptor in tinnitus.

Importance Unexplained chronic cough is common and has substantial negative quality-of-life implications, yet its causes are not well understood. A better understanding of how peripheral and central neural processes contribute to chronic cough is essential for treatment design. Objective To determine if people with chronic cough exhibit signs of abnormal neural processing over laryngeal sensorimotor cortex during voluntary laryngeal motor activity such as vocalization. Design, Setting, and Participants This was a cross-sectional study of a convenience sample of participants with chronic cough and healthy participants. Testing was performed in an acoustically and electromagnetically shielded chamber. In a single visit, electroencephalographic (EEG) signals were recorded from participants with chronic cough and healthy participants during voice production. The chronic cough group participants presented with unexplained cough of 8 weeks or longer duration with prior medical evaluation including negative results of chest imaging. None of the participants had a history of any neurologic disease known to impair vocalization or swallowing. Data collection for the healthy control group occurred from February 2 to June 28, 2018, and for the chronic cough group, from November 22, 2021, to June 21, 2022. Data analysis was performed from May 1 to October 30, 2022. Exposure Participants with or without chronic cough. Main Outcome Measures Event-related spectral perturbation over the laryngeal area of somatosensory-motor cortex from 0 to 30 Hz (ie, θ, α, and β bands) and event-related coherence as a measure of synchronous activity between somatosensory and motor cortical regions. Results The chronic cough group comprised 13 participants with chronic cough (mean [SD] age, 63.5 [7.8] years; 9 women and 4 men) and the control group, 10 healthy age-matched individuals (mean [SD] age, 60.3 [13.9] years; 6 women and 4 men). In the chronic cough group, the typical movement-related desynchronization over somatosensory-motor cortex during vocalization was significantly reduced across θ, α, and β frequency bands when compared with the control group. Conclusions and Relevance This cross-sectional study found that the typical movement-related suppression of brain oscillatory activity during vocalization is weak or absent in people with chronic cough. Thus, chronic cough affects sensorimotor cortical activity during the asymptomatic voluntary activation of laryngeal muscles.

Context modulates how information is processed in the mammalian brain. For example, brain responses are amplified to contextually unusual stimuli. This phenomenon, known as "deviance detection,"1,2 is well documented in early, primary sensory cortex, where large responses are generated to simple stimuli that deviate from their context in low-order properties, such as line orientation, size, or pitch.2,3,4,5 However, the extent to which neural deviance detection manifests (1) in broader cortical networks and (2) to simple versus complex stimuli, which deviate only in their higher-order, multisensory properties, is not known. Consistent with a predictive processing framework,6,7 we hypothesized that deviance detection manifests in a hierarchical manner across cortical networks,8,9 emerging later and further downstream when stimulus deviance is complex. To test this, we examined brain responses of awake mice to simple unisensory deviants (e.g., visual line gratings, deviating from context in their orientation alone) versus complex multisensory deviants (i.e., audiovisual pairs, deviating from context only in their audiovisual pairing but not visual or auditory content alone). We find that mouse parietal associative area-a higher cortical region-displays robust multisensory deviance detection. In contrast, primary visual cortex exhibits strong unisensory visual deviance detection but weaker multisensory deviance detection. These results suggest that deviance detection signals in the cortex may be conceptualized as "prediction errors," which are primarily fed forward-or downstream-in cortical networks.6,7.

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.

Visual stimuli that deviate from the current context elicit augmented responses in the primary visual cortex (V1). These heightened responses, known as "deviance detection," require local inhibition in the V1 and top-down input from the anterior cingulate area (ACa). Here, we investigated the mechanisms by which the ACa and V1 interact to support deviance detection. Local field potential recordings in mice during an oddball paradigm showed that ACa-V1 synchrony peaks in the theta/alpha band (≈10 Hz). Two-photon imaging in the V1 revealed that mainly pyramidal neurons exhibited deviance detection, while contextually redundant stimuli increased vasoactive intestinal peptide (VIP)-positive interneuron (VIP) activity and decreased somatostatin-positive interneuron (SST) activity. Optogenetic drive of ACa-V1 inputs at 10 Hz activated V1-VIPs but inhibited V1-SSTs, mirroring the dynamics present during the oddball paradigm. Chemogenetic inhibition of V1-VIPs disrupted Aca-V1 synchrony and deviance detection in the V1. These results outline temporal and interneuron-specific mechanisms of top-down modulation that support visual context processing.

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.

Introduction The therapeutic efficacy of deep brain stimulation (DBS) of the subthalamic nucleus (STN) for Parkinson’s disease (PD) may be limited for some patients by the presence of stimulation-related side effects. Such effects are most often attributed to electrical current spread beyond the target region. Prior computational modeling studies have suggested that changing the degree of asymmetry of the individual phases of the biphasic, stimulus pulse may allow for more selective activation of neural elements in the target region. To the extent that different neural elements contribute to the therapeutic vs. side-effect inducing effects of DBS, such improved selectivity may provide a new parameter for optimizing DBS to increase the therapeutic window. Methods We investigated the effect of six different pulse geometries on cortical and myogenic evoked potentials in eight patients with PD whose leads were temporarily externalized following STN DBS implant surgery. DBS-cortical evoked potentials were quantified using peak to peak measurements and wavelets and myogenic potentials were quantified using RMS. Results We found that the slope of the recruitment curves differed significantly as a function of pulse geometry for both the cortical- and myogenic responses. Notably, this effect was observed most frequently when stimulation was delivered using a monopolar, as opposed to a bipolar, configuration. Discussion Manipulating pulse geometry results in differential physiological effects at both the cortical and neuromuscular level. Exploiting these differences may help to expand DBS’ therapeutic window and support the potential for incorporating pulse geometry as an additional parameter for optimizing therapeutic benefit.