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

Fast ripple band high-frequency activity associated with thalamic sleep spindles in pediatric epilepsy

OBJECTIVE

To investigate high-frequency activities (HFA) associated with thalamic sleep spindles.

METHODS

We studied a cohort of ten pediatric patients with medication resistant epilepsy who were identified as potential candidates for thalamic neuromodulation. These patients had thalamic sampling as well as presumed epileptogenic zones, using stereotactic EEG (SEEG) with a sampling frequency of 2,000 Hz. We quantified the summated high-frequency activity (HFA) in the fast ripple band associated with sleep spindles using 20-minute scalp EEG and SEEG recordings during non-REM sleep and analyzed its correlation with spindle characteristics.

RESULTS

HFA, with a median peak frequency of 330 Hz, was distinctively observed in the thalamus and temporally correlated with thalamic sleep spindles. Such HFA demonstrated significant coupling with the sleep spindle range of 11-16 Hz. The duration of HFA positively correlated with higher density and longer duration of accompanying thalamic spindles. Thalamic HFA's duration negatively correlated with the presence of cortical interictal epileptiform discharges. Thalamic spindles generated in channels with HFA often coincided with sleep spindles in various brain regions.

CONCLUSION

Fast ripple band HFA associated with sleep spindles was observed exclusively in the thalamus.

SIGNIFICANCE

Thalamic HFA associated with thalamic spindles may represent a thalamus-specific physiological phenomenon.

The Phase-Amplitude Coupling Changes Induced by Smoking Cue After 12-H Abstinence in Young Smokers

Tobacco use causes more than 8 million deaths globally each year, and the number of younger smokers is growing. It is of great practical importance to explore the underlying neural mechanisms behind the behaviour of young smokers. During cue-induced craving, reward system in the brain generates neural oscillations at specific frequencies. The phase-amplitude coupling (PAC) can capture interactions between these frequencies and may be a more sensitive quantitative indicator for characterizing abnormal neural oscillations in smokers. We monitored the electroencephalography (EEG) data of 30 young smokers during a cue task after 12 h of abstinence, dividing the data into the neutral and smoking-related groups based on different experimental stimuli to analyse the relationship between PAC and craving. In addition, we computed the functional connectivity (FC) under the PAC mechanism. The results showed that the young smokers exposed to smoking-related cues under short-term abstinence conditions had significantly lower PAC values and reduced FC strength in the right prefrontal cortex. In contrast, there was a significant increase in PAC values in the parietal cortex and enhanced FC strength. The correlation analysis showed significant correlations between PAC values and craving. These findings demonstrate for the first time that PAC abnormalities in young smokers exposed to smoking-related cues under short-term abstinence conditions may be related to craving and inhibitory control.

Event-related dynamic phase-amplitude coupling analysis reveals facial emotional processing deficits in patients with major depressive disorder: a cross-sectional study

BACKGROUND

Phase-amplitude coupling (PAC) measures the interaction between neural oscillations in different frequency bands and reflects brain functional network coordination in psychiatric patients. The event-related dynamic changes in PAC characteristics and their association with the neural physiological mechanisms under emotional stimulation in major depressive disorder (MDD) remain poorly understood.

METHODS

We proposed a cross-sectional study that investigated three PAC methods using simulated data and selected the Gaussian-Copula Event-Related PAC (GC-ERPAC) method for dynamic analysis of 128-channel electroencephalogram data from 53 participants, including 24 patients with MDD and 29 healthy controls (HCs). Participants were exposed to three emotional stimuli (fearful, happy, and sad). The correlation between PAC strengths and clinical scales was then analyzed in each condition.

RESULTS

The MDD group exhibited abnormal PAC patterns. With happy stimuli, the strengths of delta-gamma coupling (DGC), theta-gamma coupling (TGC), and alpha-gamma coupling (AGC) in the frontal-parietal regions of the MDD group were lower compared to HCs. With fearful stimuli, the strength of AGC in the occipital region was higher in the MDD group. The correlation between TGC and AGC was weaker for couplings among different frequencies in the MDD group. Additionally, AGC was negatively correlated with the clinical scale in MDD but positively correlated with HCs.

CLINICAL TRIAL NUMBER

Not applicable.

CONCLUSIONS

This study confirmed that GC-ERPAC was an effective method for revealing emotion recognition features in MDD. We provided novel evidence of abnormal spatiotemporal PAC patterns linked to emotional processing deficits. Reduced DGC and TGC, along with increased AGC, suggest potential biomarkers in MDD.

Decoding Single-Pellet Retrieval Task From Local Field Potentials in Pre- and Post-Stroke Motor Areas: Insights Into Interhemispheric Connectivity Difference

OBJECTIVE

Intracortical brain-machine interfaces (iBMIs) hold promise for restoring communication and movement in stroke-paralyzed individuals. Recent studies have demonstrated the potential of using local field potentials (LFPs) for decoding single-pellet retrieval (SPR) tasks in iBMIs. However, most research has relied on LFPs from healthy rats rather than those affected by stroke. This study aimed to investigate the feasibility of utilizing LFPs from both the right and left (stroke) cortical forelimb areas (CFAs) for the SPR tasks decoding under both pre- and post-stroke conditions.

METHODS

LFPs were recorded via microelectrode arrays implanted into CFAs of eight rats trained to perform the SPR tasks. The relative spectral power method was used to represent frequency information, and random forest classification differentiated SPR tasks from resting states. We also assessed interhemispheric connectivity, including correlation, coherence, and phase-amplitude coupling (PAC), to compare differences between the SPR tasks and the resting states under both pre- and post-stroke conditions.

RESULTS

Our findings indicated that the relative PS method with LFPs achieves 87.10% 9.2% accuracy in post-stoke SPR decoding, where high gamma is crucial. Additionally, we observed changes in PACs from the right to the left sensorimotor cortex post-stroke during the SPR tasks compared to the resting states.

SIGNIFICANCE

Our work provides a comprehensive insight into the role of different frequency band from LFPs in motor function recovery mechanisms, highlighting the importance of the high gamma in motor function. This research lays the foundation for developing post-stoke SPR-related BMIs.

Phase-Amplitude Coupling of NREM Sleep Oscillations Shows Between-Night Stability and is Related to Overnight Memory Gains

There is growing evidence in humans linking the temporal coupling between spindles and slow oscillations during NREM sleep with the overnight stabilization of memories encoded from daytime experiences in humans. However, whether the type and strength of learning influence that relationship is still unknown. Here we tested whether the amount or type of verbal word-pair learning prior to sleep affects subsequent phase-amplitude coupling (PAC) between spindles and slow oscillations (SO). We measured the strength and preferred timing of such coupling in the EEG of 41 healthy human participants over a post-learning and control night to compare intra-individual changes with inter-individual differences. We leveraged learning paradigms of varying word-pair (WP) load: 40 WP learned to a minimum criterion of 60% correct (n = 11); 40 WP presented twice (n = 15); 120 WP presented twice (n = 15). There were no significant differences in the preferred phase or strength between the control and post-learning nights, in all learning conditions. We observed an overnight consolidation effect (improved performance at delayed recall) for the criterion learning condition only, and only in this condition was the overnight change in memory performance significantly positively correlated with the phase of SO-spindle coupling. These results suggest that the coupling of brain oscillations during human NREM sleep is stable traits that are not modulated by the amount of pre-sleep learning, yet are implicated in the sleep-dependent consolidation of memory-especially when overnight gains in memory are observed.

Classifying mental motor tasks from chronic ECoG-BCI recordings using phase-amplitude coupling features

Introduction

Phase-amplitude coupling (PAC), the modulation of high-frequency neural oscillations by the phase of slower oscillations, is increasingly recognized as a marker of goal-directed motor behavior. Despite this interest, its specific role and potential value in decoding attempted motor movements remain unclear.

Methods

This study investigates whether PAC-derived features can be leveraged to classify different motor behaviors from ECoG signals within Brain-Computer Interface (BCI) systems. ECoG data were collected using the WIMAGINE implant during BCI experiments with a tetraplegic patient performing mental motor tasks. The data underwent preprocessing to extract complex neural oscillation features (amplitude, phase) through spectral decomposition techniques. These features were then used to quantify PAC by calculating different coupling indices. PAC metrics served as input features in a machine learning pipeline to evaluate their effectiveness in predicting mental tasks (idle state, right-hand movement, left-hand movement) in both offline and pseudo-online modes.

Results

The PAC features demonstrated high accuracy in distinguishing among motor tasks, with key classification features highlighting the coupling of theta/low-gamma and beta/high-gamma frequency bands.

Discussion

These preliminary findings hold significant potential for advancing our understanding of motor behavior and for developing optimized BCI systems.

Phase-amplitude coupling during auditory steady-state stimulation: a methodological review

Auditory steady-state response (ASSR) is a robust method to probe gamma (>30 Hz) activity in a controlled manner. While typically the magnitude and the phase synchronization over stimulus repetitions of ASSR is assessed, other measures are being investigated. One of them is phase-amplitude coupling (PAC), which reflects the interactions between lower frequency phase and higher frequency amplitude. Considering that the number of studies assessing PAC during auditory steady-state stimulation has grown recently, in the present work, we aimed to perform a comprehensive overview of PAC methodological approaches in ASSR studies. We sought to evaluate the studies according to PAC analysis issues emphasized in empirical and theoretical PAC studies. Our work showed considerable variability in the methodology among the reviewed studies. Furthermore, the reviewed works address methodological issues and confounding factors of PAC relatively poorly and are characterized by insufficient descriptions of the applied approaches. Our review shows that systematic research of PAC in the context of ASSR is imperative in order to properly evaluate the presence of PAC during the auditory steady-state stimulation.

The structure and statistics of language jointly shape cross-frequency neural dynamics during spoken language comprehension

Humans excel at extracting structurally-determined meaning from speech despite inherent physical variability. This study explores the brain's ability to predict and understand spoken language robustly. It investigates the relationship between structural and statistical language knowledge in brain dynamics, focusing on phase and amplitude modulation. Using syntactic features from constituent hierarchies and surface statistics from a transformer model as predictors of forward encoding models, we reconstructed cross-frequency neural dynamics from MEG data during audiobook listening. Our findings challenge a strict separation of linguistic structure and statistics in the brain, with both aiding neural signal reconstruction. Syntactic features have a more temporally spread impact, and both word entropy and the number of closing syntactic constituents are linked to the phase-amplitude coupling of neural dynamics, implying a role in temporal prediction and cortical oscillation alignment during speech processing. Our results indicate that structured and statistical information jointly shape neural dynamics during spoken language comprehension and suggest an integration process via a cross-frequency coupling mechanism.

Unconscious classification of quantitative electroencephalogram features from propofol versus propofol combined with etomidate anesthesia using one-dimensional convolutional neural network

Objective

Establishing a convolutional neural network model for the recognition of characteristic raw electroencephalogram (EEG) signals is crucial for monitoring consciousness levels and guiding anesthetic drug administration.

Methods

This trial was conducted from December 2023 to March 2024. A total of 40 surgery patients were randomly divided into either a propofol group (1% propofol injection, 10 mL: 100 mg) (P group) or a propofol-etomidate combination group (1% propofol injection, 10 mL: 100 mg, and 0.2% etomidate injection, 10 mL: 20 mg, mixed at a 2:1 volume ratio) (EP group). In the P group, target-controlled infusion (TCI) was employed for sedation induction, with an initial effect site concentration set at 5-6 μg/mL. The EP group received an intravenous push with a dosage of 0.2 mL/kg. Six consciousness-related EEG features were extracted from both groups and analyzed using four prediction models: support vector machine (SVM), Gaussian Naive Bayes (GNB), artificial neural network (ANN), and one-dimensional convolutional neural network (1D CNN). The performance of the models was evaluated based on accuracy, precision, recall, and F1-score.

Results

The power spectral density (94%) and alpha/beta ratio (72%) demonstrated higher accuracy as indicators for assessing consciousness. The classification accuracy of the 1D CNN model for anesthesia-induced unconsciousness (97%) surpassed that of the SVM (83%), GNB (81%), and ANN (83%) models, with a significance level of p < 0.05. Furthermore, the mean and mean difference ± standard error of the primary power values for the EP and P groups during the induced period were as follows: delta (23.85 and 16.79, 7.055 ± 0.817, p < 0.001), theta (10.74 and 8.743, 1.995 ± 0.7045, p < 0.02), and total power (24.31 and 19.72, 4.588 ± 0.7107, p < 0.001).

Conclusion

Large slow-wave oscillations, power spectral density, and the alpha/beta ratio are effective indicators of changes in consciousness during intravenous anesthesia with a propofol-etomidate combination. These indicators can aid anesthesiologists in evaluating the depth of anesthesia and adjusting dosages accordingly. The 1D CNN model, which incorporates consciousness-related EEG features, represents a promising tool for assessing the depth of anesthesia.

Clinical Trial Registration

https://www.chictr.org.cn/index.html.

Resting-state EEG biomarkers of accelerated intermittent theta burst stimulation treatment for depression: a pilot study

Accelerated intermittent theta burst stimulation (aiTBS) is a novel and effective treatment for drug-resistant depression. While past studies have identified encephalography (EEG) features predicting repetitive transcranial magnetic stimulation (rTMS) outcomes, EEG biomarkers specifically for aiTBS in depression patients have not been explored. In this pilot trial on 5 depression patients undergoing aiTBS, we assessed clinical outcome using the Montgomery-Asberg Depression Rating Scale (MADRS) and collected resting-state EEG pre and post-treatment. All patients showed an improvement in MADRS, with 3 having at least 50% improvement. We found significant correlations between MADRS change and pre-treatment frontal beta power, midline frontal Lempel-Ziv Complexity (LZC) and alpha connectivity. We also observed a trend of increased frontal theta power post-treatment. However, no significant correlations emerged between MADRS change and change in EEG feature post-treatment. This preliminary trial highlights the potential for investigating aiTBS-specific EEG biomarkers, paving the way for larger studies to enhance personalized neurostimulation and predict treatment outcomes in drug-resistant depression patients.

Uncovering spatiotemporal dynamics of the corticothalamic network at ictal onset

OBJECTIVE

Although the clinical efficacy of deep brain stimulation targeting the anterior nucleus (AN) and centromedian nucleus (CM) of the thalamus has been actively investigated for the treatment of medication-resistant epilepsy, few studies have investigated dynamic ictal changes in corticothalamic connectivity in human electroencephalographic (EEG) recording. This study aims to establish the complex spatiotemporal dynamics of the ictal corticothalamic network associated with various seizure foci.

METHODS

We analyzed 10 patients (aged 2.7-28.1 years) with medication-resistant focal epilepsy who underwent stereotactic EEG evaluation with thalamic sampling. We examined both undirected and directed connectivity, incorporating coherence and spectral Granger causality analysis (GCA) between the diverse seizure foci and thalamic nuclei (AN and CM) at ictal onset.

RESULTS

In our analysis of 36 seizures, coherence between seizure onset and thalamic nuclei increased across all frequencies, especially in slower bands (delta, theta, alpha). GCA showed increased information flow from seizure onset to the thalamus across all frequency bands, but outflows from the thalamus were mainly in slower frequencies, particularly delta. In the subgroup analysis based on various seizure foci, the delta coherence showed a more pronounced increase at CM than at AN during frontal lobe seizures. Conversely, in limbic seizures, the delta coherence increase was greater at AN compared to CM.

SIGNIFICANCE

It appears that the delta frequency plays a pivotal role in modulating the corticothalamic network during seizures. Our results underscore the significance of comprehending the spatiotemporal dynamics of the corticothalamic network at ictal onset, and this knowledge could guide personalized responsive neuromodulation treatment strategies.

Extracting Transient Phase-Amplitude Coupling from Resting-State EEG Signals

Electroencephalography (EEG) allows us to observe brain activity through electrical signals. Phase-amplitude coupling (PAC) is a way to analyze EEG data by focusing on the interaction between the low- and high-frequency components of these signals. However, PAC analyses are often challenged by various methodological issues. We here propose a novel approach which alleviates these issues. Our method has the following features: (i) it addresses the transient nature of coupling through data epoching; (ii) it ensures the presence of low-frequency oscillations through peak detection in the power spectrum; (iii) it applies adaptive high-frequency filtering; and (iv) it performs statistical validation using surrogate data. The efficiency of our method is demonstrated through both a simulation study and the analysis of experimental EEG data, offering new insights into the intricate workings of brain signal interactions.

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.

Gait-related beta-gamma phase amplitude coupling in the subthalamic nucleus of parkinsonian patients

Analysis of coupling between the phases and amplitudes of neural oscillations has gained increasing attention as an important mechanism for large-scale brain network dynamics. In Parkinson's disease (PD), preliminary evidence indicates abnormal beta-phase coupling to gamma-amplitude in different brain areas, including the subthalamic nucleus (STN). We analyzed bilateral STN local field potentials (LFPs) in eight subjects with PD chronically implanted with deep brain stimulation electrodes during upright quiet standing and unperturbed walking. Phase-amplitude coupling (PAC) was computed using the Kullback-Liebler method, based on the modulation index. Neurophysiological recordings were correlated with clinical and kinematic measurements and individual molecular brain imaging studies ([123I]FP-CIT and single-photon emission computed tomography). We showed a dopamine-related increase in subthalamic beta-gamma PAC from standing to walking. Patients with poor PAC modulation and low PAC during walking spent significantly more time in the stance and double support phase of the gait cycle. Our results provide new insights into the subthalamic contribution to human gait and suggest cross-frequency coupling as a gateway mechanism to convey patient-specific information of motor control for human locomotion.

Citizen neuroscience: Wearable technology and open software to study the human brain in its natural habitat

Citizen science allows the public to participate in various stages of scientific research, including study design, data acquisition, and data analysis. Citizen science has a long history in several fields of the natural sciences, and with recent developments in wearable technology, neuroscience has also become more accessible to citizen scientists. This development was largely driven by the influx of minimal sensing systems in the consumer market, allowing more do-it-yourself (DIY) and quantified-self (QS) investigations of the human brain. While most subfields of neuroscience require sophisticated monitoring devices and laboratories, the study of sleep characteristics can be performed at home with relevant noninvasive consumer devices. The strong influence of sleep quality on waking life and the accessibility of devices to measure sleep are two primary reasons citizen scientists have widely embraced sleep research. Their involvement has evolved from solely contributing to data collection to engaging in more collaborative or autonomous approaches, such as instigating ideas, formulating research inquiries, designing research protocols and methodology, acting upon their findings, and disseminating results. In this article, we introduce the emerging field of citizen neuroscience, illustrating examples of such projects in sleep research. We then provide overviews of the wearable technologies for tracking human neurophysiology and various open-source software used to analyse them. Finally, we discuss the opportunities and challenges in citizen neuroscience projects and suggest how to improve the study of the human brain outside the laboratory.

A dynamical computational model of theta generation in hippocampal circuits to study theta-gamma oscillations during neurostimulation

Neurostimulation of the hippocampal formation has shown promising results for modulating memory but the underlying mechanisms remain unclear. In particular, the effects on hippocampal theta-nested gamma oscillations and theta phase reset, which are both crucial for memory processes, are unknown. Moreover, these effects cannot be investigated using current computational models, which consider theta oscillations with a fixed amplitude and phase velocity. Here, we developed a novel computational model that includes the medial septum, represented as a set of abstract Kuramoto oscillators producing a dynamical theta rhythm with phase reset, and the hippocampal formation, composed of biophysically realistic neurons and able to generate theta-nested gamma oscillations under theta drive. We showed that, for theta inputs just below the threshold to induce self-sustained theta-nested gamma oscillations, a single stimulation pulse could switch the network behavior from non-oscillatory to a state producing sustained oscillations. Next, we demonstrated that, for a weaker theta input, pulse train stimulation at the theta frequency could transiently restore seemingly physiological oscillations. Importantly, the presence of phase reset influenced whether these two effects depended on the phase at which stimulation onset was delivered, which has practical implications for designing neurostimulation protocols that are triggered by the phase of ongoing theta oscillations. This novel model opens new avenues for studying the effects of neurostimulation on the hippocampal formation. Furthermore, our hybrid approach that combines different levels of abstraction could be extended in future work to other neural circuits that produce dynamical brain rhythms.

Abnormally low sensorimotor α band nonlinearity serves as an effective EEG biomarker of Parkinson's disease

Biomarkers obtained from the neurophysiological signals of patients with Parkinson's disease (PD) have objective value in assessing their motor condition for effective diagnosis, monitoring, and clinical intervention. Prominent cortical biomarkers of PD have typically been derived from various β band wave features. This study approached the topic from an alternative perspective and attempted to estimate a recently suggested measure representing α band nonlinear autocorrelative memory from a publicly available EEG dataset that involves 15 patients with earlier-stage PD (dopaminergic medication OFF and ON states) and 16 age-matched healthy controls. The cortical nonlinearity was elevated for the PD ON state compared with the OFF state for bilateral sensorimotor channels C3 and C4 (n = 26; P = 0.003). A similar statistical difference was also identified between PD OFF state and healthy subjects (n = 26; P = 0.049). Analysis over all channels revealed that the α band nonlinearity induced upon medication was constrained to sensorimotor regions. The α nonlinearity measure was compared with a well-accepted cortical biomarker of β-γ phase-amplitude coupling (PAC). They were in moderate negative correlation (r = -0.412; P = 0.036) for only healthy subjects, but not for the patients. The nonlinearity measure was found to be insusceptible to the nonstationary variations within the particular data. Our study provides further evidence that the α band nonlinearity measure can serve as a promising cortical biomarker of PD. The suggested measure can be estimated from a noninvasive low-resolution single scalp EEG channel of patients with relatively early-stage PD, who did not yet need to undergo deep brain stimulation operation.NEW & NOTEWORTHY This study suggests a nonlinearity measure that differentiates Parkinson's disease (PD) dopamine OFF-state scalp EEG data from those of dopamine ON-state patients and healthy subjects. Unlike typical PD cortical biomarkers based on β band activity, this metric shows elevation upon dopaminergic medication in the α band. We provide evidence supporting its potential as an early-stage promising PD biomarker that can be estimated from noninvasive EEG recordings with low resolution and SNR.

Targeted memory reactivation during post-learning sleep does not enhance motor memory consolidation in older adults

Targeted memory reactivation (TMR) during sleep enhances memory consolidation in young adults by modulating electrophysiological markers of neuroplasticity. Interestingly, older adults exhibit deficits in motor memory consolidation, an impairment that has been linked to age-related degradations in the same sleep features sensitive to TMR. We hypothesised that TMR would enhance consolidation in older adults via the modulation of these markers. A total of 17 older participants were trained on a motor task involving two auditory-cued sequences. During a post-learning nap, two auditory cues were played: one associated to a learned (i.e., reactivated) sequence and one control. Performance during two delayed re-tests did not differ between reactivated and non-reactivated sequences. Moreover, both associated and control sounds modulated brain responses, yet there were no consistent differences between the auditory cue types. Our results collectively demonstrate that older adults do not benefit from specific reactivation of a motor memory trace by an associated auditory cue during post-learning sleep. Based on previous research, it is possible that auditory stimulation during post-learning sleep could have boosted motor memory consolidation in a non-specific manner.

Towards Optimization of Oscillatory Stimulation During Sleep

BACKGROUND

Oscillatory rhythms during sleep, such as slow oscillations (SOs) and spindles and, most importantly, their coupling, are thought to underlie processes of memory consolidation. External slow oscillatory transcranial direct current stimulation (so-tDCS) with a frequency of 0.75 Hz has been shown to improve this coupling and memory consolidation; however, effects varied quite markedly between individuals, studies, and species. In this study, we aimed to determine how precisely the frequency of stimulation must match the naturally occurring SO frequency in individuals to best improve SO-spindle coupling. Moreover, we systematically tested stimulation durations necessary to induce changes.

MATERIALS AND METHODS

We addressed these questions by comparing so-tDCS with individualized frequency to standardized frequency of 0.75 Hz in a within-subject design with 28 older participants during napping while stimulation train durations were systematically varied between 30 seconds, 2 minutes, and 5 minutes.

RESULTS

Stimulation trains as short as 30 seconds were sufficient to modulate the coupling between SOs and spindle activity. Contrary to our expectations, so-tDCS with standardized frequency indicated stronger aftereffects regarding SO-spindle coupling than individualized frequency. Angle and variance of spindle maxima occurrence during the SO cycle were similarly modulated.

CONCLUSIONS

In sum, short stimulation trains were sufficient to induce significant changes in sleep physiology, allowing for more trains of stimulation, which provides methodological advantages and possibly even larger behavioral effects in future studies. Regarding individualized stimulation frequency, further options of optimization need to be investigated, such as closed-loop stimulation, to calibrate stimulation frequency to the SO frequency at the time of stimulation onset.

CLINICAL TRIAL REGISTRATION

The Clinicaltrials.gov registration number for the study is NCT04714879.

Coordinated human sleeping brainwaves map peripheral body glucose homeostasis

Insufficient sleep impairs glucose regulation, increasing the risk of diabetes. However, what it is about the human sleeping brain that regulates blood sugar remains unknown. In an examination of over 600 humans, we demonstrate that the coupling of non-rapid eye movement (NREM) sleep spindles and slow oscillations the night before is associated with improved next-day peripheral glucose control. We further show that this sleep-associated glucose pathway may influence glycemic status through altered insulin sensitivity, rather than through altered pancreatic beta cell function. Moreover, we replicate these associations in an independent dataset of over 1,900 adults. Of therapeutic significance, the coupling between slow oscillations and spindles was the most significant sleep predictor of next-day fasting glucose, even more so than traditional sleep markers, relevant to the possibility of an electroencephalogram (EEG) index of hyperglycemia. Taken together, these findings describe a sleeping-brain-body framework of optimal human glucose homeostasis, offering a potential prognostic sleep signature of glycemic control.

Cross-Frequency Coupling and Intelligent Neuromodulation

Cross-frequency coupling (CFC) reflects (nonlinear) interactions between signals of different frequencies. Evidence from both patient and healthy participant studies suggests that CFC plays an essential role in neuronal computation, interregional interaction, and disease pathophysiology. The present review discusses methodological advances and challenges in the computation of CFC with particular emphasis on potential solutions to spurious coupling, inferring intrinsic rhythms in a targeted frequency band, and causal interferences. We specifically focus on the literature exploring CFC in the context of cognition/memory tasks, sleep, and neurological disorders, such as Alzheimer's disease, epilepsy, and Parkinson's disease. Furthermore, we highlight the implication of CFC in the context and for the optimization of invasive and noninvasive neuromodulation and rehabilitation. Mainly, CFC could support advancing the understanding of the neurophysiology of cognition and motor control, serve as a biomarker for disease symptoms, and leverage the optimization of therapeutic interventions, e.g., closed-loop brain stimulation. Despite the evident advantages of CFC as an investigative and translational tool in neuroscience, further methodological improvements are required to facilitate practical and correct use in cyborg and bionic systems in the field.

BOARD-FTD-PACC: a graphical user interface for the synaptic and cross-frequency analysis derived from neural signals

In order to understand the link between brain functional states and behavioral/cognitive processes, the information carried in neural oscillations can be retrieved using different analytic techniques. Processing these different bio-signals is a complex, time-consuming, and often non-automatized process that requires customization, due to the type of signal acquired, acquisition method implemented, and the objectives of each individual research group. To this end, a new graphical user interface (GUI), named BOARD-FTD-PACC, was developed and designed to facilitate the visualization, quantification, and analysis of neurophysiological recordings. BOARD-FTD-PACC provides different and customizable tools that facilitate the task of analyzing post-synaptic activity and complex neural oscillatory data, mainly cross-frequency analysis. It is a flexible and user-friendly software that can be used by a wide range of users to extract valuable information from neurophysiological signals such as phase-amplitude coupling and relative power spectral density, among others. BOARD-FTD-PACC allows researchers to select, in the same open-source GUI, different approaches and techniques that will help promote a better understanding of synaptic and oscillatory activity in specific brain structures with or without stimulation.

Toward the Identification of Neurophysiological Biomarkers for Alzheimer's Disease in Down Syndrome: A Potential Role for Cross-Frequency Phase-Amplitude Coupling Analysis

Cross-frequency coupling (CFC) mechanisms play a central role in brain activity. Pathophysiological mechanisms leading to many brain disorders, such as Alzheimer's disease (AD), may produce unique patterns of brain activity detectable by electroencephalography (EEG). Identifying biomarkers for AD diagnosis is also an ambition among research teams working in Down syndrome (DS), given the increased susceptibility of people with DS to develop early-onset AD (DS-AD). Here, we review accumulating evidence that altered theta-gamma phase-amplitude coupling (PAC) may be one of the earliest EEG signatures of AD, and therefore may serve as an adjuvant tool for detecting cognitive decline in DS-AD. We suggest that this field of research could potentially provide clues to the biophysical mechanisms underlying cognitive dysfunction in DS-AD and generate opportunities for identifying EEG-based biomarkers with diagnostic and prognostic utility in DS-AD.

Effects on resting-state EEG phase-amplitude coupling in insomnia disorder patients following 1 Hz left dorsolateral prefrontal cortex rTMS

Despite burgeoning evidence for cortical hyperarousal in insomnia disorder, the existing results on electroencephalography spectral features are highly heterogeneous. Phase-amplitude coupling, which refers to the modulation of the low-frequency phase to a high-frequency amplitude, is probably a more sensitive quantitative measure for characterizing abnormal neural oscillations and explaining the therapeutic effect of repetitive transcranial magnetic stimulation in the treatment of patients with insomnia disorder. Sixty insomnia disorder patients were randomly divided into the active and sham treatment groups to receive 4 weeks of repetitive transcranial magnetic stimulation treatment. Behavioral assessments, resting-state electroencephalography recordings, and sleep polysomnography recordings were performed before and after repetitive transcranial magnetic stimulation treatment. Forty good sleeper controls underwent the same assessment. We demonstrated that phase-amplitude coupling values in the frontal and temporal lobes were weaker in Insomnia disorder patients than in those with good sleeper controls at baseline and that phase-amplitude coupling values near the intervention area were significantly enhanced after active repetitive transcranial magnetic stimulation treatment. Furthermore, the enhancement of phase-amplitude coupling values was significantly correlated with the improvement of sleep quality. This study revealed the potential of phase-amplitude coupling in assessing the severity of insomnia disorder and the efficacy of repetitive transcranial magnetic stimulation treatment, providing new insights on the abnormal physiological mechanisms and future treatments for insomnia disorder.

Beyond linear neural envelope tracking: a mutual information approach

Objective.The human brain tracks the temporal envelope of speech, which contains essential cues for speech understanding. Linear models are the most common tool to study neural envelope tracking. However, information on how speech is processed can be lost since nonlinear relations are precluded. Analysis based on mutual information (MI), on the other hand, can detect both linear and nonlinear relations and is gradually becoming more popular in the field of neural envelope tracking. Yet, several different approaches to calculating MI are applied with no consensus on which approach to use. Furthermore, the added value of nonlinear techniques remains a subject of debate in the field. The present paper aims to resolve these open questions.Approach.We analyzed electroencephalography (EEG) data of participants listening to continuous speech and applied MI analyses and linear models.Main results.Comparing the different MI approaches, we conclude that results are most reliable and robust using the Gaussian copula approach, which first transforms the data to standard Gaussians. With this approach, the MI analysis is a valid technique for studying neural envelope tracking. Like linear models, it allows spatial and temporal interpretations of speech processing, peak latency analyses, and applications to multiple EEG channels combined. In a final analysis, we tested whether nonlinear components were present in the neural response to the envelope by first removing all linear components in the data. We robustly detected nonlinear components on the single-subject level using the MI analysis.Significance.We demonstrate that the human brain processes speech in a nonlinear way. Unlike linear models, the MI analysis detects such nonlinear relations, proving its added value to neural envelope tracking. In addition, the MI analysis retains spatial and temporal characteristics of speech processing, an advantage lost when using more complex (nonlinear) deep neural networks.

Effects of Contralateral Deep Brain Stimulation and Levodopa on Subthalamic Nucleus Oscillatory Activity and Phase-Amplitude Coupling

BACKGROUND

The modulatory effects of medication and deep brain stimulation (DBS) on subthalamic nucleus (STN) neural activity in Parkinson's disease have been widely studied. However, effects on the contralateral side to the stimulated STN, in particular, changes in local field potential (LFP) oscillatory activity and phase-amplitude coupling (PAC), have not yet been reported.

OBJECTIVE

The aim of this study was to examine changes in STN LFP activity across a range of frequency bands and STN PAC for different combinations of DBS and medication on/off on the side contralateral to the applied stimulation.

MATERIALS AND METHODS

We examined STN LFPs that were recorded using externalized leads from eight parkinsonian patients during unilateral DBS from the side contralateral to the stimulation. LFP spectral power in alpha (5 to ∼13 Hz), low beta (13 to ∼20 Hz), high beta (20-30 Hz), and high gamma plus high-frequency oscillation (high gamma+HFO) (100-400 Hz) bands were estimated for different combinations of medication and unilateral stimulation (off/on). PAC between beta and high gamma+HFO in the STN LFPs was also investigated. The effect of the condition was examined using linear mixed models.

RESULTS

PAC in the STN LFP was reduced by DBS when compared to the baseline condition (no medication and stimulation). Medication had no significant effect on PAC. Alpha power decreased with DBS, both alone and when combined with medication. Beta power decreased with DBS, medication, and DBS and medication combined. High gamma+HFO power increased during the application of contralateral DBS and was unaltered by medication.

CONCLUSIONS

The results provide new insights into the effects of DBS and levodopa on STN LFP PAC and oscillatory activity on the side contralateral to stimulation. These may have important implications in understanding mechanisms underlying motor improvements with DBS, including changes on both contralateral and ipsilateral sides, while suggesting a possible role for contralateral sensing during unilateral DBS.

Profiling the low-beta characteristics of the subthalamic nucleus in early- and late-onset Parkinson's disease

Objectives

Low-beta oscillation (13-20 Hz) has rarely been studied in patients with early-onset Parkinson's disease (EOPD, age of onset ≤50 years). We aimed to explore the characteristics of low-beta oscillation in the subthalamic nucleus (STN) of patients with EOPD and investigate the differences between EOPD and late-onset Parkinson's disease (LOPD).

Methods

We enrolled 31 EOPD and 31 LOPD patients, who were matched using propensity score matching. Patients underwent bilateral STN deep brain stimulation (DBS). Local field potentials were recorded using intraoperative microelectrode recording. We analyzed the low-beta band parameters, including aperiodic/periodic components, beta burst, and phase-amplitude coupling. We compared low-beta band activity between EOPD and LOPD. Correlation analyses were performed between the low-beta parameters and clinical assessment results for each group.

Results

We found that the EOPD group had lower aperiodic parameters, including offset (p = 0.010) and exponent (p = 0.047). Low-beta burst analysis showed that EOPD patients had significantly higher average burst amplitude (p = 0.016) and longer average burst duration (p = 0.011). Furthermore, EOPD had higher proportion of long burst (500-650 ms, p = 0.008), while LOPD had higher proportion of short burst (200-350 ms, p = 0.007). There was a significant difference in phase-amplitude coupling values between low-beta phase and fast high frequency oscillation (300-460 Hz) amplitude (p = 0.019).

Conclusion

We found that low-beta activity in the STN of patients with EOPD had characteristics that varied when compared with LOPD, and provided electrophysiological evidence for different pathological mechanisms between the two types of PD. These differences need to be considered when applying adaptive DBS on patients of different ages.

Phase-amplitude coupling-based adaptive filters for neural signal decoding

Bandpass filters play a core role in ECoG signal processing. Commonly used frequency bands such as alpha, beta, and gamma bands can reflect the normal rhythm of the brain. However, the universally predefined bands might not be optimal for a specific task. Especially the gamma band usually covers a wide frequency span (i.e., 30-200 Hz) which can be too coarse to capture features that appear in narrow bands. An ideal option is to find the optimal frequency bands for specific tasks in real-time and dynamically. To tackle this problem, we propose an adaptive band filter that selects the useful frequency band in a data-driven way. Specifically, we leverage the phase-amplitude coupling (PAC) of the coupled working mechanism of synchronizing neuron and pyramidal neurons in neuronal oscillations, in which the phase of slower oscillations modulates the amplitude of faster ones, to help locate the fine frequency bands from the gamma range, in a task-specific and individual-specific way. Thus, the information can be more precisely extracted from ECoG signals to improve neural decoding performance. Based on this, an end-to-end decoder (PACNet) is proposed to construct a neural decoding application with adaptive filter banks in a uniform framework. Experiments show that PACNet can improve neural decoding performance universally with different tasks.

Sigma oscillations protect or reinstate motor memory depending on their temporal coordination with slow waves

Targeted memory reactivation (TMR) during post-learning sleep is known to enhance motor memory consolidation but the underlying neurophysiological processes remain unclear. Here, we confirm the beneficial effect of auditory TMR on motor performance. At the neural level, TMR enhanced slow wave (SW) characteristics. Additionally, greater TMR-related phase-amplitude coupling between slow (0.5–2 Hz) and sigma (12–16 Hz) oscillations after the SW peak was related to higher TMR effect on performance. Importantly, sounds that were not associated to learning strengthened SW-sigma coupling at the SW trough. Moreover, the increase in sigma power nested in the trough of the potential evoked by the unassociated sounds was related to the TMR benefit. Altogether, our data suggest that, depending on their precise temporal coordination during post learning sleep, slow and sigma oscillations play a crucial role in either memory reinstatement or protection against irrelevant information; two processes that critically contribute to motor memory consolidation.

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.

OUP accepted manuscript

Freezing of gait is a debilitating symptom in advanced Parkinson’s disease and responds heterogeneously to treatments such as deep brain stimulation. Recent studies indicated that cortical dysfunction is involved in the development of freezing, while evidence depicting the specific role of the primary motor cortex in the multi-circuit pathology of freezing is lacking. Since abnormal beta-gamma phase-amplitude coupling recorded from the primary motor cortex in patients with Parkinson’s disease indicates parkinsonian state and responses to therapeutic deep brain stimulation, we hypothesized this metric might reveal unique information on understanding and improving therapy for freezing of gait.

Here, we directly recorded potentials in the primary motor cortex using subdural electrocorticography and synchronously captured gait freezing using optoelectronic motion-tracking systems in 16 freely-walking patients with Parkinson’s disease who received subthalamic nucleus deep brain stimulation surgery. Overall, we recorded 451 timed up-and-go walking trials and quantified 7073 s of stable walking and 3384 s of gait freezing in conditions of on/off-stimulation and with/without dual-tasking.

We found that (i) high beta-gamma phase-amplitude coupling in the primary motor cortex was detected in freezing trials (i.e. walking trials that contained freezing), but not non-freezing trials, and the high coupling in freezing trials was not caused by dual-tasking or the lack of movement; (ii) non-freezing episodes within freezing trials also demonstrated abnormally high couplings, which predicted freezing severity; (iii) deep brain stimulation of subthalamic nucleus reduced these abnormal couplings and simultaneously improved freezing; and (iv) in trials that were at similar coupling levels, stimulation trials still demonstrated lower freezing severity than no-stimulation trials.

These findings suggest that elevated phase-amplitude coupling in the primary motor cortex indicates higher probabilities of freezing. Therapeutic deep brain stimulation alleviates freezing by both decoupling cortical oscillations and enhancing cortical resistance to abnormal coupling. We formalized these findings to a novel ‘bandwidth model,’ which specifies the role of cortical dysfunction, cognitive burden and therapeutic stimulation on the emergence of freezing. By targeting key elements in the model, we may develop next-generation deep brain stimulation approaches for freezing of gait.

The Overnight Retention of Novel Metaphors Associates With Slow Oscillation-Spindle Coupling but Not With Respiratory Phase at Encoding

Accumulating evidence emphasizes the relevance of oscillatory synchrony in memory consolidation during sleep. Sleep spindles promote memory retention, especially when occurring in the depolarized upstate of slow oscillation (SO). A less studied topic is the inter-spindle synchrony, i.e. the temporal overlap and phasic coherence between spindles perceived in different electroencephalography channels. In this study, we examined how synchrony between SOs and spindles, as well as between simultaneous spindles, is associated with the retention of novel verbal metaphors. Moreover, we combined the encoding of the metaphors with respiratory phase (inhalation/exhalation) with the aim of modulating the strength of memorized items, as previous studies have shown that inhalation entrains neural activity, thereby benefiting memory in a waking condition. In the current study, 27 young adults underwent a two-night mixed-design study with a 12-h delayed memory task during both sleep and waking conditions. As expected, we found better retention over the delay containing sleep, and this outcome was strongly associated with the timing of SO–spindle coupling. However, no associations were observed regarding inter-spindle synchrony or respiratory phase. These findings contribute to a better understanding of the importance of SO–spindle coupling for memory. In contrast, the observed lack of association with inter-spindle synchrony may emphasize the local nature of spindle-related plasticity.

Multitaper estimates of phase-amplitude coupling

Phase-amplitude coupling (PAC) is the association of the amplitude of a high-frequency oscillation with the phase of a low-frequency oscillation. In neuroscience, this relationship provides a mechanism by which neural activity might be coordinated between distant regions. The dangers and pitfalls of assessing PAC with commonly used statistical measures have been well-documented. The limitations of these measures include: (1) response to non-oscillatory, high-frequency, broad-band activity, (2) response to high-frequency components of the low-frequency oscillation, (3) adhoc selection of analysis frequency-intervals, and (4) reliance upon data shuffling to assess statistical significance.Objective.To address issues (1)-(4) by introducing a nonparametric multitaper estimator of PAC.Approach.In this work, a multitaper PAC estimator is proposed that addresses these issues. Specifically, issue (1) is addressed by replacing the analytic signal envelope estimator computed using the Hilbert transform with a multitaper estimator that down-weights non-sinusoidal activity using a classical, multitaper super-resolution technique. Issue (2) is addressed by replacing coherence between the low-frequency and high-frequency components in a standard PAC estimator with multitaper partial coherence, while issue (3) is addressed with a physical argument regarding meaningful neural oscillation. Finally, asymptotic statistical assessment of the multitaper estimator is introduced to address issue (4).Main results.Multitaper estimates of PAC are introduced. Their efficacy is demonstrated in simulation and on human intracranial recordings obtained from epileptic patients.Significance.This work facilitates a more informative statistical assessment of PAC, a phenomena exhibited by many neural systems, and provides a basis upon which further nonparametric multitaper-related methods can be developed.

Measuring Phase-Amplitude Coupling Based on the Jensen-Shannon Divergence and Correlation Matrix

Phase-amplitude coupling (PAC) measures the relationship between the phase of low-frequency oscillations (LFO) and the amplitude of high-frequency oscillations (HFO). It plays an important functional role in neural information processing and cognition. Thus, we propose a novel method based on the Jensen-Shannon (JS) divergence and correlation matrix. The method takes the amplitude distributions of the HFO located in the corresponding phase bins of the LFO as multichannel inputs to construct a correlation matrix, where the elements are calculated based on the JS divergence between pairs of amplitude distributions. Then, the omega complexity extracted from the correlation matrix is used to estimate the PAC strength. The simulation results demonstrate that the proposed method can effectively reflect the PAC strength and slightly vary with the data length. Moreover, it outperforms five frequently used algorithms in the performance against additive white Gaussian noise and spike noise and the ability of detecting PAC in wide frequency ranges. To validate our proposed method with real data, it was applied to analyze the local field potential recorded from the dorsomedial striatum in a male Sprague-Dawley rat. It can replicate previous results obtained with other PAC metrics. In conclusion, these results suggest that our proposed method is a powerful tool for measuring the PAC between neural oscillations.