Simultaneous EEG and fMRI in the macaque monkey at 4.7 Tesla

Simultaneous electroencephalography (EEG)/functional magnetic resonance imaging (fMRI) acquisition can identify the brain networks involved in generating specific EEG patterns. Yet, the combination of these methodologies is hampered by strong artifacts that arise due to electromagnetic interference during magnetic resonance (MR) image acquisition. Here, we report corrections of the gradient-induced artifact in phantom measurements and in experiments with an awake behaving macaque monkey during fMRI acquisition at a magnetic field strength of 4.7 T. Ninety-one percent of the amplitude of a 10 microV, 10 Hz phantom signal could successfully be recovered without phase distortions. Using this method, we were able to extract the monkey EEG from scalp recordings obtained during MR image acquisition. Visual evoked potentials could also be reliably identified. In conclusion, simultaneous EEG/fMRI acquisition is feasible in the macaque monkey preparation at 4.7 T and holds promise for investigating the neural processes that give rise to particular EEG patterns.     

EEG quality during simultaneous functional MRI of interictal epileptiform discharges

This article concerns the evaluation of the quality of interictal epileptiform EEG discharges recorded throughout simultaneous echo planar imaging (EPI). BOLD (blood oxygen level dependent) functional MRI (fMRI) images were acquired continuously on a patient with intractable epilepsy. EEG was sampled simultaneously, during and after imaging, with removal of pulse and imaging artifacts by subtraction of channel-specific running averages. Contiguous EEG epochs recorded with and without fMRI (fMRI+ve vs. fMRI−ve) were next randomized and presented to two blinded observers. Epileptiform discharges were identified retrospectively, and comparison was made in terms of the number of identified events, their amplitude, and spatiotemporal distribution. A spectral analysis was also performed on the EEG. In the randomized comparison of EEG segments, 80 (fMRI+ve) vs. 69 (fMRI−ve) discharges were noted with good interobserver agreement (69%). There were no significant differences in amplitude or spatio-temporal distribution. Comparison of the events detected and measured by two expert observers demonstrated that the Interictal Epileptiform Discharge (IED) characteristics were indistinguishable with and without scanning. We review briefly the existing literature on EEG recording quality for combined EEG/fMRI.     

Exploring the feasibility of simultaneous electroencephalography/functional magnetic resonance imaging at 7 T

The increased blood oxygenation level-dependent contrast available at high field makes the implementation of combined EEG/fMRI experiments at 7 T highly worthwhile from the point of view of fMRI data quality, but the higher field poses greater technical challenges for achieving good quality EEG data. A study of the feasibility of recording EEG signals from human subjects at 7 T using a commercially available, MR-compatible EEG system has therefore been carried out. This involved systematic measurement of the sources of noise in EEG recordings made in the 7 T scanner and measurement of RF heating effects on a gel phantom in the presence of a 32-electrode EEG cap. Having found no significant safety concerns and identified a set-up (involving switching off the magnet's cryo-cooler pumps and mounting the EEG amplifier on a cantilever) that limited scanner-induced noise, combined EEG/fMRI experiments employing visual stimulation were then successfully carried out on two human subjects. With the use of beamformer-based analysis of the EEG data, driven responses and alpha-band, event-related desynchronisation were identified in both subjects.     

Distributed analysis of simultaneous EEG-fMRI time-series: modeling and interpretation issues

Functional magnetic resonance imaging (fMRI) and electroencephalography (EEG) represent brain activity in terms of a reliable anatomical localization and a detailed temporal evolution of neural signals. Simultaneous EEG-fMRI recordings offer the possibility to greatly enrich the significance and the interpretation of the single modality results because the same neural processes are observed from the same brain at the same time. Nonetheless, the different physical nature of the measured signals by the two techniques renders the coupling not always straightforward, especially in cognitive experiments where spatially localized and distributed effects coexist and evolve temporally at different temporal scales.     

Continuous EEG-fMRI in patients with partial epilepsy and focal interictal slow-wave discharges on EEG

Purpose

To verify whether in patients with partial epilepsy and routine electroenecephalogram (EEG) showing focal interictal slow-wave discharges without spikes combined EEG–functional magnetic resonance imaging (fMRI) would localize the corresponding epileptogenic focus, thus providing reliable information on the epileptic source.

Methods

Eight patients with partial epileptic seizures whose routine scalp EEG recordings on presentation showed focal interictal slow-wave activity underwent EEG–fMRI. EEG data were continuously recorded for 24 min (four concatenated sessions) from 18 scalp electrodes, while fMRI scans were simultaneously acquired with a 1.5-Tesla magnetic resonance imaging (MRI) scanner. After recording sessions and MRI artefact removal, EEG data were analyzed offline. We compared blood oxygen level-dependent (BOLD) signal changes on fMRI with EEG recordings obtained at rest and during activation (with and without focal interictal slow-wave discharges).

Results

In all patients, when the EEG tracing showed the onset of focal slow-wave discharges on a few lateralized electrodes, BOLD-fMRI activation in the corresponding brain area significantly increased. We detected significant concordance between focal EEG interictal slow-wave discharges and focal BOLD activation on fMRI. In patients with lesional epilepsy, the epileptogenic area corresponded to the sites of increased focal BOLD signal.

Conclusions

Even in patients with partial epilepsy whose standard EEGs show focal interictal slow-wave discharges without spikes, EEG–fMRI can visualize related focal BOLD activation thus providing useful information for pre-surgical planning.     

Estimation of the effective and functional human cortical connectivity with structural equation modeling and directed transfer function applied to high-resolution EEG

Different brain imaging devices are presently available to provide images of the human functional cortical activity, based on hemodynamic, metabolic or electromagnetic measurements. However, static images of brain regions activated during particular tasks do not convey the information of how these regions are interconnected. The concept of brain connectivity plays a central role in the neuroscience, and different definitions of connectivity, functional and effective, have been adopted in literature. While the functional connectivity is defined as the temporal coherence among the activities of different brain areas, the effective connectivity is defined as the simplest brain circuit that would produce the same temporal relationship as observed experimentally among cortical sites. The structural equation modeling (SEM) is the most used method to estimate effective connectivity in neuroscience, and its typical application is on data related to brain hemodynamic behavior tested by functional magnetic resonance imaging (fMRI), whereas the directed transfer function (DTF) method is a frequency-domain approach based on both a multivariate autoregressive (MVAR) modeling of time series and on the concept of Granger causality.

This study presents advanced methods for the estimation of cortical connectivity by applying SEM and DTF on the cortical signals estimated from high-resolution electroencephalography (EEG) recordings, since these signals exhibit a higher spatial resolution than conventional cerebral electromagnetic measures. To estimate correctly the cortical signals, we used a subject's multicompartment head model (scalp, skull, dura mater, cortex) constructed from individual MRI, a distributed source model and a regularized linear inverse source estimates of cortical current density. Before the application of SEM and DTF methodology to the cortical waveforms estimated from high-resolution EEG data, we performed a simulation study, in which different main factors (signal-to-noise ratio, SNR, and simulated cortical activity duration, LENGTH) were systematically manipulated in the generation of test signals, and the errors in the estimated connectivity were evaluated by the analysis of variance (ANOVA). The statistical analysis returned that during simulations, both SEM and DTF estimators were able to correctly estimate the imposed connectivity patterns under reasonable operative conditions, that is, when data exhibit an SNR of at least 3 and a LENGTH of at least 75 s of nonconsecutive EEG recordings at 64 Hz of sampling rate.

Hence, effective and functional connectivity patterns of cortical activity can be effectively estimated under general conditions met in any practical EEG recordings, by combining high-resolution EEG techniques and linear inverse estimation with SEM or DTF methods. We conclude that the estimation of cortical connectivity can be performed not only with hemodynamic measurements, but also with EEG signals treated with advanced computational techniques.     

Integration of fMRI, NIROT and ERP for studies of human brain function

Different methods of assessing human brain function possess specific advantages and disadvantages compared to others, but it is believed that combining different approaches will provide greater information than can be obtained from each alone. For example, functional magnetic resonance imaging (fMRI) has good spatial resolution but poor temporal resolution, whereas the converse is true for electrophysiological recordings (event-related potentials or ERPs). In this review of recent work, we highlight a novel approach to combining these modalities in a manner designed to increase information on the origins and locations of the generators of specific ERPs and the relationship between fMRI and ERP signals. Near infrared imaging techniques have also been studied as alternatives to fMRI and can be readily integrated with simultaneous electrophysiological recordings. Each of these modalities may in principle be also used in so-called steady-state acquisitions in which the correlational structure of signals from the brain may be analyzed to provide new insights into brain function.     

The BOLD response and the gamma oscillations respond differently than evoked potentials: an interleaved EEG-fMRI study

Background  

The integration of EEG and fMRI is attractive because of their complementary precision regarding time and space. But the relationship between the indirect hemodynamic fMRI signal and the more direct EEG signal is uncertain. Event-related EEG responses can be analyzed in two different ways, reflecting two different kinds of brain activity: evoked, i.e. phase-locked to the stimulus, such as evoked potentials, or induced, i.e. non phase-locked to the stimulus such as event-related oscillations. In order to determine which kind of EEG activity was more closely related with fMRI, EEG and fMRI signals were acquired together, while subjects were presented with two kinds of rare events intermingled with frequent distractors. Target events had to be signaled by pressing a button and Novel events had to be ignored.     

Special designed RF-antenna with integrated non-invasive carbon electrodes for simultaneous magnetic resonance imaging and electroencephalography acquisition at 7T

The construction of a high quality MR RF-antenna with incorporated EEG electrodes for simultaneous MRI and EEG acquisition is presented. The antenna comprises an active decoupled surface coil for receiving the MR signal and a whole body coil for transmitting the excitation RF pulses. The surface coil offers a high signal-to-noise ratio required for fMRI application and the whole body coil has a good B(1) excitation profile, which enables the application of homogeneous RF pulses. Non-invasive carbon electrodes are used in order to minimise the magnetic susceptibility artefacts that occur upon application of conductive materials. This dedicated set-up is compared to a standard set-up being a linear birdcage coil and commercially available Ag/AgCl electrodes. As the acquired EEG signals are heavily disturbed by the gradient switching, intelligent filtering is applied to obtain a clean EEG signal.     

Nonlinear responses of cerebral blood volume, blood flow and blood oxygenation signals during visual stimulation

Hemodynamic-based functional magnetic resonance imaging (fMRI) techniques provide a great utility for noninvasive functional brain mapping. However, because the hemodynamic signals reflect underlying neural activity indirectly, characterization of these signals following brain activation is essential for experimental design and data interpretation. In this report, the linear (or nonlinear) responses to neuronal activation of three hemodynamic parameters based primarily on changes of cerebral blood volume, blood flow and blood oxygenation were investigated by testing these hemodynamic responses' additivity property. Using a recently developed fMRI technique that acquires vascular space occupancy (VASO), arterial spin labeling (ASL) perfusion and blood oxygenation level-dependent (BOLD) signals simultaneously, the additivity property of the three hemodynamic responses in human visual cortex was assessed using various visual stimulus durations. Experiments on healthy volunteers showed that all three hemodynamic-weighted signals responded nonlinearly to stimulus durations less than 4 s, with the degree of nonlinearity becoming more severe as the stimulus duration decreased. Vascular space occupancy and ASL perfusion signals showed similar nonlinearity properties, whereas the BOLD signal was the most nonlinear. These data suggest that caution should be taken in the interpretation of hemodynamic-based signals in fMRI.     

An independent component analysis-based approach on ballistocardiogram artifact removing

Interest about simultaneous electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) data acquisition has rapidly increased during the last years because of the possibility that the combined method offers to join temporal and spatial resolution, providing in this way a powerful tool to investigate spontaneous and evoked brain activities. However, several intrinsic features of MRI scanning become sources of artifacts on EEG data. Noise sources of a highly predictable nature such as those related to the pulse MRI sequence and those determined by magnetic gradient switching during scanning do not represent a major problem and can be easily removed. On the contrary, the ballistocardiogram (BCG) artifact, a large signal visible on all EEG traces and related to cardiac activity inside the magnetic field, is determined by sources that are not fully stereotyped and causing important limitations in the use of artifact-removing strategies. Recently, it has been proposed to use independent component analysis (ICA) to remove BCG artifact from EEG signals. ICA is a statistical algorithm that allows blind separation of statistically independent sources when the only available information is represented by their linear combination. An important drawback with most ICA algorithms is that they exhibit a stochastic behavior: each run yields slightly different results such that the reliability of the estimated sources is difficult to assess. In this preliminary report, we present a method based on running the FastICA algorithm many times with slightly different initial conditions. Clustering structure in the signal space of the obtained components provides us with a new way to assess the reliability of the estimated sources.     

Real-time MR artifacts filtering during continuous EEG/fMRI acquisition

The purpose of this study was the development of a real-time filtering procedure of MRI artifacts in order to monitor the EEG activity during continuous EEG/fMRI acquisition. The development of a combined EEG and fMRI technique has increased in the past few years. Preliminary “spike-triggered” applications have been possible because in this method, EEG knowledge was only necessary to identify a trigger signal to start a delayed fMRI acquisition. In this way, the two methods were used together but in an interleaved manner. In real simultaneous applications, like event-related fMRI study, artifacts induced by MRI events on EEG traces represent a substantial obstacle for a right analysis. Up until now, the methods proposed to solve this problem are mainly based on procedures to remove post-processing artifacts without the possibility to control electrophysiological behavior of the patient during fMRI scan. Moreover, these methods are not characterized by a strong “prior knowledge” of the artifact, which is an imperative condition to avoid any loss of information on the physiological signals recovered after filtering. In this work, we present a new method to perform simultaneous EEG/fMRI study with real-time artifacts filtering characterized by a procedure based on a preliminary analytical study of EPI sequence parameters-related EEG-artifact shapes. Standard EEG equipment was modified in order to work properly during ultra-fast MRI acquisitions. Changes included: high-performance acquisition device; electrodes/cap/wires/cables materials and geometric design; shielding box for EEG signal receiver; optical fiber link; and software. The effects of the RF pulse and time-varying magnetic fields were minimized by using a correct head cap wires-locked environment montage and then removed during EEG/fMRI acquisition with a subtraction algorithm that takes in account the most significant EPI sequence parameters. The on-line method also allows a further post-processing utilization.     

How the NPX data format handles EEG data acquired simultaneously with fMRI

There is a growing interest in combining EEG and (f)MRI data as they provide complementary information: EEG is characterized by a high temporal resolution but poor spatial one, while fMRI is characterized by a high spatial resolution but low temporal one. However, while a standard file format for storing EEG data is available since over a decade, it does not fulfill the needs of modern protocols and devices such as those involved in simultaneous EEG and fMRI recordings. The main reasons are the limited bit resolution, some difficulties encountered in handling and storing acquisition events or trace markers for off-line analyses and the impossibility to add some protocol-specific information that is not considered in the actual data formats. This, among others, hinders the release of free analysis software and makes it difficult to share data across different laboratories as every research unit develops its own tools according to its needs, stores data in proprietary formats and a lot of time is spent building software applications for converting data from one format to another. The NPX (NeuroPhysiological signals in eXtensible Markup Language) data format was defined to overcome these and other limitations, and here its main characteristics are reported as well as how some typical problems occurring in simultaneous EEG-fMRI recordings are also treated. Many tools based on the NPX technology can be freely downloaded, including a tool for removing artifacts occurring during simultaneous EEG-fMRI recordings.     

Brain oscillatory activity during motor imagery in EEG-fMRI coregistration

The purpose of the present work was to investigate the correlation between topographical changes in brain oscillatory activity and the blood oxygenation level-dependent (BOLD) signal during a motor imagery (MI) task using electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) coregistration.     

Evaluating gradient artifact correction of EEG data acquired simultaneously with fMRI

Simultaneous electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) has become a widely used application in spite of EEG perturbations due to electromagnetic interference in the MR environment. The most prominent and disturbing artifacts in the EEG are caused by the alternating magnetic fields (gradients) of the MR scanner. Different methods for gradient artifact correction have been developed. Here we propose an approach for the systematic evaluation and comparison of these gradient artifact correction methods. Exemplarily, we evaluate different algorithms all based on artifact template subtraction--the currently most established means of gradient artifact removal. We introduce indices for the degree of gradient artifact reduction and physiological signal preservation. The combination of both indices was used as a measure for the overall performance of gradient artifact removal and was shown to be useful in identifying problems during artifact removal. We demonstrate that the evaluation as proposed here allows to reveal frequency-band specific performance differences among the algorithms. This emphasizes the importance of carefully selecting the artifact correction method appropriate for the respective case.     

Blind source separation of multichannel electroencephalogram based on wavelet transform and ICA

Combination of the wavelet transform and independent component analysis （ICA） was employed for blind source separation （BSS） of multichannel electroencephalogram （EEG）. After denoising the original signals by discrete wavelet transform, high frequency components of some noises and artifacts were removed from the original signals. The denoised signals were reconstructed again for the purpose of ICA, such that the drawback that ICA cannot distinguish noises from source signals can be overcome effectively. The practical processing results showed that this method is an effective way to BSS of multichannel EEG. The method is actually a combination of wavelet transform with adaptive neural network, so it is also useful for BBS of other complex signals.     

Integration of EEG source imaging and fMRI during continuous viewing of natural movies

Electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) are noninvasive neuroimaging tools which can be used to measure brain activity with excellent temporal and spatial resolution, respectively. By combining the neural and hemodynamic recordings from these modalities, we can gain better insight into how and where the brain processes complex stimuli, which may be especially useful in patients with different neural diseases. However, due to their vastly different spatial and temporal resolutions, the integration of EEG and fMRI recordings is not always straightforward. One fundamental obstacle has been that paradigms used for EEG experiments usually rely on event-related paradigms, while fMRI is not limited in this regard. Therefore, here we ask whether one can reliably localize stimulus-driven EEG activity using the continuously varying feature intensities occurring in natural movie stimuli presented over relatively long periods of time. Specifically, we asked whether stimulus-driven aspects in the EEG signal would be co-localized with the corresponding stimulus-driven BOLD signal during free viewing of a movie. Secondly, we wanted to integrate the EEG signal directly with the BOLD signal, by estimating the underlying impulse response function (IRF) that relates the BOLD signal to the underlying current density in the primary visual area (V1). We made sequential fMRI and 64-channel EEG recordings in seven subjects who passively watched 2-min-long segments of a James Bond movie. To analyze EEG data in this natural setting, we developed a method based on independent component analysis (ICA) to reject EEG artifacts due to blinks, subject movement, etc., in a way unbiased by human judgment. We then calculated the EEG source strength of this artifact-free data at each time point of the movie within the entire brain volume using low-resolution electromagnetic tomography (LORETA). This provided for every voxel in the brain (i.e., in 3D space) an estimate of the current density at every time point. We then carried out a correlation between the time series of visual contrast changes in the movie with that of EEG voxels. We found the most significant correlations in visual area V1, just as seen in previous fMRI studies (Bartels A, Zeki, S, Logothetis NK. Natural vision reveals regional specialization to local motion and to contrast-invariant, global flow in the human brain. Cereb Cortex 2008;18(3):705–717), but on the time scale of milliseconds rather than of seconds. To obtain an estimate of how the EEG signal relates to the BOLD signal, we calculated the IRF between the BOLD signal and the estimated current density in area V1. We found that this IRF was very similar to that observed using combined intracortical recordings and fMRI experiments in nonhuman primates. Taken together, these findings open a new approach to noninvasive mapping of the brain. It allows, firstly, the localization of feature-selective brain areas during natural viewing conditions with the temporal resolution of EEG. Secondly, it provides a tool to assess EEG/BOLD transfer functions during processing of more natural stimuli. This is especially useful in combined EEG/fMRI experiments, where one can now potentially study neural-hemodynamic relationships across the whole brain volume in a noninvasive manner.     

Multimodal integration of EEG, MEG and fMRI data for the solution of the neuroimage puzzle

In this paper, advanced methods for the modeling of human cortical activity from combined high-resolution electroencephalography (EEG), magnetoencephalography (MEG) and functional magnetic resonance imaging (fMRI) data are presented. These methods include a subject's multicompartment head model (scalp, skull, dura mater, cortex) constructed from magnetic resonance images, multidipole source model and regularized linear inverse source estimates of cortical current density. Determination of the priors in the resolution of the linear inverse problem was performed with the use of information from the hemodynamic responses of the cortical areas as revealed by block-designed (strength of activated voxels) fMRI. Examples of the application of these methods to the estimation of the time varying cortical current density activity in selected region of interest (ROI) are presented for movement-related high-resolution EEG data.     

Comparison of functional MR-venography and EPI-BOLD fMRI at 1.5 t

Functional magnetic resonance imaging (fMRI) of the brain using blood oxygenation level dependent (BOLD) contrast relies on the changes of paramagnetic deoxyhemoglobin concentration, which affects brain parenchyma and draining venous vessels. These changes in deoxyhemoglobin concentration in venous vessels can also be monitored using a high-resolution susceptibility-based MR-venography technique. Four volunteers participated in the study in which functional MR-venograms were compared with conventional echo-planar imaging (EPI)-BOLD-fMRI. In all cases, small venous vessels could be identified close to the areas of activation detected by conventional fMRI. In the venograms, task performance (finger tapping) resulted in a loss of venous vessel contrast compared to the resting state, which is consistent with a local decrease of deoxyhemoglobin concentration. MR-venography allows a direct visualization of the BOLD-effect at high spatial resolution. In combination with conventional fMRI, this technique may help to separate the contribution of brain parenchyma and venous vessels in fMRI studies.