Low-frequency fluctuations in the cardiac rate as a source of variance in the resting-state fMRI BOLD signal
Introduction
It is important to identify and characterize the sources of physiological noise in the blood-oxygenation-level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signal so that noise reduction strategies can be developed and improved. This is especially relevant as technical advances, such as the use of higher magnetic fields or multiple receiver coils, have increased the available MRI signal. Physiological noise depends on the total signal strength and therefore becomes a larger fraction of the total noise as the signal increases (Kruger and Glover, 2001, Kruger et al., 2001, Triantafyllou et al., 2005), for example with increased magnetic field strength.
Physiological noise negatively affects conventional task-based fMRI experiments and has been found to be particularly problematic in the absence of external stimuli, in resting-state functional connectivity analysis (Cordes et al., 2001). In this type of analysis neuronal connections between brain regions are inferred by measuring the spatial correlations of low-frequency (< 0.1 Hz) BOLD fMRI signal fluctuations in the brain. The minimization of physiological noise as a confound in such studies has been a preoccupation since their inception (Biswal et al., 1996, Fukunaga et al., 2006, Lowe and Sakaie, 2006, Vogt et al., 2006).
Several physiological noise sources have been identified, including effects related to the cardiac (Bhattacharyya and Lowe, 2004, Dagli et al., 1999) and respiratory (Birn et al., 2006b, Wise et al., 2004) cycles. Variations related to arterial carbon dioxide fluctuations (Wise et al., 2004) and residual movement artifacts not accounted for by rigid body registration (Lund et al., 2005) also contribute.
Various methods for reducing physiological noise have been proposed. These may operate in k-space (Hu et al., 1995, Le and Hu, 1996, Wowk et al., 1997) or in image space (Chuang and Chen, 2001, Deckers et al., 2006, Glover et al., 2000) with the latter being the preferred method since changes made in k-space affect all the voxels in the reconstructed images. This makes spatially localized noise difficult to remove and may induce spatial correlations. Fixed bandwidth band-reject filtering was initially proposed (Biswal et al., 1996) to remove noise at the fundamental cardiac and respiratory frequencies but can only be applied successfully when variations at these frequencies are stationary and adequately sampled (as in the rare case of imaging with very short TR) so that they are not aliased to lower frequencies. Moreover, such filtering cannot be applied if task-related signals are present in the rejected frequency band.
Some noise reduction methods involve the acquisition of additional physiological data (Glover et al., 2000, Hu et al., 1995, Liston et al., 2006, Vogt et al., 2006) using a photoplethysmograph and pneumatic belt, for example. Other methods utilize the MRI data itself to estimate the noise (Le and Hu, 1996, Lowe and Sakaie, 2006, Wowk et al., 1997). Some of the methods are designed for straightforward data correction (Glover et al., 2000) but most can be extended to perform ‘nuisance variable regression’ (Birn et al., 2006a, Lund et al., 2006) in which the physiological noise measures (or models derived from them) are included as regressors in a general linear model (GLM) regression analysis. This regression method is advantageous over applying a filter or additional data correction as it does not interfere with the detection of functional activation (Deckers et al., 2006).
Many studies have demonstrated that there are fMRI signal changes associated with the pulsatile cardiac motion (both at the cardiac frequency and its harmonics) (Bhattacharyya and Lowe, 2004, Biswal et al., 1996, Dagli et al., 1999, Kruger and Glover, 2001). Most of the physiological noise reduction methods model the contribution from the cardiac motion by determining the relative timing of each image volume (or slice (Liston et al., 2006, Vogt et al., 2006)) within the cardiac cycle. They assume that the cardiac cycle follows the same inherent pattern regardless of the inter-beat time interval. This assumption is not entirely accurate as there are differences in the time interval between beats that are associated more with time variation in the diastolic portion of the cardiac cycle than the systolic (Dagli et al., 1999).
Across-beat fluctuations in the heart rate have been found to occur in several frequency bands (Akselrod et al., 1981, Cohen and Taylor, 2002, Otzenberger et al., 1998) including the low-frequency region investigated in fMRI studies of resting-state functional connectivity (< 0.1 Hz). In fact, a recent study suggests that systemic cardiovascular fluctuations (low-frequency oscillations in the heart rate and arterial blood pressure) can account for about half of the information carried with low-frequency oscillations in the cerebral hemodynamics (specifically in the oxyhemoglobin concentration change) (Katura et al., 2006). Fluctuations in the heart rate may therefore be related to those in the BOLD fMRI signal, particularly since the signal is dependent on the cerebral blood flow, oxygenation and volume. Furthermore, since respiration and the cardiac rate are intimately related (Cohen and Taylor, 2002) and the respiration volume per unit time (RVT) has recently been found to be correlated with the resting-state fMRI BOLD signal (Birn et al., 2006b), it is probable that fluctuations in the cardiac rate are related to those in the resting-state fMRI BOLD signal. The effect of cardiac rate fluctuations on the fMRI signal has not been assessed or explicitly accounted for in physiological noise models.
Therefore, the aim of this study was to investigate the hypothesis that across-beat fluctuations in the cardiac rate are correlated with low-frequency variations in the resting-state BOLD fMRI signal. Temporal and spatial patterns of fMRI signal changes related to variations in the cardiac rate were characterized. To this end, resting-state BOLD fMRI data and accompanying physiological data were acquired and analyzed using both cross-correlation and regression. Subsequent to and informed by the cross-correlation analysis, several nested regression analyses were performed to determine whether adding a set of time-shifted cardiac rate timecourses to other established physiological noise regressors would explain additional variance in the data.
Section snippets
Imaging
Resting-state BOLD fMRI data were acquired using a 3-Tesla scanner (GE Signa, Milwaukee, WI, USA) equipped with a 16-channel head coil (Nova Medical, Wakefield, MA) (de Zwart et al., 2004). The imaging sequence, developed in-house, was a single-shot gradient-echo EPI with TE = 43 ms and TR = 6 s with all the slices acquired during the first 3 s of the TR interval, leaving a quiet period in which EEG acquisition could take place (Horovitz et al., in press). In each volume, 28 axial, 3 mm-thick
Results
The mean and standard deviation of the cardiac rate for each session and averaged across all the sessions are given in Table 1.
Fig. 2 illustrates the Fourier transform of the resampled cardiac rate timecourse, averaged over all subjects. The figure shows that most of the variations in the heart rate occur at very low frequencies (< 0.1 Hz) and that there are no distinct frequency peaks in the spectrum of cardiac rate variations. This was found to be the case for all the subjects.
Discussion and conclusions
Both the correlation and regression analyses performed here revealed significant correlations between the cardiac rate and resting-state fMRI signal timecourses, particularly in the GM which showed the strongest negative correlations at time shifts of around 6–12 s and large positive correlations at time shifts of 30–42 s. Regressors consisting of time-shifted measured cardiac rate timecourses explained substantial additional MRI signal variance when included in a model with other established
Acknowledgments
This research was supported by the Intramural Research Program of the National Institute of Neurological Disorders and Stroke, National Institutes of Health.
References (37)
- et al.
Cardiac-induced physiologic noise in tissue is a direct observation of cardiac-induced fluctuations
Magn. Reson. Imaging
(2004) - et al.
Separating respiratory-variation-related fluctuations from neuronal-activity-related fluctuations in fMRI
NeuroImage
(2006) - et al.
Localization of cardiac-induced signal change in fMRI
NeuroImage
(1999) - et al.
An adaptive filter for suppression of cardiac and respiratory noise in MRI time series data
NeuroImage
(2006) - et al.
Large-amplitude, spatially correlated fluctuations in BOLD fMRI signals during extended rest and early sleep stages
Magn. Reson. Imaging
(2006) - et al.
A global optimisation method for robust affine registration of brain images
Med. Image Anal.
(2001) - et al.
Quantitative evaluation of interrelations between spontaneous low-frequency oscillations in cerebral haemodynamics and systemic cardiovascular dynamics
NeuroImage
(2006) - et al.
Modelling cardiac signal as a confound in EEG–fMRI and its application in focal epilepsy studies
NeuroImage
(2006) - et al.
Motion or activity: their role in intra- and inter-subject variation in fMRI
NeuroImage
(2005) - et al.
Non-white noise in fMRI: does modelling have an impact?
NeuroImage
(2006)
Comparison of physiological noise at 1.5 T, 3 T and 7 T and optimization of fMRI acquisition parameters
NeuroImage
Resting fluctuations in arterial carbon dioxide induce significant low frequency variations in BOLD signal
NeuroImage
A general statistical analysis for fMRI data
NeuroImage
Power spectrum analysis of heart rate fluctuation: a quantitative probe of beat-to-beat cardiovascular control
Science
The use of multiple physiologic parameter regression increases gray matter temporal signal to noise by up to 50%
Functional connectivity in the motor cortex of resting human brain using echo-planar MRI
Magn. Reson. Med.
Reduction of physiological fluctuations in fMRI using digital filters
Magn. Reson. Med.
IMPACT: image-based physiological artifacts estimation and correction technique for functional MRI
Magn. Reson. Med.
Cited by (446)
In search of a marker of altered cerebrovascular function in hypertension: Analysis of the fractional amplitude of low-frequency fluctuations in UK Biobank resting state fMRI data
2024, Cerebral Circulation - Cognition and BehaviorWhy is everyone talking about brain state?
2023, Trends in NeurosciencesHead motion and physiological effects
2023, Advances in Resting-State Functional MRI: Methods, Interpretation, and ApplicationsInteroceptive influences on resting-state fMRI
2023, Advances in Resting-State Functional MRI: Methods, Interpretation, and Applications