A novel cross-frequency coupling detection method using the generalized Morse wavelets

AML-DECODER: Advanced Machine Learning for HD-sEMG Signal Classification-Decoding Lateral Epicondylitis in Forearm Muscles

BACKGROUND

Innovative algorithms for wearable devices and garments are critical for diagnosing and monitoring disease (such as lateral epicondylitis (LE)) progression. LE affects individuals across various professions and causes daily problems.

METHODS

We analyzed signals from the forearm muscles of 14 healthy controls and 14 LE patients using high-density surface electromyography. We discerned significant differences between groups by employing phase-amplitude coupling (PAC) features. Our study leveraged PAC, Daubechies wavelet with four vanishing moments (db4), and state-of-the-art techniques to train a neural network for the subject's label prediction.

RESULTS

Remarkably, PAC features achieved 100% specificity and sensitivity in predicting unseen subjects, while state-of-the-art features lagged with only 35.71% sensitivity and 28.57% specificity, and db4 with 78.57% sensitivity and 85.71 specificity. PAC significantly outperformed the state-of-the-art features (adj. p-value < 0.001) with a large effect size. However, no significant difference was found between PAC and db4 (adj. p-value = 0.147). Also, the Jeffries-Matusita (JM) distance of the PAC was significantly higher than other features (adj. p-value < 0.001), with a large effect size, suggesting PAC features as robust predictors of neuromuscular diseases, offering a profound understanding of disease pathology and new avenues for interpretation. We evaluated the generalization ability of the PAC model using 99.9% confidence intervals and Bayesian credible intervals to quantify prediction uncertainty across subjects. Both methods demonstrated high reliability, with an expected accuracy of 89% in larger, more diverse populations.

CONCLUSIONS

This study's implications might extend beyond LE, paving the way for enhanced diagnostic tools and deeper insights into the complexities of neuromuscular disorders.

Errors elicit frontoparietal theta-gamma coupling that is modulated by endogenous estradiol levels

Cognitive control-related error monitoring is intimately involved in behavioral adaptation, learning, and individual differences in a variety of psychological traits and disorders. Accumulating evidence suggests that a focus on women's health and ovarian hormones is critical to the study of such cognitive brain functions. Here we sought to identify a novel index of error monitoring using a time-frequency based phase amplitude coupling (t-f PAC) measure and examine its modulation by endogenous levels of estradiol in females. Forty-three healthy, naturally cycling young adult females completed a flanker task while continuous electroencephalogram was recorded on four occasions across the menstrual cycle. Results revealed significant error-related t-f PAC between theta phase generated in fronto-central areas and gamma amplitude generated in parietal-occipital areas. Moreover, this error-related theta-gamma coupling was enhanced by endogenous levels of estradiol both within females across the cycle as well as between females with higher levels of average circulating estradiol. While the role of frontal midline theta in error processing is well documented, this paper extends the extant literature by illustrating that error monitoring involves the coordination between multiple distributed systems with the slow midline theta activity modulating the power of gamma-band oscillatory activity in parietal regions. They further show enhancement of inter-regional coupling by endogenous estradiol levels, consistent with research indicating modulation of cognitive control neural functions by the endocrine system in females. Together, this work identifies a novel neurophysiological marker of cognitive control-related error monitoring in females that has implications for neuroscience and women's health.

[A study on the application of cross-frequency coupling characteristics of neural oscillation in the diagnosis of mild cognitive impairment]

In order to fully explore the neural oscillatory coupling characteristics of patients with mild cognitive impairment (MCI), this paper analyzed and compared the strength of the coupling characteristics for 28 MCI patients and 21 normal subjects under six different-frequency combinations. The results showed that the difference in the global phase synchronization index of cross-frequency coupling under δ-θ rhythm combination was statistically significant in the MCI group compared with the normal control group ( P = 0.025, d = 0.398). To further validate this coupling feature, this paper proposed an optimized convolutional neural network model that incorporated a time-frequency data enhancement module and batch normalization layers to prevent overfitting while enhancing the robustness of the model. Based on this optimized model, with the phase locking value matrix of δ-θ rhythm combination as the single input feature, the diagnostic accuracy of MCI patients was (95.49 ± 4.15)%, sensitivity and specificity were (93.71 ± 7.21)% and (97.50 ± 5.34)%, respectively. The results showed that the characteristics of the phase locking value matrix under the combination of δ-θ rhythms can adequately reflect the cognitive status of MCI patients, which is helpful to assist the diagnosis of MCI.

A Data-Driven System Identification Method for Random Eigenvalue Problem Using Synchrosqueezed Energy and Phase Portrait Analysis

The primary purpose of this research is to evaluate the uncertainty associated with modal parameter estimation for an inverse dynamic problem in which the structural parameters are random. The random nature of the structure's parameters will be reflected in the modal features of the respected system. However, this may result in additive/subtractive errors in modal parameter identification, affecting the identification technique's efficiency. With this in mind, the present study aims to develop an automated modal identification algorithm for a random eigenvalue problem. This is achieved by a recently developed advanced version of the wavelet transform (i.e., synchrosqueezing), which offers better resolution. Using this technique, the measured responses are transformed into a time-frequency plane, which is further processed by unsupervised learning using K-means clustering for quantification of the modal parameters. This automated identification is repeated for an ensemble of measurements to quantify the random eigenvalues in a statistical sense. The proposed methodology is first tested using simulated time histories of a two degree-of-freedom (dof) system. It is followed by an experimental validation using a beam whose mass matrix is random. The numerical results presented in this work clearly demonstrate the performance (i.e., in terms of efficiency and accuracy) of the proposed output-only automated data-driven identification scheme for random eigenvalue problems.

Brain-wide neural co-activations in resting human

Spontaneous neural activity in human as assessed with resting-state functional magnetic resonance imaging (fMRI) exhibits brain-wide coordinated patterns in the frequency of < 0.1 Hz. However, understanding of fast brain-wide networks at the timescales of neuronal events (milliseconds to sub-seconds) and their spatial, spectral, and transitional characteristics remain limited due to the temporal constraints of hemodynamic signals. With milli-second resolution and whole-head coverage, scalp-based electroencephalography (EEG) provides a unique window into brain-wide networks with neuronal-timescale dynamics, shedding light on the organizing principles of brain functions. Using the state-of-the-art signal processing techniques, we reconstructed cortical neural tomography from resting-state EEG and extracted component-based co-activation patterns (cCAPs). These cCAPs revealed brain-wide intrinsic networks and their dynamics, indicating the configuration/reconfiguration of resting human brains into recurring and transitional functional states, which are featured with the prominent spatial phenomena of global patterns and anti-state pairs of co-(de)activations. Rich oscillational structures across a wide frequency band (i.e., 0.6 Hz, 5 Hz, and 10 Hz) were embedded in the nonstationary dynamics of these functional states. We further identified a superstructure that regulated between-state immediate and long-range transitions involving the entire set of identified cCAPs and governed a significant aspect of brain-wide network dynamics. These findings demonstrated how resting-state EEG data can be functionally decomposed using cCAPs to reveal rich dynamic structures of brain-wide human neural activations.

Abnormal EEG Signal Energy in the Elderly: A Wavelet Analysis of Event-related Potentials During a Stroop Task

BACKGROUND

Previous work showed that elderly with excess in theta activity in their resting state electroencephalogram (EEG) are at higher risk of cognitive decline than those with a normal EEG. By using event-related potentials (ERP) during a counting Stroop task, our prior work showed that elderly with theta excess have a large P300 component compared with normal EEG group. This increased activity could be related to a higher EEG signal energy used during this task.

NEW METHOD

By wavelet analysis applied to ERP obtained during a counting Stroop task we quantified the energy in the different frequency bands of a group of elderly with altered EEG.

RESULTS

In theta and alpha bands, the total energy was higher in elderly subjects with theta excess, specifically in the stimulus categorization window (258-516 ms). Both groups solved the task with similar efficiency.

COMPARISON WITH EXISTING METHODS

The traditional ERP analysis in elderly compares voltage among conditions and groups for a given time windows, while the frequency composition is not usually examined. We complemented our previous ERP analysis using a wavelet methodology. Furthermore, we showed the advantages of wavelet analysis over Short Time Fourier Transform when exploring EEG signal during this task.

CONCLUSIONS

The higher EEG signal energy in ERP might reflect undergoing neurobiological mechanisms that allow the elderly with theta excess to cope with the cognitive task with similar behavioral results as the normal EEG group. This increased energy could promote a metabolic and cellular dysregulation causing a greater decline in cognitive function.

Addressing Pitfalls in Phase-Amplitude Coupling Analysis with an Extended Modulation Index Toolbox

Phase-amplitude coupling (PAC) is proposed to play an essential role in coordinating the processing of information on local and global scales. In recent years, the methods able to reveal trustworthy PAC has gained considerable interest. However, the intrinsic features of some signals can lead to the identification of spurious or waveform-dependent coupling. This prompted us to develop an easily accessible tool that could be used to differentiate spurious from authentic PAC. Here, we propose a new tool for more reliable detection of PAC named the Extended Modulation Index (eMI) based on the classical Modulation Index measure of coupling. eMI is suitable both for continuous and epoched data and allows estimation of the statistical significance of each pair of frequencies for phase and for amplitude in the whole comodulogram in the framework of extreme value statistics. We compared eMI with the reference PAC measures-direct PAC estimator (a modification of Mean Vector Length) and standard Modulation Index. All three methods were tested using computer-simulated data and actual local field potential recordings from freely moving rats. All methods exhibited similar properties in terms of sensitivity and specificity of PAC detection. eMI proved to be more selective in the dimension of frequency for phase. One of the novelty's offered by eMI is a heuristic algorithm for classification of PAC as Reliable or Ambiguous. It relies on analysis of the relation between the spectral properties of the signal and the detected coupling. Moreover, eMI generates visualizations that support further evaluation of the coupling properties. It also introduces the concept of the polar phase-histogram to study phase relations of coupled slow and fast oscillations. We discuss the extent to which eMI addresses the known problems of interpreting PAC. The Matlab® toolbox implementing eMI framework, and the two reference PAC estimators is freely available as EEGLAB plugin at https://github.com/GabrielaJurkiewicz/ePAC .

Generative models, linguistic communication and active inference

This paper presents a biologically plausible generative model and inference scheme that is capable of simulating communication between synthetic subjects who talk to each other. Building on active inference formulations of dyadic interactions, we simulate linguistic exchange to explore generative models that support dialogues. These models employ high-order interactions among abstract (discrete) states in deep (hierarchical) models. The sequential nature of language processing mandates generative models with a particular factorial structure-necessary to accommodate the rich combinatorics of language. We illustrate linguistic communication by simulating a synthetic subject who can play the 'Twenty Questions' game. In this game, synthetic subjects take the role of the questioner or answerer, using the same generative model. This simulation setup is used to illustrate some key architectural points and demonstrate that many behavioural and neurophysiological correlates of linguistic communication emerge under variational (marginal) message passing, given the right kind of generative model. For example, we show that theta-gamma coupling is an emergent property of belief updating, when listening to another.

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.

Intact Auditory Cortical Cross-Frequency Coupling in Early and Chronic Schizophrenia

BACKGROUND

Previous work has identified a hierarchical organization of neural oscillations that supports performance of complex cognitive and perceptual tasks, and can be indexed with phase-amplitude coupling (PAC) between low- and high-frequency oscillations. Our aim was to employ enhanced source localization afforded by magnetoencephalography (MEG) to expand on earlier reports of intact auditory cortical PAC in schizophrenia and to investigate how PAC may evolve over the early and chronic phases of the illness.

METHODS

Individuals with early schizophrenia (n=12) (≤5 years of illness duration), chronic schizophrenia (n=16) (>5 years of illness duration) and healthy comparators (n = 17) performed the auditory steady state response (ASSR) to 40, 30, and 20 Hz stimuli during MEG recordings. We estimated amplitude and PAC on the MEG ASSR source localized to the auditory cortices.

RESULTS

Gamma amplitude during 40-Hz ASSR exhibited a significant group by hemisphere interaction, with both patient groups showing reduced right hemisphere amplitude and no overall lateralization in contrast to the right hemisphere lateralization demonstrated in controls. We found significant PAC in the right auditory cortex during the 40-Hz entrainment condition relative to baseline, however, PAC did not differ significantly between groups.

CONCLUSIONS

In the current study, we demonstrated an apparent sparing of ASSR-related PAC across phases of the illness, in contrast with impaired cortical gamma oscillation amplitudes. The distinction between our PAC and evoked ASSR findings supports the notion of separate but interacting circuits for the generation and maintenance of sensory gamma oscillations. The apparent sparing of PAC in both early and chronic schizophrenia patients could imply that the neuropathology of schizophrenia differentially affects these mechanisms across different stages of the disease. Future studies should investigate the distinction between PAC during passive tasks and more cognitively demanding task such as working memory so that we can begin to understand the influence of schizophrenia neuropathology on the larger framework for modulating neurocomputational capacity.

Time-Frequency Based Phase-Amplitude Coupling Measure For Neuronal Oscillations

Oscillatory activity in the brain has been associated with a wide variety of cognitive processes including decision making, feedback processing, and working memory. The high temporal resolution provided by electroencephalography (EEG) enables the study of variation of oscillatory power and coupling across time. Various forms of neural synchrony across frequency bands have been suggested as the mechanism underlying neural binding. Recently, a considerable amount of work has focused on phase-amplitude coupling (PAC)– a form of cross-frequency coupling where the amplitude of a high frequency signal is modulated by the phase of low frequency oscillations. The existing methods for assessing PAC have some limitations including limited frequency resolution and sensitivity to noise, data length and sampling rate due to the inherent dependence on bandpass filtering. In this paper, we propose a new time-frequency based PAC (t-f PAC) measure that can address these issues. The proposed method relies on a complex time-frequency distribution, known as the Reduced Interference Distribution (RID)-Rihaczek distribution, to estimate both the phase and the envelope of low and high frequency oscillations, respectively. As such, it does not rely on bandpass filtering and possesses some of the desirable properties of time-frequency distributions such as high frequency resolution. The proposed technique is first evaluated for simulated data and then applied to an EEG speeded reaction task dataset. The results illustrate that the proposed time-frequency based PAC is more robust to varying signal parameters and provides a more accurate measure of coupling strength.

Phase-Amplitude Coupling of the Electroencephalogram in the Auditory Cortex in Schizophrenia

BACKGROUND

Cross-frequency interactions may coordinate neural circuits operating at different frequencies. While neural oscillations associated with particular circuits in schizophrenia (SZ) are impaired, few studies have examined cross-frequency interactions. Here we examined phase-amplitude coupling (PAC) in the electroencephalograms of individuals with SZ and healthy control subjects (HCs). We computed PAC during the baseline period of 40-Hz auditory steady-state stimulation and rest. We hypothesized that subjects with SZ would show abnormal theta/gamma coupling during stimulation, especially in the left auditory cortex, and coupling with high frequencies would be higher during stimulation than during rest.

METHODS

We reanalyzed data from 18 subjects with SZ and 18 HCs. Auditory cortex electroencephalogram activity was estimated using dipole source localization. PAC was computed using the debiased PAC measure, calculated with the generalized Morse wavelet transform. PAC clusters were identified using cluster-corrected permutation testing and interrogated in analyses of variance with correction for multiple tests.

RESULTS

Overall, coupling of high beta and gamma amplitude was higher during the auditory steady-state response, while alpha/beta PAC was higher during rest. Theta/alpha PAC was higher in subjects with SZ than in HCs. Theta/gamma PAC was lateralized to the left hemisphere in HCs but was not lateralized in subjects with SZ.

CONCLUSIONS

PAC involving high frequencies was state dependent and not abnormal in SZ. Increased theta/alpha PAC in subjects with SZ was consistent with other evidence of increased low-frequency activity. Hemispheric lateralization of theta/gamma PAC was reduced in subjects with SZ, consistent with evidence for left hemisphere auditory cortex abnormalities in subjects with SZ. PAC may reveal new insights into neural circuitry abnormalities in SZ and other neuropsychiatric disorders.

Neural and metabolic basis of dynamic resting state fMRI

Resting state fMRI (rsfMRI) as a technique showed much initial promise for use in psychiatric and neurological diseases where diagnosis and treatment were difficult. To realize this promise, many groups have moved towards examining "dynamic rsfMRI," which relies on the assumption that rsfMRI measurements on short time scales remain relevant to the underlying neural and metabolic activity. Many dynamic rsfMRI studies have demonstrated differences between clinical or behavioral groups beyond what static rsfMRI measured, suggesting a neurometabolic basis. Correlative studies combining dynamic rsfMRI and other physiological measurements have supported this. However, they also indicate multiple mechanisms and, if using correlation alone, it is difficult to separate cause and effect. Hypothesis-driven studies are needed, a few of which have begun to illuminate the underlying neurometabolic mechanisms that shape observed differences in dynamic rsfMRI. While the number of potential noise sources, potential actual neurometabolic sources, and methodological considerations can seem overwhelming, dynamic rsfMRI provides a rich opportunity in systems neuroscience. Even an incrementally better understanding of the neurometabolic basis of dynamic rsfMRI would expand rsfMRI's research and clinical utility, and the studies described herein take the first steps on that path forward.

A neural mass model of cross frequency coupling

Electrophysiological signals of cortical activity show a range of possible frequency and amplitude modulations, both within and across regions, collectively known as cross-frequency coupling. To investigate whether these modulations could be considered as manifestations of the same underlying mechanism, we developed a neural mass model. The model provides five out of the theoretically proposed six different coupling types. Within model components, slow and fast activity engage in phase-frequency coupling in conditions of low ambient noise level and with high noise level engage in phase-amplitude coupling. Between model components, these couplings can be coordinated via slow activity, giving rise to more complex modulations. The model, thus, provides a coherent account of cross-frequency coupling, both within and between components, with which regional and cross-regional frequency and amplitude modulations could be addressed.

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

Goal-directed motor behavior is associated with changes in patterns of rhythmic neuronal activity across widely distributed brain areas. In particular, movement initiation and execution are mediated by patterns of synchronization and desynchronization that occur concurrently across distinct frequency bands and across multiple motor cortical areas. To date, motor-related local oscillatory modulations have been predominantly examined by quantifying increases or suppressions in spectral power. However, beyond signal power, spectral properties such as phase and phase-amplitude coupling (PAC) have also been shown to carry information with regards to the oscillatory dynamics underlying motor processes. Yet, the distinct functional roles of phase, amplitude and PAC across the planning and execution of goal-directed motor behavior remain largely elusive. Here, we address this question with unprecedented resolution thanks to multi-site intracerebral EEG recordings in human subjects while they performed a delayed motor task. To compare the roles of phase, amplitude and PAC, we monitored intracranial brain signals from 748 sites across six medically intractable epilepsy patients at movement execution, and during the delay period where motor intention is present but execution is withheld. In particular, we used a machine-learning framework to identify the key contributions of various neuronal responses. We found a high degree of overlap between brain network patterns observed during planning and those present during execution. Prominent amplitude increases in the delta (2-4Hz) and high gamma (60-200Hz) bands were observed during both planning and execution. In contrast, motor alpha (8-13Hz) and beta (13-30Hz) power were suppressed during execution, but enhanced during the delay period. Interestingly, single-trial classification revealed that low-frequency phase information, rather than spectral power change, was the most discriminant feature in dissociating action from intention. Additionally, despite providing weaker decoding, PAC features led to statistically significant classification of motor states, particularly in anterior cingulate cortex and premotor brain areas. These results advance our understanding of the distinct and partly overlapping involvement of phase, amplitude and the coupling between them, in the neuronal mechanisms underlying motor intentions and executions.