The use of stationarity and nonstationarity in the detection and analysis of neural oscillations

Experimental evaluation of methods for real-time EEG phase-specific transcranial magnetic stimulation

OBJECTIVE

Real-time approaches for transcranial magnetic stimulation (TMS) based on a specific EEG phase are a promising avenue for more precise neuromodulation interventions. However, optimal approaches to reliably extract the EEG phase in a frequency band of interest to inform TMS are still to be identified. Here, we implement a new real-time phase detection method for closed-loop EEG-TMS for robust phase extraction. We compare this algorithm with state-of-the-art methods and evaluate its performance both in silico and experimentally.

APPROACH

We propose a new robust algorithm (Educated Temporal Prediction) for delivering real-time EEG phase-specific stimulation based on short prerecorded EEG training data. This method estimates the interpeak period from a training period and applies a bias correction to predict future peaks. We compare the accuracy and computation speed of the ETP algorithm with two existing methods (Fourier based, Autoregressive Prediction) using prerecorded resting EEG data and real-time experiments.

MAIN RESULTS

We found that Educated Temporal Prediction performs with higher accuracy than Fourier-based or Autoregressive methods both in silico and in vivo while being computationally more efficient. Further, we document the dependency of the EEG signal-to-noise ratio (SNR) on algorithm accuracy across all algorithms.

SIGNIFICANCE

Our results give important insights for real-time EEG-TMS technical development as well as experimental design. Due to its robustness and computational efficiency, our method can find broad use in experimental research or clinical applications. Through open sharing of code for all three methods, we enable broad access of TMS-EEG real-time algorithms to the community.

The dynamics of functional connectivity in neocortical focal epilepsy

Focal epilepsy is characterised by paroxysmal events, reflecting changes in underlying local brain networks. To capture brain network activity at the maximal temporal resolution of the acquired functional magnetic resonance imaging (fMRI) data, we have previously developed a novel analysis framework called Dynamic Regional Phase Synchrony (DRePS). DRePS measures instantaneous mean phase coherence within neighbourhoods of brain voxels. We use it here to examine how the dynamics of the functional connections of regional brain networks are altered in neocortical focal epilepsy. Using task-free fMRI data from 21 subjects with focal epilepsy and 21 healthy controls, we calculated the power spectral density of DRePS, which is a measure of signal variability in local connectivity estimates. Whole-brain averaged power spectral density of DRePS, or signal variability of local connectivity, was significantly higher in epilepsy subjects compared to healthy controls. Maximal increase in DRePS spectral power was seen in bilateral inferior frontal cortices, ipsilateral mid-cingulate gyrus, superior temporal gyrus, caudate head, and contralateral cerebellum. Our results provide further evidence of common brain abnormalities across people with focal epilepsy. We postulate that dynamic changes in specific cortical brain areas may help maintain brain function in the presence of pathological epileptiform network activity in neocortical focal epilepsy.

•

Focal epilepsy is a disease defined by its inherent paroxysmal brain changes.

•

However, most fMRI connectivity methods do not capture dynamic brain network activity.

•

We use our recently developed method, DRePS, to identify abnormal dynamic brain networks in focal epilepsy.

•

Fluctuations of dynamic brain network activity were increased in the frontal-temporal cortex, cingulum and cerebellum, in focal epilepsy.

•

These regions may be involved in (re-)gaining cortical homeostasis in an otherwise abnormal brain.

Elevated amygdala response to faces and gaze aversion in autism spectrum disorder

Autism spectrum disorders (ASD) are often associated with impairments in judgment of facial expressions. This impairment is often accompanied by diminished eye contact and atypical amygdala responses to face stimuli. The current study used a within-subjects design to examine the effects of natural viewing and an experimental eye-gaze manipulation on amygdala responses to faces. Individuals with ASD showed less gaze toward the eye region of faces relative to a control group. Among individuals with ASD, reduced eye gaze was associated with higher threat ratings of neutral faces. Amygdala signal was elevated in the ASD group relative to controls. This elevated response was further potentiated by experimentally manipulating gaze to the eye region. Potentiation by the gaze manipulation was largest for those individuals who exhibited the least amount of naturally occurring gaze toward the eye region and was associated with their subjective threat ratings. Effects were largest for neutral faces, highlighting the importance of examining neutral faces in the pathophysiology of autism and questioning their use as control stimuli with this population. Overall, our findings provide support for the notion that gaze direction modulates affective response to faces in ASD.

The comodulation measure of neuronal oscillations with general harmonic wavelet bicoherence and application to sleep analysis

Brain functions are related to neuronal networks of different sizes and distribution, and neuronal networks of different sizes oscillate at different frequencies. Thus the synchronization of neuronal networks is often reflected by cross-frequency interaction. The description of this cross-frequency interaction is therefore a crucial issue in understanding the modulation mechanisms between neuronal populations. A number of different kinds of interaction between frequencies have been reported. In this paper, we develop a general harmonic wavelet transform based bicoherence using a phase randomization method. This allows us to measure the comodulation of oscillations between different frequency bands in neuronal populations. The performance of the method is evaluated by a simulation study. The results show that the improved wavelet bicoherence method can detect a reliable phase coupling value, and also identify zero bicoherence for waves that are not phase-coupled. Spurious bicoherences can be effectively eliminated through the phase randomization method. Finally, this method is applied to electrocorticogram data recorded from rats during transitions between slow-wave sleep, rapid-eye movement sleep and waking. The phase coupling in rapid-eye movement sleep is statistically lower than that during slow-wave sleep, and slightly less than those in the wakeful state. The degree of phase coupling in rapid-eye movement sleep after slow-wave sleep is greater than in rapid-eye movement sleep prior to waking. This method could be applied to investigate the cross-frequency interactions in other physiological signals.

A space-frequency analysis of MEG background processes

In MEG source localization the estimated source parameters will be more reliable when the spatiotemporal covariance of the noise and background activity is taken into account. Since this covariance is in general too large to estimate based on the data and to invert efficiently, different parametrizations have been proposed in the literature. These models can be seen as special cases of the general decomposition of the covariance into a sum of Kronecker products of spatial matrices X(n) and temporal matrices T(n) (Van Loan, 2000). In this study we investigate the assumption of the matrices T(n) being Toeplitz. If so, the covariance matrix in the space-frequency domain will have an approximate block-diagonal structure, facilitating inversion, which is a prerequisite for source localization. In this study we address the question whether the Toeplitz approximation is valid for data sets obtained in visual evoked field, auditory evoked field, somatosensory evoked field experiments and data sets containing spontaneous activity. It turns out that on average 87% is in the block-diagonal of the sample covariance, which is close to the values obtained for real Toeplitz matrices T(n). This implies that the space-frequency domain is very interesting for source localization since the major part of the entire covariance can be incorporated in that domain straightforwardly. Finally, the two major processes in the background activity are characterized in terms of their spatial and frequency patterns, yielding a focal and a non-focal pattern in 8 of 10 data sets analyzed in this study. The focal pattern represents the alpha frequency at parieto-occipital areas, whereas the non-focal pattern is more widespread both in space and in frequency.

Spatial-temporal analysis of fetal bio-magnetic signals

Non-invasive technique such as magneto-encephalography (MEG), initially pioneered to study human brain signals, has found many other applications in medicine. SQUID(1) Array for Reproductive Assessment (SARA) is a unique non-invasive scanning-device developed at the University of Arkansas for Medical Sciences (UAMS) that can detect fetal brain and other signals. The fetal magneto-encephalography (fMEG) signals often have many bio-magnetic signals mixed in. Examples include the movement of the fetus or muscle contraction of the mother. As a result, the recorded signals may show unexpected patterns, other than the target signal of interest. These "interventions" make it difficult for a physician to assess the exact fetal condition, including its response to various stimuli. We propose using intervention analysis and spatial-temporal auto-regressive moving-average (STARMA) modeling to address the problem. STARMA is a statistical method that examines the relationship between the current observations as a linear combination of past observations as well as observations at neighboring sensors. Through intervention analysis, the change in a pattern due to "interfering" signals can be accounted for. When these interferences are "removed," the end product is a "template" time series, or a typical signal from the target of interest. In this research, a "universal" template is obtained. The template is then used to detect intervention in other datasets by the method of template matching. By this method, it is possible to detect if there is an intervention in any dataset. It will assist physicians in monitoring the actual signal generated by fetal brain and other organs of interest.