Invariance of firing rate and field potential dynamics to stimulus modulation rate in human auditory cortex

Gamma amplitude-envelope correlations are strongly elevated within hyperexcitable networks in focal epilepsy

Methods to quantify cortical hyperexcitability are of enormous interest for mapping epileptic networks in patients with focal epilepsy. We hypothesize that, in the resting state, cortical hyperexcitability increases firing-rate correlations between neuronal populations within seizure onset zones (SOZs). This hypothesis predicts that in the gamma frequency band (40-200 Hz), amplitude envelope correlations (AECs), a relatively straightforward measure of functional connectivity, should be elevated within SOZs compared to other areas. To test this prediction, we analyzed archived samples of interictal electrocorticographic (ECoG) signals recorded from patients who became seizure-free after surgery targeting SOZs identified by multiday intracranial recordings. We show that in the gamma band, AECs between nodes within SOZs are markedly elevated relative to those elsewhere. AEC-based node strength, eigencentrality, and clustering coefficient are also robustly increased within the SOZ with maxima in the low-gamma band (permutation test Z-scores > 8) and yield moderate discriminability of the SOZ using ROC analysis (maximal mean AUC ~ 0.73). By contrast to AECs, phase locking values (PLVs), a measure of narrow-band phase coupling across sites, and PLV-based graph metrics discriminate the seizure onset nodes weakly. Our results suggest that gamma band AECs may provide a clinically useful marker of cortical hyperexcitability in focal epilepsy.

Local field potentials, spiking activity, and receptive fields in human visual cortex

The concept of receptive field (RF) is central to sensory neuroscience. Neuronal RF properties have been substantially studied in animals, while those in humans remain nearly unexplored. Here, we measured neuronal RFs with intracranial local field potentials (LFPs) and spiking activity in human visual cortex (V1/V2/V3). We recorded LFPs via macro-contacts and discovered that RF sizes estimated from low-frequency activity (LFA, 0.5-30 Hz) were larger than those estimated from low-gamma activity (LGA, 30-60 Hz) and high-gamma activity (HGA, 60-150 Hz). We then took a rare opportunity to record LFPs and spiking activity via microwires in V1 simultaneously. We found that RF sizes and temporal profiles measured from LGA and HGA closely matched those from spiking activity. In sum, this study reveals that spiking activity of neurons in human visual cortex could be well approximated by LGA and HGA in RF estimation and temporal profile measurement, implying the pivotal functions of LGA and HGA in early visual information processing.

Single-cell activity in human STG during perception of phonemes is organized according to manner of articulation

One of the central tasks of the human auditory system is to extract sound features from incoming acoustic signals that are most critical for speech perception. Specifically, phonological features and phonemes are the building blocks for more complex linguistic entities, such as syllables, words and sentences. Previous ECoG and EEG studies showed that various regions in the superior temporal gyrus (STG) exhibit selective responses to specific phonological features. However, electrical activity recorded by ECoG or EEG grids reflects average responses of large neuronal populations and is therefore limited in providing insights into activity patterns of single neurons. Here, we recorded spiking activity from 45 units in the STG from six neurosurgical patients who performed a listening task with phoneme stimuli. Fourteen units showed significant responsiveness to the stimuli. Using a Naïve-Bayes model, we find that single-cell responses to phonemes are governed by manner-of-articulation features and are organized according to sonority with two main clusters for sonorants and obstruents. We further find that 'neural similarity' (i.e. the similarity of evoked spiking activity between pairs of phonemes) is comparable to the 'perceptual similarity' (i.e. to what extent two phonemes are judged as sounding similar) based on perceptual confusion, assessed behaviorally in healthy subjects. Thus, phonemes that were perceptually similar also had similar neural responses. Taken together, our findings indicate that manner-of-articulation is the dominant organization dimension of phoneme representations at the single-cell level, suggesting a remarkable consistency across levels of analyses, from the single neuron level to that of large neuronal populations and behavior.

Neural tracking of speech mental imagery during rhythmic inner counting

The subjective inner experience of mental imagery is among the most ubiquitous human experiences in daily life. Elucidating the neural implementation underpinning the dynamic construction of mental imagery is critical to understanding high-order cognitive function in the human brain. Here, we applied a frequency-tagging method to isolate the top-down process of speech mental imagery from bottom-up sensory-driven activities and concurrently tracked the neural processing time scales corresponding to the two processes in human subjects. Notably, by estimating the source of the magnetoencephalography (MEG) signals, we identified isolated brain networks activated at the imagery-rate frequency. In contrast, more extensive brain regions in the auditory temporal cortex were activated at the stimulus-rate frequency. Furthermore, intracranial stereotactic electroencephalogram (sEEG) evidence confirmed the participation of the inferior frontal gyrus in generating speech mental imagery. Our results indicate that a disassociated neural network underlies the dynamic construction of speech mental imagery independent of auditory perception.

Scale-Free Amplitude Modulation of Neuronal Oscillations Tracks Comprehension of Accelerated Speech

Speech comprehension is preserved up to a threefold acceleration, but deteriorates rapidly at higher speeds. Current models posit that perceptual resilience to accelerated speech is limited by the brain's ability to parse speech into syllabic units using δ/θ oscillations. Here, we investigated whether the involvement of neuronal oscillations in processing accelerated speech also relates to their scale-free amplitude modulation as indexed by the strength of long-range temporal correlations (LRTC). We recorded MEG while 24 human subjects (12 females) listened to radio news uttered at different comprehensible rates, at a mostly unintelligible rate and at this same speed interleaved with silence gaps. δ, θ, and low-γ oscillations followed the nonlinear variation of comprehension, with LRTC rising only at the highest speed. In contrast, increasing the rate was associated with a monotonic increase in LRTC in high-γ activity. When intelligibility was restored with the insertion of silence gaps, LRTC in the δ, θ, and low-γ oscillations resumed the low levels observed for intelligible speech. Remarkably, the lower the individual subject scaling exponents of δ/θ oscillations, the greater the comprehension of the fastest speech rate. Moreover, the strength of LRTC of the speech envelope decreased at the maximal rate, suggesting an inverse relationship with the LRTC of brain dynamics when comprehension halts. Our findings show that scale-free amplitude modulation of cortical oscillations and speech signals are tightly coupled to speech uptake capacity.SIGNIFICANCE STATEMENT One may read this statement in 20-30 s, but reading it in less than five leaves us clueless. Our minds limit how much information we grasp in an instant. Understanding the neural constraints on our capacity for sensory uptake is a fundamental question in neuroscience. Here, MEG was used to investigate neuronal activity while subjects listened to radio news played faster and faster until becoming unintelligible. We found that speech comprehension is related to the scale-free dynamics of δ and θ bands, whereas this property in high-γ fluctuations mirrors speech rate. We propose that successful speech processing imposes constraints on the self-organization of synchronous cell assemblies and their scale-free dynamics adjusts to the temporal properties of spoken language.

Sparse Spectro-Temporal Receptive Fields Based on Multi-Unit and High-Gamma Responses in Human Auditory Cortex

Spectro-Temporal Receptive Fields (STRFs) were estimated from both multi-unit sorted clusters and high-gamma power responses in human auditory cortex. Intracranial electrophysiological recordings were used to measure responses to a random chord sequence of Gammatone stimuli. Traditional methods for estimating STRFs from single-unit recordings, such as spike-triggered-averages, tend to be noisy and are less robust to other response signals such as local field potentials. We present an extension to recently advanced methods for estimating STRFs from generalized linear models (GLM). A new variant of regression using regularization that penalizes non-zero coefficients is described, which results in a sparse solution. The frequency-time structure of the STRF tends toward grouping in different areas of frequency-time and we demonstrate that group sparsity-inducing penalties applied to GLM estimates of STRFs reduces the background noise while preserving the complex internal structure. The contribution of local spiking activity to the high-gamma power signal was factored out of the STRF using the GLM method, and this contribution was significant in 85 percent of the cases. Although the GLM methods have been used to estimate STRFs in animals, this study examines the detailed structure directly from auditory cortex in the awake human brain. We used this approach to identify an abrupt change in the best frequency of estimated STRFs along posteromedial-to-anterolateral recording locations along the long axis of Heschl's gyrus. This change correlates well with a proposed transition from core to non-core auditory fields previously identified using the temporal response properties of Heschl's gyrus recordings elicited by click-train stimuli.

Decoding speech perception from single cell activity in humans

Deciphering the content of continuous speech is a challenging task performed daily by the human brain. Here, we tested whether activity of single cells in auditory cortex could be used to support such a task. We recorded neural activity from auditory cortex of two neurosurgical patients while presented with a short video segment containing speech. Population spiking activity (~20 cells per patient) allowed detection of word onset and decoding the identity of perceived words with significantly high accuracy levels. Oscillation phase of local field potentials (8-12Hz) also allowed decoding word identity although with lower accuracy levels. Our results provide evidence that the spiking activity of a relatively small population of cells in human primary auditory cortex contains significant information for classification of words in ongoing speech. Given previous evidence for overlapping neural representation during speech perception and production, this may have implications for developing brain-machine interfaces for patients with deficits in speech production.

Modulation of response patterns in human auditory cortex during a target detection task: an intracranial electrophysiology study

Selective attention enhances cortical activity representing an attended sound stream in human posterolateral superior temporal gyrus (PLST). It is unclear, however, what mechanisms are associated with a target detection task that necessitates sustained attention (vigilance) to a sound stream. We compared responses elicited by target and non-target sounds, and to sounds presented in a passive-listening paradigm. Subjects were neurosurgical patients undergoing invasive monitoring for medically refractory epilepsy. Stimuli were complex tones, band-limited noise bursts and speech syllables. High gamma cortical activity (70-150 Hz) was examined in all subjects using subdural grid electrodes implanted over PLST. Additionally, responses were measured from depth electrodes implanted within Heschl's gyrus (HG) in one subject. Responses to target sounds recorded from PLST were increased when compared to responses elicited by the same sounds when they were not-targets, and when they were presented during passive listening. Increases in high gamma activity to target sounds occurred during later portions (after 250 ms) of the response. These increases were related to the task and not to detailed stimulus characteristics. In contrast, earlier activity that did not vary across conditions did represent stimulus acoustic characteristics. Effects observed on PLST were not noted in HG. No consistent effects were noted in the averaged evoked potentials in either cortical region. We conclude that task dependence modulates later activity in PLST during vigilance. Later activity may represent feedback from higher cortical areas. Study of concurrently recorded activity from frontoparietal areas is necessary to further clarify task-related modulation of activity on PLST.

Temporal scaling of neural responses to compressed and dilated natural speech

Different brain areas integrate information over different timescales, and this capacity to accumulate information increases from early sensory areas to higher order perceptual and cognitive areas. It is currently unknown whether the timescale capacity of each brain area is fixed or whether it adaptively rescales depending on the rate at which information arrives from the world. Here, using functional MRI, we measured brain responses to an auditory narrative presented at different rates. We asked whether neural responses to slowed (speeded) versions of the narrative could be compressed (stretched) to match neural responses to the original narrative. Temporal rescaling was observed in early auditory regions (which accumulate information over short timescales) as well as linguistic and extra-linguistic brain areas (which can accumulate information over long timescales). The temporal rescaling phenomenon started to break down for stimuli presented at double speed, and intelligibility was also impaired for these stimuli. These data suggest that 1) the rate of neural information processing can be rescaled according to the rate of incoming information, both in early sensory regions as well as in higher order cortexes, and 2) the rescaling of neural dynamics is confined to a range of rates that match the range of behavioral performance.

Emergence of sensory patterns during sleep highlights differential dynamics of REM and non-REM sleep stages

Despite the profound reduction in conscious awareness associated with sleep, sensory cortex remains highly active during the different sleep stages, exhibiting complex interactions between different cortical sites. The potential functional significance of such spatial patterns and how they change between different sleep stages is presently unknown. In this electrocorticography study of human patients, we examined this question by studying spatial patterns of activity (broadband gamma power) that emerge during sleep (sleep patterns) and comparing them to the functional organization of sensory cortex that is activated by naturalistic stimuli during the awake state. Our results show a high correlation (p < 10(-4), permutation test) between the sleep spatial patterns and the functional organization found during wakefulness. Examining how the sleep patterns changed through the night highlighted a stage-specific difference, whereby the repertoire of such patterns was significantly larger during rapid eye movement (REM) sleep compared with non-REM stages. These results reveal that intricate spatial patterns of sensory functional organization emerge in a stage-specific manner during sleep.

Selectivity of audiovisual ECoG responses revealed under naturalistic stimuli in the human cortex

A fundamental debate in the study of cortical sensory systems concerns the scale of functional selectivity in cortical networks. Brain imaging studies have repeatedly demonstrated functional selectivity in entire cortical areas and networks using predetermined stimuli. However, it is not clear to what extent these networks are heterogeneous, i.e., whether the selectivity profiles in subregions within each sensory network show significant dissimilarity. Here, we studied local functional selectivity in the human cortex using naturalistic movie clips shown to 12 patients implanted with intracranial electrocorticography electrodes (590 in total), providing extensive cortical coverage. We examined the similarity of response profiles (40- to 80-Hz gamma-power modulations) across electrodes using a novel data driven approach without assuming any predefined category. Our results show that the functional selectivity of each highly responsive electrode was different from that of all other electrodes across the sensory cortex. Thus most responsive electrodes showed an activation profile that was unique in each patient and was similar to that of only 0.3% (1-2) of all other electrodes across all patients. Functional similarity between electrodes was linked to anatomical proximity. While in most electrodes the source of selectivity was complex, a small subset showed the well-documented selectivity to faces and actions. Our results indicate that the human sensory cortex is organized as a mosaic of functionally unique subregions in which each site manifests its own special response profile.

Encoding of ultrasonic vocalizations in the auditory cortex

One of the central tasks of the mammalian auditory system is to represent information about acoustic communicative signals, such as vocalizations. However, the neuronal computations underlying vocalization encoding in the central auditory system are poorly understood. To learn how the rat auditory cortex encodes information about conspecific vocalizations, we presented a library of natural and temporally transformed ultrasonic vocalizations (USVs) to awake rats while recording neural activity in the primary auditory cortex (A1) with chronically implanted multielectrode probes. Many neurons reliably and selectively responded to USVs. The response strength to USVs correlated strongly with the response strength to frequency-modulated (FM) sweeps and the FM rate tuning index, suggesting that related mechanisms generate responses to USVs as to FM sweeps. The response strength further correlated with the neuron's best frequency, with the strongest responses produced by neurons whose best frequency was in the ultrasonic frequency range. For responses of each neuron to each stimulus group, we fitted a novel predictive model: a reduced generalized linear-nonlinear model (GLNM) that takes the frequency modulation and single-tone amplitude as the only two input parameters. The GLNM accurately predicted neuronal responses to previously unheard USVs, and its prediction accuracy was higher than that of an analogous spectrogram-based linear-nonlinear model. The response strength of neurons and the model prediction accuracy were higher for original, rather than temporally transformed, vocalizations. These results indicate that A1 processes original USVs differentially than transformed USVs, indicating preference for temporal statistics of the original vocalizations.

The tracking of speech envelope in the human cortex

Humans are highly adept at processing speech. Recently, it has been shown that slow temporal information in speech (i.e., the envelope of speech) is critical for speech comprehension. Furthermore, it has been found that evoked electric potentials in human cortex are correlated with the speech envelope. However, it has been unclear whether this essential linguistic feature is encoded differentially in specific regions, or whether it is represented throughout the auditory system. To answer this question, we recorded neural data with high temporal resolution directly from the cortex while human subjects listened to a spoken story. We found that the gamma activity in human auditory cortex robustly tracks the speech envelope. The effect is so marked that it is observed during a single presentation of the spoken story to each subject. The effect is stronger in regions situated relatively early in the auditory pathway (belt areas) compared to other regions involved in speech processing, including the superior temporal gyrus (STG) and the posterior inferior frontal gyrus (Broca's region). To further distinguish whether speech envelope is encoded in the auditory system as a phonological (speech-related), or instead as a more general acoustic feature, we also probed the auditory system with a melodic stimulus. We found that belt areas track melody envelope weakly, and as the only region considered. Together, our data provide the first direct electrophysiological evidence that the envelope of speech is robustly tracked in non-primary auditory cortex (belt areas in particular), and suggest that the considered higher-order regions (STG and Broca's region) partake in a more abstract linguistic analysis.

Extracellular levels of lactate, but not oxygen, reflect sleep homeostasis in the rat cerebral cortex

STUDY OBJECTIVE

It is well established that brain metabolism is higher during wake and rapid eye movement (REM) sleep than in nonrapid eye movement (NREM) sleep. Most of the brain's energy is used to maintain neuronal firing and glutamatergic transmission. Recent evidence shows that cortical firing rates, extracellular glutamate levels, and markers of excitatory synaptic strength increase with time spent awake and decline throughout NREM sleep. These data imply that the metabolic cost of each behavioral state is not fixed but may reflect sleep-wake history, a possibility that is investigated in the current report.

DESIGN

Chronic (4d) electroencephalographic (EEG) recordings in the rat cerebral cortex were coupled with fixed-potential amperometry to monitor the extracellular concentration of oxygen ([oxy]) and lactate ([lac]) on a second-by-second basis across the spontaneous sleep-wake cycle and in response to sleep deprivation.

SETTING

Basic sleep research laboratory.

PATIENTS OR PARTICIPANTS

Wistar Kyoto (WKY) adult male rats.

INTERVENTIONS

N/A.

MEASUREMENTS AND RESULTS

Within 30-60 sec [lac] and [oxy] progressively increased during wake and REM sleep and declined during NREM sleep (n = 10 rats/metabolite), but with several differences. [Oxy], but not [lac], increased more during wake with high motor activity and/or elevated EEG high-frequency power. Meanwhile, only the NREM decline of [lac] reflected sleep pressure as measured by slow-wave activity, mirroring previous results for cortical glutamate.

CONCLUSIONS

The observed state-dependent changes in cortical [lac] and [oxy] are consistent with higher brain metabolism during waking and REM sleep in comparison with NREM sleep. Moreover, these data suggest that glycolytic activity, most likely through its link with glutamatergic transmission, reflects sleep homeostasis.

Mechanisms of action selection and timing in substantia nigra neurons

The timing of actions is critical for adaptive behavior. In this study we measured neural activity in the substantia nigra as mice learned to change their action duration to earn food rewards. We observed dramatic changes in single unit activity during learning: both dopaminergic and GABAergic neurons changed their activity in relation to behavior to reflect the learned instrumental contingency and the action duration. We found the emergence of "action-on" neurons that increased firing for the duration of the lever press and mirror-image "action-off" neurons that paused at the same time. This pattern is especially common among GABAergic neurons. The activity of many neurons also reflected confidence about the just completed action and the prospect of reward. Being correlated with the relative duration of the completed action, their activity could predict the likelihood of reward collection. Compared with the GABAergic neurons, the activity of dopaminergic neurons was more commonly modulated by the discriminative stimulus signaling the start of each trial, suggesting that their phasic activity reflected sensory salience rather than any reward prediction error found in previous work. In short, these results suggest that (1) nigral activity is highly plastic and modified by the learning of the instrumental contingency; (2) GABAergic output from the substantia nigra can simultaneously inhibit and disinhibit downstream structures, while the dopaminergic output also provide bidirectional modulation of the corticostriatal circuits; (3) dopaminergic and GABAergic neurons show similar task-related activity, although DA neurons are more responsive to the trial start signal.

Diversity of neural-hemodynamic relationships associated with differences in cortical processing during bilateral somatosensory activation in rats

The neural-hemodynamic relationships may vary depending on cortical processing patterns. To investigate how cortical hemodynamics reflects neural activity involving different cortical processing patterns, we delivered electrical stimulation pulses to rat hindpaws, unilaterally or bilaterally, and simultaneously measured electrophysiological (local field potential, LFP < 100 Hz; multiunit activity, MUA>300 Hz) and optical intrinsic signals associated with changes in cerebral blood volume (CBV). Unilateral stimulation evoked neural and optical signals in bilateral primary somatosensory cortices. Ipsilateral optical responses indicating an increased CBV exhibited a peak magnitude of ~30% and mediocaudal shifts relative to contralateral responses. Correlation analyses revealed different scale factors between contralateral and ipsilateral responses in LFP-MUA and LFP-CBV relationships. Bilateral stimulation at varying time intervals evoked hemodynamic responses that were strongly suppressed at 40-ms intervals. This suppression quantitatively reflected suppressed LFP responses to contralateral testing stimulation and not linear summation, with slowly fluctuating LFP responses to ipsilateral conditioning stimulation. Consequently, in the overall responses to bilateral stimulation, CBV-related responses were more linearly correlated with MUA than with LFPs. When extracting high-frequency components (>30 Hz) from LFPs, we found similar scale factors between contralateral and ipsilateral responses in LFP-MUA and LFP-CBV relationships, resulting in significant linear relationships among these components, MUA, and cortical hemodynamics in overall responses to bilateral stimulation. The dependence of LFP-MUA-hemodynamic relationships on cortical processing patterns and the LFP temporal/spectral structure is important for interpreting hemodynamic signals in complex functional paradigms driving diverse cortical processing.

Human intracranial recordings and cognitive neuroscience

The ultimate goal of neuroscience research is to understand the operating mechanism of the human brain and to exploit this understanding to devise methods for repair when it malfunctions. A key feature of this operating mechanism is electrical activity of single brain cells and cell assemblies. For obvious ethical reasons, scientists rely mostly on animal research in the study of such signals. Research in humans is often limited to electrical signals that can be recorded at the scalp or to surrogates of electrical activity, namely magnetic source imaging and measures of regional blood flow and metabolism. Invasive brain recordings performed in patients during various clinical procedures provide a unique opportunity to record high-resolution signals in vivo from the human brain-data that are otherwise unavailable. Of special value are the rare opportunities to record in awake humans the activity of single brain cells and small cellular assemblies. These recordings provide a unique view on aspects of human cognition that are impossible to study in animals, including language, imagery, episodic memory, volition, and even consciousness. In the current review we discuss the unique contribution of invasive recordings from patients to the field of cognitive neuroscience.