A pulse artifact removal method considering artifact variations in the simultaneous recording of EEG and fMRI

Preservation of EEG spectral power features during simultaneous EEG-fMRI

Introduction

Electroencephalographic (EEG) data quality is severely compromised when recorded inside the magnetic resonance (MR) environment. Here we characterized the impact of the ballistocardiographic (BCG) artifact on resting-state EEG spectral properties and compared the effectiveness of seven common BCG correction methods to preserve EEG spectral features. We also assessed if these methods retained posterior alpha power reactivity to an eyes closure-opening (EC-EO) task and compared the results from EEG-informed fMRI analysis using different BCG correction approaches.

Method

Electroencephalographic data from 20 healthy young adults were recorded outside the MR environment and during simultaneous fMRI acquisition. The gradient artifact was effectively removed from EEG-fMRI acquisitions using Average Artifact Subtraction (AAS). The BCG artifact was corrected with seven methods: AAS, Optimal Basis Set (OBS), Independent Component Analysis (ICA), OBS followed by ICA, AAS followed by ICA, PROJIC-AAS and PROJIC-OBS. EEG signal preservation was assessed by comparing the spectral power of traditional frequency bands from the corrected rs-EEG-fMRI data with the data recorded outside the scanner. We then assessed the preservation of posterior alpha functional reactivity by computing the ratio between the EC and EO conditions during the EC-EO task. EEG-informed fMRI analysis of the EC-EO task was performed using alpha power-derived BOLD signal predictors obtained from the EEG signals corrected with different methods.

Results

The BCG artifact caused significant distortions (increased absolute power, altered relative power) across all frequency bands. Artifact residuals/signal losses were present after applying all correction methods. The EEG reactivity to the EC-EO task was better preserved with ICA-based correction approaches, particularly when using ICA feature extraction to isolate alpha power fluctuations, which allowed to accurately predict hemodynamic signal fluctuations during the EEG-informed fMRI analysis.

Discussion

Current software solutions for the BCG artifact problem offer limited efficiency to preserve the EEG spectral power properties using this particular EEG setup. The state-of-the-art approaches tested here can be further refined and should be combined with hardware implementations to better preserve EEG signal properties during simultaneous EEG-fMRI. Existing and novel BCG artifact correction methods should be validated by evaluating signal preservation of both ERPs and spontaneous EEG spectral power.

EEG-fMRI: Ballistocardiogram Artifact Reduction by Surrogate Method for Improved Source Localization

For the analysis of simultaneous EEG-fMRI recordings, it is vital to use effective artifact removal tools. This applies in particular to the ballistocardiogram (BCG) artifact which is difficult to remove without distorting signals of interest related to brain activity. Here, we documented the use of surrogate source models to separate the artifact-related signals from brain signals with minimal distortion of the brain activity of interest. The artifact topographies used for surrogate separation were created automatically using principal components analysis (PCA-S) or by manual selection of artifact components utilizing independent components analysis (ICA-S). Using real resting-state data from 55 subjects superimposed with simulated auditory evoked potentials (AEP), both approaches were compared with three established BCG artifact removal methods: Blind Source Separation (BSS), Optimal Basis Set (OBS), and a mixture of both (OBS-ICA). Each method was evaluated for its applicability for ERP and source analysis using the following criteria: the number of events surviving artifact threshold scans, signal-to-noise ratio (SNR), error of source localization, and signal variance explained by the dipolar model. Using these criteria, PCA-S and ICA-S fared best overall, with highly significant differences to the established methods, especially in source localization. The PCA-S approach was also applied to a single subject Berger experiment performed in the MRI scanner. Overall, the removal of BCG artifacts by the surrogate methods provides a substantial improvement for the analysis of simultaneous EEG-fMRI data compared to the established methods.

Artifact Reduction in Simultaneous EEG-fMRI: A Systematic Review of Methods and Contemporary Usage

Simultaneous electroencephalography-functional MRI (EEG-fMRI) is a technique that combines temporal (largely from EEG) and spatial (largely from fMRI) indicators of brain dynamics. It is useful for understanding neuronal activity during many different event types, including spontaneous epileptic discharges, the activity of sleep stages, and activity evoked by external stimuli and decision-making tasks. However, EEG recorded during fMRI is subject to imaging, pulse, environment and motion artifact, causing noise many times greater than the neuronal signals of interest. Therefore, artifact removal methods are essential to ensure that artifacts are accurately removed, and EEG of interest is retained. This paper presents a systematic review of methods for artifact reduction in simultaneous EEG-fMRI from literature published since 1998, and an additional systematic review of EEG-fMRI studies published since 2016. The aim of the first review is to distill the literature into clear guidelines for use of simultaneous EEG-fMRI artifact reduction methods, and the aim of the second review is to determine the prevalence of artifact reduction method use in contemporary studies. We find that there are many published artifact reduction techniques available, including hardware, model based, and data-driven methods, but there are few studies published that adequately compare these methods. In contrast, recent EEG-fMRI studies show overwhelming use of just one or two artifact reduction methods based on literature published 15–20 years ago, with newer methods rarely gaining use outside the group that developed them. Surprisingly, almost 15% of EEG-fMRI studies published since 2016 fail to adequately describe the methods of artifact reduction utilized. We recommend minimum standards for reporting artifact reduction techniques in simultaneous EEG-fMRI studies and suggest that more needs to be done to make new artifact reduction techniques more accessible for the researchers and clinicians using simultaneous EEG-fMRI.

Adaptive optimal basis set for BCG artifact removal in simultaneous EEG-fMRI

Electroencephalography (EEG) signals recorded during simultaneous functional magnetic resonance imaging (fMRI) are contaminated by strong artifacts. Among these, the ballistocardiographic (BCG) artifact is the most challenging, due to its complex spatio-temporal dynamics associated with ongoing cardiac activity. The presence of BCG residuals in EEG data may hide true, or generate spurious correlations between EEG and fMRI time-courses. Here, we propose an adaptive Optimal Basis Set (aOBS) method for BCG artifact removal. Our method is adaptive, as it can estimate the delay between cardiac activity and BCG occurrence on a beat-to-beat basis. The effective creation of an optimal basis set by principal component analysis (PCA) is therefore ensured by a more accurate alignment of BCG occurrences. Furthermore, aOBS can automatically estimate which components produced by PCA are likely to be BCG artifact-related and therefore need to be removed. The aOBS performance was evaluated on high-density EEG data acquired with simultaneous fMRI in healthy subjects during visual stimulation. As aOBS enables effective reduction of BCG residuals while preserving brain signals, we suggest it may find wide application in simultaneous EEG-fMRI studies.

Heart-Brain Interactions in the MR Environment: Characterization of the Ballistocardiogram in EEG Signals Collected During Simultaneous fMRI

The ballistocardiographic (BCG) artifact is linked to cardiac activity and occurs in electroencephalographic (EEG) recordings acquired inside the magnetic resonance (MR) environment. Its variability in terms of amplitude, waveform shape and spatial distribution over subject's scalp makes its attenuation a challenging task. In this study, we aimed to provide a detailed characterization of the BCG properties, including its temporal dependency on cardiac events and its spatio-temporal dynamics. To this end, we used high-density EEG data acquired during simultaneous functional MR imaging in six healthy volunteers. First, we investigated the relationship between cardiac activity and BCG occurrences in the EEG recordings. We observed large variability in the delay between ECG and subsequent BCG events (ECG-BCG delay) across subjects and non-negligible epoch-by-epoch variations at the single subject level. The inspection of spatial-temporal variations revealed a prominent non-stationarity of the BCG signal. We identified five main BCG waves, which were common across subjects. Principal component analysis revealed two spatially distinct patterns to explain most of the variance (85% in total). These components are possibly related to head rotation and pulse-driven scalp expansion, respectively. Our results may inspire the development of novel, more effective methods for the removal of the BCG, capable of isolating and attenuating artifact occurrences while preserving true neuronal activity.

Gradient Artefact Correction and Evaluation of the EEG Recorded Simultaneously with fMRI Data Using Optimised Moving-Average

Over the past years, coregistered EEG-fMRI has emerged as a powerful tool for neurocognitive research and correlated studies, mainly because of the possibility of integrating the high temporal resolution of the EEG with the high spatial resolution of fMRI. However, additional work remains to be done in order to improve the quality of the EEG signal recorded simultaneously with fMRI data, in particular regarding the occurrence of the gradient artefact. We devised and presented in this paper a novel approach for gradient artefact correction based upon optimised moving-average filtering (OMA). OMA makes use of the iterative application of a moving-average filter, which allows estimation and cancellation of the gradient artefact by integration. Additionally, OMA is capable of performing the attenuation of the periodic artefact activity without accurate information about MRI triggers. By using our proposed approach, it is possible to achieve a better balance than the slice-average subtraction as performed by the established AAS method, regarding EEG signal preservation together with effective suppression of the gradient artefact. Since the stochastic nature of the EEG signal complicates the assessment of EEG preservation after application of the gradient artefact correction, we also propose a simple and effective method to account for it.