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

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.

The Scanner as the Stimulus: Deficient Gamma-BOLD Coupling in Schizophrenia at Rest

Functional magnetic resonance imaging (fMRI) scanners are unavoidably loud and uncomfortable experimental tools that are necessary for schizophrenia (SZ) neuroscience research. The validity of fMRI paradigms might be undermined by well-known sensory processing abnormalities in SZ that could exert distinct effects on neural activity in the presence of scanner background sound. Given the ubiquity of resting-state fMRI (rs-fMRI) paradigms in SZ research, elucidating the relationship between neural, hemodynamic, and sensory processing deficits during scanning is necessary to refine the construct validity of the MR neuroimaging environment. We recorded simultaneous electroencephalography (EEG)-fMRI at rest in people with SZ (n = 57) and healthy control participants without a psychiatric diagnosis (n = 46) and identified gamma EEG activity in the same frequency range as the background sounds emitted from our scanner during a resting-state sequence. In participants with SZ, gamma coupling to the hemodynamic signal was reduced in bilateral auditory regions of the superior temporal gyri. Impaired gamma-hemodynamic coupling was associated with sensory gating deficits and worse symptom severity. Fundamental sensory-neural processing deficits in SZ are present at rest when considering scanner background sound as a "stimulus." This finding may impact the interpretation of rs-fMRI activity in studies of people with SZ. Future neuroimaging research in SZ might consider background sound as a confounding variable, potentially related to fluctuations in neural excitability and arousal.

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.

Cognition is entangled with metabolism: relevance for resting-state EEG-fMRI

The brain is a living organ with distinct metabolic constraints. However, these constraints are typically considered as secondary or supportive of information processing which is primarily performed by neurons. The default operational definition of neural information processing is that (1) it is ultimately encoded as a change in individual neuronal firing rate as this correlates with the presentation of a peripheral stimulus, motor action or cognitive task. Two additional assumptions are associated with this default interpretation: (2) that the incessant background firing activity against which changes in activity are measured plays no role in assigning significance to the extrinsically evoked change in neural firing, and (3) that the metabolic energy that sustains this background activity and which correlates with differences in neuronal firing rate is merely a response to an evoked change in neuronal activity. These assumptions underlie the design, implementation, and interpretation of neuroimaging studies, particularly fMRI, which relies on changes in blood oxygen as an indirect measure of neural activity. In this article we reconsider all three of these assumptions in light of recent evidence. We suggest that by combining EEG with fMRI, new experimental work can reconcile emerging controversies in neurovascular coupling and the significance of ongoing, background activity during resting-state paradigms. A new conceptual framework for neuroimaging paradigms is developed to investigate how ongoing neural activity is "entangled" with metabolism. That is, in addition to being recruited to support locally evoked neuronal activity (the traditional hemodynamic response), changes in metabolic support may be independently "invoked" by non-local brain regions, yielding flexible neurovascular coupling dynamics that inform the cognitive context. This framework demonstrates how multimodal neuroimaging is necessary to probe the neurometabolic foundations of cognition, with implications for the study of neuropsychiatric disorders.

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.

Segmenting electroencephalography wires reduces radiofrequency shielding artifacts in simultaneous electroencephalography and functional magnetic resonance imaging at 7 T

Purpose

Simultaneous scalp electroencephalography and functional magnetic resonance imaging (EEG‐fMRI) enable noninvasive assessment of brain function with high spatial and temporal resolution. However, at ultra‐high field, the data quality of both modalities is degraded by mutual interactions. Here, we thoroughly investigated the radiofrequency (RF) shielding artifact of a state‐of‐the‐art EEG‐fMRI setup, at 7 T, and design a practical solution to limit this issue.

Methods

Electromagnetic field simulations and MR measurements assessed the shielding effect of the EEG setup, more specifically the EEG wiring. The effectiveness of segmenting the wiring with resistors to reduce the transmit field disruption was evaluated on a wire‐only EEG model and a simulation model of the EEG cap.

Results

The EEG wiring was found to exert a dominant effect on the disruption of the transmit field, whose intensity varied periodically as a function of the wire length. Breaking the electrical continuity of the EEG wires into segments shorter than one quarter RF wavelength in air (25 cm at 7 T) reduced significantly the RF shielding artifacts. Simulations of the EEG cap with segmented wires indicated similar improvements for a moderate increase of the power deposition.

Conclusion

We demonstrated that segmenting the EEG wiring into shorter lengths using commercially available nonmagnetic resistors is effective at reducing RF shielding artifacts in simultaneous EEG‐fMRI. This prevents the formation of RF‐induced standing waves, without substantial specific absorption rate (SAR) penalties, and thereby enables benefiting from the functional sensitivity boosts achievable at ultra‐high field.

Tracking of rigid head motion during MRI using an EEG system

Purpose

To demonstrate a novel method for tracking of head movements during MRI using electroencephalography (EEG) hardware for recording signals induced by native imaging gradients.

Theory and Methods

Gradient switching during simultaneous EEG–fMRI induces distortions in EEG signals, which depend on subject head position and orientation. When EEG electrodes are interconnected with high‐impedance carbon wire loops, the induced voltages are linear combinations of the temporal gradient waveform derivatives. We introduce head tracking based on these signals (CapTrack) involving 3 steps: (1) phantom scanning is used to characterize the target sequence and a fast calibration sequence; (2) a linear relation between changes of induced signals and head pose is established using the calibration sequence; and (3) induced signals recorded during target sequence scanning are used for tracking and retrospective correction of head movement without prolonging the scan time of the target sequence. Performance of CapTrack is compared directly to interleaved navigators.

Results

Head‐pose tracking at 27.5 Hz during echo planar imaging (EPI) was demonstrated with close resemblance to rigid body alignment (mean absolute difference: [0.14 0.38 0.15]‐mm translation, [0.30 0.27 0.22]‐degree rotation). Retrospective correction of 3D gradient‐echo imaging shows an increase of average edge strength of 12%/−0.39% for instructed/uninstructed motion with CapTrack pose estimates, with a tracking interval of 1561 ms and high similarity to interleaved navigator estimates (mean absolute difference: [0.13 0.33 0.12] mm, [0.28 0.15 0.22] degrees).

Conclusion

Motion can be estimated from recordings of gradient switching with little or no sequence modification, optionally in real time at low computational burden and synchronized to image acquisition, using EEG equipment already found at many research institutions.

Simultaneous EEG-fMRI: What Have We Learned and What Does the Future Hold?

Simultaneous EEG-fMRI has developed into a mature measurement technique in the past 25 years. During this time considerable technical and analytical advances have been made, enabling valuable scientific contributions to a range of research fields. This review will begin with an introduction to the measurement principles involved in EEG and fMRI and the advantages of combining these methods. The challenges faced when combining the two techniques will then be considered. An overview of the leading application fields where EEG-fMRI has made a significant contribution to the scientific literature and emerging applications in EEG-fMRI research trends is then presented.

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.

Inferring Macroscale Brain Dynamics via Fusion of Simultaneous EEG-fMRI

Advances in the instrumentation and signal processing for simultaneously acquired electroencephalography and functional magnetic resonance imaging (EEG-fMRI) have enabled new ways to observe the spatiotemporal neural dynamics of the human brain. Central to the utility of EEG-fMRI neuroimaging systems are the methods for fusing the two data streams, with machine learning playing a key role. These methods can be dichotomized into those that are symmetric and asymmetric in terms of how the two modalities inform the fusion. Studies using these methods have shown that fusion yields new insights into brain function that are not possible when each modality is acquired separately. As technology improves and methods for fusion become more sophisticated, the future of EEG-fMRI for noninvasive measurement of brain dynamics includes mesoscale mapping at ultrahigh magnetic resonance fields, targeted perturbation-based neuroimaging, and using deep learning to uncover nonlinear representations that link the electrophysiological and hemodynamic measurements.

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.

The relationship between EEG and fMRI connectomes is reproducible across simultaneous EEG-fMRI studies from 1.5T to 7T

Both electroencephalography (EEG) and functional Magnetic Resonance Imaging (fMRI) are non-invasive methods that show complementary aspects of human brain activity. Despite measuring different proxies of brain activity, both the measured blood-oxygenation (fMRI) and neurophysiological recordings (EEG) are indirectly coupled. The electrophysiological and BOLD signal can map the underlying functional connectivity structure at the whole brain scale at different timescales. Previous work demonstrated a moderate but significant correlation between resting-state functional connectivity of both modalities, however there is a wide range of technical setups to measure simultaneous EEG-fMRI and the reliability of those measures between different setups remains unknown. This is true notably with respect to different magnetic field strengths (low and high field) and different spatial sampling of EEG (medium to high-density electrode coverage). Here, we investigated the reproducibility of the bimodal EEG-fMRI functional connectome in the most comprehensive resting-state simultaneous EEG-fMRI dataset compiled to date including a total of 72 subjects from four different imaging centers. Data was acquired from 1.5T, 3T and 7T scanners with simultaneously recorded EEG using 64 or 256 electrodes. We demonstrate that the whole-brain monomodal connectivity reproducibly correlates across different datasets and that a moderate crossmodal correlation between EEG and fMRI connectivity of r ≈ 0.3 can be reproducibly extracted in low- and high-field scanners. The crossmodal correlation was strongest in the EEG-β frequency band but exists across all frequency bands. Both homotopic and within intrinsic connectivity network (ICN) connections contributed the most to the crossmodal relationship. This study confirms, using a considerably diverse range of recording setups, that simultaneous EEG-fMRI offers a consistent estimate of multimodal functional connectomes in healthy subjects that are dominantly linked through a functional core of ICNs across spanning across the different timescales measured by EEG and fMRI. This opens new avenues for estimating the dynamics of brain function and provides a better understanding of interactions between EEG and fMRI measures. This observed level of reproducibility also defines a baseline for the study of alterations of this coupling in pathological conditions and their role as potential clinical markers.

EEG Microstates Predict Concurrent fMRI Dynamic Functional Connectivity States

Brain functional connectivity measured by resting-state fMRI varies over multiple time scales, and recurrent dynamic functional connectivity (dFC) states have been identified. These have been found to be associated with different cognitive and pathological states, with potential as disease biomarkers, but their neuronal underpinnings remain a matter of debate. A number of recurrent microstates have also been identified in resting-state EEG studies, which are thought to represent the quasi-simultaneous activity of large-scale functional networks reflecting time-varying brain states. Here, we hypothesized that fMRI-derived dFC states may be associated with these EEG microstates. To test this hypothesis, we quantitatively assessed the ability of EEG microstates to predict concurrent fMRI dFC states in simultaneous EEG-fMRI data collected from healthy subjects at rest. By training a random forests classifier, we found that the four canonical EEG microstates predicted fMRI dFC states with an accuracy of 90%, clearly outperforming alternative EEG features such as spectral power. Our results indicate that EEG microstates analysis yields robust signatures of fMRI dFC states, providing evidence of the electrophysiological underpinnings of dFC while also further supporting that EEG microstates reflect the dynamics of large-scale brain networks.

7T Epilepsy Task Force Consensus Recommendations on the Use of 7T MRI in Clinical Practice

Identifying a structural brain lesion on MRI has important implications in epilepsy and is the most important factor that correlates with seizure freedom after surgery in patients with drug-resistant focal onset epilepsy. However, at conventional magnetic field strengths (1.5 and 3T), only approximately 60%-85% of MRI examinations reveal such lesions. Over the last decade, studies have demonstrated the added value of 7T MRI in patients with and without known epileptogenic lesions from 1.5 and/or 3T. However, translation of 7T MRI to clinical practice is still challenging, particularly in centers new to 7T, and there is a need for practical recommendations on targeted use of 7T MRI in the clinical management of patients with epilepsy. The 7T Epilepsy Task Force-an international group representing 21 7T MRI centers with experience from scanning over 2,000 patients with epilepsy-would hereby like to share its experience with the neurology community regarding the appropriate clinical indications, patient selection and preparation, acquisition protocols and setup, technical challenges, and radiologic guidelines for 7T MRI in patients with epilepsy. This article mainly addresses structural imaging; in addition, it presents multiple nonstructural MRI techniques that benefit from 7T and hold promise as future directions in epilepsy. Answering to the increased availability of 7T MRI as an approved tool for diagnostic purposes, this article aims to provide guidance on clinical 7T MRI epilepsy management by giving recommendations on referral, suitable 7T MRI protocols, and image interpretation.

EEG Microstates Predict Concurrent fMRI Dynamic Functional Connectivity States

Brain functional connectivity measured by resting-state fMRI varies over multiple time scales, and recurrent dynamic functional connectivity (dFC) states have been identified. These have been found to be associated with different cognitive and pathological states, with potential as disease biomarkers, but their neuronal underpinnings remain a matter of debate. A number of recurrent microstates have also been identified in resting-state EEG studies, which are thought to represent the quasi-simultaneous activity of large-scale functional networks reflecting time-varying brain states. Here, we hypothesized that fMRI-derived dFC states may be associated with these EEG microstates. To test this hypothesis, we quantitatively assessed the ability of EEG microstates to predict concurrent fMRI dFC states in simultaneous EEG-fMRI data collected from healthy subjects at rest. By training a random forests classifier, we found that the four canonical EEG microstates predicted fMRI dFC states with an accuracy of 90%, clearly outperforming alternative EEG features such as spectral power. Our results indicate that EEG microstates analysis yields robust signatures of fMRI dFC states, providing evidence of the electrophysiological underpinnings of dFC while also further supporting that EEG microstates reflect the dynamics of large-scale brain networks.

Alpha/beta power decreases track the fidelity of stimulus-specific information

Massed synchronised neuronal firing is detrimental to information processing. When networks of task-irrelevant neurons fire in unison, they mask the signal generated by task-critical neurons. On a macroscopic level, such synchronisation can contribute to alpha/beta (8-30 Hz) oscillations. Reducing the amplitude of these oscillations, therefore, may enhance information processing. Here, we test this hypothesis. Twenty-one participants completed an associative memory task while undergoing simultaneous EEG-fMRI recordings. Using representational similarity analysis, we quantified the amount of stimulus-specific information represented within the BOLD signal on every trial. When correlating this metric with concurrently-recorded alpha/beta power, we found a significant negative correlation which indicated that as post-stimulus alpha/beta power decreased, stimulus-specific information increased. Critically, we found this effect in three unique tasks: visual perception, auditory perception, and visual memory retrieval, indicating that this phenomenon transcends both stimulus modality and cognitive task. These results indicate that alpha/beta power decreases parametrically track the fidelity of both externally-presented and internally-generated stimulus-specific information represented within the cortex.

Investigating the variability of cardiac pulse artifacts across heartbeats in simultaneous EEG-fMRI recordings: A 7T study

Electroencephalography (EEG) recordings performed in magnetic resonance imaging (MRI) scanners are affected by complex artifacts caused by heart function, often termed pulse artifacts (PAs). PAs can strongly compromise EEG data quality, and remain an open problem for EEG-fMRI. This study investigated the properties and mechanisms of PA variability across heartbeats, which has remained largely unaddressed to date, and evaluated its impact on PA correction approaches. Simultaneous EEG-fMRI was performed at 7T on healthy participants at rest or under visual stimulation, with concurrent recordings of breathing and cardiac activity. PA variability was found to contribute to EEG variance with more than 500 μV² at 7T, which extrapolates to 92 μV² at 3T. Clustering analyses revealed that PA variability not only is linked to variations in head position/orientation, as previously hypothesized, but also, and more importantly, to the respiratory cycle and to heart rate fluctuations. The latter mechanisms are associated to short-timescale variability (even across consecutive heartbeats), and their importance varied across EEG channels. In light of this PA variability, three PA correction techniques were compared: average artifact subtraction (AAS), optimal basis sets (OBS), and an approach based on K-means clustering. All methods allowed the recovery of visual evoked potentials from the EEG data; nonetheless, OBS and K-means tended to outperform AAS, likely due to the inability of the latter in modeling short-timescale variability. Altogether, these results offer novel insights into the dynamics and underlying mechanisms of the pulse artifact, with important consequences for its correction, relevant to most EEG-fMRI applications.

Quantitative Analyses Help in Choosing Between Simultaneous vs. Separate EEG and fMRI

Simultaneous registration of scalp electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) is considered an attractive approach for studying brain function non-invasively. It combines the better spatial resolution of fMRI with the better temporal resolution of EEG, but comes at the cost of increased measurement artifact and the accompanying artifact preprocessing. This paper presents a study of the residual signal quality based on temporal signal to noise ratio (TSNR) for fMRI and fast Fourier transform (FFT) for EEG, after optimized conventional signal preprocessing. Measurements outside the magnetic resonance imaging scanner and inside the scanner prior to and during image acquisition were compared. For EEG, frequency and region dependent significant effects on FFT squared amplitudes were observed between separately vs. simultaneously recorded EEG and fMRI, with larger effects during image acquisition than without image acquisition inside the scanner bore. A graphical user interface was developed to aid in quality checking these measurements. For fMRI, separately recorded EEG-fMRI revealed relatively large areas with a significantly higher TSNR in right occipital and parietal regions and in the cingulum, compared to separately recorded EEG-fMRI. Simultaneously recorded EEG-fMRI showed significantly higher TSNR in inferior occipital cortex, diencephalon and brainstem, compared to separately recorded EEG-fMRI. Quantification of EEG and fMRI signals showed significant, but sometimes subtle, changes between separate compared to simultaneous EEG-fMRI measurements. To avoid interference with the experiment of interest in a simultaneous EEG-fMRI measurement, it seems warranted to perform a quantitative evaluation to ensure that there are no such uncorrectable effects present in regions or frequencies that are of special interest to the research question at hand.

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.

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.

Examples of sub-millimeter, 7T, T1-weighted EPI datasets acquired with the T123DEPI sequence

These imaging data are examples of sub-millimeter resolution T₁-weighted EPI (Echo Planar Imaging) acquired using the T₁23DEPI (T₁-imaging with 2 3D-EPIs) sequence [1]; functional MRI data with matching resolution and distortion, and MP2RAGE (Magnetization Prepared 2 Rapid Acquisition Gradient Echoes) anatomical images [2], from the same subjects. Data from two protocols and subjects presented in the paper describing the sequence [1] are made available here:

1)

0.8 mm protocol: whole brain, axial T₁23DEPI T₁-weighted images; a 5-minute fMRI run with the same orientation and 27 mm coverage in the slice selection direction, covering the primary visual cortex. fMRI data were acquired while the volunteer viewed a flashing checkerboard stimulus; the unsmoothed GLM results of the fMRI and a 0.64 mm resolution MP2RAGE from the same subject. These data are from Experiment 3 in [1]

2)

0.7 mm protocol: partial brain T₁23DEPI T₁-weighted images with longer or shorter readouts; matching coronal echo planar images again acquired while viewing a flashing checkerboard stimulus and a 0.64 mm whole brain MP2RAGE from the same subject. These data are from Experiment 1 in  [1]

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.

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.

EEG-Informed fMRI: A Review of Data Analysis Methods

The simultaneous acquisition of electroencephalography (EEG) with functional magnetic resonance imaging (fMRI) is a very promising non-invasive technique for the study of human brain function. Despite continuous improvements, it remains a challenging technique, and a standard methodology for data analysis is yet to be established. Here we review the methodologies that are currently available to address the challenges at each step of the data analysis pipeline. We start by surveying methods for pre-processing both EEG and fMRI data. On the EEG side, we focus on the correction for several MR-induced artifacts, particularly the gradient and pulse artifacts, as well as other sources of EEG artifacts. On the fMRI side, we consider image artifacts induced by the presence of EEG hardware inside the MR scanner, and the contamination of the fMRI signal by physiological noise of non-neuronal origin, including a review of several approaches to model and remove it. We then provide an overview of the approaches specifically employed for the integration of EEG and fMRI when using EEG to predict the blood oxygenation level dependent (BOLD) fMRI signal, the so-called EEG-informed fMRI integration strategy, the most commonly used strategy in EEG-fMRI research. Finally, we systematically review methods used for the extraction of EEG features reflecting neuronal phenomena of interest.

EEG synchronization measures predict epilepsy-related BOLD-fMRI fluctuations better than commonly used univariate metrics

OBJECTIVE

We hypothesize that the hypersynchronization associated with epileptic activity is best described by EEG synchronization measures, and propose to use these as predictors of epilepsy-related BOLD fluctuations.

METHODS

We computed the phase synchronization index (PSI) and global field synchronization (GFS), within two frequency bands, a broadband (1-45 Hz) and a narrower band focused on the presence of epileptic activity (3-10 Hz). The associated epileptic networks were compared with those obtained using conventional unitary regressors and two power-weighted metrics (total power and root mean square frequency), on nine simultaneous EEG-fMRI datasets from four epilepsy patients, exhibiting inter-ictal epileptiform discharges (IEDs).

RESULTS

The average PSI within 3-10 Hz achieved the best performance across several measures reflecting reliability in all datasets. The results were cross-validated through electrical source imaging of the IEDs. The applicability of PSI when no IEDs are recorded on the EEG was evaluated on three additional patients, yielding partially plausible networks in all cases.

CONCLUSIONS

Epileptic networks can be mapped based on the EEG PSI metric within an IED-specific frequency band, performing better than commonly used EEG metrics.

SIGNIFICANCE

This is the first study to investigate EEG synchronization measures as potential predictors of epilepsy-related BOLD fluctuations.

Online Reduction of Artifacts in EEG of Simultaneous EEG-fMRI Using Reference Layer Adaptive Filtering (RLAF)

Simultaneous electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) allow us to study the active human brain from two perspectives concurrently. Signal processing based artifact reduction techniques are mandatory for this, however, to obtain reasonable EEG quality in simultaneous EEG-fMRI. Current artifact reduction techniques like average artifact subtraction (AAS), typically become less effective when artifact reduction has to be performed on-the-fly. We thus present and evaluate a new technique to improve EEG quality online. This technique adds up with online AAS and combines a prototype EEG-cap for reference recordings of artifacts, with online adaptive filtering and is named reference layer adaptive filtering (RLAF). We found online AAS + RLAF to be highly effective in improving EEG quality. Online AAS + RLAF outperformed online AAS and did so in particular online in terms of the chosen performance metrics, these being specifically alpha rhythm amplitude ratio between closed and opened eyes (3-45% improvement), signal-to-noise-ratio of visual evoked potentials (VEP) (25-63% improvement), and VEPs variability (16-44% improvement). Further, we found that EEG quality after online AAS + RLAF is occasionally even comparable with the offline variant of AAS at a 3T MRI scanner. In conclusion RLAF is a very effective add-on tool to enable high quality EEG in simultaneous EEG-fMRI experiments, even when online artifact reduction is necessary.

Unimodal Versus Bimodal EEG-fMRI Neurofeedback of a Motor Imagery Task

Neurofeedback is a promising tool for brain rehabilitation and peak performance training. Neurofeedback approaches usually rely on a single brain imaging modality such as EEG or fMRI. Combining these modalities for neurofeedback training could allow to provide richer information to the subject and could thus enable him/her to achieve faster and more specific self-regulation. Yet unimodal and multimodal neurofeedback have never been compared before. In the present work, we introduce a simultaneous EEG-fMRI experimental protocol in which participants performed a motor-imagery task in unimodal and bimodal NF conditions. With this protocol we were able to compare for the first time the effects of unimodal EEG-neurofeedback and fMRI-neurofeedback versus bimodal EEG-fMRI-neurofeedback by looking both at EEG and fMRI activations. We also propose a new feedback metaphor for bimodal EEG-fMRI-neurofeedback that integrates both EEG and fMRI signal in a single bi-dimensional feedback (a ball moving in 2D). Such a feedback is intended to relieve the cognitive load of the subject by presenting the bimodal neurofeedback task as a single regulation task instead of two. Additionally, this integrated feedback metaphor gives flexibility on defining a bimodal neurofeedback target. Participants were able to regulate activity in their motor regions in all NF conditions. Moreover, motor activations as revealed by offline fMRI analysis were stronger during EEG-fMRI-neurofeedback than during EEG-neurofeedback. This result suggests that EEG-fMRI-neurofeedback could be more specific or more engaging than EEG-neurofeedback. Our results also suggest that during EEG-fMRI-neurofeedback, participants tended to regulate more the modality that was harder to control. Taken together our results shed first light on the specific mechanisms of bimodal EEG-fMRI-neurofeedback and on its added-value as compared to unimodal EEG-neurofeedback and fMRI-neurofeedback.

How to Build a Hybrid Neurofeedback Platform Combining EEG and fMRI

Multimodal neurofeedback estimates brain activity using information acquired with more than one neurosignal measurement technology. In this paper we describe how to set up and use a hybrid platform based on simultaneous electroencephalography (EEG) and functional magnetic resonance imaging (fMRI), then we illustrate how to use it for conducting bimodal neurofeedback experiments. The paper is intended for those willing to build a multimodal neurofeedback system, to guide them through the different steps of the design, setup, and experimental applications, and help them choose a suitable hardware and software configuration. Furthermore, it reports practical information from bimodal neurofeedback experiments conducted in our lab. The platform presented here has a modular parallel processing architecture that promotes real-time signal processing performance and simple future addition and/or replacement of processing modules. Various unimodal and bimodal neurofeedback experiments conducted in our lab showed high performance and accuracy. Currently, the platform is able to provide neurofeedback based on electroencephalography and functional magnetic resonance imaging, but the architecture and the working principles described here are valid for any other combination of two or more real-time brain activity measurement technologies.

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.

EEG-fMRI Gradient Artifact Correction by Multiple Motion-Related Templates

OBJECTIVES

In simultaneous electroencephalography (EEG) and functional magnetic resonance imaging (fMRI), artifacts on the EEG arise from the switching of magnetic field gradients in the MR scanner. These artifacts depend on head position, and are, therefore, difficult to remove in the presence of subject motion. In this study, gradient artifacts are modeled by multiple templates extracted from externally recorded motion information.

METHODS

Gradient artifact correction was performed in EEG-fMRI recordings by estimating artifactual templates modulated by slowly varying splines, as well as head position information. The EEG signal quality was then compared following two common methods: averaged artifact subtraction (AAS) and optimal basis sets (OBS).

RESULTS

Artifact correction using multiple templates estimated from splines or motion time courses outperformed the existing AAS and OBS approaches, as quantified by root-mean-square power across gradient epochs. Improvements were mostly seen in posterior EEG channels, where most of the residual artifacts are seen following the AAS and OBS methods. Residual spectral power was comparable to that of EEG signals recorded without fMRI scanning.

CONCLUSION

Gradient artifacts can be well modeled by multiple templates estimated from head position information, resulting in an effective artifact removal.

SIGNIFICANCE

This method can facilitate EEG-fMRI of uncooperative subjects in whom motion is inevitable, for example, to investigate high-frequency EEG activity in which gradient artifacts are particularly prominent.

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.

Ballistocardiogram artifact correction taking into account physiological signal preservation in simultaneous EEG-fMRI

The ballistocardiogram (BCG) artifact is currently one of the most challenging in the EEG acquired concurrently with fMRI, with correction invariably yielding residual artifacts and/or deterioration of the physiological signals of interest. In this paper, we propose a family of methods whereby the EEG is decomposed using Independent Component Analysis (ICA) and a novel approach for the selection of BCG-related independent components (ICs) is used (PROJection onto Independent Components, PROJIC). Three ICA-based strategies for BCG artifact correction are then explored: 1) BCG-related ICs are removed from the back-reconstruction of the EEG (PROJIC); and 2-3) BCG-related ICs are corrected for the artifact occurrences using an Optimal Basis Set (OBS) or Average Artifact Subtraction (AAS) framework, before back-projecting all ICs onto EEG space (PROJIC-OBS and PROJIC-AAS, respectively). A novel evaluation pipeline is also proposed to assess the methods performance, which takes into account not only artifact but also physiological signal removal, allowing for a flexible weighting of the importance given to physiological signal preservation. This evaluation is used for the group-level parameter optimization of each algorithm on simultaneous EEG-fMRI data acquired using two different setups at 3T and 7T. Comparison with state-of-the-art BCG correction methods showed that PROJIC-OBS and PROJIC-AAS outperformed the others when priority was given to artifact removal or physiological signal preservation, respectively, while both PROJIC-AAS and AAS were in general the best choices for intermediate trade-offs. The impact of the BCG correction on the quality of event-related potentials (ERPs) of interest was assessed in terms of the relative reduction of the standard error (SE) across trials: 26/66%, 32/62% and 18/61% were achieved by, respectively, PROJIC, PROJIC-OBS and PROJIC-AAS, for data collected at 3T/7T. Although more significant improvements were achieved at 7T, the results were qualitatively comparable for both setups, which indicate the wide applicability of the proposed methodologies and recommendations.

Spurious correlations in simultaneous EEG-fMRI driven by in-scanner movement

Simultaneous EEG-fMRI provides an increasingly attractive research tool to investigate cognitive processes with high temporal and spatial resolution. However, artifacts in EEG data introduced by the MR scanner still remain a major obstacle. This study, employing commonly used artifact correction steps, shows that head motion, one overlooked major source of artifacts in EEG-fMRI data, can cause plausible EEG effects and EEG-BOLD correlations. Specifically, low-frequency EEG (<20Hz) is strongly correlated with in-scanner movement. Accordingly, minor head motion (<0.2mm) induces spurious effects in a twofold manner: Small differences in task-correlated motion elicit spurious low-frequency effects, and, as motion concurrently influences fMRI data, EEG-BOLD correlations closely match motion-fMRI correlations. We demonstrate these effects in a memory encoding experiment showing that obtained theta power (~3-7Hz) effects and channel-level theta-BOLD correlations reflect motion in the scanner. These findings highlight an important caveat that needs to be addressed by future EEG-fMRI studies.

Presurgical brain mapping in epilepsy using simultaneous EEG and functional MRI at ultra-high field: feasibility and first results

OBJECTIVES

The aim of this study was to demonstrate that eloquent cortex and epileptic-related hemodynamic changes can be safely and reliably detected using simultaneous electroencephalography (EEG)-functional magnetic resonance imaging (fMRI) recordings at ultra-high field (UHF) for clinical evaluation of patients with epilepsy.

MATERIALS AND METHODS

Simultaneous EEG-fMRI was acquired at 7 T using an optimized setup in nine patients with lesional epilepsy. According to the localization of the lesion, mapping of eloquent cortex (language and motor) was also performed in two patients.

RESULTS

Despite strong artifacts, efficient correction of intra-MRI EEG could be achieved with optimized artifact removal algorithms, allowing robust identification of interictal epileptiform discharges. Noise-sensitive topography-related analyses and electrical source localization were also performed successfully. Localization of epilepsy-related hemodynamic changes compatible with the lesion were detected in three patients and concordant with findings obtained at 3 T. Local loss of signal in specific regions, essentially due to B 1 inhomogeneities were found to depend on the geometric arrangement of EEG leads over the cap.

CONCLUSION

These results demonstrate that presurgical mapping of epileptic networks and eloquent cortex is both safe and feasible at UHF, with the benefits of greater spatial resolution and higher blood-oxygenation-level-dependent sensitivity compared with the more traditional field strength of 3 T.