From intentions to actions: Neural oscillations encode motor processes through phase, amplitude and phase-amplitude coupling

Human local field potentials in motor and non-motor brain areas encode upcoming movement direction

Limb movement direction can be inferred from local field potentials in motor cortex during movement execution. Yet, it remains unclear to what extent intended hand movements can be predicted from brain activity recorded during movement planning. Here, we set out to probe the directional-tuning of oscillatory features during motor planning and execution, using a machine learning framework on multi-site local field potentials (LFPs) in humans. We recorded intracranial EEG data from implanted epilepsy patients as they performed a four-direction delayed center-out motor task. Fronto-parietal LFP low-frequency power predicted hand-movement direction during planning while execution was largely mediated by higher frequency power and low-frequency phase in motor areas. By contrast, Phase-Amplitude Coupling showed uniform modulations across directions. Finally, multivariate classification led to an increase in overall decoding accuracy (>80%). The novel insights revealed here extend our understanding of the role of neural oscillations in encoding motor plans.

Dissociable default-mode subnetworks subserve childhood attention and cognitive flexibility: Evidence from deep learning and stereotactic electroencephalography

Cognitive flexibility encompasses the ability to efficiently shift focus and forms a critical component of goal-directed attention. The neural substrates of this process are incompletely understood in part due to difficulties in sampling the involved circuitry. We leverage stereotactic intracranial recordings to directly resolve local-field potentials from otherwise inaccessible structures to study moment-to-moment attentional activity in children with epilepsy performing a flexible attentional task. On an individual subject level, we employed deep learning to decode neural features predictive of task performance indexed by single-trial reaction time. These models were subsequently aggregated across participants to identify predictive brain regions based on AAL atlas and FIND functional network parcellations. Through this approach, we show that fluctuations in beta (12-30 Hz) and gamma (30-80 Hz) power reflective of increased top-down attentional control and local neuronal processing within relevant large-scale networks can accurately predict single-trial task performance. We next performed connectomic profiling of these highly predictive nodes to examine task-related engagement of distributed functional networks, revealing exclusive recruitment of the dorsal default mode network during shifts in attention. The identification of distinct substreams within the default mode system supports a key role for this network in cognitive flexibility and attention in children. Furthermore, convergence of our results onto consistent functional networks despite significant inter-subject variability in electrode implantations supports a broader role for deep learning applied to intracranial electrodes in the study of human attention.

A Dynamical Systems Approach to Characterizing Brain-Body Interactions during Movement: Challenges, Interpretations, and Recommendations

Brain-body interactions (BBIs) have been the focus of intense scrutiny since the inception of the scientific method, playing a foundational role in the earliest debates over the philosophy of science. Contemporary investigations of BBIs to elucidate the neural principles of motor control have benefited from advances in neuroimaging, device engineering, and signal processing. However, these studies generally suffer from two major limitations. First, they rely on interpretations of 'brain' activity that are behavioral in nature, rather than neuroanatomical or biophysical. Second, they employ methodological approaches that are inconsistent with a dynamical systems approach to neuromotor control. These limitations represent a fundamental challenge to the use of BBIs for answering basic and applied research questions in neuroimaging and neurorehabilitation. Thus, this review is written as a tutorial to address both limitations for those interested in studying BBIs through a dynamical systems lens. First, we outline current best practices for acquiring, interpreting, and cleaning scalp-measured electroencephalography (EEG) acquired during whole-body movement. Second, we discuss historical and current theories for modeling EEG and kinematic data as dynamical systems. Third, we provide worked examples from both canonical model systems and from empirical EEG and kinematic data collected from two subjects during an overground walking task.

Effect of low-intensity transcranial ultrasound stimulation on theta and gamma oscillations in the mouse hippocampal CA1

Previous studies have demonstrated that low-intensity transcranial ultrasound stimulation (TUS) can eliminate hippocampal neural activity. However, until now, it has remained unclear how ultrasound modulates theta and gamma oscillations in the hippocampus under different behavioral states. In this study, we used ultrasound to stimulate the CA1 in mice in anesthesia, awake and running states, and we simultaneously recorded the local field potential of the stimulation location. We analyzed the power spectrum, phase-amplitude coupling (PAC) of theta and gamma oscillations, and their relationship with ultrasound intensity. The results showed that (i) TUS significantly enhanced the absolute power of theta and gamma oscillations under anesthesia and in the awake state. (ii) The PAC strength between theta and gamma oscillations is significantly enhanced under the anesthesia and awake states but is weakened under the running state with TUS. (iii) Under anesthesia, the relative power of theta decreases and that of gamma increases as ultrasound intensity increases, and the result under the awake state is opposite that under the anesthesia state. (iv) The PAC index between theta and gamma increases as ultrasound intensity increases under the anesthesia and awake states. The above results demonstrate that TUS can modulate theta and gamma oscillations in the CA1 and that the modulation effect depends on behavioral states. Our study provides guidance for the application of ultrasound in modulating hippocampal function.

Cortical theta-gamma coupling tracks the mental workload as an indicator of mental schema development during simulated quadrotor UAV operation

Objective. In the emerging field of neuroergonomics, mental workload assessment is one of the most important problems. Previous studies have made some progress on the relationship between task difficulties and mental workload, but how the mental schema, a reflection of the understanding and mastery degree of a task, affects mental workload has not been clearly discussed.Approach. There is emerging appreciation for the role of theta-gamma coupling (TGC) in high-level cognitive functions. Here, we attempt to further our understanding of how mental schema development and task difficulty had an impact on mental workload from the perspective of TGC. Specifically, the variation of TGC coupling strength and coupling pattern was estimated with different test orders and task difficulties performed by 51 students in a ten-day simulated quadrotor unmanned aerial vehicle flight training and test tasks.Main results. During the training, TGC increased with mental schema development. For the test tasks, TGC did not change with increasing task difficulty before the operator formed a mental schema but decreased with the increasing mental workload after the formation of the mental schema.Significance. Our results suggest that TGC was a robust indicator of mental schema development and could be biased by task difficulty. In conclusion, TGC can be a promising measure of mental workload, but only for experienced operators.

Comparison of event-related modulation index and traditional methods for evaluating phase-amplitude coupling using simulated brain signals

The investigation of brain oscillations and connectivity has become an important topic in the recent decade. There are several types of interactions between neuronal oscillations, and one of the most interesting among these interactions is phase-amplitude coupling (PAC). Several methods have been proposed to measure the strength of PAC, including the phase-locking value, circular-linear correlation, and modulation index. In the current study, we compared these traditional PAC methods with simulated electroencephalogram signals. Further, to assess the PAC value at each time point, we also compared two recently established methods, event-related phase-locking value and event-related circular-linear correlation, with our newly proposed event-related modulation index (ERMI). Results indicated that the ERMI has better temporal resolution and is more tolerant to noise than the other two event-related methods, suggesting the advantages of utilizing ERMI in evaluating the strength of PAC within a brain region.

Using linear parameter varying autoregressive models to measure cross frequency couplings in EEG signals

For years now, phase-amplitude cross frequency coupling (CFC) has been observed across multiple brain regions under different physiological and pathological conditions. It has been suggested that CFC serves as a mechanism that facilitates communication and information transfer between local and spatially separated neuronal populations. In non-invasive brain computer interfaces (BCI), CFC has not been thoroughly explored. In this work, we propose a CFC estimation method based on Linear Parameter Varying Autoregressive (LPV-AR) models and we assess its performance using both synthetic data and electroencephalographic (EEG) data recorded during attempted arm/hand movements of spinal cord injured (SCI) participants. Our results corroborate the potentiality of CFC as a feature for movement attempt decoding and provide evidence of the superiority of our proposed CFC estimation approach compared to other commonly used techniques.

The phase of sensorimotor mu and beta oscillations has the opposite effect on corticospinal excitability

BACKGROUND

Neural oscillations in the primary motor cortex (M1) shape corticospinal excitability. Power and phase of ongoing mu (8-13 Hz) and beta (14-30 Hz) activity may mediate motor cortical output. However, the functional dynamics of both mu and beta phase and power relationships and their interaction, are largely unknown.

OBJECTIVE

Here, we employ recently developed real-time targeting of the mu and beta rhythm, to apply phase-specific brain stimulation and probe motor corticospinal excitability non-invasively. For this, we used instantaneous read-out and analysis of ongoing oscillations, targeting four different phases (0°, 90°, 180°, and 270°) of mu and beta rhythms with suprathreshold single-pulse transcranial magnetic stimulation (TMS) to M1. Ensuing motor evoked potentials (MEPs) in the right first dorsal interossei muscle were recorded. Twenty healthy adults took part in this double-blind randomized crossover study.

RESULTS

Mixed model regression analyses showed significant phase-dependent modulation of corticospinal output by both mu and beta rhythm. Strikingly, these modulations exhibit a double dissociation. MEPs are larger at the mu trough and rising phase and smaller at the peak and falling phase. For the beta rhythm we found the opposite behavior. Also, mu power, but not beta power, was positively correlated with corticospinal output. Power and phase effects did not interact for either rhythm, suggesting independence between these aspects of oscillations.

CONCLUSION

Our results provide insights into real-time motor cortical oscillation dynamics, which offers the opportunity to improve the effectiveness of TMS by specifically targeting different frequency bands.

Advances in human intracranial electroencephalography research, guidelines and good practices

Since the second-half of the twentieth century, intracranial electroencephalography (iEEG), including both electrocorticography (ECoG) and stereo-electroencephalography (sEEG), has provided an intimate view into the human brain. At the interface between fundamental research and the clinic, iEEG provides both high temporal resolution and high spatial specificity but comes with constraints, such as the individual's tailored sparsity of electrode sampling. Over the years, researchers in neuroscience developed their practices to make the most of the iEEG approach. Here we offer a critical review of iEEG research practices in a didactic framework for newcomers, as well addressing issues encountered by proficient researchers. The scope is threefold: (i) review common practices in iEEG research, (ii) suggest potential guidelines for working with iEEG data and answer frequently asked questions based on the most widespread practices, and (iii) based on current neurophysiological knowledge and methodologies, pave the way to good practice standards in iEEG research. The organization of this paper follows the steps of iEEG data processing. The first section contextualizes iEEG data collection. The second section focuses on localization of intracranial electrodes. The third section highlights the main pre-processing steps. The fourth section presents iEEG signal analysis methods. The fifth section discusses statistical approaches. The sixth section draws some unique perspectives on iEEG research. Finally, to ensure a consistent nomenclature throughout the manuscript and to align with other guidelines, e.g., Brain Imaging Data Structure (BIDS) and the OHBM Committee on Best Practices in Data Analysis and Sharing (COBIDAS), we provide a glossary to disambiguate terms related to iEEG research.

Decoding Intracranial EEG With Machine Learning: A Systematic Review

Advances in intracranial electroencephalography (iEEG) and neurophysiology have enabled the study of previously inaccessible brain regions with high fidelity temporal and spatial resolution. Studies of iEEG have revealed a rich neural code subserving healthy brain function and which fails in disease states. Machine learning (ML), a form of artificial intelligence, is a modern tool that may be able to better decode complex neural signals and enhance interpretation of these data. To date, a number of publications have applied ML to iEEG, but clinician awareness of these techniques and their relevance to neurosurgery, has been limited. The present work presents a review of existing applications of ML techniques in iEEG data, discusses the relative merits and limitations of the various approaches, and examines potential avenues for clinical translation in neurosurgery. One-hundred-seven articles examining artificial intelligence applications to iEEG were identified from 3 databases. Clinical applications of ML from these articles were categorized into 4 domains: i) seizure analysis, ii) motor tasks, iii) cognitive assessment, and iv) sleep staging. The review revealed that supervised algorithms were most commonly used across studies and often leveraged publicly available timeseries datasets. We conclude with recommendations for future work and potential clinical applications.

Group-level inference of information-based measures for the analyses of cognitive brain networks from neurophysiological data

The reproducibility crisis in neuroimaging and in particular in the case of underpowered studies has introduced doubts on our ability to reproduce, replicate and generalize findings. As a response, we have seen the emergence of suggested guidelines and principles for neuroscientists known as Good Scientific Practice for conducting more reliable research. Still, every study remains almost unique in its combination of analytical and statistical approaches. While it is understandable considering the diversity of designs and brain data recording, it also represents a striking point against reproducibility. Here, we propose a non-parametric permutation-based statistical framework, primarily designed for neurophysiological data, in order to perform group-level inferences on non-negative measures of information encompassing metrics from information-theory, machine-learning or measures of distances. The framework supports both fixed- and random-effect models to adapt to inter-individuals and inter-sessions variability. Using numerical simulations, we compared the accuracy in ground-truth retrieving of both group models, such as test- and cluster-wise corrections for multiple comparisons. We then reproduced and extended existing results using both spatially uniform MEG and non-uniform intracranial neurophysiological data. We showed how the framework can be used to extract stereotypical task- and behavior-related effects across the population covering scales from the local level of brain regions, inter-areal functional connectivity to measures summarizing network properties. We also present an open-source Python toolbox called Frites¹ that includes the proposed statistical pipeline using information-theoretic metrics such as single-trial functional connectivity estimations for the extraction of cognitive brain networks. Taken together, we believe that this framework deserves careful attention as its robustness and flexibility could be the starting point toward the uniformization of statistical approaches.

Theta-gamma coupling as a cortical biomarker of brain-computer interface-mediated motor recovery in chronic stroke

Chronic stroke patients with upper-limb motor disabilities are now beginning to see treatment options that were not previously available. To date, the two options recently approved by the United States Food and Drug Administration include vagus nerve stimulation and brain-computer interface therapy. While the mechanisms for vagus nerve stimulation have been well defined, the mechanisms underlying brain-computer interface-driven motor rehabilitation are largely unknown. Given that cross-frequency coupling has been associated with a wide variety of higher-order functions involved in learning and memory, we hypothesized this rhythm-specific mechanism would correlate with the functional improvements effected by a brain-computer interface. This study investigated whether the motor improvements in chronic stroke patients induced with a brain-computer interface therapy are associated with alterations in phase-amplitude coupling, a type of cross-frequency coupling. Seventeen chronic hemiparetic stroke patients used a robotic hand orthosis controlled with contralesional motor cortical signals measured with EEG. Patients regularly performed a therapeutic brain-computer interface task for 12 weeks. Resting-state EEG recordings and motor function data were acquired before initiating brain-computer interface therapy and once every 4 weeks after the therapy. Changes in phase-amplitude coupling values were assessed and correlated with motor function improvements. To establish whether coupling between two different frequency bands was more functionally important than either of those rhythms alone, we calculated power spectra as well. We found that theta-gamma coupling was enhanced bilaterally at the motor areas and showed significant correlations across brain-computer interface therapy sessions. Importantly, an increase in theta-gamma coupling positively correlated with motor recovery over the course of rehabilitation. The sources of theta-gamma coupling increase following brain-computer interface therapy were mostly located in the hand regions of the primary motor cortex on the left and right cerebral hemispheres. Beta-gamma coupling decreased bilaterally at the frontal areas following the therapy, but these effects did not correlate with motor recovery. Alpha-gamma coupling was not altered by brain-computer interface therapy. Power spectra did not change significantly over the course of the brain-computer interface therapy. The significant functional improvement in chronic stroke patients induced by brain-computer interface therapy was strongly correlated with increased theta-gamma coupling in bihemispheric motor regions. These findings support the notion that specific cross-frequency coupling dynamics in the brain likely play a mechanistic role in mediating motor recovery in the chronic phase of stroke recovery.

Cross-frequency phase-amplitude coupling in repetitive movements in patients with Parkinson's disease

Bradykinesia is a cardinal motor symptom in Parkinson’s disease (PD), the pathophysiology of which is not fully understood. We analyzed the role of cross-frequency coupling of oscillatory cortical activity in motor impairment in patients with PD and healthy controls. High-density EEG signals were recorded during various motor activities and at rest. Patients performed a repetitive finger-pressing task normally, but were slower than controls during tapping. Phase-amplitude coupling (PAC) between β (13–30 Hz) and broadband γ (50–150 Hz) was computed from individual EEG source signals in the premotor, primary motor, and primary somatosensory cortices, and the primary somatosensory complex. In all four regions, averaging the entire movement period resulted in higher PAC in patients than in controls for the resting condition and the pressing task (similar performance between groups). However, this was not the case for the tapping tasks where patients performed slower. This suggests the strength of state-related β-γ PAC does not determine Parkinsonian bradykinesia. Examination of the dynamics of oscillatory EEG signals during motor transitions revealed a distinctive motif of PAC rise and decay around press onset. This pattern was also present at press offset and slow tapping onset, linking such idiosyncratic PAC changes to transitions between different movement states. The transition-related PAC modulation in patients was similar to controls in the pressing task but flattened during slow tapping, which related to normal and abnormal performance, respectively. These findings suggest that the dysfunctional evolution of neuronal population dynamics during movement execution is an important component of the pathophysiology of Parkinsonian bradykinesia.

NEW & NOTEWORTHY Our findings using noninvasive EEG recordings provide evidence that PAC dynamics might play a role in the physiological cortical control of movement execution and may encode transitions between movement states. Results in patients with Parkinson’s disease suggest that bradykinesia is related to a deficit of the dynamic regulation of PAC during movement execution rather than its absolute strength. Our findings may contribute to the development of a new concept of the pathophysiology of bradykinesia.

Unconscious Influences on "Free Will" Movement Initiation: Slow-wave Brain Stimulation and the Readiness Potential

A central objective in the study of volition has been to identify how changes in neural activity relate to voluntary-"free will"-movement. The readiness potential (RP) is observed in the EEG as a slow-building signal that precedes action onset. Many consider the RP as a marker of an underlying preparatory process for initiating voluntary movement. However, the RP may emerge from ongoing slow-wave brain oscillations that influence the timing of movement initiation in a phase-dependent manner. Transcranial alternating current stimulation (tACS) enables brain oscillations to be entrained at the frequency of stimulation. We delivered tACS at a slow-wave frequency over frontocentral motor areas while participants (n = 30) performed a simple, self-paced button press task. During the active tACS condition, participants showed a tendency to initiate actions in the phase of the tACS cycle that corresponded to increased negative potentials across the frontocentral motor region. Comparisons of premovement EEG activity observed over frontocentral and central scalp electrodes showed earlier onset and increased amplitude of RPs from active stimulation compared with sham stimulation. This suggests that movement-related activity in the brain can be modulated by the delivery of weak, nonconsciously perceptible alternating currents over frontocentral motor regions. We present novel findings that support existing theories, which suggest the timing of voluntary movement is influenced by the phase of slow-changing oscillating brain states.

Magnetoencephalography detects phase-amplitude coupling in Parkinson's disease

To characterize Parkinson’s disease, abnormal phase-amplitude coupling is assessed in the cortico-basal circuit using invasive recordings. It is unknown whether the same phenomenon might be found in regions other than the cortico-basal ganglia circuit. We hypothesized that using magnetoencephalography to assess phase-amplitude coupling in the whole brain can characterize Parkinson’s disease. We recorded resting-state magnetoencephalographic signals in patients with Parkinson’s disease and in healthy age- and sex-matched participants. We compared whole-brain signals from the two groups, evaluating the power spectra of 3 frequency bands (alpha, 8–12 Hz; beta, 13–25 Hz; gamma, 50–100 Hz) and the coupling between gamma amplitude and alpha or beta phases. Patients with Parkinson’s disease showed significant beta–gamma phase-amplitude coupling that was widely distributed in the sensorimotor, occipital, and temporal cortices; healthy participants showed such coupling only in parts of the somatosensory and temporal cortices. Moreover, beta- and gamma-band power differed significantly between participants in the two groups (P < 0.05). Finally, beta–gamma phase-amplitude coupling in the sensorimotor cortices correlated significantly with motor symptoms of Parkinson’s disease (P < 0.05); beta- and gamma-band power did not. We thus demonstrated that beta–gamma phase-amplitude coupling in the resting state characterizes Parkinson’s disease.

Amplitude of Low-Frequency Fluctuation With Different Clinical Outcomes in Patients With Generalized Tonic-Clonic Seizures

BACKGROUND

Generalized tonic-clonic seizures (GTCS) are associated with significant disability and sudden unexpected death when they cannot be controlled. We aimed to explore the underlying neural substrate of the different responses to antiseizure drugs between the seizure-free (SF) and non-seizure-free (NSF) patients with GTCS through the amplitude of low-frequency fluctuation (ALFF) method.

METHODS

We calculated ALFF among the SF group, NSF group, and healthy controls (HCs) by collecting resting-state functional magnetic resonance imaging (rs-fMRI) data. One-way ANOVA was used to compare the ALFF of the three groups, and post-hoc analysis was done at the same time. Pearson's correlation analysis between ALFF in the discrepant brain areas and the clinical characteristics (disease course and age of onset of GTCS) was calculated after then.

RESULTS

A significant group effect was found in the right fusiform gyrus (R.FG), left fusiform gyrus (L.FG), left middle occipital gyrus (L.MOG), right inferior frontal gyrus (R.IFG), right precentral gyrus (R.PreG), right postcentral gyrus (R.PostG), and left calcarine sulcus (L.CS). The SF and NSF groups both showed increased ALFF in all discrepant brain areas compared to HCs except the R.IFG in the NSF group. Significantly higher ALFF in the bilateral FG and lower ALFF in the R.IFG were found in the NSF group compared to the SF group.

CONCLUSIONS

Higher ALFF in the bilateral FG were found in the NSF group compared to the SF and HC groups. Our findings indicate that abnormal brain activity in the FG may be one potential neural substrate to interpret the failure of seizure control in patients with GTCS.

Brain oscillatory correlates of visuomotor adaptive learning

Sensorimotor adaptation involves the recalibration of the mapping between motor command and sensory feedback in response to movement errors. Although adaptation operates within individual movements on a trial-to-trial basis, it can also undergo learning when adaptive responses improve over the course of many trials. Brain oscillatory activities related to these "adaptation" and "learning" processes remain unclear. The main reason for this is that previous studies principally focused on the beta band, which confined the outcome message to trial-to-trial adaptation. To provide a wider understanding of adaptive learning, we decoded visuomotor tasks with constant, random or no perturbation from EEG recordings in different bandwidths and brain regions using a multiple kernel learning approach. These different experimental tasks were intended to separate trial-to-trial adaptation from the formation of the new visuomotor mapping across trials. We found changes in EEG power in the post-movement period during the course of the visuomotor-constant rotation task, in particular an increased (i) theta power in prefrontal region, (ii) beta power in supplementary motor area, and (iii) gamma power in motor regions. Classifying the visuomotor task with constant rotation versus those with random or no rotation, we were able to relate power changes in beta band mainly to trial-to-trial adaptation to error while changes in theta band would relate rather to the learning of the new mapping. Altogether, this suggested that there is a tight relationship between modulation of the synchronization of low (theta) and higher (essentially beta) frequency oscillations in prefrontal and sensorimotor regions, respectively, and adaptive learning.

Behavioral and Neural Dynamics of Interpersonal Synchrony Between Performing Musicians: A Wireless EEG Hyperscanning Study

Interpersonal synchrony refers to the temporal coordination of actions between individuals and is a common feature of social behaviors, from team sport to ensemble music performance. Interpersonal synchrony of many rhythmic (periodic) behaviors displays dynamics of coupled biological oscillators. The current study addresses oscillatory dynamics on the levels of brain and behavior between music duet partners performing at spontaneous (uncued) rates. Wireless EEG was measured from N = 20 pairs of pianists as they performed a melody first in Solo performance (at their spontaneous rate of performance), and then in Duet performances at each partner's spontaneous rate. Influences of partners' spontaneous rates on interpersonal synchrony were assessed by correlating differences in partners' spontaneous rates of Solo performance with Duet tone onset asynchronies. Coupling between partners' neural oscillations was assessed by correlating amplitude envelope fluctuations of cortical oscillations at the Duet performance frequency between observed partners and between surrogate (re-paired) partners, who performed the same melody but at different times. Duet synchronization was influenced by partners' spontaneous rates in Solo performance. The size and direction of the difference in partners' spontaneous rates were mirrored in the size and direction of the Duet asynchronies. Moreover, observed Duet partners showed greater inter-brain correlations of oscillatory amplitude fluctuations than did surrogate partners, suggesting that performing in synchrony with a musical partner is reflected in coupled cortical dynamics at the performance frequency. The current study provides evidence that dynamics of oscillator coupling are reflected in both behavioral and neural measures of temporal coordination during musical joint action.

Violation of rhythmic expectancies can elicit late frontal gamma activity nested in theta oscillations

Rhythm processing involves building expectations according to the hierarchical temporal structure of auditory events. Although rhythm processing has been addressed in the context of predictive coding, the properties of the oscillatory response in different cortical areas are still not clear. We explored the oscillatory properties of the neural response to rhythmic incongruence and the cross-frequency coupling between multiple frequencies to further investigate the mechanisms underlying rhythm perception. We designed an experiment to investigate the neural response to rhythmic deviations in which the tone either arrived earlier than expected or the tone in the same metrical position was omitted. These two manipulations modulate the rhythmic structure differently, with the former creating a larger violation of the general structure of the musical stimulus than the latter. Both deviations resulted in an MMN response, whereas only the rhythmic deviant resulted in a subsequent P3a. Rhythmic deviants due to the early occurrence of a tone, but not omission deviants, seemed to elicit a late high gamma response (60-80 Hz) at the end of the P3a over the left frontal region, which, interestingly, correlated with the P3a amplitude over the same region and was also nested in theta oscillations. The timing of the elicited high-frequency gamma oscillations related to rhythmic deviation suggests that it might be related to the update of the predictive neural model, corresponding to the temporal structure of the events in higher-level cortical areas.

Riding the slow wave: Exploring the role of entrained low-frequency oscillations in memory formation

Neural oscillations are proposed to support a variety of behaviors, including long-term memory, yet their functional significance remains an active area of research. Here, we explore a potential functional role of low-frequency cortical oscillations in episodic memory formation. Recent theories suggest that low-frequency oscillations orchestrate rhythmic attentional sampling of the environment by dynamically modulating neural excitability across time. When these oscillations entrain to low-frequency rhythms present in the environment, such as speech or music, the brain can build temporal predictions about the onset of relevant events so that these events can be more efficiently processed. Building upon this literature, we propose that entrained low-frequency oscillations may similarly influence the temporal dynamics of episodic memory by rhythmically modulating encoding across time (mnemonic sampling). Central to this proposal is the phenomenon of cross-frequency phase-amplitude coupling, whereby the amplitudes of faster (higher frequency) rhythms, such as gamma oscillations, couple to the phase of slower (lower-frequency) rhythms entrained to environmental stimuli. By imposing temporal structure on higher-frequency oscillatory activity previously linked to memory formation, entrained low-frequency oscillations could dynamically orchestrate memory formation and optimize encoding at specific moments in time. We discuss prior experimental and theoretical work relevant to this proposal.

Alpha and high gamma phase amplitude coupling during motor imagery and weighted cross-frequency coupling to extract discriminative cross-frequency patterns

Motor imagery modulates specific neural oscillations like actual movement does. Representatively, suppression of the alpha power (e.g., event-related desynchronization [ERD]) is the typical pattern of motor imagery in the motor cortex. However, in addition to this amplitude-based feature, the coupling across frequencies includes important information about the brain functions and the existence of such complex information has been reported in various invasive studies. Yet, the interaction across multiple frequencies during motor imagery processing is still unclear and has not been widely studied, particularly concerning the non-invasive signals. In this study, we provide empirical evidence of the comodulation between the phase of alpha rhythm and the amplitude of high gamma rhythm during the motor imagery process. We used electroencephalography (EEG) in our investigation during the imagination of left- or right-hand movement recorded from 52 healthy subjects, and quantified the ERD of alpha and phase-amplitude coupling (PAC) which is a relative change of modulation index to the base line period (before the cue). As a result, we found that the coupling between the phase of alpha (8-12 Hz) and the amplitude of high gamma (70-120 Hz) and this PAC decreases during motor imagery and then rebounds to the baseline like alpha ERD (r = 0.29 to 0.42). This correlation between PAC and ERD was particularly stronger in the ipsilateral area. In addition, trials that demonstrated higher alpha power during the ready period (before the cue) showed a larger ERD during motor imagery and similarly, trials with higher modulation index during the ready period yielded a greater decrease in PAC during imagery. In the classification analysis, we found that the effective phase frequency that showed better decoding accuracy in left and right-hand imagery, varied across subjects. Motivated by result, we proposed a weighted cross-frequency coupling (WCFC) method that extracts the maximal discriminative feature by combining band power and CFC. In the evaluation, WCFC with only two electrodes yielded a performance comparable to the conventional algorithm with 64 electrodes in classifying left and right-hand motor imagery. These results indicate that the phase-amplitude frequency plays an important role in motor imagery, and that optimizing this frequency ranges is crucial for extracting information features to decode the motor imagery types.

Neural oscillations track the maintenance and proceduralization of novel instructions

Humans are capable of flexibly converting symbolic instructions into novel behaviors. Previous evidence and theoretical models suggest that the implementation of a novel instruction requires the reformatting of its declarative content into an action-oriented code optimized for the execution of the instructed behavior. While neuroimaging research focused on identifying the brain areas involved in such a process, the temporal and electrophysiological mechanisms remain poorly understood. These mechanisms, however, can provide information about the specific cognitive processes that characterize the proceduralization of information. In the present study, we recorded EEG activity while we asked participants to either simply maintain declaratively the content of novel S-R mappings or to proactively prepare for their implementation. By means of time-frequency analyses, we isolated the oscillatory features specific to the proceduralization of instructions. Implementation of the instructed mappings elicited stronger theta activity over frontal electrodes and suppression in mu and beta activity over central electrodes. On the contrary, activity in the alpha band, which has been shown to track the attentional deployment to task-relevant items, showed no differences between tasks. Together, these results support the idea that proceduralization of information is characterized by specific component processes such as orchestrating complex task settings and configuring the motor system that are not observed when instructions are held in a declarative format.

Decoding the neural dynamics of free choice in humans

How do we choose a particular action among equally valid alternatives? Nonhuman primate findings have shown that decision-making implicates modulations in unit firing rates and local field potentials (LFPs) across frontal and parietal cortices. Yet the electrophysiological brain mechanisms that underlie free choice in humans remain ill defined. Here, we address this question using rare intracerebral electroencephalography (EEG) recordings in surgical epilepsy patients performing a delayed oculomotor decision task. We find that the temporal dynamics of high-gamma (HG, 60-140 Hz) neural activity in distinct frontal and parietal brain areas robustly discriminate free choice from instructed saccade planning at the level of single trials. Classification analysis was applied to the LFP signals to isolate decision-related activity from sensory and motor planning processes. Compared with instructed saccades, free-choice trials exhibited delayed and longer-lasting HG activity during the delay period. The temporal dynamics of the decision-specific sustained HG activity indexed the unfolding of a deliberation process, rather than memory maintenance. Taken together, these findings provide the first direct electrophysiological evidence in humans for the role of sustained high-frequency neural activation in frontoparietal cortex in mediating the intrinsically driven process of freely choosing among competing behavioral alternatives.

Characterization of EEG Data Revealing Relationships With Cognitive and Motor Symptoms in Parkinson's Disease: A Systematic Review

Recent research regards the electroencephalogram (EEG) as a promising method to study real-time brain dynamic changes in patients with Parkinson's disease (PD), but a deeper understanding is needed to discern coincident pathophysiology, patterns of changes, and diagnosis. This review summarized recent research on EEG characterization related to the cognitive and motor functions in PD patients and discussed its potential to be used as diagnostic biomarkers. Thirty papers out of 220 published from 2010 to 2020 were reviewed. Movement abnormalities and cognitive decline are related to changes in EEG spectrum and event-related potentials (ERPs) during typical oddball paradigms and/or combined motor tasks. Abnormalities in β and δ frequency bands are, respectively the main manifestation of dyskinesia and cognitive decline in PD. The review showed that PD patients have noteworthy changes in specific EEG characterizations, however, the underlying mechanism of the interrelation between gait and cognitive is still unclear. Understanding the specific nature of the relationship is essential for development of novel invasive clinical diagnostic and therapeutic methods.

The intimate relationship between coalescent generators in very premature human newborn brains: Quantifying the coupling of nested endogenous oscillations

Temporal theta slow-wave activity (TTA-SW) in premature infants is a specific neurobiomarker of the early neurodevelopment of perisylvian networks observed as early as 24 weeks of gestational age (wGA). It is present at the turning point between non-sensory driven spontaneous networks and cortical network functioning. Despite its clinical importance, the underlying mechanisms responsible for this spontaneous nested activity and its functional role have not yet been determined. The coupling between neural oscillations at different timescales is a key feature of ongoing neural activity, the characteristics of which are determined by the network structure and dynamics. The underlying mechanisms of cross-frequency coupling (CFC) are associated with several putative functions in adults. In order to show that this generic mechanism is already in place early in the course of development, we analyzed electroencephalography recordings from sleeping preterm newborns (24-27 wGA). Employing cross-frequency phase-amplitude coupling analyses, we found that TTAs were orchestrated by the SWs defined by a precise temporal relationship. Notably, TTAs were synchronized to the SW trough, and were suppressed during the SW peak. Spontaneous endogenous TTA-SWs constitute one of the very early signatures of the developing temporal neural networks with key functions, such as language and communication. The presence of a fine-tuned relationship between the slow activity and the TTA in premature neonates emphasizes the complexity and relative maturity of the intimate mechanisms that shape the CFC, the disruption of which can have severe neurodevelopmental consequences.

Spectral Power in Marmoset Frontal Motor Cortex during Natural Locomotor Behavior

During primate arboreal locomotion, substrate orientation modifies body axis orientation and biomechanical contribution of fore- and hindlimbs. To characterize the role of cortical oscillations in integrating these locomotor demands, we recorded electrocorticographic activity from left dorsal premotor, primary motor, and supplementary motor cortices of three common marmosets moving across a branch-like small-diameter pole, fixed horizontally or vertically. Animals displayed behavioral adjustments to the task, namely, the horizontal condition mainly induced quadrupedal walk with pronated/neutral forelimb postures, whereas the vertical condition induced walk and bound gaits with supinated/neutral postures. Examination of cortical activity suggests that β (16–35 Hz) and γ (75–100 Hz) oscillations could reflect different processes in locomotor adjustments. During task, modulation of γ ERS by substrate orientation (horizontal/vertical) and epoch (preparation/execution) suggests close tuning to movement dynamics and biomechanical demands. β ERD was essentially modulated by gait (walk/bound), which could illustrate contribution to movement sequence and coordination. At rest, modulation of β power by substrate orientation underlines its role in sensorimotor processes for postural maintenance.

Tensorpac: An open-source Python toolbox for tensor-based phase-amplitude coupling measurement in electrophysiological brain signals

Despite being the focus of a thriving field of research, the biological mechanisms that underlie information integration in the brain are not yet fully understood. A theory that has gained a lot of traction in recent years suggests that multi-scale integration is regulated by a hierarchy of mutually interacting neural oscillations. In particular, there is accumulating evidence that phase-amplitude coupling (PAC), a specific form of cross-frequency interaction, plays a key role in numerous cognitive processes. Current research in the field is not only hampered by the absence of a gold standard for PAC analysis, but also by the computational costs of running exhaustive computations on large and high-dimensional electrophysiological brain signals. In addition, various signal properties and analyses parameters can lead to spurious PAC. Here, we present Tensorpac, an open-source Python toolbox dedicated to PAC analysis of neurophysiological data. The advantages of Tensorpac include (1) higher computational efficiency thanks to software design that combines tensor computations and parallel computing, (2) the implementation of all most widely used PAC methods in one package, (3) the statistical analysis of PAC measures, and (4) extended PAC visualization capabilities. Tensorpac is distributed under a BSD-3-Clause license and can be launched on any operating system (Linux, OSX and Windows). It can be installed directly via pip or downloaded from Github (https://github.com/EtienneCmb/tensorpac). By making Tensorpac available, we aim to enhance the reproducibility and quality of PAC research, and provide open tools that will accelerate future method development in neuroscience.

Reactivation of Motor-Related Gamma Activity in Human NREM Sleep

Models of memory consolidation posit a central role for reactivation of brain activity patterns during sleep, especially in non-Rapid Eye Movement (NREM) sleep. While such "replay" of recent waking experiences has been well-demonstrated in rodents, electrophysiological evidence of reactivation in human sleep is still largely lacking. In this intracranial study in patients with epilepsy (N = 9) we explored the spontaneous electroencephalographic reactivation during sleep of spatial patterns of brain activity evoked by motor learning. We first extracted the gamma-band (60-140 Hz) patterns underlying finger movements during a tapping task and underlying no-movement during a short rest period just prior to the task, and trained a binary classifier to discriminate between motor movements vs. rest. We then used the trained model on NREM sleep data immediately after the task and on NREM sleep during a control sleep period preceding the task. Compared with the control sleep period, we found, at the subject level, an increase in the detection rate of motor-related patterns during sleep following the task, but without association with performance changes. These data provide electrophysiological support for the reoccurrence in NREM sleep of the neural activity related to previous waking experience, i.e. that a basic tenet of the reactivation theory does occur in human sleep.

Induced Gamma-Band Activity during Actual and Imaginary Movements: EEG Analysis

The purpose of this paper is to record and analyze induced gamma-band activity (GBA) (30-60 Hz) in cerebral motor areas during imaginary movement and to compare it quantitatively with activity recorded in the same areas during actual movement using a simplified electroencephalogram (EEG). Brain activity (basal activity, imaginary motor task and actual motor task) is obtained from 12 healthy volunteer subjects using an EEG (Cz channel). GBA is analyzed using the mean power spectral density (PSD) value. Event-related synchronization (ERS) is calculated from the PSD values of the basal GBA (GBAb), the GBA of the imaginary movement (GBAim) and the GBA of the actual movement (GBAac). The mean GBAim and GBAac values for the right and left hands are significantly higher than the GBAb value (p = 0.007). No significant difference is detected between mean GBA values during the imaginary and actual movement (p = 0.242). The mean ERS values for the imaginary movement (ERSimM (%) = 23.52) and for the actual movement (ERSacM = 27.47) do not present any significant difference (p = 0.117). We demonstrated that ERS could provide a useful way of indirectly checking the function of neuronal motor circuits activated by voluntary movement, both imaginary and actual. These results, as a proof of concept, could be applied to physiology studies, brain-computer interfaces, and diagnosis of cognitive or motor pathologies.

Cortical Oscillations during Gait: Wouldn't Walking be so Automatic?

Gait is often considered as an automatic movement but cortical control seems necessary to adapt gait pattern with environmental constraints. In order to study cortical activity during real locomotion, electroencephalography (EEG) appears to be particularly appropriate. It is now possible to record changes in cortical neural synchronization/desynchronization during gait. Studying gait initiation is also of particular interest because it implies motor and cognitive cortical control to adequately perform a step. Time-frequency analysis enables to study induced changes in EEG activity in different frequency bands. Such analysis reflects cortical activity implied in stabilized gait control but also in more challenging tasks (obstacle crossing, changes in speed, dual tasks…). These spectral patterns are directly influenced by the walking context but, when analyzing gait with a more demanding attentional task, cortical areas other than the sensorimotor cortex (prefrontal, posterior parietal cortex, etc.) seem specifically implied. While the muscular activity of legs and cortical activity are coupled, the precise role of the motor cortex to control the level of muscular contraction according to the gait task remains debated. The decoding of this brain activity is a necessary step to build valid brain-computer interfaces able to generate gait artificially.

Phase-Amplitude Coupling of Neural Oscillations Can Be Effectively Probed with Concurrent TMS-EEG

Despite the widespread use of transcranial magnetic stimulation (TMS), knowledge of its neurophysiological mode of action is still incomplete. Recently, TMS has been proposed to synchronise neural oscillators and to thereby increase the detectability of corresponding oscillations at the population level. As oscillations in the human brain are known to interact within nested hierarchies via phase-amplitude coupling, TMS might also be able to increase the macroscopic detectability of such coupling. In a concurrent TMS-electroencephalography study, we therefore examined the technique's influence on theta-gamma, alpha-gamma, and beta-gamma phase-amplitude coupling by delivering single-pulse TMS (sTMS) and repetitive TMS (rTMS) over the left motor cortex and right visual cortex of healthy participants. The rTMS pulse trains were of 5 Hz, 11 Hz, and 23 Hz for the three coupling variations, respectively. Relative to sham stimulation, all conditions showed transient but significant increases in phase-amplitude coupling at the stimulation site. In addition, we observed enhanced coupling over various other cortical sites, with a more extensive propagation during rTMS than during sTMS. By indicating that scalp-recorded phase-amplitude coupling can be effectively probed with TMS, these findings open the door to the technique's application in manipulative dissections of such coupling during human cognition and behaviour in healthy and pathological conditions.

Visbrain: A Multi-Purpose GPU-Accelerated Open-Source Suite for Multimodal Brain Data Visualization

We present Visbrain, a Python open-source package that offers a comprehensive visualization suite for neuroimaging and electrophysiological brain data. Visbrain consists of two levels of abstraction: (1) objects which represent highly configurable neuro-oriented visual primitives (3D brain, sources connectivity, etc.) and (2) graphical user interfaces for higher level interactions. The object level offers flexible and modular tools to produce and automate the production of figures using an approach similar to that of Matplotlib with subplots. The second level visually connects these objects by controlling properties and interactions through graphical interfaces. The current release of Visbrain (version 0.4.2) contains 14 different objects and three responsive graphical user interfaces, built with PyQt: Signal, for the inspection of time-series and spectral properties, Brain for any type of visualization involving a 3D brain and Sleep for polysomnographic data visualization and sleep analysis. Each module has been developed in tight collaboration with end-users, i.e., primarily neuroscientists and domain experts, who bring their experience to make Visbrain as transparent as possible to the recording modalities (e.g., intracranial EEG, scalp-EEG, MEG, anatomical and functional MRI). Visbrain is developed on top of VisPy, a Python package providing high-performance 2D and 3D visualization by leveraging the computational power of the graphics card. Visbrain is available on Github and comes with a documentation, examples, and datasets (http://visbrain.org).

The Strength of Alpha-Beta Oscillatory Coupling Predicts Motor Timing Precision

Precise timing makes the difference between harmony and cacophony, but how the brain achieves precision during timing is unknown. In this study, human participants (7 females, 5 males) generated a time interval while being recorded with magnetoencephalography. Building on the proposal that the coupling of neural oscillations provides a temporal code for information processing in the brain, we tested whether the strength of oscillatory coupling was sensitive to self-generated temporal precision. On a per individual basis, we show the presence of alpha-beta phase-amplitude coupling whose strength was associated with the temporal precision of self-generated time intervals, not with their absolute duration. Our results provide evidence that active oscillatory coupling engages α oscillations in maintaining the precision of an endogenous temporal motor goal encoded in β power; the when of self-timed actions. We propose that oscillatory coupling indexes the variance of neuronal computations, which translates into the precision of an individual's behavioral performance.SIGNIFICANCE STATEMENT Which neural mechanisms enable precise volitional timing in the brain is unknown, yet accurate and precise timing is essential in every realm of life. In this study, we build on the hypothesis that neural oscillations, and their coupling across time scales, are essential for the coding and for the transmission of information in the brain. We show the presence of alpha-beta phase-amplitude coupling (α-β PAC) whose strength was associated with the temporal precision of self-generated time intervals, not with their absolute duration. α-β PAC indexes the temporal precision with which information is represented in an individual's brain. Our results link large-scale neuronal variability on the one hand, and individuals' timing precision, on the other.

Synchronised spiking activity underlies phase amplitude coupling in the subthalamic nucleus of Parkinson's disease patients

Both phase-amplitude coupling (PAC) and beta-bursts in the subthalamic nucleus have been significantly linked to symptom severity in Parkinson's disease (PD) in humans and emerged independently as competing biomarkers for closed-loop deep brain stimulation (DBS). However, the underlying nature of subthalamic PAC is poorly understood and its relationship with transient beta burst-events has not been investigated. To address this, we studied macro- and micro electrode recordings of local field potentials (LFPs) and single unit activity from 15 hemispheres in 10 PD patients undergoing DBS surgery. PAC between beta phase and high frequency oscillation (HFO) amplitude was compared to single unit firing rates, spike triggered averages, power spectral densities, inter spike intervals and phase-spike locking, and was studied in periods of beta-bursting. We found a significant synchronisation of spiking to HFOs and correlation of mean firing rates with HFO-amplitude when the latter was coupled to beta phase (i.e. in the presence of PAC). In the presence of PAC, single unit power spectra displayed peaks in the beta and HFO frequency range and the HFO frequency was correlated with that in the LFP. Furthermore, inter spike interval frequencies peaked in the same frequencies for which PAC was observed. Finally, PAC significantly increased with beta burst-duration. Our findings offer new insight in the pathology of Parkinson's disease by providing evidence that subthalamic PAC reflects the locking of spiking activity to network beta oscillations and that this coupling progressively increases with beta-burst duration. These findings suggest that beta-bursts capture periods of increased subthalamic input/output synchronisation in the beta frequency range and have important implications for therapeutic closed-loop DBS.

Neural Oscillations and the Initiation of Voluntary Movement

The brain processes involved in the planning and initiation of voluntary action are of great interest for understanding the relationship between conscious awareness of decisions and the neural control of movement. Voluntary motor behavior has generally been considered to occur when conscious decisions trigger movements. However, several studies now provide compelling evidence that brain states indicative of forthcoming movements take place before a person becomes aware of a conscious decision to act. While such studies have created much debate over the nature of ‘free will,’ at the very least they suggest that unconscious brain processes are predictive of forthcoming movements. Recent studies suggest that slow changes in neuroelectric potentials may play a role in the timing of movement onset by pushing brain activity above a threshold to trigger the initiation of action. Indeed, recent studies have shown relationships between the phase of low frequency oscillatory activity of the brain and the onset of voluntary action. Such studies, however, cannot determine whether this underlying neural activity plays a causal role in the initiation of movement or is only associated with the intentional behavior. Non-invasive transcranial alternating current brain stimulation can entrain neural activity at particular frequencies in order to assess whether underlying brain processes are causally related to associated behaviors. In this review, we examine the evidence for neural coding of action as well as the brain states prior to action initiation and discuss whether low frequency alternating current brain stimulation could influence the timing of a persons’ decision to act.

Mapping Brain Activity with Electrocorticography: Resolution Properties and Robustness of Inverse Solutions

Electrocorticography (ECoG) is an electrophysiological technique that records brain activity directly from the cortical surface with high temporal (ms) and spatial (mm) resolution. Its major limitations are in the high invasiveness and in the restricted field-of-view of the electrode grid, which partially covers the cortex. To infer brain activity at locations different from just below the electrodes, it is necessary to solve the electromagnetic inverse problem. Limitations in the performance of source reconstruction algorithms from ECoG have been, to date, only partially addressed in the literature, and a systematic evaluation is still lacking. The main goal of this study is to provide a quantitative evaluation of resolution properties of widely used inverse methods (eLORETA and MNE) for various ECoG grid sizes, in terms of localization error, spatial dispersion, and overall amplitude. Additionally, this study aims at evaluating how the use of simultaneous electroencephalography (EEG) affects the above properties. For these purposes, we take advantage of a unique dataset in which a monkey underwent a simultaneous recording with a 128 channel ECoG grid and an 18 channel EEG grid. Our results show that, in general conditions, the reconstruction of cortical activity located more than 1 cm away from the ECoG grid is not accurate, since the localization error increases linearly with the distance from the electrodes. This problem can be partially overcome by recording simultaneously ECoG and EEG. However, this analysis enlightens the necessity to design inverse algorithms specifically targeted at taking into account the limited field-of-view of the ECoG grid.

Sleep: An Open-Source Python Software for Visualization, Analysis, and Staging of Sleep Data

We introduce Sleep, a new Python open-source graphical user interface (GUI) dedicated to visualization, scoring and analyses of sleep data. Among its most prominent features are: (1) Dynamic display of polysomnographic data, spectrogram, hypnogram and topographic maps with several customizable parameters, (2) Implementation of several automatic detection of sleep features such as spindles, K-complexes, slow waves, and rapid eye movements (REM), (3) Implementation of practical signal processing tools such as re-referencing or filtering, and (4) Display of main descriptive statistics including publication-ready tables and figures. The software package supports loading and reading raw EEG data from standard file formats such as European Data Format, in addition to a range of commercial data formats. Most importantly, Sleep is built on top of the VisPy library, which provides GPU-based fast and high-level visualization. As a result, it is capable of efficiently handling and displaying large sleep datasets. Sleep is freely available (http://visbrain.org/sleep) and comes with sample datasets and an extensive documentation. Novel functionalities will continue to be added and open-science community efforts are expected to enhance the capacities of this module.

Propofol, Sevoflurane, and Ketamine Induce a Reversible Increase in Delta-Gamma and Theta-Gamma Phase-Amplitude Coupling in Frontal Cortex of Rat

Studies from human and non-human species have demonstrated a breakdown of functional corticocortical connectivity during general anesthesia induced by anesthetics with diverse molecular, neurophysiological, and pharmacological profiles. Recent studies have demonstrated that changes in long-range neural communication, and by corollary, functional connectivity, might be influenced by cross-frequency coupling (CFC) between the phase of slow oscillations and the amplitude of local fast oscillations. Phase-amplitude coupling (PAC) between slow oscillations and alpha rhythm during general anesthesia reveal distinct patterns depending on the anesthetic. In this study, we analyzed the effect of three clinically used anesthetics (propofol: n = 6, sevoflurane: n = 10, and ketamine: n = 8) with distinct molecular mechanisms on changes in PAC in the frontal cortex of rat. The loss of righting reflex was used as a surrogate for unconsciousness. PAC was calculated using the modulation index (MI) algorithm between delta (1–4 Hz), theta (4–10 Hz), low gamma (25–55 Hz), and high gamma (65–125 Hz) bands. A linear mixed model with fixed effects was used for statistical comparisons between waking, anesthetized, and post-anesthesia recovery epochs. All three anesthetics increased the coupling between delta and low gamma (p < 0.0001) as well as between theta and low gamma (p < 0.0001) oscillations, which returned to baseline waking levels during the post-anesthetic recovery period. In addition, a reversible reduction in high gamma power (p < 0.0001) was a consistent change during anesthesia induced by all three agents. The changes in delta-high gamma and theta-high gamma PAC as well as power spectral changes in delta, theta, and low gamma bandwidths did not show a uniform response across the three anesthetics. These results encourage the study of alternative PAC patterns as drug-invariant markers of general anesthesia in humans.