The quest for EEG power band correlation with ICA derived fMRI resting state networks

Introduction

Blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) resting state networks (RSNs) have increasingly generated interest in the neuroscientific community, but the neuronal underpinnings remain unclear so far. Early studies, which examine correlations between the electroencephalography (EEG) theta, alpha, or beta band power and BOLD signal fluctuations using EEG derived regressors (Goldman et al., 2002; Laufs et al., 2003a,b, 2006; Moosmann et al., 2003; Feige et al., 2005; Goncalves et al., 2006; Scheeringa et al., 2008), report rather mixed and inconclusive BOLD correlation maps. The discovery and further analysis of RSNs (Biswal et al., 1995; Lowe et al., 2000; Cordes et al., 2001; Greicius et al., 2003; Fox et al., 2005; Damoiseaux et al., 2006; De Luca et al., 2006; Smith et al., 2009), together with the above mentioned early combined EEG-fMRI studies gave rise to the assumption that several frequency bands might be involved in distinct functional networks (Laufs et al., 2006; Mantini et al., 2007). The replication of this finding on subject level would fundamentally improve our understanding of the link with electrophysiology.

Simultaneous recordings of EEG and fMRI during resting state (RS), enables the investigation of the electrophysiological correlates of BOLD RSNs. Using simultaneous recordings, Mantini et al. (2007) reported a specific EEG frequency band power signature for RSNs on group level in eyes closed RS. However, further studies show large inter-subject variations of distinct brain areas correlated with EEG alpha band power (Goncalves et al., 2006, 2008) in RS. In a recent study by Meyer et al. (2013) electrophysiological correlation patterns (ECPs) between RSN BOLD time courses and EEG frequency band power showed large inter-subject and within subject variability. While RSNs by themselves exhibit a high reproducibility of their spatial characteristics across subjects, these studies point to less stable temporal correlations between RSNs seen in BOLD fMRI and EEG frequency band power.

Based on this evidence we hypothesize that the relationship between EEG frequency band power and RSN BOLD time courses is not stable over time. In order to assess this temporal variance in the correlation of the EEG signal and RSNs within a subject in this study, a dataset with a long RS of 34 min was split up into 15 segments and each was analyzed using the following three analysis approaches:

(1) Global frequency power correlation (GFPC) (Meyer et al., 2013) resulting in ECPs an approach that is similar to the one used by Mantini et al. (2007) who found stable correlation patterns on group level.

(2) An extended version of this method, including an anatomically informed analysis (Dale et al., 2000; Ou et al., 2010; Janssen et al., 2012) to separate the EEG based on RSN Z-maps within a subject, to obtain source frequency power correlation (SFPC), which should reduce the effect of volume conduction in the EEG.

(3) A channel wise frequency power fit (CFPF) with minimal assumptions, using the BOLD RSN time courses as the independent variable, which further reduces methodological bias.

We then calculated the temporal variance over the 15 segments for each of the three methods to estimate the temporal stability of the correlation between the two modalities.

Materials and Methods

Data Acquisition and Pre-Processing

In this study we performed a new analysis of the data sets acquired in Meyer et al. (2013). We briefly summarize the acquisition protocol and the pre-processing steps (for details, see Meyer et al., 2013): 34 min of eyes open RS were recorded from 12 healthy subjects, using combined EEG-fMRI, with approval of the local ethical committee. MR data were acquired on a 3 T Magnetom TIM Trio system (Siemens Healthcare, Erlangen, Germany) using the product 32 channel head coil. Functional data were recorded using a multi echo EPI sequence (Poser et al., 2006) (1030 Vol., TR = 2000 ms, 3.5 mm isotropic voxel size). A T1-weighted structural scan (MPRAGE) at 1 mm isotropic voxel size was also obtained (with EEG cap), to register the functional data to Montreal neurological institute (MNI) space. Five of the subjects (subjects 1, 2, 4, 10, and 11) were invited back to acquire a second T1-weighted structural scan without the EEG cap to enable the head model based analysis.

Simultaneous EEG data were recorded with a 32 channel cap (ANT WaveGuard MRI), using a BrainAmp MR plus amplifier (250 Hz low-pass analog hardware filter, 10 s time constant, 5 kHz sampling rate, 0.5 μV resolution, reference electrode: FCz) and BrainVision Recorder (BrainVision, Gilching, Germany). Two of the subjects were recorded with a 64 channel cap (BrainVision) using two BrainAmp MR plus amplifier; the same 30 channels (10–20 system) were used for all subjects in the analysis. The subjects were asked to relax, keep their eyes open, stay awake, and not think of anything specific. The room was darkened during the scan and an infrared eye tracker was used to confirm that the subject did not fall asleep. All subjects managed to stay awake for the complete duration of the experiment.

Functional magnetic resonance imaging pre-processing was performed using functions from the SPM5 software package (Welcome Department of Imaging Neuroscience, University College London, UK). The five echoes acquired at every time point were combined after SPM5 motion correction (Poser et al., 2006).

Electroencephalography pre-processing: MR related artifacts in the EEG signal were removed using Analyzer 2 (BrainVision). Trigger based average subtraction (Allen et al., 2000), as implemented in Analyzer 2, was applied to correct for gradient artifacts. The data were filtered using a Butterworth zero phase filter, 48 dB/oct, with a low cutoff at 0.8 Hz, to remove slow fluctuations from respiration, and a high cutoff at 50 Hz. Additionally, a notch filter at 50 Hz was used to remove residual mains frequency noise. Cardiac related MR artifacts were removed using the adaptive average subtraction (AAS) method of Analyzer 2 in semiautomatic mode (Allen et al., 1998). Further, eye blink related artifacts were removed using ICA and the EEG data were re-referenced to a common average.

Analysis

As motivated in the introduction, three distinct methods (see Figure 1) were used to infer whether the relationship between EEG frequency band power and RSN BOLD time courses is temporally instable. For all these methods, the preprocessed fMRI data were spatially smoothed by 5 mm and transformed to MNI space using FMRIB’s Software Library’s (FSL) Feat (version 4.1¹; Smith et al., 2004; Woolrich et al., 2009; Jenkinson et al., 2012). Group independent component analysis (ICA, as implemented in the FSL tool Melodic version 3.1) was performed on the fMRI data to obtain 30 ICs and 12 task related RSNs were selected according to Smith et al. (2009), see Figure 2 for a depiction of the RSNs. A dual regression approach was used to derive subject specific RSN maps and time courses (Filippini et al., 2009). The further analysis is described for each method separately below.

FIGURE 1

Figure 1. Overview of the three analysis methods used in this study. The highlighted regions and arrows are labeled accordingly.

FIGURE 2

Figure 2. RSNs on group level as maximum intensity projection on the central slices and their classification according to Smith et al. (2009).

Global Frequency Power Correlation

The datasets were split into 15 sections of equal length, each still longer than 2 min. For every section the EEG signal was split into 2 s segments corresponding to the TR used in the MR-acquisition. Within each section, for every segment, the mean frequency power over all channels for four frequency bands, i.e., delta: (2–4) Hz, theta: (4–8) Hz, alpha: (8–12) Hz, and beta: (12–30) Hz, was calculated, using a fast Fourier transformation (FFT), resulting in one time series for each frequency band. Motion related artifacts in the frequency power time courses were corrected. The frequency power time series were convolved with the standard SPM5 hemodynamic response function (HRF) and correlated with the RSN time courses taking into account common variance (partial correlation) between frequency bands. The correlation values were Z-transformed, using the mean over all correlation values across subjects as global mean, which resulted in time series of 15 Z-scores for every frequency band (see Figure 3). The temporal variance for each RSN and frequency band over the 15 time points was calculated and averaged over subjects (see Table 1). To estimate the temporal stability of ECPs within and across subjects, for each RSN and frequency band the Z-scores of the 15 sections were ranked from high to low, and averaged over subjects to visualize inter-subject variance.

FIGURE 3

Figure 3. Exemplary representation of data for the three methods used, depicting the results for RSN3 (lateral visual component) of subject 1. On the upper left hand side the subject specific map of RSN3 is plotted as maximum intensity projection on the central slices. The two graphs in the middle show the ECP time courses for each frequency band for GSPC and SFPC, illustrating the temporal variance. In the rank graph below the two methods are plotted on top of each other, depicting the similar variance of SFPC (solid line) and GFPC (dot-dashed line); for both methods alpha shows more negative correlation with BOLD in this visual component. On the upper right hand side the temporal sequence of ICEPs as calculated by CFPF for the four frequency bands is shown. Below the results of the spatial correlation of the subsequent ICEPs are depicted, showing no stable temporal signature.

TABLE 1

Table 1. Group mean temporal variance values (variance across the 15 sections) for GFPC and CFPF, as well as the mean temporal variance values of the same five subjects analyzed with SFPC and GFPC, for each RSN and frequency band.

Source Frequency Power Correlation

In order to get an indication for the effect of volume conduction and obtain more specific correlation patterns, in five subjects an in-house developed fMRI-informed source model was applied. In combination with a four layer realistic head model it enables to separate the EEG according to the fMRI-RSNs. This new method was tested in a separate study that employs a simple visual stimulation and is further referred to as Integrative Head Model (IHM). It merges FSL analysis, Freesurfer mesh generation (Freesurfer image analysis suite²), and a Neuroelectromagnetic Forward Head Modeling Toolbox (NFT) based head model (VER 2.0³; Acar and Makeig, 2010), to combine fMRI and EEG in an integrative way (see Figure 1). Tissue surface meshes (TSMs) from the individual T1 images are derived using Freesurfer and NFT. The scalp, inner and outer skull as well as brain TSMs are used in the Boundary Element Method (BEM) based forward model as implemented in NFT (Brain/scalp conductivity = 0.33 S/m, Skull conductivity = 0.0132 S/m, CSF conductivity = 1.79 S/m). The source space is constructed by seeding the cortical sheet with dipoles, the location, and orientation of which is derived from the pial and white matter TSMs. Sources are defined by selecting dipoles according to fMRI RSNs, mapped to the cortical sheet. A source is defined as the weighted vector sum of its active dipoles, where the weights are equal to the Z values of the fMRI activation map at dipole location, and normalized subsequently so that the sum of all weights within one source equals one. These fMRI derived sources are fed into the forward model (in NFT) to calculate the specific lead field matrix (LFM) using an electrode template, which was manually transformed to each subjects head. This specific LFM has a low dimensionality given by the number of sources times number of channel. It is inverted using a Moore–Penrose pseudo inverse and the inverted LFM is used to transform the EEG data to source specific time courses. This results in an EEG time course for each RSN. Furthermore, the same analysis steps as described in Method 1 were applied to the transformed EEG signal. For each of the 15 sections and each frequency band the fMRI derived source frequency power time courses were convolved with the standard SPM5 HRF, partial correlated with their associated RSN time course, and the correlation values were Z-transformed. The variance over the 15 sections was calculated and the Z-scores of the sections were ranked from high to low, to obtain an estimate of the temporal stability.

Channel Wise Frequency Power Fit

After pre-processing, each channel of the EEG data was band pass filtered in four frequency bands [delta: (2–4) Hz, theta: (4–8) Hz, alpha: (8–12) Hz, and beta: (12–30) Hz] using an FFT-filter (EEGlab). Power time courses were obtained from the filtered data by applying a Hilbert transform and taking the squared magnitude of the resulting signal. To correct for movement, time points where the power estimate exceeded a threshold (seven times the mean of the time course) were set to the average of the time points immediately before and after. Power time courses were segmented in to 2 s segments, according to the TR used in the fMRI acquisition; subsequently each segment was averaged over time and the resulting frequency power time course for each channel were convolved with an HRF (SPM 5). Finally the HRF convolved frequency power time course and the RSN time courses were normalized to have zero mean and a standard deviation of one. Time courses of all ICs (including noise related components) were fitted to each frequency power time course in a separate GLM for every channel. This resulted in an estimate of signal contribution for each RSN to each electrode and EEG frequency band. Plotting these contribution estimates on a scalp plot, here termed independent component expression pattern (ICEP), gives a visual representation of the electrophysiological expression of the RSN for each frequency band. Applying this approach to each of the 15 sections resulted in 15 subsequent ICEPs representing their evolution over time.

In order to obtain a comparable estimate for this method, which gives a spatial distribution as opposed to a point estimate of the other two methods, the temporal stability of the ICEPs was assessed by calculating the spatial correlation between subsequent sections within one subject, RSN and EEG frequency band. For each of those, the correlation values were Z-transformed using bootstrap statistics and the Z-scores were averaged to obtain the mean over all combinations of sections. For the bootstrapped Z-transformation a distribution was generated by repeatedly (n = 10,000) selecting 15 ICEPs at random from the entire set of ICEPs for that subject applying the same spatial correlation analysis. For each RSN and frequency band the Z-scores of the 15 sections were ranked from high to low, and for group analysis averaged over subjects. Additionally the variance over the 15 sections as well as the average variance over subjects for each RSN and frequency band was calculated.

Results

As reported in Meyer et al. (2013) we found reproducible fMRI RSNs across subjects (see Figure 2 for a depiction of the RSNs). In this study we observed very large inter-subject and intra-subject variability in the EEG frequency power correlations across all applied analysis methods. Figure 3 depicts the output of the different methods for one network (RSN3) of subject 1. It is clearly visible that GFPC and SFPC are not stable in time regarding their EEG frequency power correlation with the RSN time courses for all frequency bands. Figure 4 shows the results of the group analysis for GFPC and Figure 5 the results for the five subjects analyzed with SFPC. In both figures the group rank plots for four different RSNs show a large temporal variance within a subject – as reflected in the variance of the ranked Z-scores – in the depicted RSNs for all frequency bands. The error bars, indicating the standard deviation across subjects, show the considerable inter-subject variability. The error bars in Figure 5 are larger compared to those in Figure 4 which cannot be explained by the smaller number of analyzed subjects as controlled by performing GFPC on the same five subjects as for SFPC. Also note that, using SFPC the overall observed Z-scores are lower compared to GFPC. Strikingly, one can see in the ranking plots that for the visual components (see RSN2 and RSN3 in Figures 4 and 5) alpha power shows a more negative correlation with the BOLD signal whereas delta power shows a more positive correlation. This is also the case for the third visual component (not shown).

FIGURE 4

Figure 4. Exemplary group results of RSN2 (occipital visual component), RSN3 (lateral visual component), RSN4 (DMN), and RSN6 (sensory motor component) for GFPC as rank graph showing large temporal variance within a subject. The error bars (standard deviation across subjects) show the considerable inter-subject variability. Clearly alpha power shows a more negative correlation with the BOLD signal whereas delta power shows a more positive correlation for the visual components. Note that the connecting lines are only for visualization purposes.

FIGURE 5

Figure 5. Exemplary results of RSN2 (occipital visual component), RSN3 (lateral visual component), RSN4 (DMN), and RSN6 (sensory motor component) for the five analyzed subjects using SFPC comparable to Figure 4. The error bars (standard deviation across subjects) show the considerable inter-subject variability, which is higher compared to GFPC; note that the Z-scores are smaller compared to GFPC. Despite these differences, also when using SFPC, alpha power shows a more negative correlation with the BOLD signal whereas delta power shows a more positive correlation for the visual components. The connecting lines are only for visualization purposes.

On the right hand side of Figure 3 the temporal sequence of ICEPs as calculated by CFPF for the four frequency bands as well as the results of the spatial correlation of the subsequent ICEPs is depicted. Clearly there is no temporally stable signature in the scalp maps for different sections of the dataset. This can be also observed in the group rank plot in Figure 6. The data shown in Figure 3 as well as the rank plots in Figures 4, 5, and 6 are typical examples for all analyzed subjects, all RSNs, and the four frequency bands examined, respectively. Table 1 summarizes the results containing the group mean temporal variance across the 15 sections, for each RSN and frequency band for GFPC and CFPF, respectively, as well as for the same five subjects analyzed with SFPC and GFPC.

FIGURE 6

Figure 6. Exemplary group results of RSN2 (occipital visual component) and RSN4 (DMN) for CFPF as rank graph showing no stable temporal signature. The connecting lines are only for visualization purposes.

Discussion

The three methods used in this study were chosen to examine the temporal variability of ECPs from different perspectives with the aim to minimize methodological bias. GFPC is a very conservative approach with fairly little assumptions, taking the global EEG frequency power as independent parameter. However, due to the mixed nature of the EEG signal, volume conduction cannot be excluded, which might cause several sources contributing to a certain correlation and might explain temporally unstable ECPs. SFPC addresses this shortcoming by separating the EEG signal according to the fMRI-RSNs to correlate with, but this also did not result in temporally stable ECPs.

To test for possible methodological bias of GFPC and SFPC as source for the observed variance, we analyzed the data sets using a third approach. CFPF uses the RSN timelines as independent parameters and only uses the HRF to model the relation between the two modalities. Also this approach did not result in temporally stable correlation patterns. While the human HRF itself shows quite complex spatial dependencies (de Munck et al., 2007, 2009), in our correlation analysis it mainly causes a constant time shift and temporal smoothing, therefore it cannot be the reason for temporally instable correlation. Together with the findings of GFPC and SFPC, this leads to the suggestion that the analyzed low dimensional RSNs do not have a temporally stable relationship with EEG frequency band power fluctuations. However, the observed negative correlation of alpha power with the BOLD time courses for the visual components reaches statistical significance within three subjects for GFPC in agreement with previous literature (Goldman et al., 2002; Laufs et al., 2003a, 2006; Goncalves et al., 2006, 2008; Meyer et al., 2013), but does not reach statistical significance for either method on the group level.

One possible explanation for our observation of temporally instable ECPs might be given by Smith et al. (2012), who applied temporal ICA on high dimensional (200 spatial IC components) fMRI-RSNs and reported vast temporal dynamics within the lower dimensional (20–30 spatial IC components) RSNs. As such, there still might be a direct relation between RSNs and EEG frequency band power, but on a smaller spatial scale. However, one would expect a certain temporal stability in the results of SFPC and CFPF even if just a subcomponent of the low dimensional RSN expresses itself in a given EEG frequency band, which was not observed in our study.

An alternative explanation of our results would be, that during RS, frequency-specific power in the lower frequency bands of the EEG is not linked to changes in neuronal activity, reflected in changed oxygen consumption as measured by BOLD fMRI. This would also be supported by recent animal studies, e.g., Schölvinck et al. (2010), that show no stable correlation for the lower frequency bands in EEG with BOLD fMRI particularly in eyes open RS.

The observation that using SFPC reduced the overall observed Z-scores compared to GFPC gives rise to the assumption that the studied RSN characteristics are not related. However, one has to consider the potential limitation of the head model as used in our study; (a) the spatial resolution of the head model is limited by the relatively low number of electrodes and (b) its reduced spatial specificity in the context of spatially extended sources like RSNs, as we assume a concurrent temporal behavior within the whole source.

We therefore conclude that the correlation between lower frequency band power in EEG and BOLD RSNs time courses is at least temporally instable or even absent in eyes open RS.

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.

Acknowledgments

Financial support from Provincie Gelderland and Overijssel, as well as from the ministry of economical affairs in the framework of the project “Pieken in de Delta – ViP Brainnetworks” is acknowledged. We gratefully acknowledge technical support from ANT.

Footnotes

1. ^www.fmrib.ox.ac.uk/fsl
2. ^http://surfer.nmr.mgh.harvard.edu/
3. ^http://sccn.ucsd.edu/nft/

References

Acar, Z., and Makeig, S. (2010). Neuroelectromagnetic forward head modeling toolbox. J. Neurosci. Methods 190, 258–270. doi:10.1016/j.jneumeth.2010.04.031

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Allen, P. J., Josephs, O., and Turner, R. (2000). A method for removing imaging artifact from continuous EEG recorded during functional MRI. Neuroimage 12, 230–239. doi:10.1006/nimg.2000.0599

CrossRef Full Text

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

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Biswal, B., Yetkin, F. Z., 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

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

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. AJNR Am. J. Neuroradiol. 22, 1326–1333.

Pubmed Abstract | Pubmed Full Text

Dale, A. M., Liu, A. K., Fischl, B. R., Buckner, R. L., Belliveau, J. W., Lewine, J. D., et al. (2000). Dynamic statistical parametric mapping: combining fMRI and MEG for high-resolution imaging of cortical activity. Neuron 26, 55–67. doi:10.1016/S0896-6273(00)81138-1

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Damoiseaux, J. S., Rombouts, S. A., Barkhof, F., Scheltens, P., Stam, C. J., Smith, S. M., et al. (2006). Consistent resting-state networks across healthy subjects. Proc. Natl. Acad. Sci. U.S.A. 103, 13848–13853. doi:10.1073/pnas.0601417103

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

De Luca, M., Beckmann, C. F., De Stefano, N., Matthews, P. M., and Smith, S. M. (2006). fMRI resting state networks define distinct modes of long-distance interactions in the human brain. Neuroimage 29, 1359–1367. doi:10.1016/j.neuroimage.2005.08.035

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

de Munck, J. C., Gonçalves, S. I., Huijboom, L., Kuijer, J. P., Pouwels, P. J., 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

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

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

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Feige, B., Scheffler, K., Esposito, F., Di Salle, F., Hennig, J., and Seifritz, E. (2005). Cortical and subcortical correlates of electroencephalographic alpha rhythm modulation. J. Neurophysiol. 93, 2864–2872. doi:10.1152/jn.00721.2004

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Filippini, N., MacIntosh, B. J., Hough, M. G., Goodwin, G. M., Frisoni, G. B., Smith, S. M., et al. (2009). Distinct patterns of brain activity in young carriers of the APOE-epsilon4 allele. Proc. Natl. Acad. Sci. U.S.A. 106, 7209–7214. doi:10.1073/pnas.0811879106

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Fox, M. D., Snyder, A. Z., Vincent, J. L., Corbetta, M., Van Essen, D. C., and Raichle, M. E. (2005). The human brain is intrinsically organized into dynamic, anticorrelated functional networks. Proc. Natl. Acad. Sci. U.S.A. 102, 9673–9678. doi:10.1073/pnas.0504136102

CrossRef Full Text

Goldman, R. I., Stern, J. M., Engel, J. Jr., and Cohen, M. S. (2002). Simultaneous EEG and fMRI of the alpha rhythm. Neuroreport 13, 2487–2492. doi:10.1097/00001756-200212200-00022

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Goncalves, S. I., Bijma, F., Pouwels, P. J., Jonker, M., Kuijer, J. P., Heethaar, R. M., et al. (2008). A data and model-driven approach to explore inter-subject variability of resting-state brain activity using EEG-fMRI. IEEE J. Sel. Top. Signal Process. 2, 944–953. doi:10.1109/JSTSP.2008.2009082

CrossRef Full Text

Goncalves, S. I., de Munck, J. C., Pouwels, P. J., Schoonhoven, R., Kuijer, J. P. A., Maurits, N. M., et al. (2006). Correlating the alpha rhythm to BOLD using simultaneous EEG/fMRI: inter-subject variability. Neuroimage 30, 203–213. doi:10.1016/j.neuroimage.2005.09.062

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Greicius, M. D., Krasnow, B., Reiss, A. L., and Menon, V. (2003). Functional connectivity in the resting brain: a network analysis of the default mode hypothesis. Proc. Natl. Acad. Sci. U.S.A. 100, 253–258. doi:10.1073/pnas.0135058100

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Janssen, R. J., Meyer, M. C., Beckmann, C. F., and Barth, M. (2012). “A framework to estimate electrophysiological dynamics of fMRI derived sources,” in Third Biennial Conference on Resting State Brain Connectivity, Magdeburg, Poster 46.

Jenkinson, M., Beckmann, C. F., Behrens, T. E., Woolrich, M. W., and Smith, S. M. (2012). FSL. Neuroimage 62, 782–790. doi:10.1016/j

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

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

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Laufs, H., Kleinschmidt, A., Beyerle, A., Eger, E., Salek-Haddadi, A., Preibisch, C., et al. (2003a). EEG-correlated fMRI of human alpha activity. Neuroimage 19, 1463–1476. doi:10.1016/S1053-8119(03)00286-6

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Laufs, H., Krakow, K., Sterzer, P., Eger, E., Beyerle, A., Salek-Haddadi, A., et al. (2003b). Electroencephalographic signatures of attentional and cognitive default modes in spontaneous brain activity fluctuations at rest. Proc. Natl. Acad. Sci. U.S.A. 100, 11053–11058. doi:10.1073/pnas.1831638100

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Lowe, M. J., Dzemidzic, M., Lurito, J. T., Mathews, V. P., and Phillips, M. D. (2000). Correlations in low-frequency BOLD fluctuations reflect cortico-cortical connections. Neuroimage 12, 582–587. doi:10.1006/nimg.2000.0654

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Mantini, D., Perrucci, M. G., Del Gratta, C., Romani, G. L., and Corbetta, M. (2007). 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

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Meyer, M. C., van Oort, E. S., and Barth, M. (2013). Electrophysiological correlation patterns of resting state networks in single subjects: a combined EEG-fMRI study. Brain Topogr. 26, 98–109. doi:10.1007/s10548-012-0235-0

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

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

CrossRef Full Text

Ou, W., Nummenmaa, A., Ahveninen, J., Belliveau, J. W., Hämäläinen, M. S., and Golland, P. (2010). Multimodal functional imaging using fMRI-informed regional EEG/MEG source estimation. Neuroimage 52, 97–108. doi:10.1016/j.neuroimage.2010.03.001

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Poser, B. A., Versluis, M. J., Hoogduin, J. M., and Norris, D. G. (2006). BOLD contrast sensitivity enhancement and artifact reduction with multiecho EPI: Parallel-acquired inhomogeneity-desensitized fMRI. Magn. Reson. Med. 55, 1227–1235. doi:10.1002/mrm.20900

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Scheeringa, R., Bastiaansen, M. C., 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

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Schölvinck, M. L., Maier, A., Ye, F. Q., Duyn, J. H., and Leopold, D. A. (2010). Neural basis of global resting-state fMRI activity. Proc. Natl. Acad. Sci. U.S.A. 107, 10238–10243. doi:10.1073/pnas.0913110107

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Smith, S. M., Fox, P. T., Miller, K. L., Glahn, D. C., Fox, P. M., Mackay, C. E., et al. (2009). Correspondence of the brain’s functional architecture during activation and rest. Proc. Natl. Acad. Sci. U.S.A. 106, 13040–13045. doi:10.1073/pnas.0905267106

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Smith, S. M., Jenkinson, M., Woolrich, M. W., Beckmann, C. F., Behrens, T. E., Johansen-Berg, H., et al. (2004). Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage 23, 208–219. doi:10.1016/j.neuroimage.2004.07.051

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Smith, S. M., Miller, K. L., Moeller, S., Xu, J., Auerbach, E. J., Woolrich, M. W., et al. (2012). Temporally-independent functional modes of spontaneous brain activity. Proc. Natl. Acad. Sci. U.S.A. 109, 3131–3136. doi:10.1073/pnas.1121329109

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Woolrich, M. W., Jbabdi, S., Patenaude, B., Chappell, M., Makni, S., Behrens, T., et al. (2009). Bayesian analysis of neuroimaging data in FSL. Neuroimage 45, 173–186. doi:10.1016/j.neuroimage.2008.10.055

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text