PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Neuroimage. Author manuscript; available in PMC 2010 December 1.
Published in final edited form as:
PMCID: PMC2995588
NIHMSID: NIHMS250655

Effects of model-based physiological noise correction on default mode network anti-correlations and correlations

Abstract

Previous studies have reported that the spontaneous, resting-state time course of the default-mode network is negatively correlated with that of the “task-positive network”, a collection of regions commonly recruited in demanding cognitive tasks. However, all studies of negative correlations between the default-mode and task-positive networks have employed some form of normalization or regression of the whole-brain average signal (“global signal”); these processing steps alter the time series of voxels in an uninterpretable manner as well as introduce spurious negative correlations. Thus, the extent of negative correlations with the default mode network without global signal removal has not been well characterized, and it is has recently been hypothesized that the apparent negative correlations in many of the task-positive regions could be artifactually induced by global signal pre-processing. The present study aimed to examine negative and positive correlations with the default-mode network when model-based corrections for respiratory and cardiac noise are applied in lieu of global signal removal. Physiological noise correction consisted of (1) removal of time-locked cardiac and respiratory artifacts using RETROICOR (Glover et al., 2000), and (2) removal of low-frequency respiratory and heart rate variations by convolving these waveforms with pre-determined transfer functions (Birn et al., 2008; Chang et al., 2009) and projecting the resulting two signals out of the data. It is demonstrated that negative correlations between the default-mode network and regions of the task-positive network are present in the majority of individual subjects both with and without physiological noise correction. Physiological noise correction increased the spatial extent and magnitude of negative correlations, yielding negative correlations within task-positive regions at the group level (p<0.05, uncorrected; no regions at the group level were significant at FDR=0.05). Furthermore, physiological noise correction caused region-specific decreases in positive correlations within the default-mode network, reducing apparent false positives. It was observed that the low-frequency respiratory volume and cardiac rate regressors used within the physiological noise correction algorithm displayed significant (but not total) shared variance with the global signal, and constitute a model-based alternative to correcting for non-neural global noise.

Keywords: fMRI, physiological noise, resting state, networks, negative correlations, functional connectivity, default mode network

Introduction

Functional magnetic resonance imaging (fMRI) studies of the brain have revealed temporal synchrony between the blood oxygen level dependent (BOLD) signals of anatomically distinct regions during the absence of an explicit task, i.e. in the resting state (for reviews, see (Buckner and Vincent, 2007; Fox and Raichle, 2007)). Several such sets of temporally correlated regions, known as resting-state networks, can be identified consistently across human subjects (Beckmann et al., 2005; Damoiseaux et al., 2006; De Luca et al., 2006) and are presumed to reflect a basic functional organization of the brain.

A particular resting-state network, known as the default-mode network (DMN), has received much recent attention (see (Buckner et al., 2008) for a review). The DMN comprises brain regions that, in addition to having coherent fluctuations in the absence of a task, are commonly deactivated (exhibit sub-baseline signal deflections) in response to a broad range of cognitive and attentional tasks (Mazoyer et al., 2001; Raichle et al., 2001; Shulman et al., 1997) and may underlie self-referential processing (Gusnard et al., 2001) (see (Buckner et al., 2008) for a review). Core regions of the DMN include the precuneus/posterior cingulate cortex (PCC), medial prefrontal cortex, ventral anterior cingulate, lateral parietal cortex, inferior temporal cortex, and parahippocampal cortex.

It has been reported that the spontaneous, resting-state time course of the DMN is negatively correlated (“anti-correlated”) with that of the task-positive network (TPN), a set of regions that are (positively) activated in response to cognitive and attentional tasks (Fox et al., 2005; Fransson, 2005). The TPN includes dorsolateral prefrontal cortex (DLPFC), supramarginal gyrus and posterior parietal cortex, insula, premotor cortex, and supplementary motor area (SMA). The finding of negative correlations between the two networks has been hypothesized to reflect a competitive relationship between internally- and externally-oriented processes, or between mind wandering and focused attention, and mirrors the opposing signal changes displayed by the two networks during controlled tasks (Fox et al., 2005; Fransson, 2005; Greicius et al., 2003). Subsequent studies have asserted that the strength of the negative correlations between the DMN and TPN mediates behavioral variability (Clare Kelly et al., 2008) and may distinguish Alzheimer's patients from healthy controls (Wang et al., 2007).

However, studies reporting significant negative correlations between the DMN and the TPN in resting state have removed a whole-brain average time series from the data, either by proportional scaling (time-point-wise division), regressing it out as a pre-processing step, or including it as a nuisance covariate when assessing functional connectivity using the general linear model (Clare Kelly et al., 2008; Fair et al., 2008; Fox et al., 2005; Fransson, 2005, 2006; Greicius et al., 2003; Tian et al., 2007; Uddin et al., 2009; Wang et al., 2007). Removing a global signal is intended to remove undesired non-neural fluctuations that affect large numbers of voxels. Unfortunately, since the global signal is derived from the data itself and is an unknown mixture of neural and non-neural fluctuations, these steps alter inter-regional correlations and complicate their interpretation. Furthermore, if the global signal has been regressed out of each voxel's time series, the resulting correlation coefficients between a seed voxel and all other voxels in the brain will have a mean value that is less than or equal to zero, necessarily introducing spurious negative correlations (Fox et al., 2008; Murphy et al., 2009). Consequently, it has been hypothesized that global signal removal may be the cause of the observed negative correlations between the DMN and TPN (Murphy et al., 2009). Studies by (Golland et al., 2007; Murphy et al., 2009) did not employ global signal removal and failed to observe negative correlations at any statistical threshold, but Fransson (Fransson, 2005) noted that when reanalyzing data without performing global signal regression, negative correlations in some TPN regions were present when the statistical threshold was lowered from p<0.01 to p<0.07 (corrected), and Uddin et al. (2009) showed that bilateral clusters in the insula were negatively correlated with the PCC when data were reanalyzed without global signal regression. The question of whether negative correlation between the two networks is a real phenomenon, or simply an artifact of global signal removal, is currently under debate (Buckner et al., 2008).

The present study aimed to determine the spatial extent and magnitude of negative and positive correlations with the DMN when model-based physiological noise corrections, instead of global signal removal, are applied. Respiration and cardiac processes are known to modulate the BOLD signal (Biswal et al., 1993; Dagli et al., 1999; Jezzard et al., 1993; Shmueli et al., 2007; Weisskoff et al., 1993; Wise et al., 2004), and their effects can be reduced by monitoring respiratory and cardiac cycles throughout the scan and filtering them out of the fMRI data using a priori models (Birn et al., 2006; Birn et al., 2008; Chang et al., 2009; de Munck et al., 2008; Glover et al., 2000). In this way, known sources of noise that affect functional connectivity can be reduced in an unbiased manner. A previous study, employing the same physiological noise corrections applied here, demonstrated increased spatial specificity of positive correlations within the DMN after correction (Chang et al., 2009). It was thus hypothesized that such corrections might increase the sensitivity of detecting regions that are negatively correlated with the DMN as well.

It is important to note, however, that the determination of sensitivity and specificity in functional connectivity analysis is not straightforward. In ROI-based functional connectivity analysis, where the time series of both the reference (seed ROI) and each voxel in the brain are considered to contain both a neuronal component and a noise component, the null hypothesis (H0) is that there is no correlation between the respective neuronal components of the seed ROI and a given voxel. Accordingly, a Type I error (false positive) occurs when there is an absence of neuronal correlation, and yet the noise components correlate strongly enough so that the overall correlation coefficient exceeds a threshold; a Type II error (false negative) occurs when the neuronal components of the seed ROI and a given voxel are correlated, but the correlation is not detected because of noise in the seed ROI and/or the queried voxel. In contrast with a standard analysis of task activation, wherein the reference time series is the expected noise-less neuronal time course (often modeled by a binary stimulus waveform convolved with a hemodynamic response function), the reference time series in functional connectivity is derived from a region of the brain itself, and it is not known which part of the signal is due to neural activity and which is due to noise. Nevertheless, physiological noise correction – by reducing the noise component of BOLD signal time series – aims to improve the estimation of neuronally-driven correlations between regions of the brain. Here, after correction, it was hypothesized that: (1) positive correlations with the DMN seed time series will diminish, as correlations will be reduced in voxels that had been falsely coupled due to the common influence of respiratory and cardiac noise, and that (2) negative correlations with the DMN seed time series, if they exist, will increase due to the reduction of respiratory and cardiac noise that can obscure these dynamics.

Methods

Subjects

Participants included 15 healthy adults (7 female, aged 29 ± 11.5 years). All subjects provided written, informed consent, and all protocols were approved by the Stanford Institutional Review Board.

Imaging parameters

Magnetic resonance imaging was performed at 3.0 T using a GE whole-body scanner (GE Healthcare Systems, Milwaukee, WI). Seven of the 15 subjects were scanned on a GE Signa HDX (rev. 12M5) using a custom quadrature birdcage head coil, while the remaining 8 subjects were scanned on a GE Signa 750 (rev. 20) using an 8-channel head coil. Head movement was minimized with a bite bar. Thirty oblique axial slices were obtained parallel to the AC-PC with 4-mm slice thickness, 1-mm skip. T2-weighted fast spin echo structural images (TR = 3000 ms, TE = 68 ms, ETL = 12, FOV = 22 cm, matrix 192 × 256) were acquired for anatomical reference. A T2*-sensitive gradient echo spiral-in/out pulse sequence (Glover and Lai, 1998; Glover and Law, 2001) was used for functional imaging (TR = 2000 ms, TE = 30 ms, flip angle = 77°, matrix 64 × 64, FOV = 22 cm, same slice prescription as the anatomic images). A high-order shimming procedure was used to reduce B0 heterogeneity prior to the functional scans (Kim et al., 2002). Importantly, a frequency navigation correction was employed during reconstruction of each image to eliminate blurring from breathing-induced changes in magnetic field; no bulk mis-registration occurs from off-resonance in spiral imaging (Pfeuffer et al., 2002).

Resting state scan

All subjects underwent a resting state scan, for which they were instructed only to keep their eyes closed and remain awake. For 7 of the subjects, the duration of the scan was 8 min; for the remaining 8 subjects, the duration was 16 min, but only the first 8 min of data were used for the present analysis.

Working memory task

Seven of the 15 subjects (all scanned on the Signa HDX) performed an event-related Sternberg working memory (WM) task (Sternberg, 1969), consisting of 3 back-to-back 5 s trials (0.5 s encoding, 3 s maintenance, 1.5 s probe) followed by 45 s of fixation, repeated 10 times (total duration = 10 min). Encoding stimuli consisted of 4 randomly-selected uppercase letters in a cross-like configuration around the center of the screen, and the probe stimulus was a single lowercase letter presented at the center of the screen. This task was used only in choosing a seed coordinate for the resting-state functional connectivity analysis and for identifying regions of interest (ROIs) associated with task-induced deactivation (described below).

Physiological monitoring

Cardiac and respiratory processes were monitored using the scanner's built-in photoplethysmograph placed on a finger of the left hand and a pneumatic belt strapped around the upper abdomen, respectively. Cardiac and respiratory data were both sampled at 40 Hz on the Signa HDX, and at 100 Hz and 25 Hz, respectively, on the Signa 750. A file containing cardiac trigger times and respiratory waveforms was generated for each scan by the scanner's software. Values in the respiratory waveform were converted to a percentage of the full scale (difference between the maximum and minimum belt positions measured over the scan). Only the fractional variations in the waveform, rather than the absolute amplitude values, were of importance in the current study.

fMRI data pre-processing

Functional images were pre-processed using custom C and Matlab routines. Pre-processing included slice-timing correction using sinc interpolation, spatial smoothing with a 3D Gaussian kernel (FWHM=5 mm), and removal of linear and quadratic temporal trends. No additional temporal band-pass filtering was performed. The first 6 time frames were discarded to allow the MR signal to equilibrate. Spatial normalization to a standard template was not performed; computations were done in the original subject-space and voxel dimensions (3.4375 mm × 3.4375 mm × 4 mm).

Motion analysis

Motion parameters were calculated using methods described in (Friston et al., 1996). Image coregistration was not performed because of possible interpolation errors and unintended smoothing across voxels. Therefore, even though a bite bar was used, it was important to verify that motion was minimal. The maximum peak-to-peak excursion and root mean square (RMS) fluctuation for the 3 translational and 3 rotational motion parameter time series were calculated for the segment of resting state data used in the analysis. Rotations were converted to worst-case translations by multiplying by 65 mm, an average head radius (Thomason and Glover, 2008). Summary statistics are reported as the maximum of these values over the 6 axes of motion.

Physiological noise correction

Physiological noise correction consisted of (1) removal of time-locked cardiac and respiratory artifacts using RETROICOR (Glover et al., 2000), and (2) removal of low-frequency respiratory and heart rate effects using methods described in (Chang et al., 2009), which will be referred to as RVHRCOR. Briefly, RVHRCOR utilizes the following steps: First, a “respiratory variations” (RV) time series is computed as the root-mean-square amplitude of the respiration waveform across a 6 s sliding window, and is then convolved with an impulse response known as the “respiration response function” (Birn et al., 2008). Secondly, a heart rate (HR) time series is computed as the inverse of the average beat-to-beat interval in a 6 s sliding window, and convolved with the “cardiac response function” (Chang et al., 2009). The convolved RV and HR waveforms (RVx and HRx, respectively) are then simultaneously projected out of every voxel's time series using multiple regression (see Appendix A for further details). It was determined that reversing the order in which RETROICOR and RVHRCOR were applied made little difference. Unless otherwise specified, “physiological noise correction” will henceforth refer to the application of the 2 corrections described here.

WM (de)activation analysis

WM data were corrected for physiological noise. A binary (“boxcar”) WM task design waveform was convolved with a standard hemodynamic response function (Glover, 1999), and linear regression was applied to produce a negative (deactivation) contrast map for this covariate. The resulting individual-subject contrast maps were normalized to the SPM5 EPI template and entered into a group-level random-effects analysis using SPM5 (http://www.fil.ion.ucl.ac.uk/spm). Although only a subset of subjects performed the WM task, this task was used only to localize a seed ROI for functional connectivity and to define regions exhibiting task-induced deactivation for post-hoc ROI analysis of functional connectivity maps (described below).

Resting-state functional connectivity analysis

The DMN and its anti-correlated network were quantified using a standard ROI-based functional connectivity analysis (Biswal et al., 1995). An 8mm-radius sphere centered at (-4,-58,30) mm (Talairach coordinate = (-4,-55,30)) in the precuneus/PCC was chosen as the seed ROI. This center coordinate was a local maximum in the group-level WM deactivation map derived as above, and the resulting spherical ROI was observed to overlap with those of several previously-published studies (Fox et al., 2005; Fransson, 2005; Greicius et al., 2003). The ROI was defined in MNI space using MarsBar (http://marsbar.sourceforge.net) and reverse-normalized to each subject's mean functional image using SPM5. For each subject, the mean time series across all voxels in the ROI was extracted, and the Pearson correlation coefficient (r) was computed with the time series of every voxel in the brain. Positive and negative correlation maps were thresholded at |r|>0.15 (T>2.34, corresponding to p<0.01 uncorrected for N=240 time points).

The process of extracting the mean precuneus/PCC ROI time series and correlating it against all brain voxels was performed both before and after physiological noise correction was applied to the data. For comparison, these steps were also performed after a global signal removal (and no physiological noise correction) was applied to the data. Global signal removal was carried out by extracting the whole-brain average time series and projecting it out of every voxel's time series using linear regression. In other words, if y is the (mean-centered) time series of a voxel and g is the whole-brain average time series, then the voxel time series after global signal removal (yglob) would be yglob = yg(gTg)−1 gTy.

An additional comparison was performed, in which only the time series of white matter and cerebrospinal fluid (CSF) ROIs were regressed out of the data prior to correlation analysis. The assumption underlying this method is that white matter and CSF time series may carry physiological fluctuations that are similar to those affecting gray matter, while containing little contribution from neural activity (Behzadi et al., 2007; Fox et al., 2005). But, as previous studies of resting-state networks have removed white matter and CSF signals only in conjunction with removing the global average brain signal (e.g. Fox et al., 2005; Fair et al., 2008; Clare Kelly et al., 2008), the effects of regressing out white matter and CSF signals alone are unknown. This method is of particular interest because it can be applied without separately monitoring physiological responses, and – although derived from the fMRI data itself – does not pose the same problems as global signal removal since the white matter and CSF ROIs should not overlap with resting-state networks. In the present study, white matter and CSF (right lateral ventricle) ROIs were first defined as 6mm-diameter spheres in MNI space, where the choice of center coordinates for each ((26,-12,35) and (19, -33, 18), respectively) was guided by the SPM5 a priori masks. These spherical ROIs were then reverse-normalized to each subject's native space, and intersected with subject-specific white matter and CSF partial volume images (thresholded at a partial volume fraction >90%) to minimize overlap with gray matter. Subject-specific partial volume images were obtained by segmenting the T2-weighted anatomical image using the FAST software implemented in FSL(Smith et al., 2004; Zhang et al., 2001). Time series were extracted and projected out of the data via linear regression.

To examine negative and positive correlations with the DMN at the group level, individual subjects' correlation maps were normalized into a common space (the SPM EPI template). To perform group-level statistical analyses, individual subjects' correlation maps were first converted into linear regression weights (β = r · σyx, where σy is the standard deviation of the voxel time series and σx is the standard deviation of the seed ROI time series), and were subsequently normalized to the SPM5 EPI template, spatially resampled at a voxel size of 2 mm isotropic, and entered into a group-level random-effects analysis using SPM5. A nuisance factor of “scanner” was included since subjects had been scanned on 1 of 2 different scanners; it was furthermore verified that separate analyses of the 2 within-scanner groups of subjects yielded results consistent with the analysis of the entire group. Group-level statistical maps appear superimposed on the ch2 template brain from the MRIcro software (http://www.mricro.com), and anatomical regions/Brodmann areas were determined using MRIcro built-in atlases. All maps (individual and group-level) are shown in radiological convention.

An ROI analysis was used to quantify region-specific changes in DMN connectivity. ROIs were defined by taking the intersection of (1) the group-level random-effects map of positive correlation with the precuneus/PCC ROI (p<0.001) when no corrections were applied to the data, and (2) the working memory deactivation contrast (p<0.05). Six separate ROIs were defined by intersecting the resulting map with six anatomically-defined (AAL) ROIs (anterior cingulate cortex, bilateral angular gyri, dorso-medial prefrontal cortex (DMPFC), and bilateral parahippocampal gyri (PHG)) using the MarsBar software.

Results

Motion

Across subjects, the peak and RMS excursion in head motion were 0.67±0.42 mm and 0.10±0.07mm, respectively (mean ± SD).

Comparison of physiological noise regressors, white/CSF signals, and the global signal

Correlation coefficients between the global signal and both respiratory and cardiac RVHRCOR regressors are shown in Table 1. Correlations with the global signal ranged from -0.10 to 0.69 for RVx, and from -0.15 to 0.49 for HRx. The white matter and CSF ROIs also demonstrated significant correlations with the global signal (Table 1), ranging from -0.02 to 0.57 for white matter and -0.36 to 0.46 for CSF. Further analysis revealed that both white matter and CSF signals tended to correlate more strongly with HRx than RVx: across subjects, white matter was correlated with HRx at r = 0.24±0.17, and with RVx at r = 0.00±0.16; CSF was correlated with HRx at r = 0.17±0.21, and with RVx at 0.01±0.15 (mean ± SD).

Table 1
Correlations between the global signal and RVHRCOR physiological regressors, white matter ROI signal, and CSF ROI signal

Negative correlations with the DMN

Negative correlations (r<-0.15) with the precuneus/PCC seed region are shown for 3 different subjects in Fig. 1 (blue color scale). The spatial extent of negatively-correlated voxels before correction and after physiological noise correction, white/CSF signal removal, and global signal removal are quantified for all subjects in Fig. 2A.

Figure 1
Positive (red) and negative (blue) correlations with the precuneus/PCC ROI, thresholded at |r|>0.15. Maps for 3 different subjects (A, B, and C) depict correlation coefficients before correction (“none”), after physiological noise ...
Figure 2
Spatial extent of (A) negative (r<-0.15) and (B) positive (r>0.15) correlations for each subject before correction (blue), after white/CSF signal removal (cyan), after physiological noise correction (yellow), and after global signal removal ...

Negative correlations were present within individual subjects even without physiological or global signal removal, and Fig. 1 indicates that they were primarily localized to regions of the TPN. In 11 out of 15 subjects (including the subjects in A and C of Fig. 1), physiological noise correction increased the extent of negative correlations, though to a much lesser degree than global signal removal. Across subjects, the percent change in the number of negatively correlated (r<-0.15) voxels was 28.5% ± 9.24% after physiological noise correction, and 382% ± 106% after global signal removal (mean ± SE). White/CSF signal removal also slightly increased the extent of negative correlations, though to a lesser degree than physiological noise correction (6.9% ± 4.0%),.

To visualize the anatomical locations of negative correlations across subjects, a map was constructed that depicts – for each voxel – the number of subjects for which that voxel was negatively correlated at r<-0.15. Separate maps were constructed for each condition (no correction, physiological noise correction, white/CSF signal removal, and global signal removal); see Fig. 3. Figure 3 shows that physiological noise correction not only increased the extent of negative correlations, but tended to increase them in brain areas that were common across subjects. Anatomical regions in which at least 4 subjects exhibited suprathreshold negative correlations before and after physiological noise correction are presented in Table 2. Global signal removal greatly increased the negative correlations over more widespread regions of the brain, as well as inter-subject variability (Figs. 1--3);3); however, it is apparent that regions with the greatest consistency across subjects after global signal removal were aligned with those having the greatest consistency after physiological noise correction.

Figure 3
Anatomical locations where individual subjects had negative correlations with the precuneus/PCC ROI. The value at each pixel depicts the number of subjects that had a negative correlation (r<-0.15) at that pixel before any corrections were performed ...
Table 2
Negative correlations with precuneus/PCC in individual subjects before and after physiological noise correction

Results of a group-level random-effects analysis of negative correlations with the precuneus/PCC are shown in Fig. 4, and clusters attaining a threshold of p<0.05 (uncorrected) are reported in Table 3. Before physiological noise correction, small clusters were present in bilateral inferior parietal cortices (supramarginal gyri; BA 40) and right DLPFC (BA 46). After physiological noise correction, cluster sizes increased and included the bilateral inferior parietal cortices (BA 40), bilateral insula/frontal operculum, and right DLPFC (BA 45/46). No regions, either before or after physiological noise correction, were significant at a corrected threshold of FDR=0.05. For comparison, results after global signal removal and white/CSF signal removal are also presented in Fig. 4. After global signal removal, negatively-correlated regions included bilateral superior and inferior parietal cortices (BA 7/40), bilateral middle/inferior temporal cortices (BA 37), bilateral DLPFC and ventrolateral preforntal cortex (BA 9/45/46/47), mid cingulate/SMA (BA 24/32), and premotor cortex (BA 6). White/CSF signal removal increased the cluster sizes of regions in the uncorrected (“none”) map, and introduced an additional cluster in the right insula, but the increases were smaller than those obtained with physiological correction.

Figure 4
Regions having negative correlation (p<0.05 uncorrected) with the precuneus/PCC in a group-level (N=15), random-effects analysis. Results are shown before corrections were applied (“none”), after white/CSF signal removal (“wcsf”), ...
Table 3
Group-level negative correlations with the precuneus/PCC before and after physiological noise correction

After physiological noise correction, an increase in the (negative) value of the correlation coefficient between the seed ROI and a given voxel would arise from a reduction in the noise variance of the voxel and/or the seed time series, as well as from an increase in the covariance between a voxel and the seed. To further understand the sources of increased negative correlation, the percent increase in the magnitude of the numerator (covariance between the voxel and seed time series) and percent decrease in the denominator (product of the voxel and seed temporal standard deviations) of the correlation coefficient were separately examined, across all voxels demonstrating increased negative correlation after noise correction. From a linearization analysis, it can be seen that for small changes, the fractional change in a correlation coefficient can be approximated by the sum of the fractional changes of its numerator and denominator (Appendix B). It was observed that in 9 out of the 15 subjects, the percent increase in the numerator exceeded the percent decrease in denominator within the majority of voxels queried, while the reverse was true for the remaining 6 subjects. Hence, it appears that both the increased covariance, as well as decreased noise variance, contributed to the increased negative correlations obtained after physiological noise correction.

Positive correlations within the DMN

Positive correlations with the precuneus/PCC are shown for 3 different subjects in Fig. 1 (red color scale). Results of the group-level random-effects analysis of positive correlations are shown in Fig. 5, and the extent of positively-correlated voxels is quantified in Fig. 2B. The spatial extent of significant positive correlations is diminished after physiological noise correction is applied, and even further diminished after global signal removal. White/CSF signal removal also introduced slight decreases in the extent of positive correlations. Across subjects, the percent change in the number of positively correlated (r>0.15) voxels was -8.5% ± 1.1% after white/CSF signal removal, -19.1% ± 0.53% after physiological noise correction, -64.0% ± 1.0% after global signal removal (mean ± SE).

Figure 5
Group analysis (N=15) of regions having significant positive correlations with the precuneus/PCC before correction (“none”), after white/CSF signal removal (“wcsf”), after physiological noise correction (“phys”), ...

Figure 6 shows the mean T-scores from 6 default-mode ROIs, computed by averaging over voxels in the group-level T-maps for each of the 4 conditions (Fig. 5). The magnitude of positive correlations decreased after physiological noise correction, but decreased even more after global signal removal. Voxel-wise reductions in T-score due to physiological noise correction, however, were significant only in the PHG, while reductions due to global signal removal were significant in all 6 ROIs (p<0.01, paired t-test).

Figure 6
ROI analysis of DMN group-level statistics. For each of the 4 statistical maps of the DMN (shown in Figure 5), average T-scores were computed across voxels within 6 core regions of the DMN. ACC, anterior cingulate cortex; DMPFC, dorso-medial prefrontal ...

Magnitude reductions varied across regions; for example, the left angular gyrus and DMPFC experienced only minor reductions due to physiological noise correction (ΔT = 0.1 and 0.02, respectively), whereas decrements in the PHG were more pronounced (ΔT =1.0 for the left hemisphere, and 1.3 for the right). Global signal reduction also had differential effects across ROIs, with average T-score decrements of 1.0 in the left angular gyrus and 3.7 in the left PHG.

Discussion

In the present study, negative correlations between the DMN and multiple key regions of the TPN were observed in the resting state, despite the fact that no global signal removal or scaling steps were performed in pre-processing. While regions of negative correlation were observed in some individual subjects prior to any noise correction, the use of model-based corrections for respiratory and cardiac noise substantially increased their magnitude and extent. The present study indicates that DMN and TPN regions may exhibit, on average, opposing (or nearly 180° out of phase) fluctuations in the absence of external stimuli.

Physiological noise correction enhanced the spatial extent and inter-subject consistency of negative BOLD signal correlations with the DMN. This suggests that cardiac and respiratory processes can obscure the underlying dynamics between brain regions, and illustrates the efficacy of procedures applied herein to remove such noise. The significant correspondence between the global average brain signal and both the RVx and HRx regressors (Table 1) indicates that much of the common signal expressed by voxels is due to these physiological processes, and argues that their removal should become a standard part of fMRI analyses.

The overlap with the global signal was not complete, however, indicating that the global signal contains sources that are not modeled by RETROICOR and RVHRCOR. Some of these sources could be of neural origin, such as regions involved in resting state networks that exhibit coherent neuronal fluctuations over widespread regions of the brain. In addition, noise signals from voxels with high temporal standard deviation due to motion or pulsation, such as the eye regions, edges of the brain, and CSF, are also weighted strongly in the global signal. Lastly, unmodeled physiological noise, such as autoregulatory processes and nonlinear interactions between respiratory and cardiac effects, may also contribute to noise in the brain and possibly the global signal.

After correction for physiological noise, negative correlations at the group level were observed in the DLPFC and bilaterally in the supramarginal gyri and insula/frontal operculum, areas that are commonly recruited in tasks involving attentional control and working memory (Corbetta and Shulman, 2002; Dosenbach et al., 2006). The fronto-insular region has also been implicated in emotional arousal (Critchley et al., 2002) and interoceptive awareness (Critchley et al., 2004). Additional regions of negative correlation, observed in at least 4 out of 15 individual subjects, included bilateral middle temporal cortices, premotor cortex, SMA, and a greater amount of prefrontal cortex (extending into BA 9 and BA 47)). The existence of negative correlations between these task-positive regions and the DMN may support existing theories that resting-state brain signals, on the level of large-scale networks, reflect shifts in state between an introspective, self-reflective mode and an outwardly-oriented, attentive mode (Buckner et al., 2008; Fox et al., 2005; Fransson, 2005).

Regions having negative correlation after physiological noise correction overlapped substantially with those having the greatest negative correlation magnitudes after global signal removal, both here and in previous studies that applied global signal removal (e.g. (Fox et al., 2005; Fransson, 2005)). Although a number of negatively-correlated voxels produced by global signal removal are artificial (Fox et al., 2008; Murphy et al., 2009), it is possible that regions that emerge consistently across subjects (e.g. the brightest regions in Fig. 3, bottom row) are indeed real, especially when they are significant after only physiological noise correction. Global signal removal may thus be useful for identifying candidate regions that are anti-correlated with the DMN; however, one must be careful about interpreting the precise magnitudes of these negative correlation coefficients, since (1) global signal regression is known to skew the distribution of correlation coefficients in ROI-based functional connectivity toward zero-mean (Murphy et al., 2009), and (2) the global signal, while containing some respiratory and cardiac components (Table 3), is not based in hypothesis and – since it is derived from the data itself – may contain important neural signals in addition to noise (see (Murphy et al., 2009) for a thorough discussion of the impact of global signal removal on fMRI time series).

One recent study (Murphy et al., 2009) also performed physiological noise corrections, but negative correlations with the DMN in task-positive regions could not be found at any threshold. One potential reason for this disparity is the difference in physiological correction algorithms: although RETROICOR was applied, just as in the present study, correction for low-frequency heart rate effects was not performed and correction for respiratory variations was accomplished by computing a respiratory volume per time (RVT) signal (defined as the breath-to-breath volume divided by breath-to-breath time interval) and regressing it out of the fMRI data at time shifts ranging from -10 to 40 s. In contrast, the present study convolved the respiratory variation time series with an impulse response (respiration response function (Birn et al., 2008))) before regressing it out of the fMRI time series, and used a sliding-window measure of breathing volume (RV) rather than a breath-to-breath measure. While RV and RVT tend to be correlated during normal breathing (Chang et al., 2009), the convolution step is likely to introduce differences since the respiration response function employed here does not resemble a delta function.

However, given that the present study found negative correlations with the DMN within single subjects even before corrections were applied (e.g. Fig. 1), it is possible that other factors, such as image acquisition methods, were critical for obtaining the present results. While studies failing to observe negative correlations have used conventional echo-planar sequences (Golland et al., 2007; Murphy et al., 2009), the current study used a spiral-in/out sequence that has been shown to improve the signal-to-noise ratio (SNR) in uniform brain regions as well as to reduce susceptibility dropout in frontal and parietal regions (Glover and Law, 2001). In addition, a navigator correction was performed to correct for breathing-induced modulation of the magnetic field (Pfeuffer et al., 2002). The additional SNR provided by these techniques may have enabled the observation of anti-correlations with the DMN in the absence of any noise correction.

Most previous ROI-based functional connectivity studies have applied low-pass temporal filtering as a pre-processing step to increase the contrast-to-noise ratio (e.g. (Fox et al., 2005; Fransson, 2005; Murphy et al., 2009)), since prior work has indicated that resting-state correlations are based primarily in the low-frequency (<0.1 Hz) spectrum (Cordes et al., 2001). While no temporal filtering was originally performed in the present study, the use of the full frequency bandwidth may account for some of the observed differences with previous work – particularly, the presence of negative correlations. To understand the impact of low-pass filtering, data were re-analyzed after applying a low-pass filter (<0.1 Hz) prior to calculating positive and negative correlations with the PCC. Results after low-pass filtering were similar to those of the unfiltered data; in fact, low-pass filtering served to slightly increase the extent and magnitude of negative correlations in both the uncorrected and physiologically-corrected data (Fig. 7a), and decrease the extent and magnitude of positive correlations (Fig. 7b). Thus, low-pass filtering is not believed to be a source of difference with previous studies, and appears to exert an influence in the same direction as physiological noise correction.

Figure 7
Effect of temporal low-pass filtering on a group-level, random-effects analysis of (A) negative and (B) positive correlations with the PCC. Maps are thresholded identically to Figs. 4 and and5:5: negative correlations are shown at 1.77<T<6, ...

While physiological noise correction increased the extent of regions negatively correlated with the precuneus/PCC, the extent of positive correlations was reduced. This is consistent with previous accounts (Birn et al., 2006; Chang et al., 2009), and with observations that physiological noise exerts a nearly global influence on the BOLD signal that can artificially couple widespread regions of the brain (Birn et al., 2006; Birn et al., 2008; Chang et al., 2009; de Munck et al., 2008; Glover et al., 2000; Shmueli et al., 2007; Wise et al., 2004). While the overall extent was reduced, Fig. 6 indicates that mean T-scores in all queried ROIs except the PHG were not significantly reduced by physiological noise correction. On the other hand, global signal removal significantly diminished DMN correlations within these areas, which further cautions against its use.

Regression of white matter and CSF ROIs also elicited increases in negative correlations with the PCC (Figs. 3 and and4),4), as well as reductions in positive correlations (Fig. 5), though to a lesser extent on average than physiological noise correction. Differences between the two correction methods might arise because white matter and CSF signals may not adequately capture respiratory-related signals, since previous work has shown that respiration-induced modulation of arterial CO2 primarily affects gray matter (Birn et al., 2006; Birn et al., 2008; Chang et al., 2009; Glover et al., 2000), and the present study also observed a lack of correlation between the white/CSF signals and the respiratory variations regressor (RVx; data not shown). In addition, it is possible that varying the size and location of white/CSF ROIs can alter the outcome of the correction; future work will be necessary to resolve this question.

By contrasting Figs. 4 and and5,5, it is evident that negative correlations with the PCC were weaker and less consistent across subjects than positive correlations. One possible explanation is that residual physiological noise, as well as other unmodeled, common noise sources, skews correlation coefficients in the positive direction. Another possibility is that the anti-correlated behavior between the DMN and TPN networks is a transient, rather than persistent, phenomenon that is dependent upon particular states of cognition and/or arousal. Figure 8 depicts, for one subject, time series averaged across regions positively (red) and negatively (blue) correlated with the PCC ROI at |r|>0.15. These traces indicate that phase differences between the networks may vary over time. Further work will be necessary to illuminate and understand the dynamic behavior of these networks.

Figure 8
Mean time series of voxels positively (r>0.15; red line) and negatively (r<-0.15; blue line) correlated with the precuneus/PCC ROI, for one subject.

Conclusions

The current findings indicate that negative correlations between DMN and task-positive regions can be observed in the absence of global signal removal. A group-level random effects analysis revealed clusters at an uncorrected (p<0.05) threshold (though none were significant at FDR=0.05). Physiological noise correction, based on respiratory and cardiac monitoring in conjunction with a priori models, both increased the extent of negative correlations and decreased the extent of positive correlations. Regression of white matter and CSF ROI signals also increased and decreased the extent of negative and positive correlations, respectively, though to a lesser degree than physiological noise correction. Results suggest that modeling physiological noise will be crucial for accurately quantifying dynamic relationships within and across brain networks.

Acknowledgments

The authors gratefully acknowledge support from NIH grants F31-AG032168 (CC) and P41-RR09784 (GHG), and thank Moriah Thomason for helpful consultation regarding brain anatomy. We are also grateful to two anonymous reviewers for insightful comments.

Appendix A: RVHRCOR

Let xr and xh denote the respiratory variations and heart rate time series, respectively, and hr and hh denote the respiration response function (Birn et al., 2008) and cardiac response function (Chang et al., 2009), respectively:

equation M1
equation M2

Let zr [equivalent] xr [multiply sign in circle] hr and zh [equivalent] xh [multiply sign in circle] hh, where “[multiply sign in circle]” indicates convolution. The convolved signals, zr and zh (referred to as RVx and HRx in the text), represent the modeled BOLD signal response to low-frequency respiratory and cardiac processes. In the model proposed in (Chang et al., 2009), these 2 sources of physiological noise contribute additively to the time series of a voxel (y):

equation M3

where ε contains the desired (neural) component of the time series as well as unmodeled noise, and is assumed to be small compared to the physiological noise components. The coefficients βr and βh can be estimated using least squares:

equation M4

where Z = [zr, zh], and the corrected voxel time series is taken to be the residual:

equation M5

Appendix B

Let equation M6, where n = σxy, d = σxσy, and r is the correlation coefficient between the time series of a seed region (x) and a voxel (y). Let Δn denote the increase in the numerator, and let Δd denote the decrease in the denominator. Then equation M7, and therefore equation M8.

References

  • Beckmann CF, DeLuca M, Devlin JT, Smith SM. Investigations into resting-state connectivity using independent component analysis. Philos Trans R Soc Lond B Biol Sci. 2005;360:1001–1013. [PMC free article] [PubMed]
  • Behzadi Y, Restom K, Liau J, Liu TT. A component based noise correction method (CompCor) for BOLD and perfusion based fMRI. Neuroimage. 2007;37:90–101. [PMC free article] [PubMed]
  • Birn RM, Diamond JB, Smith MA, Bandettini PA. Separating respiratory-variation-related fluctuations from neuronal-activity-related fluctuations in fMRI. Neuroimage. 2006;31:1536–1548. [PubMed]
  • Birn RM, Smith MA, Jones TB, Bandettini PA. The respiration response function: the temporal dynamics of fMRI signal fluctuations related to changes in respiration. Neuroimage. 2008;40:644–654. [PMC free article] [PubMed]
  • Biswal B, Bandettini PA, Jesmanowicz A, Hyde JS. Time-frequency analysis of functional EPI time-course series. Proc., SMRM, 12th Annual Meeting; New York. 1992; 1993. p. 722.
  • Biswal B, Yetkin FZ, Haughton VM, Hyde JS. Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magn Reson Med. 1995;34:537–541. [PubMed]
  • Buckner RL, Andrews-Hanna JR, Schacter DL. The brain's default network: anatomy, function, and relevance to disease. Ann N Y Acad Sci. 2008;1124:1–38. [PubMed]
  • Buckner RL, Vincent JL. Unrest at rest: default activity and spontaneous network correlations. Neuroimage. 2007;37:1091–1096. discussion 1097-1099. [PubMed]
  • Chang C, Cunningham JP, Glover GH. Influence of heart rate on the BOLD signal: the cardiac response function. Neuroimage. 2009;44:857–869. [PMC free article] [PubMed]
  • Clare Kelly AM, Uddin LQ, Biswal BB, Castellanos FX, Milham MP. Competition between functional brain networks mediates behavioral variability. Neuroimage. 2008;39:527–537. [PubMed]
  • Corbetta M, Shulman GL. Control of goal-directed and stimulus-driven attention in the brain. Nat Rev Neurosci. 2002;3:201–215. [PubMed]
  • Cordes D, Haughton VM, Arfanakis K, Carew JD, Turski PA, Moritz CH, Quigley MA, Meyerand ME. Frequencies contributing to functional connectivity in the cerebral cortex in “resting-state” data. AJNR Am J Neuroradiol. 2001;22:1326–1333. [PubMed]
  • Critchley HD, Mathias CJ, Dolan RJ. Fear conditioning in humans: the influence of awareness and autonomic arousal on functional neuroanatomy. Neuron. 2002;33:653–663. [PubMed]
  • Critchley HD, Wiens S, Rotshtein P, Ohman A, Dolan RJ. Neural systems supporting interoceptive awareness. Nat Neurosci. 2004;7:189–195. [PubMed]
  • Dagli MS, Ingeholm JE, Haxby JV. Localization of cardiac-induced signal change in fMRI. Neuroimage. 1999;9:407–415. [PubMed]
  • Damoiseaux JS, Rombouts SA, Barkhof F, Scheltens P, Stam CJ, Smith SM, Beckmann CF. Consistent resting-state networks across healthy subjects. Proc Natl Acad Sci U S A. 2006;103:13848–13853. [PubMed]
  • De Luca M, Beckmann CF, De Stefano N, Matthews PM, Smith SM. fMRI resting state networks define distinct modes of long-distance interactions in the human brain. Neuroimage. 2006;29:1359–1367. [PubMed]
  • de Munck JC, Goncalves SI, Faes TJ, Kuijer JP, Pouwels PJ, Heethaar RM, Lopes da Silva FH. A study of the brain's resting state based on alpha band power, heart rate and fMRI. Neuroimage. 2008;42:112–121. [PubMed]
  • Dosenbach NU, Visscher KM, Palmer ED, Miezin FM, Wenger KK, Kang HC, Burgund ED, Grimes AL, Schlaggar BL, Petersen SE. A core system for the implementation of task sets. Neuron. 2006;50:799–812. [PMC free article] [PubMed]
  • Fair DA, Cohen AL, Dosenbach NU, Church JA, Miezin FM, Barch DM, Raichle ME, Petersen SE, Schlaggar BL. The maturing architecture of the brain's default network. Proc Natl Acad Sci U S A. 2008;105:4028–4032. [PubMed]
  • Fox MD, Raichle ME. Spontaneous fluctuations in brain activity observed with functional magnetic resonance imaging. Nat Rev Neurosci. 2007;8:700–711. [PubMed]
  • Fox MD, Snyder AZ, Raichle ME. Global Signal Regression and Anticorrelations in Resting State fMRI Data. Proc., Organization for Human Brain Mapping, 14th Annual Meeting; Melbourne. 2008; 2008. p. 575 W-AM.
  • Fox MD, Snyder AZ, Vincent JL, Corbetta M, Van Essen DC, Raichle ME. The human brain is intrinsically organized into dynamic, anticorrelated functional networks. Proc Natl Acad Sci U S A. 2005;102:9673–9678. [PubMed]
  • Fransson P. Spontaneous low-frequency BOLD signal fluctuations: an fMRI investigation of the resting-state default mode of brain function hypothesis. Hum Brain Mapp. 2005;26:15–29. [PubMed]
  • Fransson P. How default is the default mode of brain function? Further evidence from intrinsic BOLD signal fluctuations. Neuropsychologia. 2006;44:2836–2845. [PubMed]
  • Friston KJ, Williams S, Howard R, Frackowiak RS, Turner R. Movement-related effects in fMRI time-series. Magn Reson Med. 1996;35:346–355. [PubMed]
  • Glover GH. Deconvolution of impulse response in event-related BOLD fMRI. Neuroimage. 1999;9:416–429. [PubMed]
  • Glover GH, Lai S. Self-navigated spiral fMRI: interleaved versus single-shot. Magnetic Resonance in Medicine. 1998;39:361–368. [PubMed]
  • Glover GH, Law CS. Spiral-in/out BOLD fMRI for increased SNR and reduced susceptibility artifacts. Magn Reson Med. 2001;46:515–522. [PubMed]
  • Glover GH, Li TQ, Ress D. Image-based method for retrospective correction of physiological motion effects in fMRI: RETROICOR. Magn Reson Med. 2000;44:162–167. [PubMed]
  • Golland Y, Bentin S, Gelbard H, Benjamini Y, Heller R, Nir Y, Hasson U, Malach R. Extrinsic and intrinsic systems in the posterior cortex of the human brain revealed during natural sensory stimulation. Cereb Cortex. 2007;17:766–777. [PubMed]
  • Greicius MD, Krasnow B, Reiss AL, Menon V. Functional connectivity in the resting brain: a network analysis of the default mode hypothesis. Proc Natl Acad Sci U S A. 2003;100:253–258. [PubMed]
  • Gusnard DA, Raichle ME, Raichle ME. Searching for a baseline: functional imaging and the resting human brain. Nat Rev Neurosci. 2001;2:685–694. [PubMed]
  • Jezzard P, LeBihan D, Cuenod D, Pannier L, Prinster A, Turner R. An investigation of the contributiin of physiological noise in human functional MRI studies at 1.5 tesla and 4 tesla. Proc., SMRM, 12th Annual Meeting; New York. 1992; 1993. p. 1392.
  • Kim DH, Adalsteinsson E, Glover GH, Spielman DM. Regularized higher-order in vivo shimming. Magn Reson Med. 2002;48:715–722. [PubMed]
  • Mazoyer B, Zago L, Mellet E, Bricogne S, Etard O, Houde O, Crivello F, Joliot M, Petit L, Tzourio-Mazoyer N. Cortical networks for working memory and executive functions sustain the conscious resting state in man. Brain Res Bull. 2001;54:287–298. [PubMed]
  • Murphy K, Birn RM, Handwerker DA, Jones TB, Bandettini PA. The impact of global signal regression on resting state correlations: are anti-correlated networks introduced? Neuroimage. 2009;44:893–905. [PMC free article] [PubMed]
  • Pfeuffer J, Van de Moortele PF, Ugurbil K, Hu X, Glover GH. Correction of physiologically induced global off-resonance effects in dynamic echo-planar and spiral functional imaging. Magn Reson Med. 2002;47:344–353. [PubMed]
  • Raichle ME, MacLeod AM, Snyder AZ, Powers WJ, Gusnard DA, Shulman GL. A default mode of brain function. Proc Natl Acad Sci U S A. 2001;98:676–682. [PubMed]
  • Shmueli K, van Gelderen P, de Zwart JA, Horovitz SG, Fukunaga M, Jansma JM, Duyn JH. Low-frequency fluctuations in the cardiac rate as a source of variance in the resting-state fMRI BOLD signal. Neuroimage. 2007;38:306–320. [PMC free article] [PubMed]
  • Shulman GL, Fiez JA, Corbetta M, Buckner RL, Miezin FM, Raichle ME, Petersen SE. Common blood flow changes across visual tasks. II. Decreases in cerebral cortex. J Cogn Neurosci. 1997;9:648–663. [PubMed]
  • Smith SM, Jenkinson M, Woolrich MW, Beckmann CF, Behrens TE, Johansen-Berg H, Bannister PR, De Luca M, Drobnjak I, Flitney DE, Niazy RK, Saunders J, Vickers J, Zhang Y, De Stefano N, Brady JM, Matthews PM. Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage. 2004;231 1:S208–219. [PubMed]
  • Sternberg S. Memory-scanning: mental processes revealed by reaction-time experiments. Am Sci. 1969;57:421–457. [PubMed]
  • Thomason ME, Glover GH. Controlled inspiration depth reduces variance in breath-holding-induced BOLD signal. Neuroimage. 2008;39:206–214. [PMC free article] [PubMed]
  • Tian L, Jiang T, Liu Y, Yu C, Wang K, Zhou Y, Song M, Li K. The relationship within and between the extrinsic and intrinsic systems indicated by resting state correlational patterns of sensory cortices. Neuroimage. 2007;36:684–690. [PubMed]
  • Uddin LQ, Clare Kelly AM, Biswal BB, Xavier Castellanos F, Milham MP. Functional connectivity of default mode network components: Correlation, anticorrelation, and causality. Hum Brain Mapp. 2009;30:625–637. [PMC free article] [PubMed]
  • Wang K, Liang M, Wang L, Tian L, Zhang X, Li K, Jiang T. Altered functional connectivity in early Alzheimer's disease: a resting-state fMRI study. Hum Brain Mapp. 2007;28:967–978. [PubMed]
  • Weisskoff RM, Baker J, Belliveau J, Davis TL, Kwong KK, Cohen MS, Rosen BR. Power spectrum analysis of functionally-weighted MR data: what's in the noise?. Proc., SMRM, 12th Annual Meeting; New York. 1992; 1993. p. 7.
  • Wise RG, Ide K, Poulin MJ, Tracey I. Resting fluctuations in arterial carbon dioxide induce significant low frequency variations in BOLD signal. Neuroimage. 2004;21:1652–1664. [PubMed]
  • Zhang Y, Brady M, Smith S. Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm. IEEE Trans Med Imaging. 2001;20:45–57. [PubMed]