Analysis of time-variant quadratic phase couplings in the tracé alternant EEG by recursive estimation of 3rd-order time–frequency distributions

Time-frequency Techniques in Biomedical Signal Analysis

Objectives: This review outlines the method -ological fundamentals of the most frequently used non-parametric time-frequency analysis techniques in biomedicine and their main properties, as well as providing decision aids concerning their applications.

Methods: The short-term Fourier transform (STFT), the Gabor transform (GT), the S-transform (ST), the continuous Morlet wavelet transform (CMWT), and the Hilbert transform (HT) are introduced as linear transforms by using a unified concept of the time-frequency representation which is based on a standardized analytic signal. The Wigner-Ville dis -tribution (WVD) serves as an example of the ‘quadratic transforms’ class. The combination of WVD and GT with the matching pursuit (MP) decomposition and that of the HT with the empirical mode decomposition (EMD) are explained; these belong to the class of signal-adaptive approaches.

Results: Similarities between linear transforms are demonstrated and differences with regard to the time-frequency resolution and interference (cross) terms are presented in detail. By means of simulated signals the effects of different time-frequency resolutions of the GT, CMWT, and WVD as well as the resolution-related properties of the inter -ference (cross) terms are shown. The method-inherent drawbacks and their consequences for the application of the time-frequency techniques are demonstrated by instantaneous amplitude, frequency and phase measures and related time-frequency representations (spectrogram, scalogram, time-frequency distribution, phase-locking maps) of measured magnetoencephalographic (MEG) signals.

Conclusions: The appropriate selection of a method and its parameter settings will ensure readability of the time-frequency representations and reliability of results. When the time-frequency characteristics of a signal strongly correspond with the time-frequency resolution of the analysis then a method may be considered ‘optimal’. The MP-based signal-adaptive approaches are preferred as these provide an appropriate time-frequency resolution for all frequencies while simultaneously reducing interference (cross) terms.

Non-linear Analysis of Scalp EEG by Using Bispectra: The Effect of the Reference Choice

Bispectral analysis is a signal processing technique that makes it possible to capture the non-linear and non-Gaussian properties of the EEG signals. It has found various applications in EEG research and clinical practice, including the assessment of anesthetic depth, the identification of epileptic seizures, and more recently, the evaluation of non-linear cross-frequency brain functional connectivity. However, the validity and reliability of the indices drawn from bispectral analysis of EEG signals are potentially biased by the use of a non-neutral EEG reference. The present study aims at investigating the effects of the reference choice on the analysis of the non-linear features of EEG signals through bicoherence, as well as on the estimation of cross-frequency EEG connectivity through two different non-linear measures, i.e., the cross-bicoherence and the antisymmetric cross-bicoherence. To this end, four commonly used reference schemes were considered: the vertex electrode (Cz), the digitally linked mastoids, the average reference, and the Reference Electrode Standardization Technique (REST). The reference effects were assessed both in simulations and in a real EEG experiment. The simulations allowed to investigated: (i) the effects of the electrode density on the performance of the above references in the estimation of bispectral measures; and (ii) the effects of the head model accuracy in the performance of the REST. For real data, the EEG signals recorded from 10 subjects during eyes open resting state were examined, and the distortions induced by the reference choice in the patterns of alpha-beta bicoherence, cross-bicoherence, and antisymmetric cross-bicoherence were assessed. The results showed significant differences in the findings depending on the chosen reference, with the REST providing superior performance than all the other references in approximating the ideal neutral reference. In conclusion, this study highlights the importance of considering the effects of the reference choice in the interpretation and comparison of the results of bispectral analysis of scalp EEG.

Bispectral pairwise interacting source analysis for identifying systems of cross-frequency interacting brain sources from electroencephalographic or magnetoencephalographic signals

Brain cognitive functions arise through the coordinated activity of several brain regions, which actually form complex dynamical systems operating at multiple frequencies. These systems often consist of interacting subsystems, whose characterization is of importance for a complete understanding of the brain interaction processes. To address this issue, we present a technique, namely the bispectral pairwise interacting source analysis (biPISA), for analyzing systems of cross-frequency interacting brain sources when multichannel electroencephalographic (EEG) or magnetoencephalographic (MEG) data are available. Specifically, the biPISA makes it possible to identify one or many subsystems of cross-frequency interacting sources by decomposing the antisymmetric components of the cross-bispectra between EEG or MEG signals, based on the assumption that interactions are pairwise. Thanks to the properties of the antisymmetric components of the cross-bispectra, biPISA is also robust to spurious interactions arising from mixing artifacts, i.e., volume conduction or field spread, which always affect EEG or MEG functional connectivity estimates. This method is an extension of the pairwise interacting source analysis (PISA), which was originally introduced for investigating interactions at the same frequency, to the study of cross-frequency interactions. The effectiveness of this approach is demonstrated in simulations for up to three interacting source pairs and for real MEG recordings of spontaneous brain activity. Simulations show that the performances of biPISA in estimating the phase difference between the interacting sources are affected by the increasing level of noise rather than by the number of the interacting subsystems. The analysis of real MEG data reveals an interaction between two pairs of sources of central mu and beta rhythms, localizing in the proximity of the left and right central sulci.

Third order spectral analysis robust to mixing artifacts for mapping cross-frequency interactions in EEG/MEG

We present a novel approach to the third order spectral analysis, commonly called bispectral analysis, of electroencephalographic (EEG) and magnetoencephalographic (MEG) data for studying cross-frequency functional brain connectivity. The main obstacle in estimating functional connectivity from EEG and MEG measurements lies in the signals being a largely unknown mixture of the activities of the underlying brain sources. This often constitutes a severe confounder and heavily affects the detection of brain source interactions. To overcome this problem, we previously developed metrics based on the properties of the imaginary part of coherency. Here, we generalize these properties from the linear to the nonlinear case. Specifically, we propose a metric based on an antisymmetric combination of cross-bispectra, which we demonstrate to be robust to mixing artifacts. Moreover, our metric provides complex-valued quantities that give the opportunity to study phase relationships between brain sources. The effectiveness of the method is first demonstrated on simulated EEG data. The proposed approach shows a reduced sensitivity to mixing artifacts when compared with a traditional bispectral metric. It also exhibits a better performance in extracting phase relationships between sources than the imaginary part of the cross-spectrum for delayed interactions. The method is then applied to real EEG data recorded during resting state. A cross-frequency interaction is observed between brain sources at 10Hz and 20Hz, i.e., for alpha and beta rhythms. This interaction is then projected from signal to source level by using a fit-based procedure. This approach highlights a 10-20Hz dominant interaction localized in an occipito-parieto-central network.

Time-variant analysis of linear and non-linear phase couplings of and between frequency components of EEG burst patterns in full-term newborns

Time-variant (tv) phase-locking and synchronization characteristics of and between low-frequency (≤ 1.5 Hz) and high-frequency EEG oscillations (≥ 3.5 Hz) of the tracé alternant (TA) pattern in full-term newborns have been quantified to explore the origin of quadratic phase coupling (QPC, as non-linear phase coupling measure) between the frequency ranges 1 - 1.5 Hz ⇔ 3.5 - 4.5 Hz, which characterize the specific interactions of oscillations during the TA's burst activity. Using the Gabor transformation two measures of linear phase coupling, the phase-locking index (PLI) and the n:m phase synchronization index (PSI) have been determined. Phase-locking within the frequency ranges 1 -1.5 Hz and 3.5 - 4.5 Hz, and synchronization between both frequency ranges exists. These phase characteristics are significant 2 sec after burst onset and are associated with the maximum-values of the QPC(1 - 1.5 Hz ⇔ 3.5 - 4.5 Hz) which demonstrates that a specific neuronal coordination between the dynamics of phases and of amplitude-frequency dependencies must be underlying.

Coordination of the EEG and the heart rate of preterm neonates during quiet sleep

Aim of this study was to confirm that EEG bursts are associated with heart rate (HR) accelerations, and to investigate the synchronicity between quadratic phase couplings (QPC) courses of the EEG and HR before and during burst activity during quiet sleep in preterm newborns. The time-courses of QPC between frequency components of the EEG ([0.25-1.0 Hz]<-->[4.0-6.0 Hz]) as well as between the Mayer-Traube-Hering (MTH) wave and the frequency component of the HR associated to the respiratory sinus arrhythmia (RSA) ([0.02-0.15 Hz]<-->[0.4-1.5 Hz]) were investigated in five preterm neonates. During quiet sleep, the EEG alternates between burst and interburst activity. The burst onsets were used to trigger an averaging procedure for the EEG, HR, and QPC courses. It can be demonstrated that the envelopes of the EEG rise after the burst onset accompanied by an acceleration of HR before or at the burst maximum. The QPC courses show that the HR's QPC increases before or at the burst onset whereas the increase of the EEG's QPC is delayed. The synchronous changes of EEG and HR as well as of the corresponding QPC courses indicate a coupling between cortical, thalamocortical and neurovegetative brain structures. Such a coupling might be mediated by the MTH waves in the blood pressure.