EEG-Informed fMRI: A Review of Data Analysis Methods
- 1ISR-Lisboa/LARSyS and Department of Bioengineering, Instituto Superior Técnico – Universidade de Lisboa, Lisbon, Portugal
- 2Department of Neurophysiology, Centro Hospitalar Psiquiátrico de Lisboa, Lisbon, Portugal
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.
Electroencephalography (EEG) is by far the most commonly used technique to study brain function. Its millisecond temporal resolution allows the adequate sampling of the rapidly changing electrical dynamics of neuronal populations. The power spectrum of scalp EEG signals follows approximately 1/fβ power law distributions, characteristic of most scale-free dynamics found in nature. In addition, power peaks reflecting rhythmic activity, or brain oscillations at specific frequencies, can be superimposed. In general, the most relevant brain oscillations are found in the following conventional frequency bands: delta (0.5–4 Hz), theta (4–8 Hz), alpha (8–13 Hz), beta (13–30 Hz), and gamma (above 30 Hz) (Niedermeyer and Lopes Da Silva, 2005; Akay, 2006). Both rhythmic and arrhythmic ongoing wide-band EEG activities can be modulated by external stimuli or tasks (e.g., cognitive, somatosensory, motor, auditory, visual), inducing time-locked and/or phase-locked activation of specific neuronal populations. In both cases, stimuli- or task-specific electrical potentials – event-related potentials (ERPs) – are generated.
Scalp EEG signals result from the mixture of propagating electric potential fluctuations (the so-called local field potentials, LFPs), mainly reflecting the postsynaptic activity of large populations of cortical pyramidal cells. These cells display a geometric configuration and an orientation in relation to the skull that are favorable to the constructive summation of the associated electrical current sources (Michel et al., 2004). Modest contributions from multi-unit activity (MUA) associated with action potentials, as well as from glial cells, have also been reported (Nunez and Silberstein, 2000; Dale and Halgren, 2001). Unfortunately, EEG source reconstruction is a notably ill-posed inverse problem, yielding a non-unique solution for all admissible potential distributions, because the number of sources is typically much larger than the number of sensors (Michel and Murray, 2012). Available strategies for the practical resolution of this problem have been proposed, but generally provide limited spatial resolution and a heterogeneous spatial sensitivity particularly for superficial versus deep sources.
Alternatively, blood oxygenation level dependent (BOLD) functional magnetic resonance imaging (fMRI) (Ogawa et al., 1990; Belliveau et al., 1991; Kwong et al., 1992) can be used to map brain activity with excellent spatial localization power, based on neurovascular coupling mechanisms producing hemodynamic changes associated with neuronal activity. Although these are not yet completely understood, both feedforward and feedback pathways have been identified, which attempt to respond to the increased demands for oxygen and glucose of brain cells. The net increase in blood oxygenation upon brain activation leads to a bulk increase in the BOLD signal (Jezzard and Toosy, 2005), due to magnetic susceptibility differences induced by the varying concentration of paramagnetic deoxyhemoglobin relative to diamagnetic oxyhemoglobin. The temporal resolution of fMRI is limited by the relatively slow hemodynamic response, with BOLD changes being delayed by several seconds relative to the onset of neuronal activity (Logothetis et al., 2001; Lewin, 2003). Usually, whole-brain fMRI measurements are performed with a spatial resolution of a few millimeters and a temporal resolution of a few seconds. Nevertheless, a trade-off between spatial and temporal resolution is possible by manipulating the image acquisition parameters; this is ultimately limited only by the available signal-to-noise ratio (SNR), which increases with the magnetic field strength. The use of ultrahigh field strengths (such as 7 Tesla) and highly accelerated image acquisition sequences (such as simultaneous multi-slice techniques) now allow sub-millimeter spatial resolution and whole-brain coverage in under half a second (Feinberg and Setsompop, 2013; van der Zwaag et al., 2016).
Since EEG and fMRI are the two most commonly used noninvasive functional neuroimaging techniques, and because they exhibit highly complementary characteristics, their multimodal integration has been actively sought (Laufs, 2012; Jorge et al., 2014; Murta et al., 2015). It was originally motivated by the need to accurately and non-invasively map epileptic networks in patients with drug-resistant focal epilepsy undergoing pre-surgical evaluation (Ives et al., 1993; Lemieux et al., 2001; Gotman et al., 2006; LeVan and Gotman, 2009; Gotman and Pittau, 2011; Murta et al., 2012), and it was soon extended into studies of normal brain function as well. Subject safety was addressed in early studies (Lemieux et al., 1997) and the consequent hardware modifications of the EEG recording apparatus (Goldman et al., 2000; Vasios et al., 2006) were very effective at preventing any relevant side effects in the large number of recordings performed up to this date. Moreover, increasingly efficient signal processing tools have been developed for removing the MR-induced EEG artifacts (Allen et al., 1998, 2000; Niazy et al., 2005). Although to a lesser extent, attention has also been devoted to the BOLD signal distortions caused by the presence of EEG materials (Krakow et al., 2000; Stevens et al., 2007; Mullinger et al., 2008a). Modality-specific artifacts are also present, which may further confound simultaneous EEG-fMRI analyses: while eye movements and blinks, muscle activity and bad channels are typically captured by EEG (Chaumon et al., 2015), BOLD-fMRI signals are contaminated by fluctuations of non-neuronal origin (Birn, 2012; Murphy et al., 2013). Motivated by the greater sensitivity and spatial resolution/specificity of intracranial EEG (icEEG), simultaneous icEEG-fMRI recordings can now also be performed (Vulliemoz et al., 2011), after both safety (Carmichael et al., 2010) and data quality (Carmichael et al., 2012) concerns have been addressed. Compared with scalp EEG, icEEG can capture more subtle and local features of electrophysiological activity, and may therefore offer novel insights into the relationship of such features with concurrent BOLD signal changes.
When combining EEG and fMRI, an integration strategy must be chosen. Symmetrical approaches are in principle ideal to make the most of the multimodal information, since they do not constrain any of the specific modalities, potentially creating biased estimations. These can be roughly divided into model-based and data-driven EEG-fMRI fusion techniques. While the former present great challenges namely regarding model inversion (Valdes-Sosa et al., 2009), the latter have been mostly based on independent component analysis (ICA) and canonical correlation analysis (CCA) [for reviews, refer to Rosa et al. (2010a), Lei et al. (2012)]. Information theory approaches have also been employed for the integration of EEG-fMRI data, based on the study of neuronal population codes explicitly taking into account the experimentally observed stimulus-response signal probability distributions [for reviews, refer to Panzeri et al. (2008), Ostwald and Bagshaw (2011)]. In the context of EEG-fMRI, such approaches allow the quantitative evaluation of the amount of information embodied in EEG and fMRI features (separately and jointly), to determine which features are more discriminative of the brain activity under study, and the extent to which such information overlaps across the two modalities (Ostwald et al., 2010, 2011, 2012). Because of their relative conceptual and methodological simplicity, asymmetrical approaches are by far the most common (Rosa et al., 2010a). These comprise: (1) fMRI-driven EEG, whereby the activated brain regions identified with fMRI are used as spatial constraints for the EEG source reconstruction problem (for a recent review on this topic, refer to Lei et al., 2015); and (2) EEG-informed fMRI, whereby the brain activity recorded with EEG is used to predict hemodynamic changes measured with fMRI (Gotman et al., 2006; Mulert and Lemieux, 2009; Gotman and Pittau, 2011; Laufs, 2012; Jorge et al., 2014; Murta et al., 2015). The main processing pipeline steps of EEG-correlated fMRI analyses in general are illustrated in Figure 1. Regardless of the integration strategy, the process of going from the first step (data acquisition) to the final result (a brain activity/connectivity map) is always confronted with several challenges.
FIGURE 1. Main processing pipeline steps in EEG-informed fMRI analysis. Data quality is typically addressed first, taking into account modality-specific artifacts, as well as those that are induced by one technique on the other. The EEG phenomenon of interest is then identified and appropriate features are extracted, from which a BOLD signal predictor is derived for the localization of EEG-related BOLD-fMRI changes.
In this review, we cover the most important data analysis steps in EEG-informed fMRI. We start by considering the issues regarding the data quality of both the EEG and the BOLD recordings, together with the respective artifact correction techniques. We then overview the multiple approaches used for EEG-informed fMRI analyses, which rely on extracting appropriate features from the EEG data in order to derive a predictor of the BOLD signal associated with the brain activity under study.
EEG Data Quality
In this section, we provide a comprehensive characterization of multiple MR-induced EEG artifacts, as well as the respective correction techniques. The most important EEG artifacts are the gradient artifact (GA) and pulse artifact (PA), as well as head motion artifacts; other artifact sources have also been reported in the MR-environment, including the Helium cooling pump, the patient ventilation system and the room lights. In general, artifact correction techniques can be subdivided into three main types. First, the most commonly used approach consists of the time-domain subtraction of artifact templates. Typically, a template is generated for each artifact occurrence by averaging across neighboring occurrences, assuming that the artifact changes slowly over time. Second, blind source separation (BSS) techniques, particularly ICA, have also been used to separate EEG artifact sources from neuronal sources. In temporal ICA, an N × M EEG dataset, with N channels and M time points, is decomposed into a linear combination of L independent components (ICs) with N × L weights (Bell and Sejnowski, 1995; Lee et al., 1999). By removing artifact sources when back-reconstructing the EEG data to its original space, an artifact-free EEG is obtained, which makes the identification of artifact-related ICs crucial for accurate EEG cleaning. Third, some approaches rely on dedicated hardware to measure the artifact waveforms directly, followed by their subtraction from the artifact-contaminated EEG signal.
In the next two sections, we focus on the two most important MR-induced EEG artifacts, the GA and the PA, by first characterizing them and then surveying the methods used for their correction; these are summarized in Tables 1 and 2, respectively. In the third section, other MR-environment-related artifacts are also presented and discussed, as well as EEG-specific artifacts originating from sources unrelated with the concurrent acquisition of MR images.
TABLE 1. List of GA correction methods. These methods can be roughly divided into four main approaches: filtering, AAS-based (AAS and its several variations), ICA-based (as well as IVA, an extension of ICA to multiple datasets) and hardware-based.
TABLE 2. List of PA correction methods. These methods can be roughly divided into four main approaches: AAS-based (AAS and OBS), ICA-based (with selection and removal of PA-related ICs or combination with AAS/OBS), hardware-based and others.
During fMRI acquisitions, the magnetic field inside the scanner changes over time due to the application of time-varying magnetic field gradients (Allen et al., 2000; Niazy et al., 2005). According to Faraday’s law of induction, these will induce an electromotive force within the conducting loop formed by the subject’s head and the EEG hardware (electrodes, wires and amplification system). A spurious voltage is hence generated on the EEG electrodes which is usually called GA or imaging artifact (Grouiller et al., 2007). The GA waveform generated by a commonly used 2D multi-slice echo-planar imaging (EPI) sequence for fMRI acquisition at 3 Tesla is displayed in Figure 2. The amplitude of such artifact can be one hundred times greater than that of the physiological EEG signal. More importantly, its spectral content usually overlaps with frequency bands of interest of the EEG, making its removal resorting to basic filtering strategies generally inappropriate. It should be noted, however, that in a few cases, MR sequences can be designed to induce artifacts at nonessential EEG frequencies, hence allowing a Fourier domain correction (Hoffmann et al., 2000; Grouiller et al., 2007). One study employed electromagnetic theory to develop a physical model for the GA (Yan et al., 2009). This model allowed the derivation of the optimal head orientation and position inside the MRI scanner in subsequent experiments, which was shown to significantly minimize the impact of MR gradients on the EEG recordings (Mullinger et al., 2011).
FIGURE 2. Illustration of the gradient artifact (GA) generated by a 2D multi-slice EPI sequence. (Top) 5 s traces of raw EEG data from 10 channels. At approximately 15 s, the fMRI acquisition starts, completely obscuring any neuronal activity being recorded. (Bottom) The zoomed red box shows the high-amplitude electrical potentials generated by the time-varying gradients applied during the acquisition of four image slices using 2D multi-slice EPI; due to their clear periodicity and precise timing, these artifacts can be accurately corrected using channel-specific average template subtraction techniques.
Due to its strong deterministic component, most of the methods found in the literature for GA removal follow an average artifact template subtraction (AAS) approach known as the image artifact reduction (IAR) method (Allen et al., 2000), whereby artifact templates are derived for each occurrence by averaging across multiple fMRI slice or volume epochs. Prior to averaging, a high-pass filtering is typically applied (< 1 Hz) to remove slow artifact amplitude modulations; alternatively, these can be modeled as a low-frequency sinusoidal wave (Grouiller et al., 2007). An adaptive noise cancellation (ANC) filter is then usually applied to remove residual artifacts. Additionally, principal component analysis (PCA) can be performed prior to ANC, in order to capture the residual variance within a limited number of principal components (PCs), which are then linearly fitted to, and subtracted from, each residual occurrence. Alternatively, PCA can also be applied directly to the EEG data (Negishi et al., 2004), followed by a similar procedure to IAR to correct for residual artifacts (Niazy et al., 2005). Although artifact templates can be defined based on the periods of either the imaging slice or volume, a combination of both slice and volume artifact template subtraction has been found to yield superior artifact correction, particularly at frequencies in the gamma band (Gonçalves et al., 2007).
Despite the remarkable precision and reproducibility of the GA, the accuracy of artifact template computation, and consequently the efficiency of average artifact subtraction, are heavily dependent on a precise synchronization between the EEG and fMRI acquisitions (Allen et al., 2000; Niazy et al., 2005). If not synchronized, a high-precision alignment of each artifact occurrence to a given reference period may still be obtained by up-sampling the EEG signal (Allen et al., 2000; Ritter et al., 2007). Ultimately, MR sequence timing parameters, namely the volume and slice repetition times as well as the time delay between the acquisition of the first volume and respective slice, can be estimated at the original EEG sampling rate prior to artifact correction (Gonçalves et al., 2007; Koskinen and Vartiainen, 2009). However, this time alignment model assumes that artifact occurrences are highly reproducible, which is only guaranteed if electrode motion is negligible.
Electrode motion is mainly due to the subjects’ head movement and system vibrations (Yan et al., 2010), and it introduces additional variance in the GA over time (Mullinger and Bowtell, 2011). For the artifact template to capture such variance, several variations of the IAR method were proposed in the literature. In particular, the artifact template can be obtained by weighting each occurrence based on its time difference from (Ritter et al., 2007), or spectral similarity to (Freyer et al., 2009), the occurrence to be corrected for. Alternatively, artifact occurrences can be clustered based on similarity measures, and an artifact template is then computed for each cluster (de Munck et al., 2013). Artifact occurrences can also be modeled as a linear combination of an artifact template and its derivatives in a truncated Taylor’s series expansion (Wan et al., 2006a). Alternatively, head motion can be estimated from the fMRI time series, and the resulting realignment parameters used to group artifact occurrences based on the degree of motion contamination; group-specific artifact templates are then derived (Moosmann et al., 2009).
The use of ICA for GA correction poses the challenge of selecting the associated ICs. For example, ICs can be probed for artifact-related waveforms by appropriately thresholding the correlation coefficients between each IC and an artifact template (Grouiller et al., 2007). Importantly, ICA can be used even in the presence of timing errors. In this case, an average artifact template is first computed for each occurrence and randomly time-jittered within a reasonable range. The resulting templates are then subjected to ICA, and the timing errors are expected to be reflected in the resulting ICs (Ryali et al., 2009). An extension of ICA for multiple datasets – independent vector analysis (IVA) – has also been proposed for GA correction (Acharjee et al., 2015). In this method, EEG data are first segmented using the fMRI scanning triggers as time-locking events, and the multiple segmented channels are then entered as multiple datasets for IVA. The resulting IVA components are maximally dependent across channels, and maximally independent across segments for a given channel. GA sources are thus separated from those of different origins, as in ICA, but also taking into account spatial dependencies of GA waveforms across electrodes, and hence allowing for a more accurate, channel-wise GA source estimation (Acharjee et al., 2015).
New EEG caps have been developed, incorporating a second set of EEG electrodes that overlay those in contact with the scalp, separated by a reference layer (Chowdhury et al., 2014). While the electrodes located below the reference layer are physically attached to the scalp as in a standard EEG cap, capturing a mixture of artifacts and brain signals, the electrodes located above the reference layer are electrically isolated from the scalp so that they do not pick up brain signals and hence measure only MR-induced artifacts. These include not only the GA, but also the PA, as well as motion-driven artifacts. By taking the difference between the signals measured by the scalp electrodes and the corresponding reference electrodes, an artifact-free EEG signal can be obtained.
In a different approach, an MR-compatible camera-tracker device is attached to the subject’s head in order to monitor head motion. Prospective motion correction (PMC) of the EEG signal can then be performed based on the head translation and rotation parameters estimated along the three main axes. These estimates are typically used to improve MR image quality, by updating the specifications of the RF pulses and MR gradients during the image acquisition in real-time (recent reviews on these approaches can be found in Maclaren et al., 2013; Zaitsev et al., 2016). In the work by Maziero et al. (2016), these motion parameters are additionally used to model, and subsequently regress out, the motion-induced voltages on the concurrently acquired EEG.
At least two comparative studies of multiple GA correction methods can be found in the literature (Grouiller et al., 2007; Ritter et al., 2007). By taking into account not only the extent of artifact removal, but also the associated degradation of the EEG signal of interest, both studies arrived at the conclusion that the performance of each method is highly dependent on the type of EEG activity of interest, namely in terms of its frequency band. Appropriate balancing between artifact and physiological signal removal can be achieved by directly comparing performance measures between scan and non-scan periods (Freyer et al., 2009). This requires the acquisition of artifact-free EEG data during MR-silent periods, which can be obtained using a dedicated EPI sequence characterized by a long repetition time (TR; 4070 ms) (Anami et al., 2003). Ritter et al. (2007) investigated the impact on the GA correction algorithms of using different parameters (namely the number of artifact windows to build the template, and whether PCA should be employed to remove residual artifacts or not), and found sets of optimal parameter values for each GA correction method.
In our experience, the deterministic component of GA renders the use of the original IAR method suitable for most applications, avoiding the challenging selection of GA-related ICs in ICA-based methods, or the introduction of additional hardware to the typically intricate EEG-fMRI setups. Whenever facing unsatisfactory results, one can start by varying the number of artifact occurrences used to generate the artifact template, and thus differently weighting the amount of artifact reduction (lower number of occurrences) relative to the preservation of EEG physiological signal (higher number of occurrences). Variations of the original IAR can also be tested; however, in our experience these variations were only successful at capturing GA variability at the slice level and lack sensitivity at the volume level.
The PA, also commonly referred to as the ballistocardiogram artifact, is currently one of the most challenging artifacts in the EEG acquired concurrently with fMRI, mainly due to its non-stationary nature determined by the temporal variability of the cardiac pulse. As discussed in Yan et al. (2010) and Mullinger et al. (2013), several mechanisms contribute to the PA by inducing voltage changes on the scalp EEG electrodes due to interactions with the strong, static magnetic field of the MRI scanner, namely: (1) bulk head rotation due to cardiac blood ejection (Bonmassar et al., 2002); (2) scalp expansion due to arterial pulsation (Debener et al., 2008); and (3) Hall effect caused by the pulsatile flow of the blood, which is an electrically conductive fluid (Tenforde et al., 1983). One study found that most of the artifact variance was explained by flow-induced Hall voltage and pulse-driven head rotation (Mullinger et al., 2013). Importantly, the PA amplitude significant increases with the static magnetic field strength, B0 (Debener et al., 2008; Neuner et al., 2013), such that it can severely hamper the visual inspection of typical EEG behaviors at 7 Tesla (Jorge et al., 2015a,b).
Because the PA is roughly periodic, with artifact occurrences being approximately time-locked with the cardiac cycle, an AAS algorithm can be employed, similarly to the case of the GA. In this case, an artifact template is extracted from the EEG signal by averaging across multiple cardiac cycles, followed by a time-domain subtraction procedure (illustrated in Figure 3) (Allen et al., 1998). Unlike the GA, however, the much greater variability over time must be taken into account for effectively reducing the PA. This can be achieved by computing the temporal PCA over all time-locked occurrences of the artifact in order to build an optimal basis set (OBS), comprising a given number of PCs that explain the PA variance to some extent (Niazy et al., 2005); this basis set is then fitted to, and subtracted from, each artifact occurrence.
FIGURE 3. Illustration of the AAS technique to correct for the pulse artifact (PA). The EEG traces are shown for three channels, before and after correction, as well as for the ECG channel, over a time period of 10 s including 11 artifact occurrences. The segmented windows for each artifact occurrence (blue and red boxes) were averaged to compute the artifact template (red trace); this was then subtracted from the corresponding artifact occurrence (red box), yielding the artifact-corrected signal (green trace).
If ICA can accurately separate the PA sources from other sources contributing to the EEG signal, a PA-free EEG signal may be obtained by reconstructing the EEG without the PA-related ICs (Bénar et al., 2003; Srivastava et al., 2005; Mantini et al., 2007a). However, the objective and accurate classification of PA-related ICs remains a major concern, and several criteria can be found in the literature for that purpose; these are reviewed in Vanderperren et al. (2010) and their performance is compared in Abreu et al. (2016b). Briefly, such criteria can be roughly divided into five main types, whereby PA-related ICs are identified by: (1) thresholding the correlation coefficients between each IC and the simultaneously acquired ECG signal (Mantini et al., 2007a), or a PA template (Srivastava et al., 2005); (2) probing the auto-correlation function of each IC for peaks located at the distance between two consecutive QRS peaks (Deburchgraeve et al., 2008); (3) identifying ICs that exhibit spectral peaks at cardiac-related frequencies, exploring the periodic nature of the PA (Vanderperren et al., 2007); (4) thresholding the peak-to-peak (PTP) values of the back-projected QRS-triggered EEG data when using each IC individually, under the hypothesis that PA-related ICs should exhibit higher PTP values (Vanderperren et al., 2010); and (5) thresholding the amount of artifact variance that is explained by each source (Debener et al., 2008). Following the rationale of the latter, we have recently proposed the PROJIC (PROJection onto Independent Components) algorithm, whereby the average PA waveform of each channel is first projected onto the IC space by means of the corresponding un-mixing matrix estimated using ICA. The power of each projection, computed as the squared sum for each time instant, is then clustered by k-means; ICs assigned to high-powered clusters are deemed PA-related (Abreu et al., 2016a).
Due to the relative advantages and disadvantages of AAS- and ICA-based methods, their combination has been proposed (Debener et al., 2005, 2007). In one method, the OBS method is first used to remove most of the PA contribution from the EEG data and ICA is then employed to further remove residual artifacts by excluding PA-related ICs from the back-reconstruction of the EEG signal. The use of a modified version of OBS in the IC space instead has also been proposed, whereby non-PA-related ICs are corrected for the PA, while PA-related ICs are removed (Liu et al., 2012). Such procedure was motivated by the hypothesis that the artifact contributes to all ICs in varying degrees due to its non-stationary nature. Following a similar rationale, we have recently proposed the application of AAS and OBS in the IC space, but only to PA-related ICs, which are subsequently kept in the EEG signal reconstruction together with the original non-PA-related ICs (Abreu et al., 2016a). This approach aims to preserve the physiological signal as much as possible without compromising artifact removal efficiency.
Motion sensors can be used to monitor pulse-driven as well as voluntary head movements during concurrent EEG-fMRI acquisitions; the sensor signals can be used to estimate head motion and associated EEG artifacts, which can then be subtracted from the EEG data (Bonmassar et al., 2002). Different kinds of motion sensors have been used for reducing the PA in simultaneous EEG-fMRI recordings. PA waveforms can be obtained by means of loops of carbon-fiber wire that are physically attached to, but electrically isolated from, the subject’s head (Masterton et al., 2007). More recently, a simple modification of a standard EEG cap has been proposed, whereby four of the existing electrodes are insulated from the scalp and directly connected to the reference electrode so that they can be used as motion sensors (Jorge et al., 2015a). PA waveforms can also be extracted from a layer of electrodes added to a standard EEG cap and separated from the scalp by two insulating layers (Xia et al., 2014a,b). A practical approach in this case consists of defining a set of insulated electrodes spatially surrounding each uninsulated electrode, in order to build an accurate estimate of the local PA (Xia et al., 2014b). Following this rationale, complex methods to determine the optimal minimum number of electrodes to insulate have been developed based on the spatial redundancy of PA measured from neighboring electrodes (Xia et al., 2014a). In an alternative approach, an MR-compatible camera-tracker device has also been shown to improve PA correction, by first converting the measured six motion parameters into velocities, and then modeling the PA-induced EEG voltages as a linear combination of the low-pass filtered velocities (LeVan et al., 2013). In a similar approach, which, however, avoids additional hardware in order to minimize the complexity of the simultaneous EEG-fMRI setup, the PA is estimated directly from the EEG as a linear combination of several cardiac-related harmonics (Krishnaswamy et al., 2016).
Simultaneous EEG-fMRI measurements of event-related activity are typically used for assessing the performance of a given correction method, based on metrics computed from the ERPs such as: the inter-trial variability (Vanderperren et al., 2010), the SNR (Debener et al., 2007), and the difference between the ERPs extracted from the inside-MR EEG datasets and those that are obtained from the PA-free outside-MR EEG data (Mantini et al., 2007a). When the frequency content of an event-related EEG signal is known, the power within that frequency band can also be computed before and after PA correction (Xia et al., 2014a). In resting-state EEG-fMRI (rs-fMRI) studies, the PA occurrences can be compared before and after correction based on their root mean square (RMS) or PTP values (Chowdhury et al., 2014). Additionally, the total spectral power within windows around the cardiac fundamental frequency and its first harmonics can be computed, in order to quantify the amount of PA that is removed (Liu et al., 2012). Similarly to the evaluation pipeline described in Freyer et al. (2009) for the GA, we have proposed to assess the trade-off between PA and physiological signal reductions by computing ratios of the power over specific frequency bands and then linearly combining them using a weighting factor that describes the importance given to the preservation of physiological signal relative to artifact correction (Abreu et al., 2016a).
Grouiller et al. (2007) found that AAS was the method of choice if highly accurate QRS detection was achieved. Additionally, and as discussed in Vanderperren et al. (2010), OBS and ICA-based approaches only yielded comparable results if the ICA parameters were fine-tuned. In fact, optimizing the algorithms’ parameters has been shown to critically affect the efficiency and reliability of the subsequent analyses (Vanderperren et al., 2010; Abreu et al., 2016a), particularly at high magnetic field strengths (Debener et al., 2008). According to our own study (Abreu et al., 2016a), AAS was the method exhibiting the second best results in terms of accurately removing the artifact while preserving the physiological signal of interest. In contrast, purely ICA-based methods either resulted in substantial residual artifacts, or significantly distorted the physiological signal. The best results were obtained by combining ICA to separate the PA sources, with AAS to correct the artifact occurrences in the IC space. We believe that this may be a simple but effective solution for PA correction, which does not require additional hardware.
Other Sources of EEG Artifacts
Several artifacts are induced on the EEG in the MR environment, even without an ongoing fMRI acquisition. Typically, these artifacts are caused by electrode motion as a result of MR scanner vibrations associated with the Helium compression pumps used for cooling down MR components, the patient ventilation system, and the room lights (Mullinger et al., 2008a; Mulert and Lemieux, 2009; Nierhaus et al., 2013; Neuner et al., 2014; Rothlübbers et al., 2014). In particular, the Helium pump artifact has been characterized in some systems by prominent peaks in the EEG spectrum around frequencies of 50 and 100 Hz (Rothlübbers et al., 2014). The ventilation system and room lights are reflected in spectral peaks at other specific frequencies, with the former depending on the ventilation level (Nierhaus et al., 2013). While turning off the room lights does not present any clear compromise, switching off the ventilation system may cause patient discomfort. In principle, the Helium pump cooling system can also be turned off; however, this may carry the associated risk of Helium boil-off in certain systems (Mullinger et al., 2008a) and it is not permitted in some clinical sites for safety and procedural reasons. Due to its repetitive nature, the Helium pump artifact can be adequately removed by employing AAS-based approaches (Rothlübbers et al., 2014). Alternatively, PCA has been recursively applied to EEG segments in order to separate the components associated with the Helium pump artifact; components exhibiting a single peak within a frequency range typically spanning the artifact were removed from the data (Kim et al., 2015).
Although not caused by the MR environment, other sources of artifacts that contaminate the EEG must also be accounted for in simultaneous EEG-fMRI studies, namely: eye movements, saccades and blinks, muscle activity, and bad channels (Chaumon et al., 2015). Eye movements are usually picked up by frontal electrodes, although artifactual waveforms may also be observed on distant electrodes (Urigüen and Garcia-Zapirain, 2015). The degree of contamination by this type of artifact is determined by the proximity of the electrodes to the eyes, as well as the direction of the movement (Croft and Barry, 2000). Eye saccades can be either horizontal or vertical, both mainly captured by frontal electrodes, and are characterized by abrupt changes in the EEG amplitude. These changes can be mistaken with EEG activity in the gamma band (Yuval-Greenberg et al., 2008). Eye blinking is also more prominent in frontal electrodes, usually inducing high-frequency artifacts due to its abrupt nature (Croft and Barry, 2000; Urigüen and Garcia-Zapirain, 2015). The EEG topographies associated with eye blinks resemble those of vertical eye saccades (Chaumon et al., 2015). Muscle artifacts result from the myogenic activity caused by contracting muscles, particularly those surrounding the mandible, and their effects depend on the degree of contraction and the type of muscle (Goncharova et al., 2003). Muscle artifacts usually span a wide frequency band of the EEG, although they considerably overlap with beta activity (∼15–30 Hz) (Goncharova et al., 2003; McMenamin et al., 2010). This type of artifacts can also affect a large cortical surface area due to volume conduction of myogenic activity from different head muscles (Urigüen and Garcia-Zapirain, 2015). The so-called bad channels are usually associated with high-impedance electrodes and are typically characterized by strong fluctuations that are uncorrelated with the remaining electrodes (Chaumon et al., 2015).
Some of these types of EEG artifacts are expected to be mainly stationary, making the use of ICA particularly suitable. For this purpose, techniques have been proposed for the automatic selection of the artifactual ICs, which are then removed from the back-reconstruction of the EEG data (for reviews, refer to Vanderperren et al., 2010; Chaumon et al., 2015; Urigüen and Garcia-Zapirain, 2015). Most of those techniques, however, rely on artifact-specific a priori information, rendering them unsuitable for other types of EEG activity (Campos Viola et al., 2009; Nolan et al., 2010; Mognon et al., 2011; Daly et al., 2014; Chaumon et al., 2015).
MR Data Quality
In this section, we start by briefly characterizing MR image artifacts that are induced by the presence of EEG hardware inside the MR scanner. Although unrelated with EEG, signal fluctuations of non-neuronal origin are known to contaminate BOLD-fMRI data, especially in resting-state studies of spontaneous activity. For this reason, such fluctuations are also considered here, including their characterization and an overview of the methods used for their modeling and removal.
EEG-Induced Image Quality Degradation
The presence of EEG hardware during MR image acquisitions is known to degrade image quality, although to a substantially lesser degree when compared with the effects of such acquisitions on the EEG signal (Mulert and Lemieux, 2009; Mullinger and Bowtell, 2011). Since MR image quality is closely related with the amplitude and homogeneity of both static (B0) and oscillating (B1) magnetic fields, and because the EEG system will directly interfere with those magnetic fields, appreciable effects on MR data quality are expected.
Firstly, magnetic susceptibility differences between the EEG system materials and the head tissues induce perturbations on B0, which in turn cause geometric distortions and signal loss in the MR images. These effects scale linearly with both the magnetic susceptibility difference and the field strength (Mullinger et al., 2008b; Jorge et al., 2015b), and depend on the spatial orientation of the material relative to B0 (Krakow et al., 2000). Secondly, the conductive materials of the EEG system, particularly leads outside the EEG cap because of the longer wires connected to them, also cause perturbations on the transmitted B1 field. These alter the effective flip angle of imaging sequences and hence induce local image intensity variations, which are exacerbated if B0 homogeneity is also compromised (Mullinger et al., 2008b). Additionally, B1 field shielding also occurs due to loop currents induced on the electrically conductive materials of the EEG electrodes and leads (Stevens et al., 2007), which is aggravated when using high-density EEG caps (Mullinger and Bowtell, 2011). Overall, these effects of the EEG on B1 lead to image SNR losses, and also to local changes in the specific absorption rate (SAR) of MR image acquisition sequences (Vasios et al., 2006; Stevens et al., 2007). The latter scale linearly with the number of electrodes and the B0 field strength (Angelone et al., 2006). The use of specifically designed low-SAR sequences is therefore crucial to ensure subject safety and comfort in simultaneous EEG-fMRI recordings at ultra-high fields (Nöth et al., 2012).
The electromagnetic noise generated by the EEG recording components can be minimized through appropriate shielding and the use of suitable materials, which should be as diamagnetic as possible (Krakow et al., 2000; Stevens et al., 2007). However, modifications of the electrodes are limited: on the one hand, they must not compromise the overall functionality of the EEG system; on the other hand, they must ensure subject safety, namely by including current-limiting resistors (Lemieux et al., 1997). Moreover, only a minimal quantity of conductive gel should be used, so that it does not induce appreciable image artifacts while still providing acceptable electrode impedance (Krakow et al., 2000). Novel EEG caps have been developed using different technologies in order to minimize the impact of simultaneous EEG-fMRI recordings on subject safety and MR data quality (Vasios et al., 2006), and ultimately to make feasible the recording of high-density EEG data at ultra-high magnetic field strengths (Poulsen et al., 2017).
In addition, B0 and B1 inhomogeneities can be accounted for through post-processing based on B0 and B1 field maps obtained using dedicated MR sequences (Ericsson et al., 1995; Yarnykh, 2007). Although increasingly faster field mapping sequences are being developed (Eggenschwiler et al., 2012), the associated additional scanning time may not be available in conventional imaging settings. Nevertheless, current literature suggests that negligible effects are observed on image quality in humans when using commercially available, low-density EEG caps at field strengths of 3 Tesla (Mullinger et al., 2008b). At 7 Tesla, however, important B1 field distortions have been reported (Jorge et al., 2015b). Most importantly, although the spatial SNR is significantly affected, the temporal SNR (tSNR) of fMRI seems to be relatively preserved, in part because physiological noise is also reduced with the overall signal loss; in fact, no significant differences have been reported when comparing BOLD-fMRI detection sensitivity with or without the EEG cap in place (Luo and Glover, 2012; Klein et al., 2015).
BOLD-fMRI Physiological Noise
BOLD-fMRI signal changes result from contributions from both neuronal and non-neuronal origins. The latter include cardiac and respiratory sources commonly referred to as physiological noise (Hutton et al., 2011; Brooks et al., 2013; Caballero-Gaudes and Reynolds, 2017). The arterial pulsation associated with cardiac function produces brain tissue movements, as well as changes in cerebral blood volume (CBV) and cerebral blood flow (CBF), across the cardiac cycle (Purdon and Weisskoff, 1998; Krüger and Glover, 2001; Greitz et al., 2010). As for respiration, the thoracic modulation within each respiratory cycle produces bulk head motion, as well as changes in B0 (Raj et al., 2001) and in the arterial CO2 partial pressure (Wise et al., 2004). These changes associated with the cardiac and respiratory cycles induce correlated, quasi-periodic BOLD fluctuations, located predominantly near and within large blood vessels and more generally across the brain, respectively (Birn, 2012). The flow of cerebrospinal fluid (CSF) is also modulated by both cardiac and respiratory cycles, resulting in associated BOLD signal changes in CSF-filled regions (Klose et al., 2000). Bulk motion related with the cardiac and respiratory cycles leads to confounds similar to the ones produced by voluntary head motion, even with the use of head restraints (Murphy et al., 2013). Typically, they manifest as correlated signal changes at the edges of the brain and in regions with large spatial variations in image contrast. Additionally, non-periodic BOLD signal fluctuations are also produced due to changes in cardiac rate (Shmueli et al., 2007), as well as in breathing depth and rate leading to changes in the CO2 arterial partial pressure (Birn et al., 2006).
These physiologically driven effects, if left uncorrected, may compromise the analysis of fMRI data, particularly when studying spontaneous activity (Biswal et al., 1995; Cordes et al., 2001; Birn, 2012; Murphy et al., 2013), and they are therefore a concern also in simultaneous EEG-fMRI (Liston et al., 2006; van Houdt et al., 2009; Abreu et al., 2017b) (Figure 4). Since a typical fMRI acquisition sequence uses a TR of ∼2–3 s, aliasing of both cardiac (≈ 1 Hz) and respiratory (≈ 0.3 Hz) fundamental frequencies will inevitably occur (Bhattacharyya and Lowe, 2004), making the use of temporal filtering strategies unsuitable unless a very short TR < 0.4 s is used (Biswal et al., 1996). Such a short TR would allow only partial brain coverage with conventional sequences. However, faster sequences have recently been developed, allowing whole-brain coverage with such short TR values (Feinberg and Setsompop, 2013), which opens up new possibilities for physiological noise characterization and correction. A number of physiological noise correction approaches have been proposed, including both model-based (relying on external physiological recordings and/or the fMRI data itself) and data-driven techniques; these are reviewed next, and summarized in Table 3.
FIGURE 4. Illustration of the presence of physiological-related fluctuations in the BOLD signal. (Top) Structural brain image highlighting the brainstem (red dashed circle), a brain structure located close to major arteries and CSF-filled spaces, thus particularly susceptible to physiological fluctuations. (Middle) Average physiological noise-related BOLD time-course and respective power spectrum computed from a region near the brainstem. (Bottom) Average BOLD time-course across GM; in contrast with the brainstem, the GM time-course presents a clear spectrum, with most of its power located at frequencies below 0.1 Hz.
TABLE 3. List of BOLD-fMRI physiological noise correction methods. These methods can be roughly divided into four main approaches: filtering, physiological recordings-based, image-based and data-driven.
Physiological Recordings-Based Approaches
An important class of physiological noise correction methods relies on recording, concurrently with the fMRI acquisition, external cardiac and respiratory signals by means of appropriate sensors, usually a plethysmograph and a respiratory belt, respectively. A retrospective correction method that works directly in k-space (RETROKCOR) was first introduced, aiming to remove the contribution of physiological noise prior to image reconstruction (Hu et al., 1995). However, most methods operate in the image space, and are generally extensions of the retrospective image-based correction (RETROICOR) method (Glover et al., 2000). The contribution of periodic cardiac and respiratory processes is described by a low-order Fourier expansion in terms of the phases of the cardiac and respiratory signals in relation to the fMRI acquisition time; several model orders have been tested, as well as their interactions (Harvey et al., 2008). These terms (cardiac, respiratory and interactions) are then estimated and regressed out from the fMRI data using a general linear model (GLM) framework. Since RETROICOR relies critically on the timing of image acquisition relative to each cardiac and respiratory cycle, and because such timing can be confounded by head motion, the incorporation of timing errors introduced by volume registration into the Fourier expansion of the cardiac phase at each voxel has been proposed (Jones et al., 2008).
Non-periodic, respiratory-related contributions to the BOLD signal were first tackled by measuring the end-tidal CO2 partial pressure (PETCO2) using a capnograph, and including its time-course as a confounding regressor in the GLM (Wise et al., 2004). Alternatively, the respiratory volume per unit time (RVT) can be obtained directly from the respiratory recordings, and is commonly used as a surrogate of PETCO2 (Birn et al., 2006). Similarly, changes in heart rate (HR) have also been found to induce confounding BOLD fluctuations (Shmueli et al., 2007). In order to describe the contributions of RVT and HR to the BOLD signal in a linear systems framework, a respiratory response function (RRF, Birn et al., 2008) and a cardiac response function (CRF, Chang et al., 2009) have been estimated, respectively. These are then convolved with the RVT and HR time courses, respectively, following a similar rationale to the one underlying the use of the hemodynamic response function (HRF) to describe the BOLD response to neuronal activity (Friston et al., 1998). Alternatively, RVT and HR regressors can be shifted by an appropriate time lag in order to maximize the BOLD signal variance explained (Bianciardi et al., 2009a; Jorge et al., 2013). In order to account for the large inter- and intra-subject variability, both RRF/CRF estimation and time lag optimization can be performed specifically for each subject or even brain region (Falahpour et al., 2013; Pinto et al., 2017). Our own (Abreu et al., 2017b) and other studies (van Houdt et al., 2009) have shown that, if only cardiac data is available, it is possible to derive the corresponding respiratory data from this, and to successfully denoise BOLD-fMRI data using the resulting surrogate of the respiratory signal.
Another method of physiological noise correction is based on the extraction of confounding regressors from the fMRI data itself (Murphy et al., 2013). Because BOLD signal fluctuations of neuronal origin should be mainly located in gray matter (GM), CSF and white matter (WM) fluctuations are likely to reflect predominantly physiological noise contributions. Typically, the average BOLD time-courses within CSF and/or WM masks are computed and regressed out from the fMRI data in a GLM framework. Variations of this approach can be found in the literature in terms of the definition of the noise-related regions, as well as the features to be extracted from them. As for the former, thresholding of voxel-specific temporal standard deviation (tSTD) measures can be performed, supported by the observation of a positive correlation between the variance explained by the respective RETROICOR regressors and the tSTD value for a given voxel (Behzadi et al., 2007). A biophysically inspired measure of robust tSNR has also been proposed, whereby a mixture of Gaussians is fitted to this metric in each voxel using an expectation-maximization approach (Tierney et al., 2015). Regarding the extraction of features, PCA can be applied to the BOLD signals from noise-related regions, so that the variability of the physiological fluctuations can be taken into account (Behzadi et al., 2007; Bianciardi et al., 2009b; Jorge et al., 2013; Tierney et al., 2015). The identification of the optimal number of PCs to keep is then crucial, in order to avoid under-/over-estimation of the contribution of physiological fluctuations if an excessively low/high number of PCs is included in the model, respectively (Behzadi et al., 2007).
Furthermore, head motion parameters are also commonly regressed out from the BOLD signal. These are usually extracted from the fMRI time-series in the motion correction pre-processing step, in which each brain volume in the series is aligned to a reference volume by estimating a rigid body transformation characterized by three translation and three rotation parameters, yielding a total of six motion regressors (Murphy et al., 2013). Because large, abrupt head movements are not accurately estimated using such affine transformations, metrics identifying fMRI volumes affected by this type of motion can be computed, so that these volumes may be removed from subsequent analyses (Power et al., 2012; Tierney et al., 2015). In order to improve the temporal resolution of the motion regressors, several approaches have been proposed in the literature, mostly based on the use of navigator echoes (Thesen et al., 2000; Welch et al., 2002), camera-tracker devices (Zaitsev et al., 2016), and active markers (Ooi et al., 2011). Although they do not need additional hardware unlike tracker devices and active markers, the use of navigator echoes requires longer acquisition times (a recent review on this topic can be found in Godenschweger et al., 2016). Alternatively, in simultaneous EEG-fMRI sessions, the EEG data may be used to derive such highly sampled regressors, allowing the correction of motion at the slice level retrospectively (Zotev et al., 2012; Wong et al., 2016). Regardless of the estimation procedure, special attention must be devoted to cases where head motion is correlated with task parameters of interest, seriously confounding subsequent analyses, because plausible, but yet spurious, effects may be observed (Fellner et al., 2016).
Independent component analysis can be used for fMRI de-noising by separating sources of scanner artifact, physiological noise and brain activation (Beckmann and Smith, 2004; Brooks et al., 2008). By removing the contribution of ICs reflecting non-neuronal fluctuations from the back-reconstruction of the data, a noise-free BOLD signal is obtained. Such identification procedure is critical and it can be done manually (Griffanti et al., 2017) or by resorting to automatic classification tools, based on the temporal and spatial characteristics of the expected artifacts (De Martino et al., 2007; Tohka et al., 2008; Churchill et al., 2012b; Smith et al., 2013; Salimi-Khorshidi et al., 2014).
Data Pre-processing and Physiological Noise Correction
There is currently no consensus regarding the order by which the various steps in the processing pipeline of fMRI data should be taken to include physiological noise correction. In particular, the question whether RETROICOR should precede slice timing correction or not, and what the order between this and motion correction should be, remains open. One study found that RETROICOR should be performed prior to slice timing correction (Jones et al., 2008). Another recent study investigated the optimal processing pipeline more systematically, by testing the impact of 48 different combinations of the main processing steps on task-based activation maps (Churchill et al., 2012a). The combination of motion correction with a second order polynomial detrending yielded the highest performance on average, but a large variability across subjects was found for the optimal order of the processing steps, indicating that a subject-specific optimization should preferentially be carried out. In general, inter-subject variability and the need for subject-specific optimization of physiological noise correction have been reported in several studies (Falahpour et al., 2013; Golestani et al., 2015; Nunes et al., 2015), including our own (Abreu et al., 2017b).
In fact, in our study we showed the importance of removing not only cardiac and respiratory fluctuations, as well as those of WM, CSF and head motion origins, although the former explained BOLD signal variance to a much lesser extent than the latter. This claim is in line with the results and recommendations made by Jo et al. (2010), indicating that the removal of cardiac and respiratory fluctuations is crucial, particularly when subjects exhibit a high variability in HR and breathing depth and rate during the acquisitions. Another important question relates to the inclusion of additional terms in the physiological noise model. Typically, statistical testing is applied in a nested model approach to ascertain if the additional variance explained by a given term relative to that explained by those already included in the model is significant. Although an order of two is most commonly used for RETROICOR, with the frequent addition of RV/RVT, HR, WM, CSF and motion parameters, in the so-called extended RETROICOR approaches, it may be important to investigate the impact of specifically identifying which regressors to include in the model in each study.
EEG-Informed fMRI Integration Methods
In this section, we review the data integration methods used for mapping brain networks using an EEG-informed fMRI approach. These are largely sub-divided into univariate and multivariate methods and they are summarized in Table 4. In univariate methods, a limited number of EEG time-courses (often a single time-course) representative of the phenomena of interest is selected, and temporal or spectral features are then extracted and used to predict BOLD changes. In contrast, multivariate methods consider multiple EEG channels in this feature extraction step. The rationale underlying the choice of method for the extraction of such features mainly depends on the type of activity of interest. Epileptic activity is particularly relevant in the scope of simultaneous EEG-fMRI studies, given the suitability of this technique for the localization of brain networks associated with epileptic discharges, and therefore extensive literature on the extraction of epilepsy-related EEG features predictive of BOLD fluctuations is available (Gotman et al., 2006; Marques et al., 2009; Mulert and Lemieux, 2009; Gotman and Pittau, 2011; Laufs, 2012; Jorge et al., 2014; Murta et al., 2015; Abreu et al., 2018).
TABLE 4. List of EEG features predictive of BOLD signal fluctuations of interest. The methods used to derive such EEG features can be roughly divided into univariate (temporal, spectral and intra-cranial features) and multivariate (spatial correlation features, functional connectivity methods and others).
In univariate methods, a time-course (or a limited number of time-courses) representative of the phenomenon of interest must be obtained prior to the extraction of features predictive of the BOLD signal. Both temporal and spectral features are commonly extracted and will be specifically described here.
EEG Time-Courses Representative of Phenomena of Interest
If the phenomenon of interest of brain activity is reflected within a limited set of electrodes (e.g., posterior alpha rhythms measured by occipital electrodes), subsequent analyses can be focused on those electrodes alone (Goldman et al., 2002). However, in general, this dimensionality reduction problem is not trivial, and spatial filtering strategies must be applied to the EEG data. Assuming that fluctuations have a specific temporal structure, these can be extracted from the EEG using linear prediction algorithms that incorporate such prior knowledge about the frequency and location of the sources of interest to be isolated (Ferdowsi et al., 2015). This prior knowledge can also be incorporated in a semi blind source separation (s-BSS) technique named Functional Source Separation (FSS, Porcaro et al., 2010), biasing the estimation procedure toward sources comprising the physiological aspects of interest. If such information is not available, conventional BSS techniques must be employed.
Temporal ICA is the most commonly used spatial filtering strategy in this context (Marques et al., 2009; Formaggio et al., 2011; Leite et al., 2013). However, the selection of the ICs that best reflect the phenomenon of interest is crucial, and different selection criteria can be found in the literature. In a visual inspection approach, the IC topographies can be probed for spatial patterns resembling that expected for a given activity of interest (Debener et al., 2005). A more objective approach consists of identifying the ICs that exhibit the highest weights within a pre-specified set of electrodes expected to reflect the activity of interest (Scheeringa et al., 2008; Ke and Shen, 2010). The temporal dynamics of the IC time-courses can also be inspected for events of interest (Jann et al., 2008), and a number of quantitative methods have been proposed based on spectral criteria, particularly the power within a given frequency band (Formaggio et al., 2011). In studies where more than one EEG acquisition is performed, the reproducibility of the ICs across runs can also be used as a selection criterion (Leite et al., 2013).
A few studies have developed automatic methods for selecting ICs representing specific activities of interest based on the use of respective templates: component assessment (COMPASS) (Wessel and Ullsperger, 2011), spatiotemporal templates for independent component selection (STTICS) (Abreu et al., 2015), and PROJIC (Abreu et al., 2016b). The COMPASS and STTICS algorithms are based on the explicit similarity of the ICs with spatiotemporal templates of the activity of interest. Both methods assume that a single spatial map describes such activity, which is, however, not true if a given EEG source has a certain degree of non-stationarity. In contrast, and as briefly described above, PROJIC first projects a temporal template (an average EEG event) onto the IC space, and then clusters the resulting projections based on their power. In our experience, PROJIC has been found to accurately identify both epilepsy-related ICs as well as PA-related ICs (Abreu et al., 2016a,b). Such versatility arises from the fact that PROJIC is not based on the explicit similarity to spatiotemporal templates, which renders it potentially suitable in many other applications.
In the simplest approach for combining EEG and fMRI, the EEG signal is visually inspected for events of interest. This is essentially used in epilepsy: while inter-ictal epileptiform discharges (IEDs) are treated as zero-duration events and modeled as stick functions, ictal activity is represented by boxcar functions (both usually referred to as unitary regressors) (Lemieux et al., 2001; Bagshaw et al., 2005; Jacobs et al., 2009; LeVan and Gotman, 2009; Leal et al., 2016). A finer representation of the seizure dynamics can be achieved by dividing these epileptic events into a succession of stages (e.g., early ictal, clinical seizure onset and late ictal), and modeling each stage by a separate boxcar function (Thornton et al., 2010). In the case of IEDs, the amplitude of the stick functions can be modulated by different IED features, namely their amplitude, energy or width, which has been found to improve the correlation with BOLD signal changes (Bénar et al., 2002; LeVan et al., 2010). When recording icEEG, the amplitude, width, slope of the rising phase, energy and spatial extent of IEDs are possible predictors of epilepsy-related BOLD changes; in a recent study, only the width was found to explain additional variance to that of unitary regressors (Murta et al., 2016).
More generally, in stimuli/task-based EEG-fMRI studies, features extracted from the associated ERPs on a trial-by-trial basis can also be used to predict BOLD fluctuations, namely the trial-specific amplitude and response latency (Debener et al., 2005; Bénar et al., 2007; Fuglø et al., 2012; Nguyen and Cunnington, 2014; Wirsich et al., 2014).
In order to account for the rich temporal and spectral profiles of the EEG, more complex transfer functions between the EEG and BOLD signals have been proposed based on time-frequency decompositions of the EEG signal (Goldman et al., 2002; Moosmann et al., 2003; Laufs et al., 2006; Scheeringa et al., 2008). One of the first studies using spectral features mapped EEG alpha fluctuations by extracting the average EEG power within the alpha band across four occipital channels, at epochs of the same duration as the repetition time of the fMRI acquisition sequence (Goldman et al., 2002). A similar procedure can be applied to other EEG rhythms. By including the EEG power over multiple frequency bands as regressors in a GLM analysis of the fMRI data, their individual contributions to the BOLD signal, as well as their interactions, can be investigated (Mantini et al., 2007b; de Munck et al., 2009). Besides the EEG power across specific frequency bands, several other features of the spectrogram have also been proposed to explain the BOLD signal, namely: total power (Wan et al., 2006b), linear combination of band-specific power values (Goense and Logothetis, 2008), mean frequency (Rosa et al., 2010b), and root mean squared frequency (RMSF) (Kilner et al., 2005; Rosa et al., 2010b). Studies comparing the predictive power of these features in healthy subjects (Rosa et al., 2010b), and in an epilepsy case study (Leite et al., 2013), found that RMSF outperformed other power-weighted metrics, in agreement with the heuristic model proposed by Kilner et al. (2005). In icEEG recordings, the single-trial phase-amplitude coupling strength when performing a motor task has been found to explain BOLD variance in addition to that from a combination of EEG power across several frequency bands (Murta et al., 2017).
Multivariate methods use data from multiple EEG channels in order to capture spatial information, which cannot be assessed by univariate approaches. They include methods based on the spatial correlation of the EEG with reference spatial maps, functional connectivity measures across different EEG channels, as well as other multiway decomposition methods.
Spatial Correlation Methods
In resting-state studies of healthy volunteers, BOLD correlates of EEG microstates have been investigated, based on the hypothesis that resting-state networks are reflected in both signals. Predictors of spontaneous BOLD fluctuations occurring during rest have been obtained, by spatially correlating the concurrent EEG topographies at each time point with the previously identified EEG microstates (Britz et al., 2010; Yuan et al., 2012). Such EEG microstates can be derived by clustering (Britz et al., 2010), ICA (Yuan et al., 2012), or topographic time-frequency decomposition (Schwab et al., 2015). The latter combines knowledge from techniques for time-domain spatial analysis of EEG and time-frequency decomposition of single time-courses (Koenig et al., 2001).
In a similar approach, one study proposed to use epilepsy-specific spatial templates derived from separate, long-term EEG recordings of epileptic activity (Grouiller et al., 2011). They found that the spatial correlation between these templates and the EEG scalp topographies measured at each time point during simultaneous EEG-fMRI recordings provided a good BOLD predictor. This approach may be advantageous in cases where it is not possible to detect epileptic events on the EEG recorded simultaneously with fMRI.
In a somewhat related approach, local estimates of electrical activity can be obtained by electrical source imaging (ESI). This tool estimates the location of EEG sources in the brain responsible for generating a given topography measured at the scalp. Considering a topography representative of epileptic activity, Vulliemoz et al. (2010) determined the averaged current density within the ESI solution for the whole EEG (continuous ESI), and used it to predict the associated local BOLD changes. One study showed that this approach yielded a more accurate epileptic network mapping when compared to that described in Grouiller et al. (2011), providing concordant electro-clinical localization of the epileptic focus in all investigated patients (Elshoff et al., 2012).
Functional Connectivity Methods
Other multivariate EEG measures have been employed, particularly with the purpose of reflecting functional connectivity across the brain. In one study, the partial directed coherence (a directed measure of functional connectivity) across different frequency bands has been correlated with BOLD, in order to map the intra- and inter-hemispheric flow of information measured with EEG (Biazoli et al., 2013). Most interestingly, EEG phase synchronization measures can also be used, with the advantage that they do not depend on the amplitude of the EEG signal, in contrast with most temporal and spectral features, which renders them less susceptible to artifacts. However, such measures have been scarcely used as predictors of BOLD. In particular, Mizuhara et al. (2005) successfully mapped task-dependent BOLD signal changes using the phase synchronization index (PSI) computed for a specific frequency and channel pair. A global connectivity measure has also been employed, global field synchronization (GFS), which quantifies the overall EEG synchrony across the scalp. In a task-based EEG-fMRI study, Kottlow et al. (2012) used GFS to predict BOLD changes associated with face integration. In a resting-state study, GFS measures within the lower (8.5–10.5 Hz) and upper (10.5–12.5 Hz) alpha band were found to be positively correlated with the BOLD-fMRI-derived dorsal attention network and the default mode network, respectively (Jann et al., 2009). More recently, a study by our group applied both PSI and GFS in an EEG-fMRI study of epilepsy (Abreu et al., 2018). We showed that PSI within a frequency band of interest outperformed power-weighted metrics, as well as GFS, in predicting BOLD changes associated with epileptic activity. More importantly, we showed that PSI more specifically reflects epileptic activity, rather than motion-related spurious signal changes. We believe that the potential use of EEG synchronization measures for predicting BOLD fluctuations of interest should thus be further explored in different applications.
Other Multivariate Methods
By applying multiway decomposition methods, the multi-dimensional EEG spectrum can be blindly decomposed into patterns characterized by spatial, spectral and temporal signatures. The use of these methods allowed the mapping of thalamic substructures associated with scalp EEG signals (Schwab et al., 2015), as well as task-related BOLD signal changes more accurately than commonly used power-weighted EEG features (Marecek et al., 2016).
In a different approach, Meir-Hasson et al. (2014) tested the accuracy of different combinations of features (frequency band and time delay) extracted from each channel to predict BOLD fluctuations of interest in the visual cortex. In each iteration, the prediction error is used as feedback to guide the next combination of features to be tested; upon convergence, an optimal EEG-derived model of BOLD is obtained for each channel, and thus avoiding the spatial filtering step.
Other Integration Approaches
Although a GLM framework is by far the most commonly used for the integration of EEG and fMRI, alternative methods have also been proposed in order to overcome assumptions regarding the shape of the HRF, the linearity between EEG activity and the BOLD signal, and the probability distribution of noise in the data (e.g., Caballero-Gaudes et al., 2013). Even within the GLM framework, several approaches have been proposed to account for variable HRF shapes, and hence improve the sensitivity of detecting BOLD changes. A popular method consists of allowing variations around the canonical HRF, by adding its temporal and dispersion derivatives to the GLM, in order to account for deviations in both time-to-peak and time of onset (Friston et al., 1998). More flexible methods can be used by considering Fourier or finite basis sets (Goutte et al., 2000; Thornton et al., 2010). Ultimately, the shape of the HRF can be estimated freely as a finite impulse response for each voxel separately (Glover, 1999; Lu et al., 2006; de Munck et al., 2007; Storti et al., 2013). One study employed a set of gamma-based HRFs in which the two main parameters (time-to-peak and time of onset) were systematically varied, in order to choose the best suited to model the HRF in epilepsy patients (Grouiller et al., 2010). Subject- and pathology-specific HRFs can also be defined by estimating combinations of multiple HRFs peaking at different latencies (Bagshaw et al., 2005). Nonetheless, studies assessing the minimal degree of spatial specificity (e.g., voxel versus brain regions) that is required in order to accurately map BOLD fluctuations of interest without the risk of overfitting are still missing.
Most EEG-informed fMRI studies focus on predicting the BOLD signal measured at each voxel based on the EEG. However, new opportunities for EEG-informed fMRI have recently been created by the growing interest in the study of the temporal fluctuations of BOLD functional connectivity across the brain, the so-called dynamic functional connectivity (dFC) (Hutchison et al., 2013; Calhoun et al., 2014; Preti et al., 2016). In fact, a number of studies have already attempted to incorporate EEG data in dFC analyses, which are reviewed in Tagliazucchi and Laufs (2015). Normal brain function was first investigated by correlating EEG fluctuations of interest (mainly band-specific EEG power) with dFC fluctuations, in order to identify electrophysiological correlates of functional connectivity patterns (Scheeringa et al., 2012; Tagliazucchi et al., 2012; Chang et al., 2013; Allen et al., 2017). The relevance of using EEG to inform analyses of BOLD dFC fluctuations in epilepsy is also starting to be investigated, with promising results regarding the identification of epileptic networks that are no longer assumed to be static over time, as in standard EEG-informed fMRI studies as those described in the previous sections (Preti et al., 2014; Abreu et al., 2017a). Although the field of dFC is relatively recent and further studies are needed to better understand the physiological meaning of functional connectivity fluctuations, the inclusion of EEG in dFC analyses appears to be a promising avenue to achieve such goal, opening new lines of EEG-fMRI research.
In this review, we overviewed the several challenges associated with each step of the data analysis pipeline in EEG-informed fMRI, and provided a comprehensive description and discussion of the plethora of methods available to address each of those challenges. The motivation underlying the concurrent, multimodal acquisition of EEG and fMRI was first highlighted, including a brief description of the fundamentals of each neuroimaging modality. Special attention was then given to the critical problems concerning EEG and MRI data quality, by characterizing the artifacts induced by each modality on the other, as well as the most important modality-specific artifacts, and describing the respective artifact reduction techniques. Finally, we focused on multimodal data integration in the context of the EEG-informed fMRI approach, surveying both univariate and multivariate methods used to extract EEG features that may predict BOLD signal changes. This review may help the identification of the processing pipeline that best fits each study, in order to optimize data quality as well as the sensitivity and specificity of the brain networks obtained by EEG-informed fMRI analysis. The optimal method for the integration of data from the two modalities remains an open question, mainly because a deeper understanding about the substrates of each modality and the extent to which these substrates overlap is still needed. Furthermore, more extensive, critical and independent validation studies are needed to guide the interpretation of the findings obtained using the methods described above.
RA was responsible for reviewing the literature on simultaneous EEG-fMRI methods, and writing the review itself. AL provided crucial feedback on the more clinical aspects of the review. PF supervised the literature review and the organization of the manuscript, and participated in writing the manuscript.
We acknowledge the Portuguese Science Foundation (FCT) for financial support through Project PTDC/SAUENB/112294/2009, Project PTDC/EEIELC/3246/2012, Grant UID/EEA/50009/2013, and the Doctoral Grant PD/BD/105777/2014.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We acknowledge Prof. Fernando Lopes da Silva for his feedback and help in structuring this paper.
Abreu, R., Leal, A., and Figueiredo, P. (2017a). “Improved mapping of epileptic networks based on the correlation of BOLD-fMRI dynamic functional connectivity components with simultaneous EEG,” in 2017 25th International Society for Magnetic Resonance in Medicine (ISMRM) Congress, Concord, CA.
Abreu, R., Leal, A., Lopes da Silva, F., and Figueiredo, P. (2018). EEG synchronization measures predict epilepsy-related BOLD-fMRI fluctuations better than commonly used univariate metrics. Clin. Neurophysiol. (in press). doi: 10.1016/j.clinph.2017.12.038
Abreu, R., Leite, M., Jorge, J., Grouiller, F., van der Zwaag, W., Leal, A., et al. (2016a). Ballistocardiogram artifact correction taking into account physiological signal preservation in simultaneous EEG-fMRI. Neuroimage 135, 45–63. doi: 10.1016/j.neuroimage.2016.03.034
Abreu, R., Leite, M., Leal, A., and Figueiredo, P. (2015). “STTICS: a template-based algorithm for the objective selection of epilepsy-related EEG ICA components,” in 2015 IEEE 12th International Symposium on Biomedical Imaging (ISBI) (IEEE), Brooklyn, NY, 343–346. doi: 10.1109/ISBI.2015.7163883
Abreu, R., Leite, M., Leal, A., and Figueiredo, P. (2016b). Objective selection of epilepsy-related independent components from EEG data. J. Neurosci. Methods 258, 67–78. doi: 10.1016/j.jneumeth.2015.10.003
Abreu, R., Nunes, S., Leal, A., and Figueiredo, P. (2017b). Physiological noise correction using ECG-derived respiratory signals for enhanced mapping of spontaneous neuronal activity with simultaneous EEG-fMRI. Neuroimage 154, 115–127. doi: 10.1016/j.neuroimage.2016.08.008
Acharjee, P. P., Phlypo, R., Wu, L., Calhoun, V. D., and Adali, T. (2015). Independent vector analysis for gradient artifact removal in concurrent EEG-fMRI Data. IEEE Trans. Biomed. Eng. 62, 1750–1758. doi: 10.1109/TBME.2015.2403298
Allen, P. J., Polizzi, G., Krakow, K., Fish, D. R., and Lemieux, L. (1998). Identification of EEG events in the MR scanner: the problem of pulse artifact and a method for its subtraction. Neuroimage 8, 229–239. doi: 10.1006/nimg.1998.0361
Anami, K., Mori, T., Tanaka, F., Kawagoe, Y., Okamoto, J., Yarita, M., et al. (2003). Stepping stone sampling for retrieving artifact-free electroencephalogram during functional magnetic resonance imaging. Neuroimage 19, 281–295. doi: 10.1016/S1053-8119(03)00048-X
Angelone, L. M., Vasios, C. E., Wiggins, G., Purdon, P. L., and Bonmassar, G. (2006). On the effect of resistive EEG electrodes and leads during 7 T MRI: simulation and temperature measurement studies. Magn. Reson. Imaging 24, 801–812. doi: 10.1016/j.mri.2006.01.006
Bagshaw, A. P., Hawco, C., Bénar, C.-G., Kobayashi, E., Aghakhani, Y., Dubeau, F., et al. (2005). Analysis of the EEG-fMRI response to prolonged bursts of interictal epileptiform activity. Neuroimage 24, 1099–1112. doi: 10.1016/j.neuroimage.2004.10.010
Behzadi, Y., Restom, K., Liau, J., and Liu, T. T. (2007). A component based noise correction method (CompCor) for BOLD and perfusion based fMRI. Neuroimage 37, 90–101. doi: 10.1016/j.neuroimage.2007.04.042
Belliveau, J. W., Kennedy, D. N., McKinstry, R. C., Buchbinder, B. R., Weisskoff, R. M., Cohen, M. S., et al. (1991). Functional mapping of the human visual cortex by magnetic resonance imaging. Science 254, 716–719. doi: 10.2307/2879286
Bénar, C.-G., Aghakhani, Y., Wang, Y., Izenberg, A., Al-Asmi, A., Dubeau, F., et al. (2003). Quality of EEG in simultaneous EEG-fMRI for epilepsy. Clin. Neurophysiol. 114, 569–580. doi: 10.1016/S1388-2457(02)00383-8
Bénar, C.-G., Schön, D., Grimault, S., Nazarian, B., Burle, B., Roth, M., et al. (2007). Single-trial analysis of oddball event-related potentials in simultaneous EEG-fMRI. Hum. Brain Mapp. 28, 602–613. doi: 10.1002/hbm.20289
Bhattacharyya, P. K., and Lowe, M. J. (2004). Cardiac-induced physiologic noise in tissue is a direct observation of cardiac-induced fluctuations. Magn. Reson. Imaging 22, 9–13. doi: 10.1016/j.mri.2003.08.003
Bianciardi, M., Fukunaga, M., van Gelderen, P., Horovitz, S. G., de Zwart, J. A., Shmueli, K., et al. (2009a). Sources of functional magnetic resonance imaging signal fluctuations in the human brain at rest: a 7 T study. Magn. Reson. Imaging 27, 1019–1029. doi: 10.1016/j.mri.2009.02.004
Bianciardi, M., van Gelderen, P., Duyn, J. H., Fukunaga, M., and de Zwart, J. A. (2009b). Making the most of fMRI at 7 T by suppressing spontaneous signal fluctuations. Neuroimage 44, 448–454. doi: 10.1016/j.neuroimage.2008.08.037
Biazoli, C. E., Sturzbecher, M., White, T. P., dos Santos Onias, H. H., Andrade, K. C., de Araujo, D. B., et al. (2013). Application of partial directed coherence to the analysis of resting-state EEG-fMRI data. Brain Connect. 3, 563–568. doi: 10.1089/brain.2012.0135
Birn, R. M., Diamond, J. B., Smith, M. A., and Bandettini, P. A. (2006). Separating respiratory-variation-related fluctuations from neuronal-activity-related fluctuations in fMRI. Neuroimage 31, 1536–1548. doi: 10.1016/j.neuroimage.2006.02.048
Birn, R. M., Smith, M. A., Jones, T. B., and Bandettini, P. A. (2008). The respiration response function: the temporal dynamics of fMRI signal fluctuations related to changes in respiration. Neuroimage 40, 644–654. doi: 10.1016/j.neuroimage.2007.11.059
Biswal, B., Zerrin Yetkin, F., Haughton, V. M., and Hyde, J. S. (1995). Functional connectivity in the motor cortex of resting human brain using echo-planar mri. Magn. Reson. Med. 34, 537–541. doi: 10.1002/mrm.1910340409
Bonmassar, G., Purdon, P. L., Jääskeläinen, I. P., Chiappa, K., Solo, V., Brown, E. N., et al. (2002). Motion and Ballistocardiogram artifact removal for interleaved recording of EEG and EPs during MRI. Neuroimage 16, 1127–1141. doi: 10.1006/nimg.2002.1125
Brooks, J. C. W., Beckmann, C. F., Miller, K. L., Wise, R. G., Porro, C. A., Tracey, I., et al. (2008). Physiological noise modelling for spinal functional magnetic resonance imaging studies. Neuroimage 39, 680–692. doi: 10.1016/j.neuroimage.2007.09.018
Caballero-Gaudes, C., Van de Ville, D., Grouiller, F., Thornton, R., Lemieux, L., Seeck, M., et al. (2013). Mapping interictal epileptic discharges using mutual information between concurrent EEG and fMRI. Neuroimage 68, 248–262. doi: 10.1016/j.neuroimage.2012.12.011
Calhoun, V. D., Miller, R., Pearlson, G., and Adalı, T. (2014). The chronnectome: time-varying connectivity networks as the next frontier in fMRI data discovery. Neuron 84, 262–274. doi: 10.1016/j.neuron.2014.10.015
Campos Viola, F., Thorne, J., Edmonds, B., Schneider, T., Eichele, T., and Debener, S. (2009). Semi-automatic identification of independent components representing EEG artifact. Clin. Neurophysiol. 120, 868–877. doi: 10.1016/j.clinph.2009.01.015
Carmichael, D. W., Thornton, J. S., Rodionov, R., Thornton, R., McEvoy, A. W., Ordidge, R. J., et al. (2010). Feasibility of simultaneous intracranial EEG-fMRI in humans: a safety study. Neuroimage 49, 379–390. doi: 10.1016/j.neuroimage.2009.07.062
Carmichael, D. W., Vulliemoz, S., Rodionov, R., Thornton, J. S., McEvoy, A. W., and Lemieux, L. (2012). Simultaneous intracranial EEG-fMRI in humans: protocol considerations and data quality. Neuroimage 63, 301–309. doi: 10.1016/j.neuroimage.2012.05.056
Chaumon, M., Bishop, D. V. M., and Busch, N. A. (2015). A practical guide to the selection of independent components of the electroencephalogram for artifact correction. J. Neurosci. Methods 250, 47–63. doi: 10.1016/j.jneumeth.2015.02.025
Chowdhury, M. E. H., Mullinger, K. J., Glover, P., and Bowtell, R. (2014). Reference layer artefact subtraction (RLAS): A novel method of minimizing EEG artefacts during simultaneous fMRI. Neuroimage 84, 307–319. doi: 10.1016/j.neuroimage.2013.08.039
Churchill, N. W., Oder, A., Abdi, H., Tam, F., Lee, W., Thomas, C., et al. (2012a). Optimizing preprocessing and analysis pipelines for single-subject fMRI. I. Standard temporal motion and physiological noise correction methods. Hum. Brain Mapp. 33, 609–627. doi: 10.1002/hbm.21238
Churchill, N. W., Yourganov, G., Spring, R., Rasmussen, P. M., Lee, W., Ween, J. E., et al. (2012b). PHYCAA: Data-driven measurement and removal of physiological noise in BOLD fMRI. Neuroimage 59, 1299–1314. doi: 10.1016/j.neuroimage.2011.08.021
Cordes, D., Haughton, V. M., Arfanakis, K., Carew, J. D., Turski, P. A., Moritz, C. H., et al. (2001). Frequencies contributing to functional connectivity in the cerebral cortex in “resting-state” data. Am. J. Neuroradiol. 22, 1326–1333.
Daly, I., Scherer, R., Billinger, M., and Muller-Putz, G. (2014). FORCe: fully online and automated artifact removal for brain-computer interfacing. IEEE Trans. Neural Syst. Rehabil. Eng. 23, 725–736. doi: 10.1109/TNSRE.2014.2346621
De Martino, F., Gentile, F., Esposito, F., Balsi, M., Di Salle, F., Goebel, R., et al. (2007). Classification of fMRI independent components using IC-fingerprints and support vector machine classifiers. Neuroimage 34, 177–194. doi: 10.1016/j.neuroimage.2006.08.041
de Munck, J. C., Gonçalves, S. I., Huijboom, L., Kuijer, J. P. A., Pouwels, P. J. W., Heethaar, R. M., et al. (2007). The hemodynamic response of the alpha rhythm: An EEG/fMRI study. Neuroimage 35, 1142–1151. doi: 10.1016/j.neuroimage.2007.01.022
de Munck, J. C., Gonçalves, S. I., Mammoliti, R., Heethaar, R. M., and Lopes da Silva, F. H. (2009). Interactions between different EEG frequency bands and their effect on alpha-fMRI correlations. Neuroimage 47, 69–76. doi: 10.1016/j.neuroimage.2009.04.029
de Munck, J. C., van Houdt, P. J., Gonçalves, S. I., van Wegen, E., and Ossenblok, P. P. W. (2013). Novel artefact removal algorithms for co-registered EEG/fMRI based on selective averaging and subtraction. Neuroimage 64, 407–415. doi: 10.1016/j.neuroimage.2012.09.022
Debener, S., Mullinger, K. J., Niazy, R. K., and Bowtell, R. W. (2008). Properties of the ballistocardiogram artefact as revealed by EEG recordings at 1.5, 3 and 7 T static magnetic field strength. Int. J. Psychophysiol. 67, 189–199. doi: 10.1016/j.ijpsycho.2007.05.015
Debener, S., Strobel, A., Sorger, B., Peters, J., Kranczioch, C., Engel, A. K., et al. (2007). Improved quality of auditory event-related potentials recorded simultaneously with 3-T fMRI: removal of the ballistocardiogram artefact. Neuroimage 34, 587–597. doi: 10.1016/j.neuroimage.2006.09.031
Debener, S., Ullsperger, M., Siegel, M., Fiehler, K., von Cramon, D. Y., and Engel, A. K. (2005). Trial-by-trial coupling of concurrent electroencephalogram and functional magnetic resonance imaging identifies the dynamics of performance monitoring. J. Neurosci. 25, 11730–11737. doi: 10.1523/JNEUROSCI.3286-05.2005
Deburchgraeve, W., Cherian, P. J., De Vos, M., Swarte, R. M., Blok, J. H., Visser, G. H., et al. (2008). Automated neonatal seizure detection mimicking a human observer reading EEG. Clin. Neurophysiol. 119, 2447–2454. doi: 10.1016/j.clinph.2008.07.281
Elshoff, L., Groening, K., Grouiller, F., Wiegand, G., Wolff, S., Michel, C., et al. (2012). The value of EEG-fMRI and EEG source analysis in the presurgical setup of children with refractory focal epilepsy. Epilepsia 53, 1597–1606. doi: 10.1111/j.1528-1167.2012.03587.x
Ericsson, A., Weis, J., Hemmingsson, A., Wikström, M., and Sperber, G. O. (1995). Measurements of magnetic field variations in the human brain using a 3D-FT multiple gradient echo technique. Magn. Reson. Med. 33, 171–177. doi: 10.1002/mrm.1910330205
Falahpour, M., Refai, H., and Bodurka, J. (2013). Subject specific BOLD fMRI respiratory and cardiac response functions obtained from global signal. Neuroimage 72, 252–264. doi: 10.1016/j.neuroimage.2013.01.050
Fellner, M.-C., Volberg, G., Mullinger, K. J., Goldhacker, M., Wimber, M., Greenlee, M. W., et al. (2016). Spurious correlations in simultaneous EEG-fMRI driven by in-scanner movement. Neuroimage 133, 354–366. doi: 10.1016/j.neuroimage.2016.03.031
Formaggio, E., Storti, S. F., Bertoldo, A., Manganotti, P., Fiaschi, A., and Toffolo, G. M. (2011). Integrating EEG and fMRI in epilepsy. Neuroimage 54, 2719–2731. doi: 10.1016/j.neuroimage.2010.11.038
Freyer, F., Becker, R., Anami, K., Curio, G., Villringer, A., and Ritter, P. (2009). Ultrahigh-frequency EEG during fMRI: Pushing the limits of imaging-artifact correction. Neuroimage 48, 94–108. doi: 10.1016/j.neuroimage.2009.06.022
Fuglø, D., Pedersen, H., Rostrup, E., Hansen, A. E., and Larsson, H. B. W. (2012). Correlation between single-trial visual evoked potentials and the blood oxygenation level dependent response in simultaneously recorded electroencephalography-functional magnetic resonance imaging. Magn. Reson. Med. 68, 252–260. doi: 10.1002/mrm.23227
Glover, G. H., Li, T. Q., and Ress, D. (2000). Image-based method for retrospective correction of physiological motion effects in fMRI: RETROICOR. Magn. Reson. Med. 44, 162–167. doi: 10.1002/1522-2594(200007)44:1<162::AID-MRM23>3.0.CO;2-E
Golestani, A. M., Chang, C., Kwinta, J. B., Khatamian, Y. B., and Jean Chen, J. (2015). Mapping the end-tidal CO2 response function in the resting-state BOLD fMRI signal: spatial specificity, test-retest reliability and effect of fMRI sampling rate. Neuroimage 104, 266–277. doi: 10.1016/j.neuroimage.2014.10.031
Gonçalves, S. I., Pouwels, P. J. W., Kuijer, J. P. A., Heethaar, R. M., and de Munck, J. C. (2007). Artifact removal in co-registered EEG/fMRI by selective average subtraction. Clin. Neurophysiol. 118, 2437–2450. doi: 10.1016/j.clinph.2007.08.017
Goncharova, I., McFarland, D., Vaughan, T., and Wolpaw, J. (2003). EMG contamination of EEG: spectral and topographical characteristics. Clin. Neurophysiol. 114, 1580–1593. doi: 10.1016/S1388-2457(03)00093-2
Gotman, J., Kobayashi, E., Bagshaw, A. P., Bénar, C. G., and Dubeau, F. (2006). Combining EEG and fMRI: a multimodal tool for epilepsy research. J. Magn. Reson. Imaging 23, 906–920. doi: 10.1002/jmri.20577
Greitz, D., Franck, A., and Nordell, B. (2010). On the pulsatile nature of intracranial and spinal CSF-circulation demonstrated by MR imaging. Acta Radiol. 34, 321–328. doi: 10.1177/028418519303400403
Griffanti, L., Douaud, G., Bijsterbosch, J., Evangelisti, S., Alfaro-Almagro, F., Glasser, M. F., et al. (2017). Hand classification of fMRI ICA noise components. Neuroimage 154, 188–205. doi: 10.1016/j.neuroimage.2016.12.036
Grouiller, F., Thornton, R. C., Groening, K., Spinelli, L., Duncan, J. S., Schaller, K., et al. (2011). With or without spikes: localization of focal epileptic activity by simultaneous electroencephalography and functional magnetic resonance imaging. Brain 134, 2867–2886. doi: 10.1093/brain/awr156
Grouiller, F., Vercueil, L., Krainik, A., Segebarth, C., Kahane, P., and David, O. (2007). A comparative study of different artefact removal algorithms for EEG signals acquired during functional MRI. Neuroimage 38, 124–137. doi: 10.1016/j.neuroimage.2007.07.025
Grouiller, F., Vercueil, L., Krainik, A., Segebarth, C., Kahane, P., and David, O. (2010). Characterization of the hemodynamic modes associated with interictal epileptic activity using a deformable model-based analysis of combined EEG and functional MRI recordings. Hum. Brain Mapp. 31, 1157–1173. doi: 10.1002/hbm.20925
Harvey, A. K., Pattinson, K. T. S., Brooks, J. C. W., Mayhew, S. D., Jenkinson, M., and Wise, R. G. (2008). Brainstem functional magnetic resonance imaging: disentangling signal from physiological noise. J. Magn. Reson. Imaging 28, 1337–1344. doi: 10.1002/jmri.21623
Hoffmann, A., Jäger, L., Werhahn, K. J., Jaschke, M., Noachtar, S., and Reiser, M. (2000). Electroencephalography during functional echo-planar imaging: detection of epileptic spikes using post-processing methods. Magn. Reson. Med. 44, 791–798. doi: 10.1002/1522-2594(200011)44:5<791::AID-MRM17>3.0.CO;2-2
Hu, X., Le, T. H., Parrish, T., and Erhard, P. (1995). Retrospective estimation and correction of physiological fluctuation in functional MRI. Magn. Reson. Med. 34, 201–212. doi: 10.1002/mrm.1910340211
Hutchison, R. M., Womelsdorf, T., Allen, E. A., Bandettini, P. A., Calhoun, V. D., Corbetta, M., et al. (2013). Dynamic functional connectivity: promise, issues, and interpretations. Neuroimage 80, 360–378. doi: 10.1016/j.neuroimage.2013.05.079
Hutton, C., Josephs, O., Stadler, J., Featherstone, E., Reid, A., Speck, O., et al. (2011). The impact of physiological noise correction on fMRI at 7T. Neuroimage 57, 101–112. doi: 10.1016/j.neuroimage.2011.04.018
Ives, J. R., Warach, S., Schmitt, F., Edelman, R. R., and Schomer, D. L. (1993). Monitoring the patient’s EEG during echo planar MRI. Electroencephalogr. Clin. Neurophysiol. 87, 417–420. doi: 10.1016/0013-4694(93)90156-P
Jacobs, J., Levan, P., Moeller, F., Boor, R., Stephani, U., Gotman, J., et al. (2009). Hemodynamic changes preceding the interictal EEG spike in patients with focal epilepsy investigated using simultaneous EEG-fMRI. Neuroimage 45, 1220–1231. doi: 10.1016/j.neuroimage.2009.01.014
Jann, K., Dierks, T., Boesch, C., Kottlow, M., Strik, W., and Koenig, T. (2009). BOLD correlates of EEG alpha phase-locking and the fMRI default mode network. Neuroimage 45, 903–916. doi: 10.1016/j.neuroimage.2009.01.001
Jann, K., Wiest, R., Hauf, M., Meyer, K., Boesch, C., Mathis, J., et al. (2008). BOLD correlates of continuously fluctuating epileptic activity isolated by independent component analysis. Neuroimage 42, 635–648. doi: 10.1016/j.neuroimage.2008.05.001
Jezzard, P., and Toosy, A. (2005). “Functional MRI,” in MR Imaging in White Matter Diseases of the Brain and Spinal Cord Medical Radiology Diagnostic Imaging, eds M. Filippi, N. Stefano, V. Dousset, and J. C. McGowan (Berlin: Springer-Verlag), 93–110. doi: 10.1007/3-540-27644-0
Jo, H. J., Saad, Z. S., Simmons, W. K., Milbury, L. A., and Cox, R. W. (2010). Mapping sources of correlation in resting state FMRI, with artifact detection and removal. Neuroimage 52, 571–582. doi: 10.1016/j.neuroimage.2010.04.246
Jorge, J., Figueiredo, P., van der Zwaag, W., and Marques, J. P. (2013). Signal fluctuations in fMRI data acquired with 2D-EPI and 3D-EPI at 7 Tesla. Magn. Reson. Imaging 31, 212–220. doi: 10.1016/j.mri.2012.07.001
Jorge, J., Grouiller, F., Gruetter, R., van der Zwaag, W., and Figueiredo, P. (2015a). Towards high-quality simultaneous EEG-fMRI at 7T: Detection and reduction of EEG artifacts due to head motion. Neuroimage 120, 143–153. doi: 10.1016/j.neuroimage.2015.07.020
Jorge, J., Grouiller, F., Ipek,Ö., Stoermer, R., Michel, C. M., Figueiredo, P., et al. (2015b). Simultaneous EEG-fMRI at ultra-high field: artifact prevention and safety assessment. Neuroimage 105, 132–144. doi: 10.1016/j.neuroimage.2014.10.055
Ke, L., and Shen, J. (2010). “Classification of EEG signals by ICA and OVR-CSP,” in Proceedings of the 3rd International Congress Image and Signal Processing CISP, Yantai, 2980–2984. doi: 10.1109/CISP.2010.5647534
Kim, H.-C., Yoo, S.-S., and Lee, J.-H. (2015). Recursive approach of EEG-segment-based principal component analysis substantially reduces cryogenic pump artifacts in simultaneous EEG–fMRI data. Neuroimage 104, 437–451. doi: 10.1016/j.neuroimage.2014.09.049
Klose, U., Strik, C., Kiefer, C., and Grodd, W. (2000). Detection of a relation between respiration and CSF pulsation with an echoplanar technique. J. Magn. Reson. Imaging 11, 438–444. doi: 10.1002/(SICI)1522-2586(200004)11:4<438::AID-JMRI12>3.0.CO;2-O
Koskinen, M., and Vartiainen, N. (2009). Removal of imaging artifacts in EEG during simultaneous EEG/fMRI recording: reconstruction of a high-precision artifact template. Neuroimage 46, 160–167. doi: 10.1016/j.neuroimage.2009.01.061
Kottlow, M., Jann, K., Dierks, T., and Koenig, T. (2012). Increased phase synchronization during continuous face integration measured simultaneously with EEG and fMRI. Clin. Neurophysiol. 123, 1536–1548. doi: 10.1016/j.clinph.2011.12.019
Krakow, K., Allen, P. J., Symms, M. R., Lemieux, L., Josephs, O., and Fish, D. R. (2000). EEG recording during fMRI experiments: image quality. Hum. Brain Mapp. 10, 10–15. doi: 10.1002/(SICI)1097-0193(200005)10:1<10::AID-HBM20>3.0.CO;2-T
Krishnaswamy, P., Bonmassar, G., Poulsen, C., Pierce, E. T., Purdon, P. L., and Brown, E. N. (2016). Reference-free removal of EEG-fMRI ballistocardiogram artifacts with harmonic regression. Neuroimage 128, 398–412. doi: 10.1016/j.neuroimage.2015.06.088
Kwong, K. K., Belliveau, J. W., Chesler, D. A., Goldberg, I. E., Weisskoff, R. M., Poncelet, B. P., et al. (1992). Dynamic magnetic resonance imaging of human brain activity during primary sensory stimulation. Proc. Natl. Acad. Sci. U.S.A. 89, 5675–5679. doi: 10.1073/pnas.89.12.5675
Laufs, H., Holt, J. L., Elfont, R., Krams, M., Paul, J. S., Krakow, K., et al. (2006). Where the BOLD signal goes when alpha EEG leaves. Neuroimage 31, 1408–1418. doi: 10.1016/j.neuroimage.2006.02.002
Leal, A., Vieira, J. P., Lopes, R., Nunes, R. G., Gonçalves, S. I., Lopes da Silva, F., et al. (2016). Dynamics of epileptic activity in a peculiar case of childhood absence epilepsy and correlation with thalamic levels of GABA. Epilepsy Behav. Case Rep. 5, 57–65. doi: 10.1016/j.ebcr.2016.03.004
Lee, T.-W., Girolami, M., and Sejnowski, T. J. (1999). Independent component analysis using an extended infomax algorithm for mixed subgaussian and supergaussian sources. Neural Comput. 11, 417–441. doi: 10.1162/089976699300016719
Lei, X., Valdes-Sosa, P. A., and Yao, D. (2012). EEG/fMRI fusion based on independent component analysis: integration of data-driven and model-driven methods. J. Integr. Neurosci. 11, 313–337. doi: 10.1142/S0219635212500203
Lemieux, L., Salek-Haddadi, A., Josephs, O., Allen, P., Toms, N., Scott, C., et al. (2001). Event-related fMRI with simultaneous and continuous EEG: description of the method and initial case report. Neuroimage 14, 780–787. doi: 10.1006/nimg.2001.0853
LeVan, P., and Gotman, J. (2009). Independent component analysis as a model-free approach for the detection of BOLD changes related to epileptic spikes: a simulation study. Hum. Brain Mapp. 30, 2021–2031. doi: 10.1002/hbm.20647
LeVan, P., Maclaren, J., Herbst, M., Sostheim, R., Zaitsev, M., and Hennig, J. (2013). Ballistocardiographic artifact removal from simultaneous EEG-fMRI using an optical motion-tracking system. Neuroimage 75, 1–11. doi: 10.1016/j.neuroimage.2013.02.039
Liston, A., Lund, T. E., Salek-Haddadi, A., Hamandi, K., Friston, K. J., and Lemieux, L. (2006). Modelling cardiac signal as a confound in EEG-fMRI and its application in focal epilepsy studies. Neuroimage 30, 827–834. doi: 10.1016/j.neuroimage.2005.10.025
Liu, Z., de Zwart, J. A., van Gelderen, P., Kuo, L. W., and Duyn, J. H. (2012). Statistical feature extraction for artifact removal from concurrent fMRI-EEG recordings. Neuroimage 59, 2073–2087. doi: 10.1016/j.neuroimage.2011.10.042
Lu, Y., Bagshaw, A. P., Grova, C., Kobayashi, E., Dubeau, F., and Gotman, J. (2006). Using voxel-specific hemodynamic response function in EEG-fMRI data analysis. Neuroimage 32, 238–247. doi: 10.1016/j.neuroimage.2005.11.040
Mantini, D., Perrucci, M. G., Cugini, S., Ferretti, A., Romani, G. L., and Del Gratta, C. (2007a). Complete artifact removal for EEG recorded during continuous fMRI using independent component analysis. Neuroimage 34, 598–607. doi: 10.1016/j.neuroimage.2006.09.037
Mantini, D., Perrucci, M. G., Del Gratta, C., Romani, G. L., and Corbetta, M. (2007b). Electrophysiological signatures of resting state networks in the human brain. Proc. Natl. Acad. Sci. U.S.A. 104, 13170–13175. doi: 10.1073/pnas.0700668104
Marecek, R., Lamos, M., Mikl, M., Barton, M., Fajkus, J., Rektor, R., et al. (2016). What can be found in scalp EEG spectrum beyond common frequency bands. EEG–fMRI study. J. Neural Eng. 13:46026. doi: 10.1088/1741-2560/13/4/046026
Marques, J. P., Rebola, J., Figueiredo, P., Pinto, A., Sales, F., and Castelo-Branco, M. (2009). ICA decomposition of EEG signal for fMRI processing in epilepsy. Hum. Brain Mapp. 30, 2986–2996. doi: 10.1002/hbm.20723
Masterton, R. A. J., Abbott, D. F., Fleming, S. W., and Jackson, G. D. (2007). Measurement and reduction of motion and ballistocardiogram artefacts from simultaneous EEG and fMRI recordings. Neuroimage 37, 202–211. doi: 10.1016/j.neuroimage.2007.02.060
Maziero, D., Velasco, T. R., Hunt, N., Payne, E., Lemieux, L., Salmon, C. E. G., et al. (2016). Towards motion insensitive EEG-fMRI: correcting motion-induced voltages and gradient artefact instability in EEG using an fMRI prospective motion correction (PMC) system. Neuroimage 138, 13–27. doi: 10.1016/j.neuroimage.2016.05.003
McMenamin, B. W., Shackman, A. J., Maxwell, J. S., Bachhuber, D. R. W., Koppenhaver, A. M., Greischar, L. L., et al. (2010). Validation of ICA-based myogenic artifact correction for scalp and source-localized EEG. Neuroimage 49, 2416–2432. doi: 10.1016/j.neuroimage.2009.10.010
Meir-Hasson, Y., Kinreich, S., Podlipsky, I., Hendler, T., and Intrator, N. (2014). An EEG Finger-Print of fMRI deep regional activation. Neuroimage 102, 128–141. doi: 10.1016/j.neuroimage.2013.11.004
Mizuhara, H., Wang, L. Q., Kobayashi, K., and Yamaguchi, Y. (2005). Long-range EEG phase synchronization during an arithmetic task indexes a coherent cortical network simultaneously measured by fMRI. Neuroimage 27, 553–563. doi: 10.1016/j.neuroimage.2005.04.030
Mognon, A., Jovicich, J., Bruzzone, L., and Buiatti, M. (2011). ADJUST: an automatic EEG artifact detector based on the joint use of spatial and temporal features. Psychophysiology 48, 229–240. doi: 10.1111/j.1469-8986.2010.01061.x
Moosmann, M., Ritter, P., Krastel, I., Brink, A., Thees, S., Blankenburg, F., et al. (2003). Correlates of alpha rhythm in functional magnetic resonance imaging and near infrared spectroscopy. Neuroimage 20, 145–158. doi: 10.1016/S1053-8119(03)00344-6
Moosmann, M., Schönfelder, V. H., Specht, K., Scheeringa, R., Nordby, H., and Hugdahl, K. (2009). Realignment parameter-informed artefact correction for simultaneous EEG-fMRI recordings. Neuroimage 45, 1144–1150. doi: 10.1016/j.neuroimage.2009.01.024
Mullinger, K., Brookes, M., Stevenson, C., Morgan, P., and Bowtell, R. (2008a). Exploring the feasibility of simultaneous electroencephalography/functional magnetic resonance imaging at 7 T. Magn. Reson. Imaging 26, 968–977. doi: 10.1016/j.mri.2008.02.014
Mullinger, K., Debener, S., Coxon, R., and Bowtell, R. (2008b). Effects of simultaneous EEG recording on MRI data quality at 1.5, 3 and 7 tesla. Int. J. Psychophysiol. 67, 178–188. doi: 10.1016/j.ijpsycho.2007.06.008
Mullinger, K. J., Havenhand, J., and Bowtell, R. (2013). Identifying the sources of the pulse artefact in EEG recordings made inside an MR scanner. Neuroimage 71, 75–83. doi: 10.1016/j.neuroimage.2012.12.070
Mullinger, K. J., Yan, W. X., and Bowtell, R. (2011). Reducing the gradient artefact in simultaneous EEG-fMRI by adjusting the subject’s axial position. Neuroimage 54, 1942–1950. doi: 10.1016/j.neuroimage.2010.09.079
Murta, T., Chaudhary, U. J., Tierney, T. M., Dias, A., Leite, M., Carmichael, D. W., et al. (2017). Phase–amplitude coupling and the BOLD signal: a simultaneous intracranial EEG (icEEG) - fMRI study in humans performing a finger-tapping task. Neuroimage 146, 438–451. doi: 10.1016/j.neuroimage.2016.08.036
Murta, T., Hu, L., Tierney, T. M., Chaudhary, U. J., Walker, M. C., Carmichael, D. W., et al. (2016). A study of the electro-haemodynamic coupling using simultaneously acquired intracranial EEG and fMRI data in humans. Neuroimage 142, 371–380. doi: 10.1016/j.neuroimage.2016.08.001
Murta, T., Leal, A., Garrido, M. I., and Figueiredo, P. (2012). Dynamic causal modelling of epileptic seizure propagation pathways: a combined EEG–fMRI study. Neuroimage 62, 1634–1642. doi: 10.1016/j.neuroimage.2012.05.053
Murta, T., Leite, M., Carmichael, D. W., Figueiredo, P., and Lemieux, L. (2015). Electrophysiological correlates of the BOLD signal for EEG-informed fMRI. Hum. Brain Mapp. 36, 391–414. doi: 10.1002/hbm.22623
Negishi, M., Abildgaard, M., Nixon, T., and Constable, R. T. (2004). Removal of time-varying gradient artifacts from EEG data acquired during continuous fMRI. Clin. Neurophysiol. 115, 2181–2192. doi: 10.1016/j.clinph.2004.04.005
Neuner, I., Arrubla, J., Felder, J., and Shah, N. J. (2014). Simultaneous EEG-fMRI acquisition at low, high and ultra-high magnetic fields up to 9.4 T: perspectives and challenges. Neuroimage 102(Pt 1), 71–79. doi: 10.1016/j.neuroimage.2013.06.048
Neuner, I., Warbrick, T., Arrubla, J., Felder, J., Celik, A., Reske, M., et al. (2013). EEG acquisition in ultra-high static magnetic fields up to 9.4 T. Neuroimage 68, 214–220. doi: 10.1016/j.neuroimage.2012.11.064
Nguyen, V. T., and Cunnington, R. (2014). The superior temporal sulcus and the N170 during face processing: single trial analysis of concurrent EEG–fMRI. Neuroimage 86, 492–502. doi: 10.1016/j.neuroimage.2013.10.047
Niazy, R. K., Beckmann, C. F., Iannetti, G. D., Brady, J. M., and Smith, S. M. (2005). Removal of FMRI environment artifacts from EEG data using optimal basis sets. Neuroimage 28, 720–737. doi: 10.1016/j.neuroimage.2005.06.067
Nierhaus, T., Gundlach, C., Goltz, D., Thiel, S. D., Pleger, B., and Villringer, A. (2013). Internal ventilation system of MR scanners induces specific EEG artifact during simultaneous EEG-fMRI. Neuroimage 74, 70–76. doi: 10.1016/j.neuroimage.2013.02.016
Nöth, U., Laufs, H., Stoermer, R., and Deichmann, R. (2012). Simultaneous electroencephalography-functional MRI at 3 T: an analysis of safety risks imposed by performing anatomical reference scans with the EEG equipment in place. J. Magn. Reson. Imaging 35, 561–571. doi: 10.1002/jmri.22843
Nunes, S., Bianciardi, M., Dias, A., Abreu, R., Rodrigues, J., Silveira, L. M., et al. (2015). “Subject-specific modeling of physiological noise in resting-state fMRI at 7T,” in Proceedings of the 23th International Society for Magnetic Resonance in Medicine (ISMRM) Congress, Toronto, ON.
Nunez, P. L., and Silberstein, R. B. (2000). On the relationship of synaptic activity to macroscopic measurements: does co-registration of EEG with fMRI Make Sense? Brain Topogr. 13, 79–96. doi: 10.1023/A:1026683200895
Ogawa, S., Lee, T.-M., Nayak, A. S., and Glynn, P. (1990). Oxygenation-sensitive contrast in magnetic resonance image of rodent brain at high magnetic fields. Magn. Reson. Med. 14, 68–78. doi: 10.1002/mrm.1910140108
Ooi, M. B., Krueger, S., Muraskin, J., Thomas, W. J., and Brown, T. R. (2011). Echo-planar imaging with prospective slice-by-slice motion correction using active markers. Magn. Reson. Med. 66, 73–81. doi: 10.1002/mrm.22780
Ostwald, D., Porcaro, C., and Bagshaw, A. P. (2010). An information theoretic approach to EEG-fMRI integration of visually evoked responses. Neuroimage 49, 498–516. doi: 10.1016/j.neuroimage.2009.07.038
Ostwald, D., Porcaro, C., Mayhew, S. D., and Bagshaw, A. P. (2012). EEG-fMRI based information theoretic characterization of the human perceptual decision system. PLOS ONE 7:e33896. doi: 10.1371/journal.pone.0033896
Panzeri, S., Magri, C., and Logothetis, N. K. (2008). On the use of information theory for the analysis of the relationship between neural and imaging signals. Magn. Reson. Imaging 26, 1015–1025. doi: 10.1016/j.mri.2008.02.019
Pinto, J., Nunes, S., Bianciardi, M., Dias, A., Silveira, L. M., Wald, L. L., et al. (2017). Improved 7 Tesla resting-state fMRI connectivity measurements by cluster-based modeling of respiratory volume and heart rate effects. Neuroimage 153, 262–272. doi: 10.1016/j.neuroimage.2017.04.009
Porcaro, C., Ostwald, D., and Bagshaw, A. P. (2010). Functional source separation improves the quality of single trial visual evoked potentials recorded during concurrent EEG-fMRI. Neuroimage 50, 112–123. doi: 10.1016/j.neuroimage.2009.12.002
Poulsen, C., Wakeman, D. G., Atefi, S. R., Luu, P., Konyn, A., and Bonmassar, G. (2017). Polymer thick film technology for improved simultaneous dEEG/MRI recording: Safety and MRI data quality. Magn. Reson. Med. 77, 895–903. doi: 10.1002/mrm.26116
Power, J. D., Barnes, K. A., Snyder, A. Z., Schlaggar, B. L., and Petersen, S. E. (2012). Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. Neuroimage 59, 2142–2154. doi: 10.1016/j.neuroimage.2011.10.018
Preti, M. G., Leonardi, N., Karahanoglu, F. I., Grouiller, F., Genetti, M., Seeck, M., et al. (2014). “Epileptic network activity revealed by dynamic functional connectivity in simultaneous EEG-fMRI,” in 2014 IEEE 11th International Symposium on Biomedical Imaging (ISBI) (IEEE), Beijing, 9–12. doi: 10.1109/ISBI.2014.6867796
Purdon, P. L., and Weisskoff, R. M. (1998). Effect of temporal autocorrelation due to physiological noise and stimulus paradigm on voxel-level false-positive rates in fMRI. Hum. Brain Mapp. 6, 239–249. doi: 10.1002/(SICI)1097-0193(1998)6:4<239::AID-HBM4>3.0.CO;2-4
Raj, D., Anderson, A. W., and Gore, J. C. (2001). Respiratory effects in human functional magnetic resonance imaging due to bulk susceptibility changes. Phys. Med. Biol. 46, 3331–3340. doi: 10.1088/0031-9155/46/12/318
Ritter, P., Becker, R., Graefe, C., and Villringer, A. (2007). Evaluating gradient artifact correction of EEG data acquired simultaneously with fMRI. Magn. Reson. Imaging 25, 923–932. doi: 10.1016/j.mri.2007.03.005
Rosa, M. J., Daunizeau, J., and Friston, K. J. (2010a). EEG-fMRI integration: a critical review of biophysical modeling and data analysis approaches. J. Integr. Neurosci. 9, 453–476. doi: 10.1142/S0219635210002512
Rosa, M. J., Kilner, J., Blankenburg, F., Josephs, O., and Penny, W. (2010b). Estimating the transfer function from neuronal activity to BOLD using simultaneous EEG-fMRI. Neuroimage 49, 1496–1509. doi: 10.1016/j.neuroimage.2009.09.011
Rothlübbers, S., Relvas, V., Leal, A., Murta, T., Lemieux, L., and Figueiredo, P. (2014). Characterisation and reduction of the EEG artefact caused by the helium cooling pump in the MR environment: validation in epilepsy patient data. Brain Topogr. 28, 208–220. doi: 10.1007/s10548-014-0408-0
Ryali, S., Glover, G. H., Chang, C., and Menon, V. (2009). Development, validation, and comparison of ICA-based gradient artifact reduction algorithms for simultaneous EEG-spiral in/out and echo-planar fMRI recordings. Neuroimage 48, 348–361. doi: 10.1016/j.neuroimage.2009.06.072
Salimi-Khorshidi, G., Douaud, G., Beckmann, C. F., Glasser, M. F., Griffanti, L., and Smith, S. M. (2014). Automatic denoising of functional MRI data: combining independent component analysis and hierarchical fusion of classifiers. Neuroimage 90, 449–468. doi: 10.1016/j.neuroimage.2013.11.046
Scheeringa, R., Bastiaansen, M. C. M., Petersson, K. M., Oostenveld, R., Norris, D. G., and Hagoort, P. (2008). Frontal theta EEG activity correlates negatively with the default mode network in resting state. Int. J. Psychophysiol. 67, 242–251. doi: 10.1016/j.ijpsycho.2007.05.017
Scheeringa, R., Petersson, K. M., Kleinschmidt, A., Jensen, O., and Bastiaansen, M. C. M. (2012). EEG alpha power modulation of fMRI resting-state connectivity. Brain Connect. 2, 254–264. doi: 10.1089/brain.2012.0088
Schwab, S., Koenig, T., Morishima, Y., Dierks, T., Federspiel, A., and Jann, K. (2015). Discovering frequency sensitive thalamic nuclei from EEG microstate informed resting state fMRI. Neuroimage 118, 368–375. doi: 10.1016/j.neuroimage.2015.06.001
Shmueli, K., van Gelderen, P., de Zwart, J. A., Horovitz, S. G., Fukunaga, M., Jansma, J. M., et al. (2007). Low-frequency fluctuations in the cardiac rate as a source of variance in the resting-state fMRI BOLD signal. Neuroimage 38, 306–320. doi: 10.1016/j.neuroimage.2007.07.037
Smith, S. M., Beckmann, C. F., Andersson, J., Auerbach, E. J., Bijsterbosch, J., Douaud, G., et al. (2013). Resting-state fMRI in the human connectome project. Neuroimage 80, 144–168. doi: 10.1016/j.neuroimage.2013.05.039
Srivastava, G., Crottaz-Herbette, S., Lau, K. M., Glover, G. H., and Menon, V. (2005). ICA-based procedures for removing ballistocardiogram artifacts from EEG data acquired in the MRI scanner. Neuroimage 24, 50–60. doi: 10.1016/j.neuroimage.2004.09.041
Storti, S. F., Formaggio, E., Bertoldo, A., Manganotti, P., Fiaschi, A., and Toffolo, G. M. (2013). Modelling hemodynamic response function in epilepsy. Clin. Neurophysiol. 124, 2108–2118. doi: 10.1016/j.clinph.2013.05.024
Tagliazucchi, E., von Wegner, F., Morzelewski, A., Brodbeck, V., and Laufs, H. (2012). Dynamic BOLD functional connectivity in humans and its electrophysiological correlates. Front. Hum. Neurosci. 6:339. doi: 10.3389/fnhum.2012.00339
Tenforde, T. S., Gaffey, C. T., Moyer, B. R., and Budlnger, T. F. (1983). Cardiovascular alterations in Macaca monkeys exposed to stationary magnetic fields: experimental observations and theoretical analysis. Bioelectromagnetics 4, 1–9. doi: 10.1002/bem.2250040102
Thesen, S., Heid, O., Mueller, E., and Schad, L. R. (2000). Prospective acquisition correction for head motion with image-based tracking for real-time fMRI. Magn. Reson. Med. 44, 457–465. doi: 10.1002/1522-2594(200009)44:3<457::AID-MRM17>3.0.CO;2-R
Thornton, R. C., Rodionov, R., Laufs, H., Vulliemoz, S., Vaudano, A., Carmichael, D., et al. (2010). Imaging haemodynamic changes related to seizures: comparison of EEG-based general linear model, independent component analysis of fMRI and intracranial EEG. Neuroimage 53, 196–205. doi: 10.1016/j.neuroimage.2010.05.064
Tierney, T. M., Croft, L. J., Centeno, M., Shamshiri, E. A., Perani, S., Baldeweg, T., et al. (2015). FIACH: A biophysical model for automatic retrospective noise control in fMRI. Neuroimage 124, 1009–1020. doi: 10.1016/j.neuroimage.2015.09.034
Tohka, J., Foerde, K., Aron, A. R., Tom, S. M., Toga, A. W., and Poldrack, R. A. (2008). Automatic independent component labeling for artifact removal in fMRI. Neuroimage 39, 1227–1245. doi: 10.1016/j.neuroimage.2007.10.013
Valdes-Sosa, P. A., Sanchez-Bornot, J. M., Sotero, R. C., Iturria-Medina, Y., Aleman-Gomez, Y., Bosch-Bayard, J., et al. (2009). Model driven EEG/fMRI fusion of brain oscillations. Hum. Brain Mapp. 30, 2701–2721. doi: 10.1002/hbm.20704
van der Zwaag, W., Schäfer, A., Marques, J. P., Turner, R., and Trampel, R. (2016). Recent applications of UHF-MRI in the study of human brain function and structure: a review. NMR Biomed. 29, 1274–1288. doi: 10.1002/nbm.3275
van Houdt, P. J., Ossenblok, P. P. W., Boon, P. A. J. M., Leijten, F. S. S., Velis, D. N., Stam, C. J., et al. (2009). Correction for pulse height variability reduces physiological noise in functional MRI when studying spontaneous brain activity. Hum. Brain Mapp. 31, 311–325. doi: 10.1002/hbm.20866
Vanderperren, K., De Vos, M., Ramautar, J. R., Novitskiy, N., Mennes, M., Assecondi, S., et al. (2010). Removal of BCG artifacts from EEG recordings inside the MR scanner: a comparison of methodological and validation-related aspects. Neuroimage 50, 920–934. doi: 10.1016/j.neuroimage.2010.01.010
Vasios, C. E., Angelone, L. M., Purdon, P. L., Ahveninen, J., Belliveau, J. W., and Bonmassar, G. (2006). EEG/(f)MRI measurements at 7 Tesla using a new EEG cap (“InkCap”). Neuroimage 33, 1082–1092. doi: 10.1016/j.neuroimage.2006.07.038
Vulliemoz, S., Carmichael, D. W., Rosenkranz, K., Diehl, B., Rodionov, R., Walker, M. C., et al. (2011). Simultaneous intracranial EEG and fMRI of interictal epileptic discharges in humans. Neuroimage 54, 182–190. doi: 10.1016/j.neuroimage.2010.08.004
Vulliemoz, S., Rodionov, R., Carmichael, D. W., Thornton, R., Guye, M., Lhatoo, S. D., et al. (2010). Continuous EEG source imaging enhances analysis of EEG-fMRI in focal epilepsy. Neuroimage 49, 3219–3229. doi: 10.1016/j.neuroimage.2009.11.055
Wan, X., Iwata, K., Riera, J., Kitamura, M., and Kawashima, R. (2006a). Artifact reduction for simultaneous EEG/fMRI recording: adaptive FIR reduction of imaging artifacts. Clin. Neurophysiol. 117, 681–692. doi: 10.1016/j.clinph.2005.07.025
Wan, X., Riera, J., Iwata, K., Takahashi, M., Wakabayashi, T., and Kawashima, R. (2006b). The neural basis of the hemodynamic response nonlinearity in human primary visual cortex: implications for neurovascular coupling mechanism. Neuroimage 32, 616–625. doi: 10.1016/j.neuroimage.2006.03.040
Welch, E. B., Manduca, A., Grimm, R. C., Ward, H. A., and Jack, C. R. Jr. (2002). Spherical navigator echoes for full 3D rigid body motion measurement in MRI. Magn. Reson. Med. 47, 32–41. doi: 10.1002/mrm.10012
Wessel, J. R., and Ullsperger, M. (2011). Selection of independent components representing event-related brain potentials: a data-driven approach for greater objectivity. Neuroimage 54, 2105–2115. doi: 10.1016/j.neuroimage.2010.10.033
Wirsich, J., Bénar, C., Ranjeva, J.-P., Descoins, M., Soulier, E., Le Troter, A., et al. (2014). Single-trial EEG-informed fMRI reveals spatial dependency of BOLD signal on early and late IC-ERP amplitudes during face recognition. Neuroimage 100, 325–336. doi: 10.1016/j.neuroimage.2014.05.075
Wise, R. G., Ide, K., Poulin, M. J., and Tracey, I. (2004). Resting fluctuations in arterial carbon dioxide induce significant low frequency variations in BOLD signal. Neuroimage 21, 1652–1664. doi: 10.1016/j.neuroimage.2003.11.025
Wong, C.-K., Zotev, V., Misaki, M., Phillips, R., Luo, Q., and Bodurka, J. (2016). Automatic EEG-assisted retrospective motion correction for fMRI (aE-REMCOR). Neuroimage 129, 133–147. doi: 10.1016/J.NEUROIMAGE.2016.01.042
Xia, H., Ruan, D., and Cohen, M. S. (2014a). Removing ballistocardiogram (BCG) artifact from full-scalp EEG acquired inside the MR scanner with Orthogonal Matching Pursuit (OMP). Front. Neurosci. 8:218. doi: 10.3389/fnins.2014.00218
Yarnykh, V. L. (2007). Actual flip-angle imaging in the pulsed steady state: a method for rapid three-dimensional mapping of the transmitted radiofrequency field. Magn. Reson. Med. 57, 192–200. doi: 10.1002/mrm.21120
Yuan, H., Zotev, V., Phillips, R., Drevets, W. C., and Bodurka, J. (2012). Spatiotemporal dynamics of the brain at rest - Exploring EEG microstates as electrophysiological signatures of BOLD resting state networks. Neuroimage 60, 2062–2072. doi: 10.1016/j.neuroimage.2012.02.031
Yuval-Greenberg, S., Tomer, O., Keren, A. S., Nelken, I., and Deouell, L. Y. (2008). Transient induced gamma-band response in EEG as a manifestation of miniature saccades. Neuron 58, 429–441. doi: 10.1016/j.neuron.2008.03.027
Keywords: simultaneous EEG-fMRI, data quality, neurovascular coupling
Citation: Abreu R, Leal A and Figueiredo P (2018) EEG-Informed fMRI: A Review of Data Analysis Methods. Front. Hum. Neurosci. 12:29. doi: 10.3389/fnhum.2018.00029
Received: 03 August 2017; Accepted: 18 January 2018;
Published: 06 February 2018.
Edited by:Iiro P. Jääskeläinen, Aalto University School of Science, Finland
Reviewed by:Xu Lei, Southwest University, China
Camillo Porcaro, Istituto di Scienze e Tecnologie della Cognizione (ISTC) – CNR, Italy
Copyright © 2018 Abreu, Leal and Figueiredo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rodolfo Abreu, email@example.com