Multivariate cross-frequency coupling via generalized eigendecomposition

Multivariate sharp-wave ripples in schizophrenia during awake state

AIMS

Schizophrenia (SZ) is a brain disorder characterized by psychotic symptoms and cognitive dysfunction. Recently, irregularities in sharp-wave ripples (SPW-Rs) have been reported in SZ. As SPW-Rs play a critical role in memory, their irregularities can cause psychotic symptoms and cognitive dysfunction in patients with SZ. In this study, we investigated the SPW-Rs in human SZ.

METHODS

We measured whole-brain activity using magnetoencephalography (MEG) in patients with SZ (n = 20) and sex- and age-matched healthy participants (n = 20) during open-eye rest. We identified SPW-Rs and analyzed their occurrence and time-frequency traits. Furthermore, we developed a novel multivariate analysis method, termed "ripple-gedMEG" to extract the global features of SPW-Rs. We also examined the association between SPW-Rs and brain state transitions. The outcomes of these analyses were modeled to predict the positive and negative syndrome scale (PANSS) scores of SZ.

RESULTS

We found that SPW-Rs in the SZ (1) occurred more frequently, (2) the delay of the coupling phase (3) appeared in different brain areas, (4) consisted of a less organized spatiotemporal pattern, and (5) were less involved in brain state transitions. Finally, some of the neural features associated with the SPW-Rs were found to be PANSS-positive, a pathological indicator of SZ. These results suggest that widespread but disorganized SPW-Rs underlies the symptoms of SZ.

CONCLUSION

We identified irregularities in SPW-Rs in SZ and confirmed that their alternations were strongly associated with SZ neuropathology. These results suggest a new direction for human SZ research.

A dynamic generative model can extract interpretable oscillatory components from multichannel neurophysiological recordings

Modern neurophysiological recordings are performed using multichannel sensor arrays that are able to record activity in an increasingly high number of channels numbering in the 100s to 1000s. Often, underlying lower-dimensional patterns of activity are responsible for the observed dynamics, but these representations are difficult to reliably identify using existing methods that attempt to summarize multivariate relationships in a post hoc manner from univariate analyses or using current blind source separation methods. While such methods can reveal appealing patterns of activity, determining the number of components to include, assessing their statistical significance, and interpreting them requires extensive manual intervention and subjective judgment in practice. These difficulties with component selection and interpretation occur in large part because these methods lack a generative model for the underlying spatio-temporal dynamics. Here, we describe a novel component analysis method anchored by a generative model where each source is described by a bio-physically inspired state-space representation. The parameters governing this representation readily capture the oscillatory temporal dynamics of the components, so we refer to it as oscillation component analysis. These parameters - the oscillatory properties, the component mixing weights at the sensors, and the number of oscillations - all are inferred in a data-driven fashion within a Bayesian framework employing an instance of the expectation maximization algorithm. We analyze high-dimensional electroencephalography and magnetoencephalography recordings from human studies to illustrate the potential utility of this method for neuroscience data.

Cross-Frequency Coupling as a Biomarker for Early Stroke Recovery

BACKGROUND

The application of neuroimaging-based biomarkers in stroke has enriched our understanding of post-stroke recovery mechanisms, including alterations in functional connectivity based on synchronous oscillatory activity across various cortical regions. Phase-amplitude coupling, a type of cross-frequency coupling, may provide additional mechanistic insight.

OBJECTIVE

To determine how the phase of prefrontal cortex delta (1-3 Hz) oscillatory activity mediates the amplitude of motor cortex beta (13-20 Hz) oscillations in individual's early post-stroke.

METHODS

Participants admitted to an inpatient rehabilitation facility completed resting and task-based EEG recordings and motor assessments around the time of admission and discharge along with structural neuroimaging. Unimpaired controls completed EEG procedures during a single visit. Mixed-effects linear models were performed to assess within- and between-group differences in delta-beta prefrontomotor coupling. Associations between coupling and motor status and injury were also determined.

RESULTS

Thirty individuals with stroke and 17 unimpaired controls participated. Coupling was greater during task versus rest conditions for all participants. Though coupling during affected extremity task performance decreased during hospitalization, coupling remained elevated at discharge compared to controls. Greater baseline coupling was associated with better motor status at admission and discharge and positively related to motor recovery. Coupling demonstrated both positive and negative associations with injury involving measures of lesion volume and overlap injury to anterior thalamic radiation, respectively.

CONCLUSIONS

This work highlights the utility of prefrontomotor cross-frequency coupling as a potential motor status and recovery biomarker in stroke. The frequency- and region-specific neurocircuitry featured in this work may also facilitate novel treatment strategies in stroke.

A dynamic generative model can extract interpretable oscillatory components from multichannel neurophysiological recordings

Modern neurophysiological recordings are performed using multichannel sensor arrays that are able to record activity in an increasingly high number of channels numbering in the 100's to 1000's. Often, underlying lower-dimensional patterns of activity are responsible for the observed dynamics, but these representations are difficult to reliably identify using existing methods that attempt to summarize multivariate relationships in a post-hoc manner from univariate analyses, or using current blind source separation methods. While such methods can reveal appealing patterns of activity, determining the number of components to include, assessing their statistical significance, and interpreting them requires extensive manual intervention and subjective judgement in practice. These difficulties with component selection and interpretation occur in large part because these methods lack a generative model for the underlying spatio-temporal dynamics. Here we describe a novel component analysis method anchored by a generative model where each source is described by a bio-physically inspired state space representation. The parameters governing this representation readily capture the oscillatory temporal dynamics of the components, so we refer to it as Oscillation Component Analysis (OCA). These parameters - the oscillatory properties, the component mixing weights at the sensors, and the number of oscillations - all are inferred in a data-driven fashion within a Bayesian framework employing an instance of the expectation maximization algorithm. We analyze high-dimensional electroencephalography and magnetoencephalography recordings from human studies to illustrate the potential utility of this method for neuroscience data.

Adaptive adjustment after conflict with group opinion: evidence from neural electrophysiology

Individuals inherently seek social consensus when making decisions or judgments. Previous studies have consistently indicated that dissenting group opinions are perceived as social conflict that demands attitude adjustment. However, the neurocognitive processes of attitude adjustment are unclear. In this electrophysiological study, participants were recruited to perform a face attractiveness judgment task. After forming their own judgment of a face, participants were informed of a purported group judgment (either consistent or inconsistent with their judgment), and then, critically, the same face was presented again. The neural responses to the second presented faces were measured. The second presented faces evoked a larger late positive potential after conflict with group opinions than those that did not conflict, suggesting that more motivated attention was allocated to stimulus. Moreover, faces elicited greater midfrontal theta (4-7 Hz) power after conflict with group opinions than after consistency with group opinions, suggesting that cognitive control was initiated to support attitude adjustment. Furthermore, the mixed-effects model revealed that single-trial theta power predicted behavioral change in the Conflict condition, but not in the No-Conflict condition. These findings provide novel insights into the neurocognitive processes underlying attitude adjustment, which is crucial to behavioral change during conformity.

A recruitment through coherence theory of working memory

The interactions between prefrontal cortex and other areas during working memory have been studied for decades. Here we outline a conceptual framework describing interactions between these areas during working memory, and review evidence for key elements of this model. We specifically suggest that a top-down signal sent from prefrontal to sensory areas drives oscillations in these areas. Spike timing within sensory areas becomes locked to these working-memory-driven oscillations, and the phase of spiking conveys information about the representation available within these areas. Downstream areas receiving these phase-locked spikes from sensory areas can recover this information via a combination of coherent oscillations and gating of input efficacy based on the phase of their local oscillations. Although the conceptual framework is based on prefrontal interactions with sensory areas during working memory, we also discuss the broader implications of this framework for flexible communication between brain areas in general.

Beyond rhythm - a framework for understanding the frequency spectrum of neural activity

Cognitive and behavioral processes are often accompanied by changes within well-defined frequency bands of the local field potential (LFP i.e., the voltage induced by neuronal activity). These changes are detectable in the frequency domain using the Fourier transform and are often interpreted as neuronal oscillations. However, aside some well-known exceptions, the processes underlying such changes are difficult to track in time, making their oscillatory nature hard to verify. In addition, many non-periodic neural processes can also have spectra that emphasize specific frequencies. Thus, the notion that spectral changes reflect oscillations is likely too restrictive. In this study, we use a simple yet versatile framework to understand the frequency spectra of neural recordings. Using simulations, we derive the Fourier spectra of periodic, quasi-periodic and non-periodic neural processes having diverse waveforms, illustrating how these attributes shape their spectral signatures. We then show how neural processes sum their energy in the local field potential in simulated and real-world recording scenarios. We find that the spectral power of neural processes is essentially determined by two aspects: (1) the distribution of neural events in time and (2) the waveform of the voltage induced by single neural events. Taken together, this work guides the interpretation of the Fourier spectrum of neural recordings and indicates that power increases in specific frequency bands do not necessarily reflect periodic neural activity.

Neurophysiological alterations in mice and humans carrying mutations in APP and PSEN1 genes

BACKGROUND

Studies in animal models of Alzheimer's disease (AD) have provided valuable insights into the molecular and cellular processes underlying neuronal network dysfunction. Whether and how AD-related neurophysiological alterations translate between mice and humans remains however uncertain.

METHODS

We characterized neurophysiological alterations in mice and humans carrying AD mutations in the APP and/or PSEN1 genes, focusing on early pre-symptomatic changes. Longitudinal local field potential recordings were performed in APP/PS1 mice and cross-sectional magnetoencephalography recordings in human APP and/or PSEN1 mutation carriers. All recordings were acquired in the left frontal cortex, parietal cortex, and hippocampus. Spectral power and functional connectivity were analyzed and compared with wildtype control mice and healthy age-matched human subjects.

RESULTS

APP/PS1 mice showed increased absolute power, especially at higher frequencies (beta and gamma) and predominantly between 3 and 6 moa. Relative power showed an overall shift from lower to higher frequencies over almost the entire recording period and across all three brain regions. Human mutation carriers, on the other hand, did not show changes in power except for an increase in relative theta power in the hippocampus. Mouse parietal cortex and hippocampal power spectra showed a characteristic peak at around 8 Hz which was not significantly altered in transgenic mice. Human power spectra showed a characteristic peak at around 9 Hz, the frequency of which was significantly reduced in mutation carriers. Significant alterations in functional connectivity were detected in theta, alpha, beta, and gamma frequency bands, but the exact frequency range and direction of change differed for APP/PS1 mice and human mutation carriers.

CONCLUSIONS

Both mice and humans carrying APP and/or PSEN1 mutations show abnormal neurophysiological activity, but several measures do not translate one-to-one between species. Alterations in absolute and relative power in mice should be interpreted with care and may be due to overexpression of amyloid in combination with the absence of tau pathology and cholinergic degeneration. Future studies should explore whether changes in brain activity in other AD mouse models, for instance, those also including tau pathology, provide better translation to the human AD continuum.

Midfrontal theta dynamics index the monitoring of postural stability

Stepping is a common strategy to recover postural stability and maintain upright balance. Postural perturbations have been linked to neuroelectrical markers such as the N1 potential and theta frequency dynamics. Here, we investigated the role of cortical midfrontal theta dynamics of balance monitoring, driven by balance perturbations at different initial standing postures. We recorded electroencephalography, electromyography, and motion tracking of human participants while they stood on a platform that delivered a range of forward and backward whole-body balance perturbations. The participants' postural threat was manipulated prior to the balance perturbation by instructing them to lean forward or backward while keeping their feet-in-place in response to the perturbation. We hypothesized that midfrontal theta dynamics index the engagement of a behavioral monitoring system and, therefore, that perturbation-induced theta power would be modulated by the initial leaning posture and perturbation intensity. Targeted spatial filtering in combination with mixed-effects modeling confirmed our hypothesis and revealed distinct modulations of theta power according to postural threat. Our results provide novel evidence that midfrontal theta dynamics subserve action monitoring of human postural balance. Understanding of cortical mechanisms of balance control is crucial for studying balance impairments related to aging and neurological conditions (e.g. stroke).

Ultra-high-resolution time-frequency analysis of EEG to characterise brain functional connectivity with the application in Alzheimer's disease

OBJECTIVE

This study aims to explore the potential of high-resolution brain functional connectivity based on electroencephalogram (EEG), a non-invasive low-cost technique, to be translated into a long-overdue biomarker and a diagnostic method for Alzheimer's disease (AD).

APPROACH

The paper proposes a novel ultra-high-resolution time-frequency nonlinear cross-spectrum method to construct a promising biomarker of AD pathophysiology. Specifically, using the peak frequency estimated from a Revised Hilbert-Huang Transformation cross-spectrum as a biomarker, the Support Vector Machine classifier is used to distinguish AD from healthy controls (HC).

MAIN RESULTS

With the combinations of the proposed biomarker and machine learning, we achieved a promising accuracy of 89%. The proposed method performs better than the wavelet cross-spectrum and other functional connectivity measures in the temporal or frequency domain, particularly in the Full, Delta and Alpha bands. Besides, a novel visualisation approach developed from topography is introduced to represent the brain functional connectivity, with which the difference between AD and HCs can be clearly displayed. The interconnections between posterior and other brain regions are obviously affected in AD.

SIGNIFICANCE

Those findings imply that the proposed RHHT approach could better track dynamic and nonlinear functional connectivity information, paving the way for the development of a novel diagnostic approach.

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.

Measuring Phase-Amplitude Coupling between Neural Oscillations of Different Frequencies via the Wasserstein Distance

BACKGROUND

Phase-amplitude coupling (PAC) is a key neuronal mechanism. Here, a novel method for quantifying PAC via the Wasserstein distance is presented.

NEW METHOD

The Wasserstein distance is an optimization algorithm for minimizing transportation cost and distance. For the first time, the author has applied this distance function to quantify PAC and named the Wasserstein Modulation Index (wMI). As the wMI accommodates the product of the amplitude value in each phase position and the coupling phase position, it allows for extraction of more detailed PAC features from the data.

RESULTS

The validity of the wMI calculations was examined using various simulation data, including sinusoidal and non-sinusoidal waves and empirical data sets. The current findings showed that the wMI is a more robust and stable index for quantifying PAC under various measuring conditions. Specifically, it can better reflect the timing of coupling and distinguish the shape of the coupling distribution than other measurements, both of which are the most significant parameters related to the functionality of PAC. Furthermore, the wMI is also suitable for many applications, such as more data-driven approaches and direct comparisons.

COMPARISON WITH EXISTING METHOD(S)

Compared with Euler-based PAC methods and the MI, the wMI is not easily affected by the non-sinusoidal nature of neural oscillation and the short data length and enables better reflection of the natures of PAC, such as the timing of coupling and the amplitude distribution in the phase plane, than the MI.

CONCLUSION

The wMI is expected to extract more detailed PAC characteristics, which could considerably contribute to the neuroscience field.

A tutorial on generalized eigendecomposition for denoising, contrast enhancement, and dimension reduction in multichannel electrophysiology

The goal of this paper is to present a theoretical and practical introduction to generalized eigendecomposition (GED), which is a robust and flexible framework used for dimension reduction and source separation in multichannel signal processing. In cognitive electrophysiology, GED is used to create spatial filters that maximize a researcher-specified contrast. For example, one may wish to exploit an assumption that different sources have different frequency content, or that sources vary in magnitude across experimental conditions. GED is fast and easy to compute, performs well in simulated and real data, and is easily adaptable to a variety of specific research goals. This paper introduces GED in a way that ties together myriad individual publications and applications of GED in electrophysiology, and provides sample MATLAB and Python code that can be tested and adapted. Practical considerations and issues that often arise in applications are discussed.

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.

Multivariate Analysis of Bivariate Phase-Amplitude Coupling in EEG Data Using Tensor Robust PCA

Cross-frequency coupling is emerging as a crucial mechanism that coordinates the integration of spectrally and spatially distributed neuronal oscillations. Recently, phase-amplitude coupling, a form of cross-frequency coupling, where the phase of a slow oscillation modulates the amplitude of a fast oscillation, has gained attention. Existing phase-amplitude coupling measures are mostly confined to either coupling within a region or between pairs of brain regions. Given the availability of multi-channel electroencephalography recordings, a multivariate analysis of phase amplitude coupling is needed to accurately quantify the coupling across multiple frequencies and brain regions. In the present work, we propose a tensor based approach, i.e., higher order robust principal component analysis, to identify response-evoked phase-amplitude coupling across multiple frequency bands and brain regions. Our experiments on both simulated and electroencephalography data demonstrate that the proposed multivariate phase-amplitude coupling method can capture the spatial and spectral dynamics of phase-amplitude coupling more accurately compared to existing methods. Accordingly, we posit that the proposed higher order robust principal component analysis based approach filters out the background phase-amplitude coupling activity and predominantly captures the event-related phase-amplitude coupling dynamics to provide insight into the spatially distributed brain networks across different frequency bands.

A Revised Hilbert-Huang Transformation to Track Non-Stationary Association of Electroencephalography Signals

The time-varying cross-spectrum method has been used to effectively study transient and dynamic brain functional connectivity between non-stationary electroencephalography (EEG) signals. Wavelet-based cross-spectrum is one of the most widely implemented methods, but it is limited by the spectral leakage caused by the finite length of the basic function that impacts the time and frequency resolutions. This paper proposes a new time-frequency brain functional connectivity analysis framework to track the non-stationary association of two EEG signals based on a Revised Hilbert-Huang Transform (RHHT). The framework can estimate the cross-spectrum of decomposed components of EEG, followed by a surrogate significance test. The results of two simulation examples demonstrate that, within a certain statistical confidence level, the proposed framework outperforms the wavelet-based method in terms of accuracy and time-frequency resolution. A case study on classifying epileptic patients and healthy controls using interictal seizure-free EEG data is also presented. The result suggests that the proposed method has the potential to better differentiate these two groups benefiting from the enhanced measure of dynamic time-frequency association.

The metastable brain associated with autistic-like traits of typically developing individuals

Metastability in the brain is thought to be a mechanism involved in the dynamic organization of cognitive and behavioral functions across multiple spatiotemporal scales. However, it is not clear how such organization is realized in underlying neural oscillations in a high-dimensional state space. It was shown that macroscopic oscillations often form phase-phase coupling (PPC) and phase-amplitude coupling (PAC), which result in synchronization and amplitude modulation, respectively, even without external stimuli. These oscillations can also make spontaneous transitions across synchronous states at rest. Using resting-state electroencephalographic signals and the autism-spectrum quotient scores acquired from healthy humans, we show experimental evidence that the PAC combined with PPC allows amplitude modulation to be transient, and that the metastable dynamics with this transient modulation is associated with autistic-like traits. In individuals with a longer attention span, such dynamics tended to show fewer transitions between states by forming delta-alpha PAC. We identified these states as two-dimensional metastable states that could share consistent patterns across individuals. Our findings suggest that the human brain dynamically organizes inter-individual differences in a hierarchy of macroscopic oscillations with multiple timescales by utilizing metastability.

The human brain organizes cognitive and behavioral functions dynamically. For decades, the dynamic organization of underlying neural oscillations has been a fundamental topic in neuroscience research. Even without external stimuli, macroscopic oscillations often form phase-phase coupling and phase-amplitude coupling (PAC) that result in synchronization and amplitude modulation, respectively, and can make spontaneous transitions across synchronous states at rest. Using resting-state electroencephalography signals acquired from healthy humans, we show evidence that these two neural couplings enable amplitude modulation to be transient, and that this transient modulation can be viewed as the transition among oscillatory states with different PAC strengths. We also demonstrate that such transition dynamics are associated with the ability to maintain attention to detail and to switch attention, as measured by autism-spectrum quotient scores. These individual dynamics were visualized as a trajectory among states with attracting tendencies, and involved consistent brain states across individuals. Our findings have significant implications for unraveling the variability in the individual brains showing typical and atypical development.

Inverse effects of time-on-task in task-related and task-unrelated theta activity

The phenomenon of mental fatigue has recently been investigated extensively by means of the EEG. Studies deploying spectral analysis consistently reported an increase of spectral power in the lower frequencies with increasing time-on-task, whereas event-related studies observed decreases in various measures related to task engagement and attentional resources. The results from these two lines of research cannot be aligned easily. (Frontal) theta power has been linked to cognitive control and was found to increase with time-on-task. In contrast, theoretical frameworks on mental fatigue suggest a decline in task-engagement as causal for the performance decline observed in mental fatigue. The goal of the present study was to investigate mental fatigue in time-frequency space using linear regression on single-trial data in order to obtain a better understanding about how time-on-task affects theta oscillatory activity. A data-driven analysis approach indicated an increase of alpha and theta power during the intertrial interval. In contrast, task-related theta activity declined. This reduction of stimulus-locked theta power may be interpreted as a reduction of task engagement with increasing mental fatigue. The increase of theta spectral power in the intertrial interval, moreover, could possibly be explained by an increased idling of cognitive control networks. Alternatively, it might be the case that the increase of theta power with time-on-task is a by-product an alpha power increase. As alpha peak frequency systematically decreases with time-on-task, the theta band might be affected as well.

Beta2 Oscillations in Hippocampal-Cortical Circuits During Novelty Detection

Novelty detection is a core feature of behavioral adaptation and involves cascades of neuronal responses-from initial evaluation of the stimulus to the encoding of new representations-resulting in the behavioral ability to respond to unexpected inputs. In the past decade, a new important novelty detection feature, beta2 (~20-30 Hz) oscillations, has been described in the hippocampus (HC). However, the interactions between beta2 and the hippocampal network are unknown, as well as the role-or even the presence-of beta2 in other areas involved with novelty detection. In this work, we combined multisite local field potential (LFP) recordings with novelty-related behavioral tasks in mice to describe the oscillatory dynamics associated with novelty detection in the CA1 region of the HC, parietal cortex, and mid-prefrontal cortex. We found that transient beta2 power increases were observed only during interaction with novel contexts and objects, but not with familiar contexts and objects. Also, robust theta-gamma phase-amplitude coupling was observed during the exploration of novel environments. Surprisingly, bursts of beta2 power had strong coupling with the phase of delta-range oscillations. Finally, the parietal and mid-frontal cortices had strong coherence with the HC in both theta and beta2. These results highlight the importance of beta2 oscillations in a larger hippocampal-cortical circuit, suggesting that beta2 plays a role in the mechanism for detecting and modulating behavioral adaptation to novelty.

Synchronization between Keyboard Typing and Neural Oscillations

Rhythmic neural activity synchronizes with certain rhythmic behaviors, such as breathing, sniffing, saccades, and speech. The extent to which neural oscillations synchronize with higher-level and more complex behaviors is largely unknown. Here, we investigated electrophysiological synchronization with keyboard typing, which is an omnipresent behavior daily engaged by an uncountably large number of people. Keyboard typing is rhythmic, with frequency characteristics roughly the same as neural oscillatory dynamics associated with cognitive control, notably through midfrontal theta (4-7 Hz) oscillations. We tested the hypothesis that synchronization occurs between typing and midfrontal theta and breaks down when errors are committed. Thirty healthy participants typed words and sentences on a keyboard without visual feedback, while EEG was recorded. Typing rhythmicity was investigated by interkeystroke interval analyses and by a kernel density estimation method. We used a multivariate spatial filtering technique to investigate frequency-specific synchronization between typing and neuronal oscillations. Our results demonstrate theta rhythmicity in typing (around 6.5 Hz) through the two different behavioral analyses. Synchronization between typing and neuronal oscillations occurred at frequencies ranging from 4 to 15 Hz, but to a larger extent for lower frequencies. However, peak synchronization frequency was idiosyncratic across participants, therefore not specific to theta nor to midfrontal regions, and correlated somewhat with peak typing frequency. Errors and trials associated with stronger cognitive control were not associated with changes in synchronization at any frequency. As a whole, this study shows that brain-behavior synchronization does occur during keyboard typing but is not specific to midfrontal theta.

Predicting Grating Orientations With Cross-Frequency Coupling and Least Absolute Shrinkage and Selection Operator in V1 and V4 of Rhesus Monkeys

Orientation selectivity, as an emergent property of neurons in the visual cortex, is of critical importance in the processing of visual information. Characterizing the orientation selectivity based on neuronal firing activities or local field potentials (LFPs) is a hot topic of current research. In this paper, we used cross-frequency coupling and least absolute shrinkage and selection operator (LASSO) to predict the grating orientations in V1 and V4 of two rhesus monkeys. The experimental data were recorded by utilizing two chronically implanted multi-electrode arrays, which were placed, respectively, in V1 and V4 of two rhesus monkeys performing a selective visual attention task. The phase-amplitude coupling (PAC) and amplitude-amplitude coupling (AAC) were employed to characterize the cross-frequency coupling of LFPs under sinusoidal grating stimuli with different orientations. Then, a LASSO logistic regression model was constructed to predict the grating orientation based on the strength of PAC and AAC. Moreover, the cross-validation method was used to evaluate the performance of the model. It was found that the average accuracy of the prediction based on the combination of PAC and AAC was 73.9%, which was higher than the predicting accuracy with PAC or AAC separately. In conclusion, a LASSO logistic regression model was introduced in this study, which can predict the grating orientations with relatively high accuracy by using PAC and AAC together. Our results suggest that the principle behind the LASSO model is probably an alternative direction to explore the mechanism for generating orientation selectivity.

Narrowband multivariate source separation for semi-blind discovery of experiment contrasts

BACKGROUND

Electrophysiological recordings contain mixtures of signals from distinct neural sources, impeding a straightforward interpretation of the sensor-level data. This mixing is particularly detrimental when distinct sources resonate in overlapping frequencies. Fortunately, the mixing is linear and instantaneous. Multivariate source separation methods may therefore successfully separate statistical sources, even with overlapping spatial distributions.

NEW METHOD

We demonstrate a feature-guided multivariate source separation method that is tuned to narrowband frequency content as well as binary condition differences. This method - comparison scanning generalized eigendecomposition, csGED - harnesses the covariance structure of multichannel data to find directions (i.e., eigenvectors) that maximally separate two subsets of data. To drive condition specificity and frequency specificity, our data subsets were taken from different task conditions and narrowband-filtered prior to applying GED.

RESULTS

To validate the method, we simulated MEG data in two conditions with shared noise characteristics and unique signal. csGED outperformed the best sensor at reconstructing the ground truth signals, even in the presence of large amounts of noise. We next applied csGED to a published empirical MEG dataset on visual perception vs. imagery. csGED identified sources in alpha, beta, and gamma bands, and successfully separated distinct networks in the same frequency band.

COMPARISON WITH EXISTING METHOD(S)

GED is a flexible feature-guided decomposition method that has previously successfully been applied. Our combined frequency- and condition-tuning is a novel adaptation that extends the power of GED in cognitive electrophysiology.

CONCLUSIONS

We demonstrate successful condition-specific source separation by applying csGED to simulated and empirical data.

Complex interplay between spectral harmonicity and different types of cross-frequency couplings in nonlinear oscillators and biologically plausible neural network models

Cross-frequency coupling (CFC) refers to the nonlinear interaction between oscillations in different frequency bands, and it is a rather ubiquitous phenomenon that has been observed in a variety of physical and biophysical systems. In particular, the coupling between the phase of slow oscillations and the amplitude of fast oscillations, referred as phase-amplitude coupling (PAC), has been intensively explored in the brain activity recorded from animals and humans. However, the interpretation of these CFC patterns remains challenging since harmonic spectral correlations characterizing nonsinusoidal oscillatory dynamics can act as a confounding factor. Specialized signal processing techniques are proposed to address the complex interplay between spectral harmonicity and different types of CFC, not restricted only to PAC. For this, we provide an in-depth characterization of the time locked index (TLI) as a tool aimed to efficiently quantify the harmonic content of noisy time series. It is shown that the proposed TLI measure is more robust and outperforms traditional phase coherence metrics (e.g., phase locking value, pairwise phase consistency) in several aspects. We found that a nonlinear oscillator under the effect of additive noise can produce spurious CFC with low spectral harmonic content. On the other hand, two coupled oscillatory dynamics with independent fundamental frequencies can produce true CFC with high spectral harmonic content via a rectification mechanism or other post-interaction nonlinear processing mechanisms. These results reveal a complex interplay between CFC and harmonicity emerging in the dynamics of biologically plausible neural network models and more generic nonlinear and parametric oscillators. We show that, contrary to what is usually assumed in the literature, the high harmonic content observed in nonsinusoidal oscillatory dynamics is neither a sufficient nor necessary condition to interpret the associated CFC patterns as epiphenomenal. There is mounting evidence suggesting that the combination of multimodal recordings, specialized signal processing techniques, and theoretical modeling is becoming a required step to completely understand CFC patterns observed in oscillatory rich dynamics of physical and biophysical systems.

A data-driven method to identify frequency boundaries in multichannel electrophysiology data

BACKGROUND

Electrophysiological recordings of the brain often exhibit neural oscillations, defined as narrowband bumps that deviate from the background power spectrum. These narrowband dynamics are grouped into frequency ranges, and the study of how activities in these ranges are related to cognition and disease is a major part of the neuroscience corpus. Frequency ranges are nearly always defined according to integer boundaries, such as 4-8 Hz for the theta band and 8-12 Hz for the alpha band.

NEW METHOD

A data-driven multivariate method is presented to identify empirical frequency boundaries based on clustering of spatiotemporal similarities across a range of frequencies. The method, termed gedBounds, identifies patterns in covariance matrices that maximally separate narrowband from broadband activity, and then identifies clusters in the correlation matrix of those spatial patterns over all frequencies, using the dbscan clustering algorithm. Those clusters are empirically derived frequency bands, from which boundaries can be extracted.

RESULTS

gedBounds recovers ground truth results in simulated data with high accuracy. The method was tested on EEG resting-state data from Parkinson's patients and control, and several features of the frequency components differed between patients and controls.

COMPARISON WITH EXISTING METHODS

The proposed method offers higher precision in defining subject-specific frequency boundaries compared to the current standard approach.

CONCLUSIONS

gedBounds can increase the precision and feature extraction of spectral dynamics in electrophysiology data.

Multiple Midfrontal Thetas Revealed by Source Separation of Simultaneous MEG and EEG

Theta-band (∼6 Hz) rhythmic activity within and over the medial PFC ("midfrontal theta") has been identified as a distinctive signature of "response conflict," the competition between multiple actions when only one action is goal-relevant. Midfrontal theta is traditionally conceptualized and analyzed under the assumption that it is a unitary signature of conflict that can be uniquely identified at one electrode (typically FCz). Here we recorded simultaneous MEG and EEG (total of 328 sensors) in 9 human subjects (7 female) and applied a feature-guided multivariate source-separation decomposition to determine whether conflict-related midfrontal theta is a unitary or multidimensional feature of the data. For each subject, a generalized eigendecomposition yielded spatial filters (components) that maximized the ratio between theta and broadband activity. Components were retained based on significance thresholding and midfrontal EEG topography. All of the subjects individually exhibited multiple (mean 5.89, SD 2.47) midfrontal components that contributed to sensor-level midfrontal theta power during the task. Component signals were temporally uncorrelated and asynchronous, suggesting that each midfrontal theta component was unique. Our findings call into question the dominant notion that midfrontal theta represents a unitary process. Instead, we suggest that midfrontal theta spans a multidimensional space, indicating multiple origins, but can manifest as a single feature at the sensor level because of signal mixing.SIGNIFICANCE STATEMENT "Midfrontal theta" is a rhythmic electrophysiological signature of the competition between multiple response options. Midfrontal theta is traditionally considered to reflect a single process. However, this assumption could be erroneous because of "mixing" (multiple sources contributing to the activity recorded at a single electrode). We investigated the dimensionality of midfrontal theta by applying advanced multivariate analysis methods to a multimodal MEG/EEG dataset. We identified multiple topographically overlapping neural sources that drove response conflict-related midfrontal theta. Midfrontal theta thus reflects multiple uncorrelated signals that manifest with similar EEG scalp projections. In addition to contributing to the cognitive control literature, we demonstrate both the feasibility and the necessity of signal demixing to understand the narrowband neural dynamics underlying cognitive processes.

Mesoscopic-scale functional networks in the primate amygdala

The primate amygdala performs multiple functions that may be related to the anatomical heterogeneity of its nuclei. Individual neurons with stimulus- and task-specific responses are not clustered in any of the nuclei, suggesting that single-units may be too-fine grained to shed light on the mesoscale organization of the amygdala. We have extracted from local field potentials recorded simultaneously from multiple locations within the primate (Macaca mulatta) amygdala spatially defined and statistically separable responses to visual, tactile, and auditory stimuli. A generalized eigendecomposition-based method of source separation isolated coactivity patterns, or components, that in neurophysiological terms correspond to putative subnetworks. Some component spatial patterns mapped onto the anatomical organization of the amygdala, while other components reflected integration across nuclei. These components differentiated between visual, tactile, and auditory stimuli suggesting the presence of functionally distinct parallel subnetworks.

Spectral fingerprints of correct vestibular discrimination of the intensity of body accelerations

Perceptual decision-making is a complex task that requires multiple processing steps performed by spatially distinct brain regions interacting in order to optimize perception and motor response. Most of our knowledge on these processes and interactions were derived from unimodal stimulations of the visual system which identified the lateral intraparietal area and the posterior parietal cortex as critical regions. Unlike the visual system, the vestibular system has no primary cortical areas and it is associated with separate multisensory areas within the temporo-parietal cortex with the parieto-insular vestibular cortex, PIVC, being the core region. The aim of the presented experiment was to investigate the transition from sensation to perception and to reveal the main structures of the cortical vestibular system involved in perceptual decision-making. Therefore, an EEG analysis was performed in 35 healthy subjects during linear whole-body accelerations of different intensities on a motor-driven motion platform (hexapod). We used a discrimination task in order to judge the intensity of the accelerations. Furthermore, we manipulated the expectation of the upcoming stimulus by indicating the probability (25%, 50%, 75%, 100%) of the motion direction. The analysis of the vestibular evoked potentials (VestEPs) showed that the decision-making process leads to a second positive peak (P2b) which was not observed in previous task-free experiments. The comparison of the estimated neural generators of the P2a and P2b components showed significant activity differences in the anterior cingulus, the parahippocampal and the middle temporal gyri. Taking into account the time courses of the P2 components, the physical properties of the stimuli, and the responses given by the subjects we conclude that the P2b likely reflects the transition from the processing of sensory information to perceptual evaluation. Analyzing the decision-uncertainty reported by the subjects, a persistent divergence of the time courses starting at 188 ​ms after the acceleration was found at electrode Pz. This finding demonstrated that meta-cognition by means of confidence estimation starts in parallel with the decision-making process itself. Further analyses in the time-frequency domain revealed that a correct classification of acceleration intensities correlated with an inter-trial phase clustering at electrode Cz and an inter-site phase clustering of theta oscillations over frontal, central, and parietal cortical areas. The sites where the phase clustering was observed corresponded to core decision-making brain areas known from neuroimaging studies in the visual domain.

Nonlinear interaction decomposition (NID): A method for separation of cross-frequency coupled sources in human brain

Cross-frequency coupling (CFC) between neuronal oscillations reflects an integration of spatially and spectrally distributed information in the brain. Here, we propose a novel framework for detecting such interactions in Magneto- and Electroencephalography (MEG/EEG), which we refer to as Nonlinear Interaction Decomposition (NID). In contrast to all previous methods for separation of cross-frequency (CF) sources in the brain, we propose that the extraction of nonlinearly interacting oscillations can be based on the statistical properties of their linear mixtures. The main idea of NID is that nonlinearly coupled brain oscillations can be mixed in such a way that the resulting linear mixture has a non-Gaussian distribution. We evaluate this argument analytically for amplitude-modulated narrow-band oscillations which are either phase-phase or amplitude-amplitude CF coupled. We validated NID extensively with simulated EEG obtained with realistic head modelling. The method extracted nonlinearly interacting components reliably even at SNRs as small as -15 dB. Additionally, we applied NID to the resting-state EEG of 81 subjects to characterize CF phase-phase coupling between alpha and beta oscillations. The extracted sources were located in temporal, parietal and frontal areas, demonstrating the existence of diverse local and distant nonlinear interactions in resting-state EEG data. All codes are available publicly via GitHub.

Cross-Channel Phase-Amplitude Transfer Entropy Conceptualize Long-Range Transmission in sleep: a case study

A causal algorithmic framework quantifying cross-channel phase-amplitude transfer entropy was proposed to measure long-range transmission dynamics between frontal and occipital brain areas during sleep. To this end, a noise-assisted multivariate empirical mode decomposition method was used to guarantee the consistent scales across multivariate signals. On the other side, transfer entropy was applied to measure information transfers from a low-frequency phase to a high-frequency amplitude across different brain regions. Our results showed δ phase may modulate either θ or α amplitude. The frontal cortex transferred information to the occipital brain area more than its inverse direction during Awake and N3 sleep stages, whereas N1 was more likely of serving as a transition state. On the other side, the information flow transferred from the occipital area to the frontal cortex surpassed its inverse flow in the N2 sleep stage. The proposed causal algorithmic framework facilitated identifying information flow and driving force across brain regions in sleep.

Midfrontal theta phase coordinates behaviorally relevant brain computations during cognitive control

Neural oscillations are thought to provide a cyclic time frame for orchestrating brain computations. Following this assumption, midfrontal theta oscillations have recently been proposed to temporally organize brain computations during conflict processing. Using a multivariate analysis approach, we show that brain-behavior relationships during conflict tasks are modulated according to the phase of ongoing endogenous midfrontal theta oscillations recorded by scalp EEG. We found reproducible results in two independent datasets, using two different conflict tasks: brain-behavior relationships (correlation between reaction time and theta power) were theta phase-dependent in a subject-specific manner, and these "behaviorally optimal" theta phases were also associated with fronto-parietal cross-frequency dynamics emerging as theta phase-locked beta power bursts. These effects were present regardless of the strength of conflict. Thus, these results provide empirical evidence that midfrontal theta oscillations are involved in cyclically orchestrating brain computations likely related to response execution during the tasks rather than purely related to conflict processing. More generally, this study supports the hypothesis that phase-based computation is an important mechanism giving rise to cognitive processing.

A statistical framework to assess cross-frequency coupling while accounting for confounding analysis effects

Cross frequency coupling (CFC) is emerging as a fundamental feature of brain activity, correlated with brain function and dysfunction. Many different types of CFC have been identified through application of numerous data analysis methods, each developed to characterize a specific CFC type. Choosing an inappropriate method weakens statistical power and introduces opportunities for confounding effects. To address this, we propose a statistical modeling framework to estimate high frequency amplitude as a function of both the low frequency amplitude and low frequency phase; the result is a measure of phase-amplitude coupling that accounts for changes in the low frequency amplitude. We show in simulations that the proposed method successfully detects CFC between the low frequency phase or amplitude and the high frequency amplitude, and outperforms an existing method in biologically-motivated examples. Applying the method to in vivo data, we illustrate examples of CFC during a seizure and in response to electrical stimuli.

Frontal theta and posterior alpha in resting EEG: A critical examination of convergent and discriminant validity

Prior research has identified two resting EEG biomarkers with potential for predicting functional outcomes in depression: theta current density in frontal brain regions (especially rostral anterior cingulate cortex) and alpha power over posterior scalp regions. As little is known about the discriminant and convergent validity of these putative biomarkers, a thorough evaluation of these psychometric properties was conducted toward the goal of improving clinical utility of these markers. Resting 71-channel EEG recorded from 35 healthy adults at two sessions (1-week retest) were used to systematically compare different quantification techniques for theta and alpha sources at scalp (surface Laplacian or current source density [CSD]) and brain (distributed inverse; exact low resolution electromagnetic tomography [eLORETA]) level. Signal quality was evaluated with signal-to-noise ratio, participant-level spectra, and frequency PCA covariance decomposition. Convergent and discriminant validity were assessed within a multitrait-multimethod framework. Posterior alpha was reliably identified as two spectral components, each with unique spatial patterns and condition effects (eyes open/closed), high signal quality, and good convergent and discriminant validity. In contrast, frontal theta was characterized by one low-variance component, low signal quality, lack of a distinct spectral peak, and mixed validity. Correlations between candidate biomarkers suggest that posterior alpha components constitute reliable, convergent, and discriminant biometrics in healthy adults. Component-based identification of spectral activity (CSD/eLORETA-fPCA) was superior to fixed, a priori frequency bands. Improved quantification and conceptualization of frontal theta is necessary to determine clinical utility.

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.

Non-linear auto-regressive models for cross-frequency coupling in neural time series

We address the issue of reliably detecting and quantifying cross-frequency coupling (CFC) in neural time series. Based on non-linear auto-regressive models, the proposed method provides a generative and parametric model of the time-varying spectral content of the signals. As this method models the entire spectrum simultaneously, it avoids the pitfalls related to incorrect filtering or the use of the Hilbert transform on wide-band signals. As the model is probabilistic, it also provides a score of the model “goodness of fit” via the likelihood, enabling easy and legitimate model selection and parameter comparison; this data-driven feature is unique to our model-based approach. Using three datasets obtained with invasive neurophysiological recordings in humans and rodents, we demonstrate that these models are able to replicate previous results obtained with other metrics, but also reveal new insights such as the influence of the amplitude of the slow oscillation. Using simulations, we demonstrate that our parametric method can reveal neural couplings with shorter signals than non-parametric methods. We also show how the likelihood can be used to find optimal filtering parameters, suggesting new properties on the spectrum of the driving signal, but also to estimate the optimal delay between the coupled signals, enabling a directionality estimation in the coupling.

Neural oscillations synchronize information across brain areas at various anatomical and temporal scales. Of particular relevance, slow fluctuations of brain activity have been shown to affect high frequency neural activity, by regulating the excitability level of neural populations. Such cross-frequency-coupling can take several forms. In the most frequently observed type, the power of high frequency activity is time-locked to a specific phase of slow frequency oscillations, yielding phase-amplitude-coupling (PAC). Even when readily observed in neural recordings, such non-linear coupling is particularly challenging to formally characterize. Typically, neuroscientists use band-pass filtering and Hilbert transforms with ad-hoc correlations. Here, we explicitly address current limitations and propose an alternative probabilistic signal modeling approach, for which statistical inference is fast and well-posed. To statistically model PAC, we propose to use non-linear auto-regressive models which estimate the spectral modulation of a signal conditionally to a driving signal. This conditional spectral analysis enables easy model selection and clear hypothesis-testing by using the likelihood of a given model. We demonstrate the advantage of the model-based approach on three datasets acquired in rats and in humans. We further provide novel neuroscientific insights on previously reported PAC phenomena, capturing two mechanisms in PAC: influence of amplitude and directionality estimation.

Disharmony in neural oscillations

Cross-frequency phase coupling (PPC) may play an important role in neural processing and cognition. However, a new study unveils a statistical bias in how PPC is detected in neural recordings, questions prior evidence for PPC in hippocampus, and shows PPC tests are dramatically flawed by their confounds with oscillation harmonics.