Ballistocardiogram artifact removal with a reference layer and standard EEG cap

Distributed fMRI dynamics predict distinct EEG rhythms across sleep and wakefulness

The brain exhibits rich oscillatory dynamics that vary across tasks and states, such as the EEG oscillations that define sleep. These oscillations play critical roles in cognition and arousal, but the brainwide mechanisms underlying them are not yet described. Using simultaneous EEG and fast fMRI in subjects drifting between sleep and wakefulness, we developed a machine learning approach to investigate which brainwide fMRI dynamics predict alpha (8-12 Hz) and delta (1-4 Hz) rhythms. We predicted moment-by-moment EEG power from fMRI activity in held-out subjects, and found that information about alpha power was represented by a remarkably small set of regions, segregated in two distinct networks linked to arousal and visual systems. Conversely, delta rhythms were diffusely represented on a large spatial scale across the cortex. These results identify distributed networks that predict delta and alpha rhythms, and establish a computational framework for investigating fMRI brainwide dynamics underlying EEG oscillations.

NeuXus open-source tool for real-time artifact reduction in simultaneous EEG-fMRI

The simultaneous acquisition of electroencephalography and functional magnetic resonance imaging (EEG-fMRI) allows the complementary study of the brain's electrophysiology and hemodynamics with high temporal and spatial resolution. One application with great potential is neurofeedback training of targeted brain activity, based on the real-time analysis of the EEG and/or fMRI signals. This depends on the ability to reduce in real time the severe artifacts affecting the EEG signal acquired with fMRI, mainly the gradient and pulse artifacts. A few methods have been proposed for this purpose, but they are either slow, hardware-dependent, publicly unavailable, or proprietary software. Here, we present a fully open-source and publicly available tool for real-time EEG artifact reduction in simultaneous EEG-fMRI recordings that is fast and applicable to any hardware. Our tool is integrated in the Python toolbox NeuXus for real-time EEG processing and adapts to a real-time scenario well-established artifact average subtraction methods combined with a long short-term memory network for R peak detection. We benchmarked NeuXus on three different datasets, in terms of artifact power reduction and background signal preservation in resting state, alpha-band power reactivity to eyes closure, and event-related desynchronization during motor imagery. We showed that NeuXus performed at least as well as the only available real-time tool for conventional hardware setups (BrainVision's RecView) and a well-established offline tool (EEGLAB's FMRIB plugin). We also demonstrated NeuXus' real-time ability by reporting execution times under 250 ms. In conclusion, we present and validate the first fully open-source and hardware-independent solution for real-time artifact reduction in simultaneous EEG-fMRI studies.

An open-access dataset of naturalistic viewing using simultaneous EEG-fMRI

In this work, we present a dataset that combines functional magnetic imaging (fMRI) and electroencephalography (EEG) to use as a resource for understanding human brain function in these two imaging modalities. The dataset can also be used for optimizing preprocessing methods for simultaneously collected imaging data. The dataset includes simultaneously collected recordings from 22 individuals (ages: 23-51) across various visual and naturalistic stimuli. In addition, physiological, eye tracking, electrocardiography, and cognitive and behavioral data were collected along with this neuroimaging data. Visual tasks include a flickering checkerboard collected outside and inside the MRI scanner (EEG-only) and simultaneous EEG-fMRI recordings. Simultaneous recordings include rest, the visual paradigm Inscapes, and several short video movies representing naturalistic stimuli. Raw and preprocessed data are openly available to download. We present this dataset as part of an effort to provide open-access data to increase the opportunity for discoveries and understanding of the human brain and evaluate the correlation between electrical brain activity and blood oxygen level-dependent (BOLD) signals.

EEG-LLAMAS: A low-latency neurofeedback platform for artifact reduction in EEG-fMRI

Simultaneous EEG-fMRI is a powerful multimodal technique for imaging the brain, but its use in neurofeedback experiments has been limited by EEG noise caused by the MRI environment. Neurofeedback studies typically require analysis of EEG in real time, but EEG acquired inside the scanner is heavily contaminated with ballistocardiogram (BCG) artifact, a high-amplitude artifact locked to the cardiac cycle. Although techniques for removing BCG artifacts do exist, they are either not suited to real-time, low-latency applications, such as neurofeedback, or have limited efficacy. We propose and validate a new open-source artifact removal software called EEG-LLAMAS (Low Latency Artifact Mitigation Acquisition Software), which adapts and advances existing artifact removal techniques for low-latency experiments. We first used simulations to validate LLAMAS in data with known ground truth. We found that LLAMAS performed better than the best publicly-available real-time BCG removal technique, optimal basis sets (OBS), in terms of its ability to recover EEG waveforms, power spectra, and slow wave phase. To determine whether LLAMAS would be effective in practice, we then used it to conduct real-time EEG-fMRI recordings in healthy adults, using a steady state visual evoked potential (SSVEP) task. We found that LLAMAS was able to recover the SSVEP in real time, and recovered the power spectra collected outside the scanner better than OBS. We also measured the latency of LLAMAS during live recordings, and found that it introduced a lag of less than 50 ms on average. The low latency of LLAMAS, coupled with its improved artifact reduction, can thus be effectively used for EEG-fMRI neurofeedback. A limitation of the method is its use of a reference layer, a piece of EEG equipment which is not commercially available, but can be assembled in-house. This platform enables closed-loop experiments which previously would have been prohibitively difficult, such as those that target short-duration EEG events, and is shared openly with the neuroscience community.

The MotoNet: A 3 Tesla MRI-Conditional EEG Net with Embedded Motion Sensors

We introduce a new electroencephalogram (EEG) net, which will allow clinicians to monitor EEG while tracking head motion. Motion during MRI limits patient scans, especially of children with epilepsy. EEG is also severely affected by motion-induced noise, predominantly ballistocardiogram (BCG) noise due to the heartbeat.

METHODS

The MotoNet was built using polymer thick film (PTF) EEG leads and motion sensors on opposite sides in the same flex circuit. EEG/motion measurements were made with a standard commercial EEG acquisition system in a 3 Tesla (T) MRI. A Kalman filtering-based BCG correction tool was used to clean the EEG in healthy volunteers.

RESULTS

MRI safety studies in 3 T confirmed the maximum heating below 1 °C. Using an MRI sequence with spatial localization gradients only, the position of the head was linearly correlated with the average motion sensor output. Kalman filtering was shown to reduce the BCG noise and recover artifact-clean EEG.

CONCLUSIONS

The MotoNet is an innovative EEG net design that co-locates 32 EEG electrodes with 32 motion sensors to improve both EEG and MRI signal quality. In combination with custom gradients, the position of the net can, in principle, be determined. In addition, the motion sensors can help reduce BCG noise.

Neural activity induced by sensory stimulation can drive large-scale cerebrospinal fluid flow during wakefulness in humans

Cerebrospinal fluid (CSF) flow maintains healthy brain homeostasis, facilitating solute transport and the exchange of brain waste products. CSF flow is thus important for brain health, but the mechanisms that control its large-scale movement through the ventricles are not well understood. While it is well established that CSF flow is modulated by respiratory and cardiovascular dynamics, recent work has also demonstrated that neural activity is coupled to large waves of CSF flow in the ventricles during sleep. To test whether the temporal coupling between neural activity and CSF flow is in part due to a causal relationship, we investigated whether CSF flow could be induced by driving neural activity with intense visual stimulation. We manipulated neural activity with a flickering checkerboard visual stimulus and found that we could drive macroscopic CSF flow in the human brain. The timing and amplitude of CSF flow was matched to the visually evoked hemodynamic responses, suggesting neural activity can modulate CSF flow via neurovascular coupling. These results demonstrate that neural activity can contribute to driving CSF flow in the human brain and that the temporal dynamics of neurovascular coupling can explain this effect.

A temporal sequence of thalamic activity unfolds at transitions in behavioral arousal state

Awakening from sleep reflects a profound transformation in neural activity and behavior. The thalamus is a key controller of arousal state, but whether its diverse nuclei exhibit coordinated or distinct activity at transitions in behavioral arousal state is unknown. Using fast fMRI at ultra-high field (7 Tesla), we measured sub-second activity across thalamocortical networks and within nine thalamic nuclei to delineate these dynamics during spontaneous transitions in behavioral arousal state. We discovered a stereotyped sequence of activity across thalamic nuclei and cingulate cortex that preceded behavioral arousal after a period of inactivity, followed by widespread deactivation. These thalamic dynamics were linked to whether participants subsequently fell back into unresponsiveness, with unified thalamic activation reflecting maintenance of behavior. These results provide an outline of the complex interactions across thalamocortical circuits that orchestrate behavioral arousal state transitions, and additionally, demonstrate that fast fMRI can resolve sub-second subcortical dynamics in the human brain.

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.

Investigating mechanisms of fast BOLD responses: The effects of stimulus intensity and of spatial heterogeneity of hemodynamics

Recent studies have demonstrated that fast fMRI can track neural activity well above the temporal limit predicted by the canonical hemodynamic response model. While these findings are promising, the biophysical mechanisms underlying these fast fMRI phenomena remain underexplored. In this study, we discuss two aspects of the hemodynamic response, complementary to several existing hypotheses, that can accommodate faster fMRI dynamics beyond those predicted by the canonical model. First, we demonstrate, using both visual and somatosensory paradigms, that the timing and shape of hemodynamic response functions (HRFs) vary across graded levels of stimulus intensity-with lower-intensity stimulation eliciting faster and narrower HRFs. Second, we show that as the spatial resolution of fMRI increases, voxel-wise HRFs begin to deviate from the canonical model, with a considerable portion of voxels exhibiting faster temporal dynamics than predicted by the canonical HRF. Collectively, both stimulus/task intensity and image resolution can affect the sensitivity of fMRI to fast brain activity, which may partly explain recent observations of fast fMRI signals. It is further noteworthy that, while the present investigations focus on fast neural responses, our findings suggest that a revised hemodynamic model may benefit the many fMRI studies using paradigms with wide ranges of contrast levels (e.g., resting or naturalistic conditions) or with modern, high-resolution MR acquisitions.

Unified Retrospective EEG Motion Educated Artefact Suppression for EEG-fMRI to Suppress Magnetic Field Gradient Artefacts During Motion

The data quality of simultaneously acquired electroencephalography and functional magnetic resonance imaging (EEG-fMRI) can be strongly affected by motion. Recent work has shown that the quality of fMRI data can be improved by using a Moiré-Phase-Tracker (MPT)-camera system for prospective motion correction. The use of the head position acquired by the MPT-camera-system has also been shown to correct motion-induced voltages, ballistocardiogram (BCG) and gradient artefact residuals separately. In this work we show the concept of an integrated framework based on the general linear model to provide a unified motion informed model of in-MRI artefacts. This model (retrospective EEG motion educated gradient artefact suppression, REEG-MEGAS) is capable of correcting voltage-induced, BCG and gradient artefact residuals of EEG data acquired simultaneously with prospective motion corrected fMRI. In our results, we have verified that applying REEG-MEGAS correction to EEG data acquired during subject motion improves the data quality in terms of motion induced voltages and also GA residuals in comparison to standard Artefact Averaging Subtraction and Retrospective EEG Motion Artefact Suppression. Besides that, we provide preliminary evidence that although adding more regressors to a model may slightly affect the power of physiological signals such as the alpha-rhythm, its application may increase the overall quality of a dataset, particularly when strongly affected by motion. This was verified by analysing the EEG traces, power spectra density and the topographic distribution from two healthy subjects. We also have verified that the correction by REEG-MEGAS improves higher frequency artefact correction by decreasing the power of Gradient Artefact harmonics. Our method showed promising results for decreasing the power of artefacts for frequencies up to 250 Hz. Additionally, REEG-MEGAS is a hybrid framework that can be implemented for real time prospective motion correction of EEG and fMRI data. Among other EEG-fMRI applications, the approach described here may benefit applications such as EEG-fMRI neurofeedback and brain computer interface, which strongly rely on the prospective acquisition and application of motion artefact removal.

Automated pipeline for EEG artifact reduction (APPEAR) recorded during fMRI

Objective.Simultaneous electroencephalography-functional magnetic resonance imaging (EEG-fMRI) recordings offer a high spatiotemporal resolution approach to study human brain and understand the underlying mechanisms mediating cognitive and behavioral processes. However, the high susceptibility of EEG to MRI-induced artifacts hinders a broad adaptation of this approach. More specifically, EEG data collected during fMRI acquisition are contaminated with MRI gradients and ballistocardiogram artifacts, in addition to artifacts of physiological origin. There have been several attempts for reducing these artifacts with manual and time-consuming pre-processing, which may result in biasing EEG data due to variations in selecting steps order, parameters, and classification of artifactual independent components. Thus, there is a strong urge to develop a fully automatic and comprehensive pipeline for reducing all major EEG artifacts. In this work, we introduced an open-access toolbox with a fully automatic pipeline for reducing artifacts from EEG data collected simultaneously with fMRI (refer to APPEAR).Approach.The pipeline integrates average template subtraction and independent component analysis to suppress both MRI-related and physiological artifacts. To validate our results, we tested APPEAR on EEG data recorded from healthy control subjects during resting-state (n= 48) and task-based (i.e. event-related-potentials (ERPs);n= 8) paradigms. The chosen gold standard is an expert manual review of the EEG database.Main results.We compared manually and automated corrected EEG data during resting-state using frequency analysis and continuous wavelet transformation and found no significant differences between the two corrections. A comparison between ERP data recorded during a so-called stop-signal task (e.g. amplitude measures and signal-to-noise ratio) also showed no differences between the manually and fully automatic fMRI-EEG-corrected data.Significance.APPEAR offers the first comprehensive open-source toolbox that can speed up advancement of EEG analysis and enhance replication by avoiding experimenters' preferences while allowing for processing large EEG-fMRI cohorts composed of hundreds of subjects with manageable researcher time and effort.

Safety and data quality of EEG recorded simultaneously with multi-band fMRI

PURPOSE

Simultaneously recorded electroencephalography and functional magnetic resonance imaging (EEG-fMRI) is highly informative yet technically challenging. Until recently, there has been little information about EEG data quality and safety when used with newer multi-band (MB) fMRI sequences. Here, we measure the relative heating of a MB protocol compared with a standard single-band (SB) protocol considered to be safe. We also evaluated EEG quality recorded concurrently with the MB protocol on humans.

MATERIALS AND METHODS

We compared radiofrequency (RF)-related heating at multiple electrodes and magnetic field magnitude, B1+RMS, of a MB fMRI sequence with whole-brain coverage (TR = 440 ms, MB factor = 4) against a previously recommended, safe SB sequence using a phantom outfitted with a 64-channel EEG cap. Next, 9 human subjects underwent eyes-closed resting state EEG-fMRI using the MB sequence. Additionally, in three of the subjects resting state EEG was recorded also during the SB sequence and in an fMRI-free condition to directly compare EEG data quality across scanning conditions. EEG data quality was assessed by the ability to remove gradient and cardioballistic artifacts along with a clean spectrogram.

RESULTS

The heating induced by the MB sequence was lower than that of the SB sequence by a factor of 0.73 ± 0.38. This is consistent with an expected heating ratio of 0.64, calculated from the square of the ratio of B1+RMS values of the sequences. In the resting state EEG data, gradient and cardioballistic artifacts were successfully removed using traditional template subtraction. All subjects showed an individual alpha peak in the spectrogram with a posterior topography characteristic of eyes-closed EEG. The success of artifact rejection for the MB sequence was comparable to that in traditional SB sequences.

CONCLUSIONS

Our study shows that B1+RMS is a useful indication of the relative heating of fMRI protocols. This observation indicates that simultaneous EEG-fMRI recordings using this MB sequence can be safe in terms of RF-related heating, and that EEG data recorded using this sequence is of acceptable quality after traditional artifact removal techniques.

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.

Ballistocardiogram Artifact Reduction in Simultaneous EEG-fMRI Using Deep Learning

OBJECTIVE

The concurrent recording of electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) is a technique that has received much attention due to its potential for combined high temporal and spatial resolution. However, the ballistocardiogram (BCG), a large-amplitude artifact caused by cardiac induced movement contaminates the EEG during EEG-fMRI recordings. Removal of BCG in software has generally made use of linear decompositions of the corrupted EEG. This is not ideal as the BCG signal propagates in a manner which is non-linearly dependent on the electrocardiogram (ECG). In this paper, we present a novel method for BCG artifact suppression using recurrent neural networks (RNNs).

METHODS

EEG signals were recovered by training RNNs on the nonlinear mappings between ECG and the BCG corrupted EEG. We evaluated our model's performance against the commonly used Optimal Basis Set (OBS) method at the level of individual subjects, and investigated generalization across subjects.

RESULTS

We show that our algorithm can generate larger average power reduction of the BCG at critical frequencies, while simultaneously improving task relevant EEG based classification.

CONCLUSION

The presented deep learning architecture can be used to reduce BCG related artifacts in EEG-fMRI recordings.

SIGNIFICANCE

We present a deep learning approach that can be used to suppress the BCG artifact in EEG-fMRI without the use of additional hardware. This method may have scope to be combined with current hardware methods, operate in real-time and be used for direct modeling of the BCG.

Coupled electrophysiological, hemodynamic, and cerebrospinal fluid oscillations in human sleep

Sleep is essential for both cognition and maintenance of healthy brain function. Slow waves in neural activity contribute to memory consolidation, whereas cerebrospinal fluid (CSF) clears metabolic waste products from the brain. Whether these two processes are related is not known. We used accelerated neuroimaging to measure physiological and neural dynamics in the human brain. We discovered a coherent pattern of oscillating electrophysiological, hemodynamic, and CSF dynamics that appears during non-rapid eye movement sleep. Neural slow waves are followed by hemodynamic oscillations, which in turn are coupled to CSF flow. These results demonstrate that the sleeping brain exhibits waves of CSF flow on a macroscopic scale, and these CSF dynamics are interlinked with neural and hemodynamic rhythms.

Direct modelling of gradient artifacts for EEG-fMRI denoising and motion tracking

OBJECTIVE

Simultaneous electroencephalography and functional magnetic resonance imaging recording (EEG-fMRI) has been widely used in neuroscientific and clinical research. The artifacts in the recorded EEG resulting from rapidly switching magnetic field gradients are usually corrected by average-artifact subtraction (AAS) due to their repetitive nature. But the performance of AAS is often disrupted by altered artifact waveforms across epochs, notably due to head motion.

APPROACH

Here, a method is proposed to make use of the known MR sequence gradient waveforms for a direct modelling of gradient artifacts. After accounting for filtering effects on the gradient artifacts, a continuous modulation of the gradient waveforms superimposed on the EEG signal is obtained.

MAIN RESULTS

Although a moving AAS template can adjust to slow drifts in gradient artifact variation, it fails to adapt to abrupt motion, resulting in residual noise. We demonstrate how this modelling approach can reduce motion-affected gradient artifacts without distorting the underlying neuronal signals. Moreover, the method provides useful head motion information highly correlated with motion tracked by an optical camera.

SIGNIFICANCE

Our work provides a novel way to improve gradient artifact removal in EEG-fMRI, and shows a potential to detect head motion without requiring additional hardware.

Simultaneous EEG-fMRI: Evaluating the Effect of the EEG Cap-Cabling Configuration on the Gradient Artifact

Electroencephalography (EEG) data recorded during simultaneous EEG-fMRI experiments are contaminated by large gradient artifacts (GA). The amplitude of the GA depends on the area of the wire loops formed by the EEG leads, as well as on the rate of switching of the magnetic field gradients, which are essential for MR imaging. Average artifact subtraction (AAS), the most commonly used method for GA correction, relies on the EEG amplifier having a large enough dynamic range to characterize the artifact voltages. Low-pass filtering (250 Hz cut-off) is generally used to attenuate the high-frequency voltage fluctuations of the GA, but even with this precaution channel saturation can occur, particularly during acquisition of high spatial resolution MRI data. Previous work has shown that the ribbon cable, used to connect the EEG cap and amplifier, makes a significant contribution to the GA, since the cable geometry produces large effective wire-loop areas. However, by appropriately connecting the wires of the ribbon cable to the EEG cap it should be possible to minimize the overall range and root mean square (RMS) amplitude of the GA by producing partial cancelation of the cap and cable contributions. Here by modifying the connections of the EEG cap to a 1 m ribbon cable we were able to reduce the range of the GA for a high-resolution coronal echo planar Imaging (EPI) acquisition by a factor of ∼ 1.6 and by a factor of ∼ 1.15 for a standard axial EPI acquisition. These changes could potentially be translated into a reduction in the required dynamic range, an increase in the EEG bandwidth or an increase in the achievable image resolution without saturation, all of which could be beneficially exploited in EEG-fMRI studies. The re-wiring could also prevent the system from saturating when small subject movements occur using the standard recording bandwidth.

Conductive gel bridge sensor for motion tracking in simultaneous EEG-fMRI recordings

EEG-fMRI allows the localization of the hemodynamic correlates of neural activity and has been shown to be useful as a diagnostic tool in pre-surgical evaluation of refractory epilepsy. However, EEG recordings may be highly contaminated by artifacts induced by movements inside the magnetic field thus rendering the scan difficult for interpretation. Existing methods for motion correction require additional equipment or hardware modification. We introduce a simple method for motion artifact detection, the conductive gel bridge sensor (CGBS), easily applicable using the standard setup. We report examples of CGBS use in two patients with epilepsy and demonstrate the method's ability to successfully differentiate between epochs of brain activity and those of movement.

Adaptive and Wireless Recordings of Electrophysiological Signals During Concurrent Magnetic Resonance Imaging

Strong electromagnetic fields that occur during functional magnetic resonance imaging (fMRI) presents a challenging environment for concurrent electrophysiological recordings. Here, we present a miniaturized, wireless platform-"MR-Link" (Multimodal Recording Link) that provides a hardware solution for simultaneous electrophysiological and fMRI signal acquisition. The device detects the changes in the electromagnetic field during fMRI to synchronize amplification and sampling of electrophysiological signals with minimal artifacts. It wirelessly transmits the recorded data at a frequency detectable by the MR-receiver coil. The transmitted data is readily separable from MRI in the frequency domain. To demonstrate its efficacy, we used this device to record electrocardiograms and somatosensory evoked potential during concurrent fMRI scans. The device minimized the fMRI-induced artifacts in electrophysiological data and wirelessly transmitted the data back to the receiver coil without compromising the fMRI signal quality. The device is compact (22 mm dia., 2 gms) and can be placed within the MRI bore to precisely synchronize with fMRI. Therefore, MR-Link offers an inexpensive system by eliminating the need for amplifiers with a high dynamic range, high-speed sampling, additional storage, or synchronization hardware for electrophysiological signal acquisition. It is expected to enable a broader range of applications of simultaneous fMRI and electrophysiology in animals and humans.

Exploring the relative efficacy of motion artefact correction techniques for EEG data acquired during simultaneous fMRI

Simultaneous EEG‐fMRI allows multiparametric characterisation of brain function, in principle enabling a more complete understanding of brain responses; unfortunately the hostile MRI environment severely reduces EEG data quality. Simply eliminating data segments containing gross motion artefacts [MAs] (generated by movement of the EEG system and head in the MRI scanner's static magnetic field) was previously believed sufficient. However recently the importance of removal of all MAs has been highlighted and new methods developed. A systematic comparison of the ability to remove MAs and retain underlying neuronal activity using different methods of MA detection and post‐processing algorithms is needed to guide the neuroscience community. Using a head phantom, we recorded MAs while simultaneously monitoring the motion using three different approaches: Reference Layer Artefact Subtraction (RLAS), Moiré Phase Tracker (MPT) markers and Wire Loop Motion Sensors (WLMS). These EEG recordings were combined with EEG responses to simple visual tasks acquired on a subject outside the MRI environment. MAs were then corrected using the motion information collected with each of the methods combined with different analysis pipelines. All tested methods retained the neuronal signal. However, often the MA was not removed sufficiently to allow accurate detection of the underlying neuronal signal. We show that the MA is best corrected using the RLAS combined with post‐processing using a multichannel, recursive least squares (M‐RLS) algorithm. This method needs to be developed further to enable practical utility; thus, WLMS combined with M‐RLS currently provides the best compromise between EEG data quality and practicalities of motion detection.

Exploring the origins of EEG motion artefacts during simultaneous fMRI acquisition: Implications for motion artefact correction

Motion artefacts (MAs) are induced within EEG data collected simultaneously with fMRI when the subject's head rotates relative to the magnetic field. The effects of these artefacts have generally been ameliorated by removing periods of data during which large artefact voltages appear in the EEG traces. However, even when combined with other standard post-processing methods, this strategy does not remove smaller MAs which can dominate the neuronal signals of interest. A number of methods are therefore being developed to characterise the MA by measuring reference signals and then using these in artefact correction. These methods generally assume that the head and EEG cap, plus any attached sensors, form a rigid body which can be characterised by a standard set of six motion parameters. Here we investigate the motion of the head/EEG cap system to provide a better understanding of MAs. We focus on the reference layer artefact subtraction (RLAS) approach, as this allows measurement of a separate reference signal for each electrode that is being used to measure brain activity.

Through a series of experiments on phantoms and subjects, we find that movement of the EEG cap relative to the phantom and skin on the forehead is relatively small and that this non-rigid body movement does not appear to cause considerable discrepancy in artefacts between the scalp and reference signals. However, differences in the amplitude of these signals is observed which may be due to differences in geometry of the system from which the reference signals are measured compared with the brain signals. In addition, we find that there is non-rigid body movement of the skull and skin which produces an additional MA component for a head shake, which is not present for a head nod. This results in a large discrepancy in the amplitude and temporal profile of the MA measured on the scalp and reference layer, reducing the efficacy of MA correction based on the reference signals.

Together our data suggest that the efficacy of the correction of MA using any reference-based system is likely to differ for different types of head movement with head shake being the hardest to correct. This provides new information to inform the development of hardware and post-processing methods for removing MAs from EEG data acquired simultaneously with fMRI data.

Clustering-Constrained ICA for Ballistocardiogram Artifacts Removal in Simultaneous EEG-fMRI

Combination of electroencephalogram (EEG) recording and functional magnetic resonance imaging (fMRI) plays a potential role in neuroimaging due to its high spatial and temporal resolution. However, EEG is easily influenced by ballistocardiogram (BCG) artifacts and may cause false identification of the related EEG features, such as epileptic spikes. There are many related methods to remove them, however, they do not consider the time-varying features of BCG artifacts. In this paper, a novel method using clustering algorithm to catch the BCG artifacts' features and together with the constrained ICA (ccICA) is proposed to remove the BCG artifacts. We first applied this method to the simulated data, which was constructed by adding the BCG artifacts to the EEG signal obtained from the conventional environment. Then, our method was tested to demonstrate the effectiveness during EEG and fMRI experiments on 10 healthy subjects. In simulated data analysis, the value of error in signal amplitude (Er) computed by ccICA method was lower than those from other methods including AAS, OBS, and cICA (p < 0.005). In vivo data analysis, the Improvement of Normalized Power Spectrum (INPS) calculated by ccICA method in all electrodes was much higher than AAS, OBS, and cICA methods (p < 0.005). We also used other evaluation index (e.g., power analysis) to compare our method with other traditional methods. In conclusion, our novel method successfully and effectively removed BCG artifacts in both simulated and vivo EEG data tests, showing the potentials of removing artifacts in EEG-fMRI applications.

Removal of BCG artefact from concurrent fMRI-EEG recordings based on EMD and PCA

BACKGROUND

Simultaneous electroencephalography (EEG) and functional magnetic resonance image (fMRI) acquisitions provide better insight into brain dynamics. Some artefacts due to simultaneous acquisition pose a threat to the quality of the data. One such problematic artefact is the ballistocardiogram (BCG) artefact.

METHODS

We developed a hybrid algorithm that combines features of empirical mode decomposition (EMD) with principal component analysis (PCA) to reduce the BCG artefact. The algorithm does not require extra electrocardiogram (ECG) or electrooculogram (EOG) recordings to extract the BCG artefact.

RESULTS

The method was tested with both simulated and real EEG data of 11 participants. From the simulated data, the similarity index between the extracted BCG and the simulated BCG showed the effectiveness of the proposed method in BCG removal. On the other hand, real data were recorded with two conditions, i.e. resting state (eyes closed dataset) and task influenced (event-related potentials (ERPs) dataset). Using qualitative (visual inspection) and quantitative (similarity index, improved normalized power spectrum (INPS) ratio, power spectrum, sample entropy (SE)) evaluation parameters, the assessment results showed that the proposed method can efficiently reduce the BCG artefact while preserving the neuronal signals.

COMPARISON WITH EXISTING METHODS

Compared with conventional methods, namely, average artefact subtraction (AAS), optimal basis set (OBS) and combined independent component analysis and principal component analysis (ICA-PCA), the statistical analyses of the results showed that the proposed method has better performance, and the differences were significant for all quantitative parameters except for the power and sample entropy.

CONCLUSIONS

The proposed method does not require any reference signal, prior information or assumption to extract the BCG artefact. It will be very useful in circumstances where the reference signal is not available.

Reference layer adaptive filtering (RLAF) for EEG artifact reduction in simultaneous EEG-fMRI

OBJECTIVE

Simultaneous electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) combines advantages of both methods, namely high temporal resolution of EEG and high spatial resolution of fMRI. However, EEG quality is limited due to severe artifacts caused by fMRI scanners.

APPROACH

To improve EEG data quality substantially, we introduce methods that use a reusable reference layer EEG cap prototype in combination with adaptive filtering. The first method, reference layer adaptive filtering (RLAF), uses adaptive filtering with reference layer artifact data to optimize artifact subtraction from EEG. In the second method, multi band reference layer adaptive filtering (MBRLAF), adaptive filtering is performed on bandwidth limited sub-bands of the EEG and the reference channels.

MAIN RESULTS

The results suggests that RLAF outperforms the baseline method, average artifact subtraction, in all settings and also its direct predecessor, reference layer artifact subtraction (RLAS), in lower (<35 Hz) frequency ranges. MBRLAF is computationally more demanding than RLAF, but highly effective in all EEG frequency ranges. Effectivity is determined by visual inspection, as well as root-mean-square voltage reduction and power reduction of EEG provided that physiological EEG components such as occipital EEG alpha power and visual evoked potentials (VEP) are preserved. We demonstrate that both, RLAF and MBRLAF, improve VEP quality. For that, we calculate the mean-squared-distance of single trial VEP to the mean VEP and estimate single trial VEP classification accuracies. We found that the average mean-squared-distance is lowest and the average classification accuracy is highest after MBLAF. RLAF was second best.

SIGNIFICANCE

In conclusion, the results suggests that RLAF and MBRLAF are potentially very effective in improving EEG quality of simultaneous EEG-fMRI. Highlights We present a new and reusable reference layer cap prototype for simultaneous EEG-fMRI We introduce new algorithms for reducing EEG artifacts due to simultaneous fMRI The algorithms combine a reference layer and adaptive filtering Several evaluation criteria suggest superior effectivity in terms of artifact reduction We demonstrate that physiological EEG components are preserved.

Effectiveness of Reference Signal-Based Methods for Removal of EEG Artifacts Due to Subtle Movements During fMRI Scanning

OBJECTIVE

Subtle motion of an epileptic patient examined with co-registered EEG and functional MRI (EEG-fMRI) may often lead to spurious fMRI activation patterns when true epileptic spikes are contaminated with motion artefacts. In recent years, methods relying on reference signals for correcting these subtle movements in the EEG have emerged. In this study, the performance of two reference-based devices are compared to the template-based method with regard to their ability to remove movement-related artifacts in EEG measured during scanning.

METHODS

Measurements were performed with a novel double layer cap consisting of 29 EEG and 29 reference electrodes, and with a current loop cap consisting of 60 electrodes and three current loop wires attached to the cap. EEG was acquired inside the scanner during resting state, as well as when the subject was performing a cued movement task. For the double layer cap recordings, newly developed artifact removal algorithms are introduced and both reference signal-based methods are compared to a template-based correction method.

RESULTS

The BCG artifacts occurring at resting state could be removed successfully by both the reference signal-based methods as well as by the template-based method. However, the reference signal-based methods were also capable of removing EEG artifacts induced by subtle movements, whereas the template-based method failed to remove these artifacts.

CONCLUSION

Reference signal-based methods enable to correct for artifacts due to subtle movements, which are not removed by commonly used template-based removal algorithms.

SIGNIFICANCE

Sensitivity of EEG-fMRI analysis in patients with focal epilepsy is improved by avoiding erroneous detections of subtle movements as epileptic spikes in the EEG.

Dual array EEG-fMRI: An approach for motion artifact suppression in EEG recorded simultaneously with fMRI

OBJECTIVE

Although simultaneous recording of EEG and MRI has gained increasing popularity in recent years, the extent of its clinical use remains limited by various technical challenges. Motion interference is one of the major challenges in EEG-fMRI. Here we present an approach which reduces its impact with the aid of an MR compatible dual-array EEG (daEEG) in which the EEG itself is used both as a brain signal recorder and a motion sensor.

METHODS

We implemented two arrays of EEG electrodes organized into two sets of nearly orthogonally intersecting wire bundles. The EEG was recorded using referential amplifiers inside a 3T MR-scanner. Virtual bipolar measurements were taken both along bundles (creating a small wire loop and therefore minimizing artifact) and across bundles (creating a large wire loop and therefore maximizing artifact). Independent component analysis (ICA) was applied. The resulting ICA components were classified into brain signal and noise using three criteria: 1) degree of two-dimensional spatial correlation between ICA coefficients along bundles and across bundles; 2) amplitude along bundles vs. across bundles; 3) correlation with ECG. The components which passed the criteria set were transformed back to the channel space. Motion artifact suppression and the ability to detect interictal epileptic spikes following daEEG and Optimal Basis Set (OBS) procedures were compared in 10 patients with epilepsy.

RESULTS

The SNR achieved by daEEG was 11.05±3.10 and by OBS was 8.25±1.01 (p<0.00001). In 9 of 10 patients, more spikes were detected after daEEG than after OBS (p<0.05).

SIGNIFICANCE

daEEG improves signal quality in EEG-fMRI recordings, expanding its clinical and research potential.

Polymer thick film technology for improved simultaneous dEEG/MRI recording: Safety and MRI data quality

PURPOSE

To develop a 256-channel dense-array electroencephalography (dEEG) sensor net (the Ink-Net) using high-resistance polymer thick film (PTF) technology to improve safety and data quality during simultaneous dEEG/MRI.

METHODS

Heating safety was assessed with temperature measurements in an anthropomorphic head phantom during a 30-min, induced-heating scan at 7T. MRI quality assessment used B1 field mapping and functional MRI (fMRI) retinotopic scans in three humans at 3T. Performance of the 256-channel PTF Ink-Net was compared with a 256-channel MR-conditional copper-wired electroencephalography (EEG) net and to scans with no sensor net. A visual evoked potential paradigm assessed EEG quality within and outside the 3T scanner.

RESULTS

Phantom temperature measurements revealed nonsignificant heating (ISO 10974) in the presence of either EEG net. In human B1 field and fMRI scans, the Ink-Net showed greatly reduced cross-modal artifact and less signal degradation than the copper-wired net, and comparable quality to MRI without sensor net. Cross-modal ballistocardiogram artifact in the EEG was comparable for both nets.

CONCLUSION

High-resistance PTF technology can be effectively implemented in a 256-channel dEEG sensor net for MR conditional use at 7T and with significantly improved structural and fMRI data quality as assessed at 3T. Magn Reson Med 77:895-903, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Carbon-wire loop based artifact correction outperforms post-processing EEG/fMRI corrections--A validation of a real-time simultaneous EEG/fMRI correction method

Simultaneous EEG-fMRI combines two powerful neuroimaging techniques, but the EEG signal suffers from severe artifacts in the MRI environment that are difficult to remove. These are the MR scanning artifact and the blood-pulsation artifact--strategies to remove them are a topic of ongoing research. Additionally large, unsystematic artifacts are produced across the full frequency spectrum by the magnet's helium pump (and ventilator) systems which are notoriously hard to remove. As a consequence, experimenters routinely deactivate the helium pump during simultaneous EEG-fMRI acquisitions which potentially risks damaging the MRI system and necessitates more frequent and expensive helium refills. We present a novel correction method addressing both helium pump and ballisto-cardiac (BCG) artifacts, consisting of carbon-wire loops (CWL) as additional sensors to accurately track unpredictable artifacts related to subtle movements in the scanner, and an EEGLAB plugin to perform artifact correction. We compare signal-to-noise metrics of EEG data, corrected with CWL and three conventional correction methods, for helium pump off and on measurements. Because the CWL setup records signals in real-time, it fits requirements of applications where immediate correction is necessary, such as neuro-feedback applications or stimulation time-locked to specific sleep oscillations. The comparison metrics in this paper relate to: (1) the EEG signal itself, (2) the "eyes open vs. eyes closed" effect, and (3) an assessment of how the artifact corrections impacts the ability to perform meaningful correlations between EEG alpha power and the BOLD signal. Results show that the CWL correction corrects for He pump artifact and also produces EEG data more comparable to EEG obtained outside the magnet than conventional post-processing methods.

Towards high-quality simultaneous EEG-fMRI at 7 T: Detection and reduction of EEG artifacts due to head motion

The enhanced functional sensitivity offered by ultra-high field imaging may significantly benefit simultaneous EEG-fMRI studies, but the concurrent increases in artifact contamination can strongly compromise EEG data quality. In the present study, we focus on EEG artifacts created by head motion in the static B0 field. A novel approach for motion artifact detection is proposed, based on a simple modification of a commercial EEG cap, in which four electrodes are non-permanently adapted to record only magnetic induction effects. Simultaneous EEG-fMRI data were acquired with this setup, at 7 T, from healthy volunteers undergoing a reversing-checkerboard visual stimulation paradigm. Data analysis assisted by the motion sensors revealed that, after gradient artifact correction, EEG signal variance was largely dominated by pulse artifacts (81-93%), but contributions from spontaneous motion (4-13%) were still comparable to or even larger than those of actual neuronal activity (3-9%). Multiple approaches were tested to determine the most effective procedure for denoising EEG data incorporating motion sensor information. Optimal results were obtained by applying an initial pulse artifact correction step (AAS-based), followed by motion artifact correction (based on the motion sensors) and ICA denoising. On average, motion artifact correction (after AAS) yielded a 61% reduction in signal power and a 62% increase in VEP trial-by-trial consistency. Combined with ICA, these improvements rose to a 74% power reduction and an 86% increase in trial consistency. Overall, the improvements achieved were well appreciable at single-subject and single-trial levels, and set an encouraging quality mark for simultaneous EEG-fMRI at ultra-high field.