Reduced Dynamic Coupling Between Spontaneous BOLD-CBF Fluctuations in Older Adults: A Dual-Echo pCASL Study

Introduction

Since the first observation that spontaneous Blood Oxygenation Level Dependent (BOLD) signal fluctuations in the left and right motor cortex are correlated in the absence of a task (Biswal et al., 1995), resting-state fMRI has witnessed an exponential growth of interest.

Since it does not require any task, resting-state fMRI is particularly attractive for studies on patients, children and elderlies, reducing problems related to participant’s compliance or intersubject variability due to task performance. Indeed, there is an increasing number of investigations using resting-state fMRI as a sensitive biomarker to study normal and pathological aging (D’Esposito et al., 1999; Fox and Greicius, 2010; Brier et al., 2014; Gardini et al., 2014; Li et al., 2015; Vecchio et al., 2015; Esposito et al., 2018).

Most fMRI studies are based on the BOLD technique (Bandettini et al., 1992; Kwong et al., 1992; Ogawa et al., 1992) which offers a large sensitivity and easy of implementation. However, due to the complex nature of the BOLD effect, the quantitative interpretation of this fMRI signal can be problematic. Indeed, the BOLD signal change is modulated by local variations in deoxyhemoglobin content stemming from changes in cerebral blood flow (CBF), cerebral blood volume (CBV) and cerebral metabolic rate of oxygen consumption (CMRO₂) induced by neuronal activation or vascular challenges (Buxton, 2012). Depending on the complex interplay of these variables, the magnitude and dynamics of the BOLD signal change may not always reflect the underlying change in neuronal activity or the variation of a specific hemodynamic/metabolic quantity (Ances et al., 2008; Buxton, 2010; Griffeth and Buxton, 2011; Moradi et al., 2012). This limitation is especially important in studies comparing populations with different neurovascular properties (Liu, 2013). In particular, this issue has been largely recognized when comparing elderly vs. young adults since significant vascular changes are known to occur during adult life (O’Rourke and Hashimoto, 2007; Samanez-Larkin and D’Esposito, 2008; Chen et al., 2011; Lu et al., 2011; Gauthier and Hoge, 2013; De Vis et al., 2015).

To overcome the problem of BOLD signal ambiguity, the acquisition of concurrent BOLD and CBF data from arterial spin labeling (ASL) sequences has been proposed in an early biophysical model that allows the calculation of fractional changes of the different hemodynamic and metabolic variables involved in brain activation (Davis et al., 1998). ASL uses magnetically labeled arterial blood water as an endogenous tracer (Detre et al., 1992) and can quantify both baseline levels and activity induced variation of regional CBF. The Davis model is currently applied with different strategies, usually requiring additional calibration measurements based on gas challenges (Davis et al., 1998; Hoge et al., 1999; Chiarelli et al., 2007) and has been expanded to allow the measurement of absolute CMRO₂ changes or baseline levels of oxygen consumption (Bulte et al., 2012; Gauthier et al., 2012; Wise et al., 2013; Germuska and Bulte, 2014; Germuska et al., 2016). Although gas challenges might be less tolerated by patients or elderlies, applications of these techniques to aging studies have been recently reported (Mohtasib et al., 2012; Hutchison et al., 2013; Liu et al., 2013; De Vis et al., 2015; Garrett et al., 2017). Alternative calibration techniques without gas administration have also been proposed based on particular MRI acquisition sequences or breath-hold tasks (Kastrup et al., 1999; Fujita et al., 2006; Bulte et al., 2009; Blockley et al., 2012, 2015). Nevertheless, the calibrated BOLD technique has been mostly applied and validated in studies using block paradigms, and its extension to more dynamic experimental designs is not straightforward (Kida et al., 2006; Simon and Buxton, 2015).

In this regard, however, the combined acquisition of BOLD and CBF data using ASL has recently attracted increasing interest to study brain function, even without calibration measurements (Chen et al., 2015; Simon and Buxton, 2015; Storti et al., 2017).

Indeed, although CBF is still an indirect measurement of cerebral metabolism and neuronal activity, it constitutes a well defined and fundamental physiological process that is altered in different pathologies and with physiological aging. Furthermore, despite ASL has lower sensitivity compared to BOLD, it can extend the study of resting-state brain function beyond that of functional connectivity, allowing quantitative CBF measurements and the investigation of BOLD-CBF coupling when the two signals are acquired simultaneously. In particular, the possibility to capture spontaneous fluctuations of cerebral blood flow with ASL has recently received an increasing attention (De Luca et al., 2006; Chuang et al., 2008; Fukunaga et al., 2008; Zou et al., 2009; Viviani et al., 2011; Liang et al., 2014; Tak et al., 2014; Chen et al., 2015; Fernández-Seara et al., 2015; Jann et al., 2015).

This is an appealing application of concurrent BOLD and CBF dynamic data acquisition, potentially able to offer insights on mechanisms underlying resting-state brain functioning and physiology. In this regard, recent evidence demonstrated that the resting-state dynamic relationship between BOLD and CBF is approximately linear across the brain (Fukunaga et al., 2008; Wu et al., 2009), with a significantly stronger coupling between spontaneous BOLD and CBF fluctuations within the major nodes of established resting-state networks (Tak et al., 2014; Cohen et al., 2017). Noteworthy, the study of BOLD-CBF dynamic coupling during resting-state could offer innovative metrics to assess brain health (Chen et al., 2015), in addition to the more commonly used functional connectivity metrics. However, to the best of our knowledge, no investigation based on this method has been performed on aging so far.

In the present study we investigated how the dynamic BOLD-CBF coupling during resting-state is affected by aging. Specifically, we compared a group of healthy elderlies with a group of young subjects, addressing between-group differences in: (i) the linear correlation between the two signals; and (ii) the effect on the calculated correlation when introducing a relative temporal shift between BOLD and CBF timecourses.

Materials and Methods

Fifteen healthy young adults (Young: mean age = 26.4, SD = 4.2) and 17 healthy elderlies (Elderly: mean age = 63.4, SD = 8.1) were included in the study. All individuals were right-handed, gave their written informed consent according to the Declaration of Helsinki (World Medical Association Declaration of Helsinki, 1997) and all procedures were approved by the Ethics Committee for Biomedical Research of the provinces of Chieti and Pescara and the “G. D’Annunzio” University of Chieti and Pescara. None of the participants reported a history of neurological or psychiatric disease, or used psychopharmacological drugs. Subjects with any drug or alcohol abuse within the previous 6 months were also excluded to avoid confounding effects on the fMRI signal. Other exclusion criteria included implanted metals, pregnancy and abnormal findings in their structural brain MRI. Mild forms of hypertension and hyperlipidemia were accepted. Six older adults used antihypertensive medication.

Elderlies were screened with Mini Mental State Examination (MMSE; Folstein et al., 1975) to evaluate the global cognitive status (reported score range of included subjects: 25.5 ÷ 28.9), Babcock story test to evaluate prose memory (reported score range: 5.2 ÷ 8.8 for immediate recall, 2.5 ÷ 8.6 for delayed recall), and the Frontal Assessment Battery (FAB) to assess global executive functions (reported score range: 14.7 ÷ 18.3). Young group was screened using the Trail Making Test to evaluate sustained visuo-spatial attention (reported score range: 5.5s ÷ 48.8), MMSE (27.7 ÷ 29.3), Babcock story test (7.2 ÷ 8.9 for immediate recall, 5.1 ÷ 8.8 for delayed recall) and FAB (15.1 ÷ 18.9). Statistical between-groups comparisons were performed with parametric or non-parametric tests, depending on normality distribution verified using the Shapiro test. A significant between-group difference was observed for MMSE (p = 0.03, Wilcoxon test), whereas Babcock and FAB tests did not reveal significant effects (p = 0.19 and p = 0.43 respectively, unpaired t-tests).

All participants were required to refrain from caffeine, alcohol and nicotine for at least 6 h before the MRI session.

MRI was performed with a 3T Philips Achieva scanner (Philips Medical Systems, Best, Netherlands), using a whole-body radiofrequency coil for signal excitation and an 8-channel phased-array head coil for signal reception.

Subjects were instructed to keep their eyes closed and not to engage in structured thoughts during acquisition.

Resting-state CBF and BOLD data were simultaneously acquired with a dual-echo pseudo-continuous ASL (pCASL) sequence (Dai et al., 2008) with the following parameters: TR/TE1/TE2: 3500/10/28 ms, FOV 230 mm × 230 mm, matrix 64 × 64, voxel size 3.6 mm × 3.6 mm × 5 mm, SENSE factor 2.3, 19 slices acquired in ascending order, 90 dynamics. The label duration was 1650 ms and the postlabel delay was 1000 ms.

Baseline perfusion was measured using a pCASL sequence optimizing the labeling parameters for a reliable quantification of CBF for both young and elderly subjects (Alsop et al., 2015): TR/TE 4269/10 ms, FOV 230 mm × 230 mm, matrix 64 × 64, voxel size 3.6 mm × 3.6 mm × 5 mm, SENSE factor 2.3, 19 slices acquired in ascending order, 60 dynamics. The label duration was 1750 ms and the postlabel delay was 1900 ms. Background suppression pulses at 2110 ms and 3260 ms after start of labeling were used. The labeling plane of the pCASL sequences was positioned 85 mm below the AC-PC line, according to recent guidelines (Aslan et al., 2010; Alsop et al., 2015). An equilibrium magnetization image (M₀) was also acquired with scan parameters identical to the pCASL sequence (same matrix and readout) but using a long TR (10,000 ms) and without labeling or background suppression pulses.

A high resolution structural volume was finally acquired via a 3D fast field echo T1-weighted sequence with the following parameters: 1 mm isotropic voxel size, TR/TE = 8.1/3.7 ms, flip angle = 8°, 160 sections, SENSE factor = 2.

During fMRI, physiological signals related to respiratory and cardiac cycles were registered using a pneumatic belt strapped around the upper abdomen and a pulse oximeter placed on a finger of the right hand, respectively. Respiratory and cardiac data were both sampled at 100 Hz and stored in a logfile for each run.

Resting state pCASL fMRI data were analyzed using AFNI (Cox et al., 2006¹) and custom-written software implemented in Python². First, the dual-echo pCASL data were split into four EPI timeseries, corresponding to label and control images for acquisitions at TE1 and TE2.

Then, initial preprocessing was performed for both echoes on the label and control ASL images separately (Restom et al., 2006; Wang et al., 2008) according to the following steps:

(i) RETROICOR was applied to remove signal fluctuations related to cardiac and respiratory cycles (Glover et al., 2000);

(ii) slice-timing correction using sinc interpolation;

(iii) motion correction was performed using rigid body registration to realign all time frames to a base image represented by the first label/control volume (dummy scans prevented the need to discard the first few volumes). A summary statistic of motion was defined as the root mean square (RMS) of the six realignment parameters (three translations and three rotations), in order to inspect for EPI timeseries affected by excessive motion. One elderly subject exceeded an RMS value of 1.5 mm for both label and control images and was discarded from further analysis. Further motion assessment was performed considering the RMS value of the differentiated EPI timeseries (DVARS) within a whole brain spatial mask (Power et al., 2012, 2014). The run-averaged DVARS metrics was then compared across groups using ANOVA in order to characterize potential motion effects not accounted for by spatial registration and regression of motion parameters;

(iv) an additional coregistration was performed to minimize effects of the spatial offset between label and control timeseries caused by possible head motion between the reference control image and the reference label image. The corresponding spatial transformation was determined using the mutual information-based approach (Studholme et al., 1998) which is much less sensitive to the labeling related signal intensity differences that could be interpreted by the standard rigid body motion correction algorithm as apparent head motion between label and control images (Friston et al., 1995; Woods et al., 1998).

After these steps, the coregistration matrix between the structural data set and the preprocessed timeseries was determined using an affine transformation.

Then, additional preprocessing was performed using ANATICOR (Jo et al., 2010). Briefly, we first obtained individual masks of large ventricles and white matter from the structural scans segmentation using FreeSurfer³. The white matter mask was slightly eroded (one functional voxel) to prevent partial volume effects that might include signal from gray matter voxels in the mask. This step was not performed for the CSF mask, since with the functional voxel size used in this study the eroded CSF mask would not contain enough voxels in most subjects. Then, for each run, a global nuisance regressor was obtained extracting the EPI average time course within the ventricle mask and local nuisance regressors were obtained calculating for each gray matter voxel the average signal time course for all white matter voxels within a 3 cm radius (Jo et al., 2010). These nuisance regressors and the six regressors derived from motion parameters were removed from the EPI timeseries using AFNI’s @ANATICOR.

At this stage, pCASL time series with interleaved control and label volumes were rebuilt using the preprocessed data. These preprocessed pCASL data were then temporal filtered to separate CBF and BOLD signals (Chuang et al., 2008). Briefly, the ASL time series can be view as the sum of a component with a rapid modulation reflecting the alternating control and label images, and a component that is not modulated by the labeling process, reflecting slower BOLD weighted signal variations. An high-pass filtering with a cut-off frequency corresponding to 1/4TR has been demonstrated to be effective in retaining the modulated component while minimizing the BOLD contamination (Chuang et al., 2008; Tak et al., 2014). Note that this high-pass filtering does not restrict the CBF timeseries to the high frequency band since the CBF signal is derived afterward from the difference between subsequent scans (i.e., control—label). Specifically, CBF timeseries were obtained from TE1 data by high-pass filtering the corresponding preprocessed pCASL signal (>0.071 Hz, corresponding to 1/4TR), multiplying it by cos[π(n − 1)] (where n is the frame number), and then summing together every two images. BOLD timeseries were obtained from TE2 data by low-pass filtering (<0.071 Hz) the corresponding preprocessed pCASL signal and then summing together every two images.

After these processing procedures, both BOLD and CBF data were spatially normalized using Advanced Normalization Tools (ANTs; Avants et al., 2009). Briefly, the individual T1-weighted anatomical images of young and elderly subjects were first bias corrected using the N4ITK algorithm (Tustison et al., 2010) and then used to create a study specific template encompassing the age range in our study (Avants et al., 2010). Non-brain removal was also performed using BET (Smith, 2002) and a nonlinear warping of the resulting template to the MNI standard space (3 mm isotropic spatial resolution) was computed using symmetric diffeomorphic image normalization (Avants et al., 2008). Finally, BOLD and CBF data were spatial smoothed (6 mm FWHM) and band-pass filtered (0.01–0.071 Hz).

After these preprocessing procedures, linear correlation (Pearson r) between BOLD and CBF timeseries was calculated for each voxel to obtain individual maps representing the BOLD-CBF coupling during resting-state. In this step no time-shift was introduced between the two timecourses and the resulting correlation coefficients were indicated as r₀. Random effects group maps (one sample t-test) were then obtained after individual r₀ to z-Fisher transform for Young and Elderly separately and then compared between groups (two-sample unpaired t-test). Group maps were thresholded at p < 0.05, FDR corrected.

In a second step, the cross-correlation between BOLD and CBF timecourses was calculated again after introducing a variable time-shift τ in the BOLD signal (eq.1), following previous work (Fukunaga et al., 2008; Tak et al., 2014). In this calculation, both BOLD (t) and CBF (t) were first upsampled to a resolution of 100 ms using sync interpolation as implemented in the AFNI program 3Dtshift. The maximum correlation value (r_max) was retained for τ ranging between ±7 s (with a 350 ms step). Then, after transforming individual r_max to z-Fisher, the random effects group maps were calculated again for Young and Elderly and compared between groups.

$r_{\max }={\max }_{\tau }{\displaystyle {\sum }_{t}BOLD(t+\tau )CBF(t)-7s<\tau <7s}(1)$

Note that previous studies mostly considered the r_max approach, in order to maximize the statistical significance of correlation by correcting for potential temporal mismatch between the two hemodynamic signals (Fukunaga et al., 2008; Tak et al., 2014). Here we investigated both r₀ and r_max and in particular how the two groups compare with respect to these metrics. Our choice was motivated by the fact that the temporal mismatch between BOLD and CBF might be group dependent (due to e.g., age-related vascular effects). A quantitative estimation of the effect of time-shift correction on the correlation coefficient was performed in selected regions of interest (ROIs) as follows.

We focused on ROIs in the default mode network (DMN) and frontoparietal network (FPN) that are the two most investigated brain networks. In order to allow the definition of ROIs independent from the previous analysis, independent component analysis (ICA) was performed, using a probabilistic algorithm (Beckmann and Smith, 2004) as implemented in the MELODIC tool of FSL (FMRIB Software Library). Briefly, individual functional datasets were first temporally concatenated across subjects, groups (Young and Elderly) and modalities (BOLD and CBF) to form a single 4D data set to be used as input for the probabilistic ICA algorithm. The number of components was fixed to 20. Dual-regression (Beckmann et al., 2009) was then used to identify individual spatial maps for each independent component to be used as input for the second-level group analysis. In this calculation random-effects group statistical maps were obtained for each component using permutation-based non-parametric testing (5000 permutations; Nichols and Holmes, 2002). Multiple comparisons correction was addressed applying a cluster-based threshold of Z > 2.3 and a family-wise-error (FWE) corrected cluster significance of p < 0.01 for the suprathreshold clusters. Using these group maps, the DMN and FPN were easily identified from their characteristic spatial pattern, including posterior cingulate cortex, bilateral angular gyrus and ventromedial prefrontal cortex for DMN, bilateral inferior parietal lobe and bilateral middle frontal gyrus for FPN. For our ROIs definition we pooled BOLD and CBF data because a high degree of spatial overlap has been shown for the two modalities in both resting state and task paradigms (Mayhew et al., 2014; Jann et al., 2015). Furthermore, in order to consider gray matter voxels only, a binary mask was defined by averaging the normalized individual gray matter masks obtained from the FreeSurfer segmentation and thresholding at 0.3 (Jann et al., 2015). The final masks representing the investigated ROIs were defined multiplying this binary mask with the ICA clusters.

Then, mean r₀ and r_max values were extracted from these ROIs and the r_max − r₀ difference was compared between groups (two-sample unpaired t-test). Correction for multiple comparisons was performed using FDR. We expect this difference to be larger for elderlies due to an increased time-shift between BOLD and CBF timecourses, possibly reflecting age related modifications of vascular response.

In addition, we computed the resting state fluctuation amplitude (RSFA; Kannurpatti and Biswal, 2008) for both BOLD and CBF timeseries in Elderly and Young, using the AFNI program 3dRSFC. The corresponding values were extracted from our ROIs and compared between groups (two-sample unpaired t-tests with correction for multiple comparisons using FDR). This analysis aimed at evaluating the impact of potential age dependent differences in the amplitude of resting state BOLD and CBF fluctuations on the observed linear correlation between the two signals, since a reduced fluctuation amplitude could lead to a decreased correlation coefficient in presence of noise (Liu, 2013).

Finally, we also extracted regional perfusion values from quantitative CBF maps provided by the background suppressed ASL data and a single compartment model (Buxton et al., 1998; Alsop et al., 2015):

$CBF=\frac{6000\lambda (S{I}_c-S{I}_L)e^{\frac{PLD}{T_{1A}}}}{2\alpha {\alpha }_{\text{inv}}{T}_{1\text{A}}{M}_0(1-e^{-\frac{\tau }{T1\text{A}}})}[\text{ml}/100g/\min ](2)$

Where λ is the blood-brain partition coefficient (0.9 ml/g), SI_C and SI_L are the means over time of the control and label images respectively, PLD is the slice-dependent post label delay (1900–2800 ms), T_1A is the longitudinal relaxation time of arterial blood (1650 ms at 3T), α is the labeling efficiency (0.85; Alsop et al., 2015), α_inv is a correction factor for the background suppression (0.83; van Osch et al., 2009), M₀ is the equilibrium magnetization signal, and τ is the label duration (1750 ms).

Results

Random effects group maps obtained using the z-Fisher transformed voxelwise correlation (r₀, see “Materials and Methods” section) between BOLD and CBF timeseries during resting state are shown in Figure 1A for the two groups. In young subjects, a significant correlation between spontaneous fluctuations of BOLD and CBF signals was observed in most cortical areas, whereas in elderly subjects this correlation was markedly reduced (mean r₀ values in significant areas were 0.24 ± 0.03 for Young and 0.18 ± 0.04 for Elderly). The BOLD-CBF coupling was stronger in medial, parietal and frontal regions, which are part of well-established resting state brain networks. No significant correlation was observed in white matter or other brain structures. The voxelwise contrast between groups showed a significantly lower correlation for elderly subjects in the left supramarginal gyrus (MNI coordinates: −59, −34, 33; Figure 1B).

FIGURE 1

Figure 1. (A) Random effect group maps showing the dynamic coupling between spontaneous blood oxygenation level dependent (BOLD) and cerebral blood flow (CBF) fluctuations for Young (Y) and Elderly (E) groups when potential time shifts between the two signals were not compensated for (r₀, see “Materials and Methods” section). (B) Between-group contrast comparing the dynamic coupling between spontaneous BOLD and CBF fluctuations for Y and E, using r₀ values. A significant age related decrease of BOLD-CBF coupling is observed in the left supramarginal gyrus (MNI coordinates: −59, −34, 33). The statistical maps were thresholded at p < 0.05 (corrected for multiple comparisons using FDR) and superimposed on the partially inflated study specific template.

Random effects group maps obtained with the z-Fisher transformed maximum voxelwise correlation between BOLD and CBF timeseries (r_max, obtained introducing a temporal shift between the two signals as described in “Materials and Methods” section) are reported in Figure 2A. As expected, the number of cortical voxels showing a significant coupling increased for both groups. Interestingly, this increase was more pronounced for Elderly (mean r_max values in significant areas were 0.32 ± 0.02 for Young and 0.30 ± 0.03 for Elderly). However, the voxelwise contrast between groups still showed a significantly lower correlation for elderly subjects in the left supramarginal gyrus (MNI coordinates: −58, −36, 34; Figure 2B), in a cluster of voxels largely overlapping with that observed in Figure 1B.

FIGURE 2

Figure 2. (A) Random effect group maps showing the dynamic coupling between spontaneous BOLD and CBF fluctuations for Y and E groups when potential time shifts between the two signals were compensated for (r_max, see “Materials and Methods” section). (B) Between-group contrast comparing the dynamic coupling between spontaneous BOLD and CBF fluctuations for Y and E, using r_max values. Note that a significant age related decrease of BOLD-CBF coupling is observed in the left supramarginal gyrus (MNI coordinates: −58, −36, 34) even after temporal shift correction, in a cluster of voxels largely overlapping with that observed in Figure 1B. The statistical maps were thresholded at p < 0.05 (corrected for multiple comparisons using FDR) and superimposed on the partially inflated study specific template.

DMN and FPN were easily identified in the ICA results, with one component showing the typical DMN pattern, one component showing the left FPN and another component the right FPN (Figure 3). The results of the ROI approach for the selected regions in DMN and FPN are reported in Figure 4A. The comparison of the r_max − r₀ difference between groups showed significantly larger values for Elderly in all ROIs, indicating a more pronounced effect of time-shift corrections on the linear correlation between BOLD and CBF timecourses in elderly subjects. An effect size estimate (Cohen’s d) was also computed for the between-group comparison of the r_max − r₀ difference in these ROIs. The observed values were quite large, ranging from 0.78 to 1.24.

FIGURE 3

Figure 3. Investigated nodes of default mode network (DMN) and frontoparietal network (FPN), defined with independent component analysis (ICA) calculated pooling the two groups and modalities (FWE corrected cluster significance of p < 0.01). Regions of interest (ROIs) were defined masking these clusters with a gray matter binary mask (L_AG/R_AG, left/right angular gyrus; PCC, posterior cingulate cortex; Med_FG, medial frontal gyrus; L_IPL/R_IPL, left/right inferior parietal lobule; L_MidFG/R_MidFG, left/right middle frontal gyrus).

FIGURE 4

Figure 4. (A) Between-group comparison of r_max −r₀ values, resting state fluctuation amplitudes (RSFAs) and baseline perfusion values extracted in ROIs of Figure 3 (two-sample unpaired t-test; *p < 0.05 FDR corrected, **p < 0.01 FDR corrected, ***p < 0.001 FDR corrected, ****p < 0.0005 FDR corrected). (B) The same values extracted from the ROI obtained pooling the two clusters in Figures 1B, 2B and masking with gray matter. Error bars are standard errors.

A larger r_max − r₀ difference for Elderly with respect to Young was also observed in the supramarginal gyrus, without reaching statistical significance (Figure 4B, values extracted from a mask obtained pooling the clusters in Figures 1B, 2A). The time-shift parameters in the different ROIs ranged from −1.12 s to −1.81 s for Elderly and from −0.4 s to −1.4 s for Young. The between-group difference was significant in the right angular gyrus (−1.8 s for Elderly and −0.5 s for Young, p < 0.03, unpaired t-test), indicating an increased delay of BOLD dynamics (with respect to CBF) in elderlies.

The baseline CBF values showed the well known age related decrease for all the investigated regions (Figure 3). However, no significant correlation of the BOLD-CBF coupling with baseline CBF values was observed across subjects, for either Elderly or Young groups.

The control analysis on the amplitude of resting state low frequency fluctuations revealed that RSFA values in these regions did not change significantly with age, except for the BOLD-RSFA in the frontal regions of FPN (Figure 4A). However, no region showed a significant across-subjects correlation of the values of BOLD-CBF coupling with BOLD-RSFA or CBF-RSFA values, for either Elderly or Young groups.

The comparison of the run-averaged DVARS metrics did not show statistically significant differences across groups for either control or label images (F_(2,120) = 1.29; p = 0.26). Additional control analysis showed no significant correlation of DVARS values with the BOLD-CBF coupling for either Elderly or Young groups.

Discussion

The present results showed that in young subjects spontaneous BOLD and CBF fluctuations are significantly synchronized in most cortical areas and especially within the major nodes of prominent resting-state networks, confirming previous evidence (Tak et al., 2014, 2015; Chen et al., 2015). As a new finding, this resting-state BOLD-CBF dynamic coupling was reduced in elderly individuals, especially in the left supramarginal gyrus. Moreover, this decrease was not related to a reduced amplitude of either BOLD or CBF spontaneous fluctuations. Furthermore, elderlies showed a larger increase in the correlation coefficient after the introduction of a relative time-shift between BOLD and CBF timecourses.

Recently, an increasing evidence that intrinsic CBF fluctuations are a major contributor to the resting state BOLD signal has been reported. Indeed, different studies showed a high level of spatial overlap between connectivity maps obtained from BOLD and CBF timecourses (Viviani et al., 2011; Jann et al., 2015). Moreover, a set of studies investigated the dynamic relationship between spontaneous CBF and BOLD fluctuations (Fukunaga et al., 2008; Wu et al., 2009; Tak et al., 2014, 2015; Chen et al., 2015; Cohen et al., 2017). Fukunaga et al. (2008) were the first to investigate the dynamic coupling between BOLD and perfusion fluctuations during resting state, showing that the two signals are correlated in most part of the cortex, with little involvement of white matter and cerebrospinal fluid. Furthermore, using the BOLD/perfusion ratio to target resting state oxidative metabolism fluctuations, they found a similar BOLD/CBF coupling with respect to a visual task induced brain activity, adding evidence to a metabolic/neuronal origin of spontaneous BOLD and CBF fluctuations. In a subsequent work, Wu et al. (2009) obtained functional connectivity maps using CMRO₂ time series derived from simultaneous BOLD and CBF time series acquired with ASL. Although the CMRO₂ time courses were estimated assuming the steady-state biophysical BOLD model, that might have some limitations in a more dynamical situation (Simon and Buxton, 2015), the strong similarity observed by Wu et al. (2009) between functional connectivity maps obtained from BOLD, CBF and CMRO₂ time courses added further evidence to a significant BOLD-CBF coupling during resting state. In a recent work the study of dynamic BOLD-CBF relationship during resting state was addressed with significant methodological improvements, including physiological noise correction in the tag and control ASL images and taking into account the influence of global cardiac fluctuations (Tak et al., 2014). This work demonstrated that the resting-state BOLD-CBF coupling strength, although varying across the brain, was stronger in the gray matter and in particular in the major nodes of well-established functional networks. Moreover, the BOLD-CBF coupling observed by Tak and collaborators was significantly reduced in voxels associated with a high macrovascular content, suggesting that the component of spontaneous BOLD signal fluctuations that is more directly driven by dynamic CBF fluctuations is more likely related to neuronal activity, which is known to modulate the microvascular response.

Keeping in line with previous evidence, our results confirm that in young individuals spontaneous BOLD and CBF fluctuations are significantly synchronized in most cortical areas and especially in regions of the major resting state networks. As a new finding, we observed a general age-related decrease of the BOLD-CBF dynamic coupling during resting state across cortical areas. The between-group voxelwise contrast showed that this decrease was statistically significant in a cluster of voxels in the left supramarginal gyrus.

Slight discrepancies between BOLD and CBF dynamics have been previously reported even in young subjects using e.g., visual or motor stimuli (Obata et al., 2004; Cavusoglu et al., 2012). These discrepancies can be expected, due to possible temporal uncoupling of the involved physiological responses (i.e., CBF, CBV and CMRO₂) that have a competing effect on the BOLD signal amplitude. Indeed, a temporal uncoupling of e.g., CBF and venous CBV has been proposed in early models trying to explain BOLD signal transients like the post-stimulus undershoot or the initial overshoot during long stimulation blocks (Buxton et al., 1998, 2004; Buxton, 2012). Although a neuronal contribution to these transients has been recognized (Sadaghiani et al., 2009; Mullinger et al., 2014, 2017), there is increasing evidence that vascular mechanisms related to delayed compliance of venous vessels play an important role as well (Chen and Pike, 2009; Havlicek et al., 2017). In particular, a recent work suggests that the venous CBV response can be an order of magnitude slower than either CBF or CMRO₂ (Simon and Buxton, 2015).

Despite the specific physiological mechanisms responsible for the decreased BOLD-CBF dynamic coupling that we observed in elderly cannot be determined using the present data alone, it could be argued that increased vessel stiffening due to age (Podlutsky et al., 2010; Trott et al., 2011; Csiszár et al., 2015; Tsvetanov et al., 2015; Chiarelli et al., 2017; Tan et al., 2017; Tarantini et al., 2017; Toth et al., 2017) can lead to increased delay in venous compliance, thus introducing further temporal discrepancies between BOLD and CBF timecourses. Animal studies also showed that modifications of vessel compliance with age decrease the ratio of CBV to CBF responses (Dubeau et al., 2011; Desjardins et al., 2014) that in turn would also affect the BOLD response. Moreover, the ratio of CBV to CBF response is also reflected in the Grubb’s parameter (Grubb et al., 1974) that is involved in calibrated BOLD models. Importantly, this parameter has been shown to vary not only with aging (Dubeau et al., 2011) but also during different phases of functional stimulation, leading to transient relationships between CBF, CBV and BOLD changes (Kida et al., 2006), suggesting a possible effect on BOLD-CBF coupling. The potential role of delayed vessel compliance is also supported by our second finding, i.e., that introducing a time-shift correction between BOLD and CBF signals has more effect for elderlies than for young subjects in increasing the calculated correlation coefficient. However, even with time-shift correction, the left supramarginal gyrus showed a significant between-group difference in the voxelwise contrast, indicating that the dynamic uncoupling of BOLD and CBF timecourses in this region is more affected by the aging process and can be only partially explained with a simple time-shift between the two signals. Interestingly, different neuroimaging and transcranial magnetic stimulation studies reported that the left supramarginal gyrus is involved in verbal working memory and episodic memory (Romero et al., 2006; Koelsch et al., 2009; Deschamps et al., 2014; Thakral et al., 2017), which are functions known to be impaired with age. Considering that our investigated population was selected among healthy elderlies that did not show a compromised performance in memory tests, our results acquire an increased interest. Indeed, the observed age related differences in the dynamic BOLD-CBF coupling could be considered a potential biomarker that might anticipate changes in neuronal function. In this regard, this type of index could be spatially more specific than e.g., age-related decrease of baseline perfusion that we and others (Parkes et al., 2004; Ambarki et al., 2015; De Vis et al., 2015) clearly observed in most cortical areas. However, we are aware that a full validation of this hypothesis would require additional data, with a larger sample size and possibly longitudinal studies.

Potential Limitations and Caveats

As in all studies using dynamic ASL to assess brain function, a major concern is minimizing BOLD contamination of the CBF signal. We addressed this issue by using the short echo time data (which have minimal sensitivity to the BOLD effect) and by high pass filtering the ASL signal followed by demodulation to derive the CBF timecourse (Chuang et al., 2008). This approach can be considered a generalization of previously proposed methods like e.g., sinc interpolation (Liu and Wong, 2005), with similar efficiency in removing BOLD contamination from the perfusion timeseries. Another important concern in resting state fMRI studies is the physiological noise due to cardiac and respiratory activity. While different methods of noise cleanup have been established for the BOLD signal (for a review see; Murphy et al., 2013; Liu, 2016), there is scarce evidence on the application of similar procedures to ASL. We chose the approach using RETROICOR to correct time-locked effects of cardiac and respiratory fluctuations on the EPI signal of label and control separately (Restom et al., 2006). Indeed, modeling the effect of the cardiac and respiratory cycles on the derived CBF weighted timeseries had minimal effect on signal quality (Restom et al., 2006). However, as observed by the same authors, future approaches using additional information such as end tidal carbon dioxide measurements could allow improved modeling of the effect of physiological noise on the perfusion timecourses, extending previous work on BOLD (Wise et al., 2004). An additional source of noise in the fMRI signals is represented by head motion, that has been shown to introduce spurious correlations in functional networks investigated in resting state BOLD studies (Power et al., 2012). This issue is particularly important when comparing different populations that could differ in the amount of head motion, like e.g healthy subjects vs. patients or young subjects vs. elderlies and censoring procedures have been proposed, in addition to regression of motion parameters, to mitigate these effects in resting state BOLD timeseries (Power et al., 2014). Again, the extension of these methods to CBF timeseries is not straightforward because the perfusion information is carried by the difference image between label and control that should be censored in pairs, thus strongly reducing the number of usable timepoints. However, motion parameters derived during preprocessing of label and control images did not differ significantly between groups in our data, suggesting a negligible effect of motion on the observed results.

In addition to noise concerns, potential between-group differences in BOLD and/or CBF signal fluctuation amplitude could also affect the calculated correlations. Indeed, while from a mathematical point of view the linear correlation of two timecourses is not affected by the amplitude of the signals, in presence of noise a reduced fluctuation amplitude could lead to a decreased correlation coefficient (Liu, 2013). However, the two groups showed comparable levels of signal fluctuations for both BOLD and CBF in the majority of investigated ROIs, suggesting no bias related to this issue. Even in the two regions showing a significant decrease of BOLD fluctuations in elderlies, any significant correlation of RSFA with BOLD-CBF correlation was observed. A reduced fluctuation amplitude in elderlies could have been expected due to gray matter atrophy. In this regard, the observation of similar fluctuation amplitudes of the functional signals in the two groups also helps to mitigate the concern of potential confounds due to the presence of atrophy in elderlies that could have biased our results.

Potential biases in the results could also be introduced by the six older adults using antihypertensive medication that would possibly alter the hemodynamic response in these subjects. In our group of subjects a control analysis revealed similar BOLD-CBF coupling when comparing this hypertensive subgroup with the other elderlies, with no significant differences. Nevertheless, future studies should address this issue using larger sample sizes.

A further concern could be raised regarding the postlabel delay (1000 ms for the first slice) used for the study of dynamic BOLD-CBF coupling that can be too short to ensure a complete delivery of the bolus to the tissue at the time of acquisition, in particular for elderlies. However, while this is an important issue for absolute CBF quantification, short postlabel delays are often used or even recommended for functional applications of ASL because offer a larger sensitivity and temporal resolution without introducing bright spots in activation maps due to large vessel contribution or introducing distortions of timecourses (Gonzalez-At et al., 2000; Zappe et al., 2008). In contrast, for the quantification of baseline CBF, we used a longer postlabel delay (1900 ms for the first slice), following current guidelines for both young and healthy elderlies (Alsop et al., 2015).

Finally, another limitation of the study could arise from the small sample size, potentially resulting in a low statistical power. Indeed, despite statistical significance was obtained for the between-group comparison of our main metrics of interest, future studies with larger sample size possibly including elderlies with impaired cognitive functions would allow further assessment of the significance of our findings.

Conclusion

In this study we observed an age-related decrease of the temporal correlation between the BOLD and CBF spontaneous fluctuations in the cortex. Interestingly, when compensating for potential time-shifts between the two signals, the increase of the correlation coefficient was larger for elderlies. However a significant between-group difference was observed in the left supramargynal gyrus even after time-shift correction, suggesting a more pronounced age-related BOLD-CBF uncoupling in this area. These results suggest that the study of dynamic coupling between spontaneous BOLD-CBF fluctuations using the simultaneous acquisition of the two signals with ASL is a promising technique to study-resting state brain function in aging and disease.

Author Contributions

PC and AF: conceived and designed the study and wrote the article. FC and GB: enrolled and tested the subjects. PC and MGP: performed the experiments. PC: implemented AFNI and Python scripts. PC, FC and AF: performed the analyses.

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.

Footnotes

1. ^ afni.nimh.nih.gov/afni
2. ^ http://www.python.org
3. ^ http://surfer.nmr.mgh.harvard.edu

References

Alsop, D. C., Detre, J. A., Golay, X., Günther, M., Hendrikse, J., Hernandez-Garcia, L., et al. (2015). Recommended implementation of arterial spin-labeled perfusion MRI for clinical applications: a consensus of the ISMRM perfusion study group and the European consortium for ASL in dementia. Magn. Reson. Med. 73, 102–116. doi: 10.1002/mrm.25197

PubMed Abstract | CrossRef Full Text | Google Scholar

Ambarki, K., Wahlin, A., Zarrinkoob, L., Wirestam, R., Petr, J., Malm, J., et al. (2015). Accuracy of parenchymal cerebral blood flow measurements using pseudocontinuous arterial spin-labeling in healthy volunteers. Am. J. Neuroradiol. 36, 1816–1821. doi: 10.3174/ajnr.A4367

PubMed Abstract | CrossRef Full Text | Google Scholar

Ances, B. M., Leontiev, O., Perthen, J. E., Liang, C., Lansing, A. E., and Buxton, R. B. (2008). Regional differences in the coupling of cerebral blood flow and oxygen metabolism changes in response to activation: implications for BOLD-fMRI. Neuroimage 39, 1510–1521. doi: 10.1016/j.neuroimage.2007.11.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Aslan, S., Xu, F., Wang, P. L., Uh, J., Yezhuvath, U. S., van Osch, M., et al. (2010). Estimation of labeling efficiency in pseudocontinuous arterial spin labeling. Magn. Reson. Med. 63, 765–771. doi: 10.1002/mrm.22245

PubMed Abstract | CrossRef Full Text | Google Scholar

Avants, B. B., Cook, P., Pluta, J., Duda, J. T., Rao, H., Giannetta, J., et al. (2009). Multivariate diffeomorphic analysis of longitudinal increase in white matter directionality and decrease in cortical thickness between ages 14 and 18. Neuroimage 47:S105. doi: 10.1016/s1053-8119(09)70918-8

CrossRef Full Text | Google Scholar

Avants, B. B., Epstein, C., Grossman, M., and Gee, J. (2008). Symmetric diffeomorphic image registration with cross-correlation: evaluating automated labeling of elderly and neurodegenerative brain. Med. Image Anal. 12, 26–41. doi: 10.1016/j.media.2007.06.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Avants, B. B., Yushkevich, P., Pluta, J., Minkoff, D., Korczykowski, M., Detre, J., et al. (2010). The optimal template effect in hippocampus studies of diseased populations. Neuroimage 49, 2457–2466. doi: 10.1016/j.neuroimage.2009.09.062

PubMed Abstract | CrossRef Full Text | Google Scholar

Bandettini, P. A., Wong, E. C., Hinks, R. S., Tikofsky, R. S., and Hyde, J. S. (1992). Time course EPI of human brain function during task activation. Magn. Reson. Med. 25, 390–397. doi: 10.1002/mrm.1910250220

PubMed Abstract | CrossRef Full Text | Google Scholar

Beckmann, C. F., Mackay, C. E., Filippini, N., and Smith, S. M. (2009). Group comparison of resting-state FMRI data using multi-subject ICA and dual regression. Neuroimage 47:S148. doi: 10.1016/s1053-8119(09)71511-3

CrossRef Full Text | Google Scholar

Beckmann, C. F., and Smith, S. M. (2004). Probabilistic independent component analysis for functional magnetic resonance imaging. IEEE Trans. Med. Imaging 23, 137–152. doi: 10.1109/tmi.2003.822821

PubMed Abstract | CrossRef Full Text | Google Scholar

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 | CrossRef Full Text | Google Scholar

Blockley, N. P., Griffeth, V. E. M., and Buxton, R. B. (2012). A general analysis of calibrated BOLD methodology for measuring CMRO₂ responses: comparison of a new approach with existing methods. Neuroimage 60, 279–289. doi: 10.1016/j.neuroimage.2011.11.081

PubMed Abstract | CrossRef Full Text | Google Scholar

Blockley, N. P., Griffeth, V. E. M., Simon, A. B., Dubowitz, D. J., and Buxton, R. B. (2015). Calibrating the BOLD response without administering gases: comparison of hypercapnia calibration with calibration using an asymmetric spin echo. Neuroimage 104, 423–429. doi: 10.1016/j.neuroimage.2014.09.061

PubMed Abstract | CrossRef Full Text | Google Scholar

Brier, M. R., Thomas, J. B., and Ances, B. M. (2014). Network dysfunction in Alzheimer’s disease: refining the disconnection hypothesis. Brain Connect. 4, 299–311. doi: 10.1089/brain.2014.0236

PubMed Abstract | CrossRef Full Text | Google Scholar

Bulte, D. P., Drescher, K., and Jezzard, P. (2009). Comparison of hypercapnia-based calibration techniques for measurement of cerebral oxygen metabolism with MRI. Magn. Reson. Med. 61, 391–398. doi: 10.1002/mrm.21862

PubMed Abstract | CrossRef Full Text | Google Scholar

Bulte, D. P., Kelly, M., Germuska, M., Xie, J., Chappell, M. A., Okell, T. W., et al. (2012). Quantitative measurement of cerebral physiology using respiratory-calibrated MRI. Neuroimage 60, 582–591. doi: 10.1016/j.neuroimage.2011.12.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Buxton, R. B. (2010). Interpreting oxygenation-based neuroimaging signals: the importance and the challenge of understanding brain oxygen metabolism. Front. Neuroenergetics 2:8. doi: 10.3389/fnene.2010.00008

PubMed Abstract | CrossRef Full Text | Google Scholar

Buxton, R. B. (2012). Dynamic models of BOLD contrast. Neuroimage 62, 953–961. doi: 10.1016/j.neuroimage.2012.01.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Buxton, R. B., Frank, L. R., Wong, E. C., Siewert, B., Warach, S., and Edelman, R. R. (1998). A general kinetic model for quantitative perfusion imaging with arterial spin labeling. Magn. Reson. Med. 40, 383–396. doi: 10.1002/mrm.1910400308

PubMed Abstract | CrossRef Full Text | Google Scholar

Buxton, R. B., Uludağ, K., Dubowitz, D. J., and Liu, T. T. (2004). Modeling the hemodynamic response to brain activation. Neuroimage 23, S220–S233. doi: 10.1016/j.neuroimage.2004.07.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Cavusoglu, M., Bartels, A., Yesilyurt, B., and Uludağ, K. (2012). Retinotopic maps and hemodynamic delays in the human visual cortex measured using arterial spin labeling. Neuroimage 59, 4044–4054. doi: 10.1016/j.neuroimage.2011.10.056

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, J. J., Jann, K., and Wang, D. J. J. (2015). Characterizing resting-state brain function using arterial spin labeling. Brain Connect. 5, 527–542. doi: 10.1089/brain.2015.0344

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, J. J., and Pike, G. B. (2009). Origins of the BOLD post-stimulus undershoot. Neuroimage 46, 559–568. doi: 10.1016/j.neuroimage.2009.03.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, J. J., Rosas, H. D., and Salat, D. H. (2011). Age-associated reductions in cerebral blood flow are independent from regional atrophy. Neuroimage 55, 468–478. doi: 10.1016/j.neuroimage.2010.12.032

PubMed Abstract | CrossRef Full Text | Google Scholar

Chiarelli, P. A., Bulte, D. P., Wise, R., Gallichan, D., and Jezzard, P. (2007). A calibration method for quantitative BOLD fMRI based on hyperoxia. Neuroimage 37, 808–820. doi: 10.1016/j.neuroimage.2007.05.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Chiarelli, A. M., Fletcher, M. A., Tan, C. H., Low, K. A., Maclin, E. L., Zimmerman, B., et al. (2017). Individual differences in regional cortical volumes across the life span are associated with regional optical measures of arterial elasticity. Neuroimage 162, 199–213. doi: 10.1016/j.neuroimage.2017.08.064

PubMed Abstract | CrossRef Full Text | Google Scholar

Chuang, K.-H., van Gelderen, P., Merkle, H., Bodurka, J., Ikonomidou, V. N., Koretsky, A. P., et al. (2008). Mapping resting-state functional connectivity using perfusion MRI. Neuroimage 40, 1595–1605. doi: 10.1016/j.neuroimage.2008.01.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Cohen, A. D., Nencka, A. S., Lebel, R. M., and Wang, Y. (2017). Multiband multi-echo imaging of simultaneous oxygenation and flow timeseries for resting state connectivity. PLoS One 12:e0169253. doi: 10.1371/journal.pone.0169253

PubMed Abstract | CrossRef Full Text | Google Scholar

Csiszár, A., Csiszar, A., Pinto, J. T., Gautam, T., Kleusch, C., Hoffmann, B., et al. (2015). Resveratrol encapsulated in novel fusogenic liposomes activates Nrf2 and attenuates oxidative stress in cerebromicrovascular endothelial cells from aged rats. J. Gerontol. A Biol. Sci. Med. Sci. 70, 303–313. doi: 10.1093/gerona/glu029

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, W., Garcia, D., de Bazelaire, C., and Alsop, D. C. (2008). Continuous flow-driven inversion for arterial spin labeling using pulsed radio frequency and gradient fields. Magn. Reson. Med. 60, 1488–1497. doi: 10.1002/mrm.21790

PubMed Abstract | CrossRef Full Text | Google Scholar

Davis, T. L., Kwong, K. K., Weisskoff, R. M., and Rosen, B. R. (1998). Calibrated functional MRI: mapping the dynamics of oxidative metabolism. Proc. Natl. Acad. Sci. U S A 95, 1834–1839. doi: 10.1073/pnas.95.4.1834

PubMed Abstract | CrossRef Full Text | Google Scholar

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 | CrossRef Full Text | Google Scholar

De Vis, J. B., Hendrikse, J., Bhogal, A., Adams, A., Kappelle, L. J., and Petersen, E. T. (2015). Age-related changes in brain hemodynamics; A calibrated MRI study. Hum. Brain Mapp. 36, 3973–3987. doi: 10.1002/hbm.22891

PubMed Abstract | CrossRef Full Text | Google Scholar

Deschamps, I., Baum, S. R., and Gracco, V. L. (2014). On the role of the supramarginal gyrus in phonological processing and verbal working memory: evidence from rTMS studies. Neuropsychologia 53, 39–46. doi: 10.1016/j.neuropsychologia.2013.10.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Desjardins, M., Berti, R., Pouliot, P., Dubeau, S., and Lesage, F. (2014). Multimodal study of the hemodynamic response to hypercapnia in anesthetized aged rats. Neurosci. Lett. 563, 33–37. doi: 10.1016/j.neulet.2014.01.027

PubMed Abstract | CrossRef Full Text | Google Scholar

D’Esposito, M., Zarahn, E., Aguirre, G. K., and Neuroimage, B. R. (1999). The effect of normal aging on the coupling of neural activity to the bold hemodynamic response. Neuroimage 10, 6–14. doi: 10.1006/nimg.1999.0444

PubMed Abstract | CrossRef Full Text | Google Scholar

Detre, J. A., Leigh, J. S., Williams, D. S., and Koretsky, A. P. (1992). Perfusion imaging. Magn. Reson. Med. 23, 37–45. doi: 10.1002/mrm.1910230106

PubMed Abstract | CrossRef Full Text | Google Scholar

Dubeau, S., Ferland, G., Gaudreau, P., Beaumont, E., and Lesage, F. (2011). Cerebrovascular hemodynamic correlates of aging in the Lou/c rat: a model of healthy aging. Neuroimage 56, 1892–1901. doi: 10.1016/j.neuroimage.2011.03.076

PubMed Abstract | CrossRef Full Text | Google Scholar

Esposito, R., Cieri, F., Chiacchiaretta, P., Cera, N., Lauriola, M., Di Giannantonio, M., et al. (2018). Modifications in resting state functional anticorrelation between default mode network and dorsal attention network: comparison among young adults, healthy elders and mild cognitive impairment patients. Brain Imaging Behav. 12, 127–141. doi: 10.1007/s11682-017-9686-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Fernández-Seara, M. A., Mengual, E., Vidorreta, M., Castellanos, G., Irigoyen, J., Erro, E., et al. (2015). Resting state functional connectivity of the subthalamic nucleus in Parkinson’s disease assessed using arterial spin-labeled perfusion fMRI. Hum. Brain Mapp. 36, 1937–1950. doi: 10.1002/hbm.22747

PubMed Abstract | CrossRef Full Text | Google Scholar

Folstein, M. F., Folstein, S. E., and McHugh, P. R. (1975). “Mini-mental state”. A practical method for grading the cognitive state of patients for the clinician. J. Psychiatr. Res. 12, 189–198. doi: 10.1016/0022-3956(75)90026-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Fox, M. D., and Greicius, M. (2010). Clinical applications of resting state functional connectivity. Front. Syst. Neurosci. 4:19. doi: 10.3389/fnsys.2010.00019

PubMed Abstract | CrossRef Full Text | Google Scholar

Friston, K. J., Ashburner, J., Frith, C. D., Poline, J. B., Heather, J. D., and Frackowiak, R. S. J. (1995). Spatial registration and normalization of images. Hum. Brain Mapp. 3, 165–189. doi: 10.1002/hbm.460030303

CrossRef Full Text | Google Scholar

Fujita, N., Matsumoto, K., Tanaka, H., Watanabe, Y., and Murase, K. (2006). Quantitative study of changes in oxidative metabolism during visual stimulation using absolute relaxation rates. NMR Biomed. 19, 60–68. doi: 10.1002/nbm.1001

PubMed Abstract | CrossRef Full Text | Google Scholar

Fukunaga, M., Horovitz, S. G., de Zwart, J. A., van Gelderen, P., Balkin, T. J., Braun, A. R., et al. (2008). Metabolic origin of bold signal fluctuations in the absence of stimuli. J. Cereb. Blood Flow Metab. 28, 1377–1387. doi: 10.1038/jcbfm.2008.25

PubMed Abstract | CrossRef Full Text | Google Scholar

Gardini, S., Sambataro, F., Cuetos, F., Fasano, F., Marchi, M., Crisi, G., et al. (2014). Increased default mode network functional connectivity in mild cognitive impairment: a detrimental brain mechanism associated with poor semantic memory performance. Alzheimers Dement. 10, P823–P824. doi: 10.1016/j.jalz.2014.05.1623

CrossRef Full Text | Google Scholar

Garrett, D. D., Lindenberger, U., Hoge, R. D., and Gauthier, C. J. (2017). Age differences in brain signal variability are robust to multiple vascular controls. Sci. Rep. 7:10149. doi: 10.1038/s41598-017-09752-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Gauthier, C. J., Desjardins-Crépeau, L., Madjar, C., Bherer, L., and Hoge, R. D. (2012). Absolute quantification of resting oxygen metabolism and metabolic reactivity during functional activation using QUO2 MRI. Neuroimage 63, 1353–1363. doi: 10.1016/j.neuroimage.2012.07.065

PubMed Abstract | CrossRef Full Text | Google Scholar

Gauthier, C. J., and Hoge, R. D. (2013). A generalized procedure for calibrated MRI incorporating hyperoxia and hypercapnia. Hum. Brain Mapp. 34, 1053–1069. doi: 10.1002/hbm.21495

PubMed Abstract | CrossRef Full Text | Google Scholar

Germuska, M., and Bulte, D. P. (2014). MRI measurement of oxygen extraction fraction, mean vessel size and cerebral blood volume using serial hyperoxia and hypercapnia. Neuroimage 92, 132–142. doi: 10.1016/j.neuroimage.2014.02.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Germuska, M., Merola, A., Murphy, K., Babic, A., Richmond, L., Khot, S., et al. (2016). A forward modelling approach for the estimation of oxygen extraction fraction by calibrated fMRI. Neuroimage 139, 313–323. doi: 10.1016/j.neuroimage.2016.06.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Glover, G. H., Li, T. Q., and Ress, D. (2000). Image-based method for retrospective correction of physiological motion effects in fMRI: RETROICOR. Magn. Reson. Med. 44, 162–167. doi: 10.1002/1522-2594(200007)44:1<162::aid-mrm23>3.3.co;2-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Gonzalez-At, J. B., Alsop, D. C., and Detre, J. A. (2000). Cerebral perfusion and arterial transit time changes during task activation determined with continuous arterial spin labeling. Magn. Reson. Med. 43, 739–746. doi: 10.1002/(sici)1522-2594(200005)43:5<739::aid-mrm17>3.0.co;2-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Griffeth, V. E. M., and Buxton, R. B. (2011). A theoretical framework for estimating cerebral oxygen metabolism changes using the calibrated-BOLD method: modeling the effects of blood volume distribution, hematocrit, oxygen extraction fraction, and tissue signal properties on the BOLD signal. Neuroimage 58, 198–212. doi: 10.1016/j.neuroimage.2011.05.077

PubMed Abstract | CrossRef Full Text | Google Scholar

Grubb, R. L. Jr., Raichle, M. E., Eichling, J. O., and Ter-Pogossian, N. M. (1974). The effects of changes in PaCO₂ on cerebral blood volume, blood flow, and vascular mean transit time. Stroke 5, 630–639. doi: 10.1161/01.str.5.5.630

PubMed Abstract | CrossRef Full Text | Google Scholar

Havlicek, M., Ivanov, D., Poser, B. A., and Uludağ, K. (2017). Echo-time dependence of the BOLD response transients—A window into brain functional physiology. Neuroimage 159, 355–370. doi: 10.1016/j.neuroimage.2017.07.034

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoge, R. D., Atkinson, J., Gill, B., Crelier, G. R. R., Marrett, S., and Pike, G. B. (1999). Investigation of BOLD signal dependence on cerebral blood flow and oxygen consumption: the deoxyhemoglobin dilution model. Magn. Reson. Med. 42, 849–863. doi: 10.1002/(sici)1522-2594(199911)42:5<849::aid-mrm4>3.0.co;2-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Hutchison, J. L., Lu, H., and Cortex, B. R. C. (2013). Neural mechanisms of age-related slowing: the ΔCBF/ΔCMRO₂ ratio mediates age-differences in BOLD signal and human performance. Cereb. Cortex 23, 2337–2346. doi: 10.1093/cercor/bhs233

PubMed Abstract | CrossRef Full Text | Google Scholar

Jann, K., Hernandez, L. M., Beck-Pancer, D., McCarron, R., Smith, R. X., Dapretto, M., et al. (2015). Altered resting perfusion and functional connectivity of default mode network in youth with autism spectrum disorder. Brain Behav. 5:e00358. doi: 10.1002/brb3.358

PubMed Abstract | CrossRef Full Text | Google Scholar

Jo, H. J., Saad, Z. S., Simmons, W. K., Milbury, L. A., and Cox, R. W. (2010). Mapping sources of correlation in resting state FMRI, with artifact detection and removal. Neuroimage 52, 571–582. doi: 10.1016/j.neuroimage.2010.04.246

PubMed Abstract | CrossRef Full Text | Google Scholar

Kannurpatti, S. S., and Biswal, B. B. (2008). Detection and scaling of task-induced fMRI-BOLD response using resting state fluctuations. Neuroimage 40, 1567–1574. doi: 10.1016/j.neuroimage.2007.09.040

PubMed Abstract | CrossRef Full Text | Google Scholar

Kastrup, A., Krüger, G., Glover, G. H., and Moseley, M. E. (1999). Assessment of cerebral oxidative metabolism with breath holding and fMRI. Magn. Reson. Med. 42, 608–611. doi: 10.1002/(sici)1522-2594(199909)42:3<608::aid-mrm26>3.3.co;2-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Kida, I., Rothman, D. L., and Hyder, F. (2006). Dynamics of changes in blood flow, volume, and oxygenation: implications for dynamic functional magnetic resonance imaging calibration. J. Cereb. Blood Flow Metab. 27, 690–696. doi: 10.1038/sj.jcbfm.9600409

PubMed Abstract | CrossRef Full Text | Google Scholar

Koelsch, S., Schulze, K., Sammler, D., Fritz, T., Müller, K., and Gruber, O. (2009). Functional architecture of verbal and tonal working memory: an FMRI study. Hum. Brain Mapp. 30, 859–873. doi: 10.1002/hbm.20550

PubMed Abstract | CrossRef Full Text | Google Scholar

Kwong, K. K., Belliveau, J. W., Chesler, D. A., Goldberg, I. E., Weisskoff, R. M., Poncelet, B. P., et al. (1992). Dynamic magnetic resonance imaging of human brain activity during primary sensory stimulation. Proc. Natl. Acad. Sci. U S A 89, 5675–5679. doi: 10.1073/pnas.89.12.5675

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H.-J., Hou, X.-H., Liu, H.-H., Yue, C.-L., He, Y., and Zuo, X.-N. (2015). Toward systems neuroscience in mild cognitive impairment and Alzheimer’s disease: a meta-analysis of 75 fMRI studies. Hum. Brain Mapp. 36, 1217–1232. doi: 10.1002/hbm.22689

PubMed Abstract | CrossRef Full Text | Google Scholar

Liang, X., Connelly, A., and Calamante, F. (2014). Graph analysis of resting-state ASL perfusion MRI data: nonlinear correlations among CBF and network metrics. Neuroimage 87, 265–275. doi: 10.1016/j.neuroimage.2013.11.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, T. T. (2013). Neurovascular factors in resting-state functional MRI. Neuroimage 80, 339–348. doi: 10.1016/j.neuroimage.2013.04.071

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, T. T. (2016). Noise contributions to the fMRI signal: an overview. Neuroimage 143, 141–151. doi: 10.1016/j.neuroimage.2016.09.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, P., Hebrank, A. C., Rodrigue, K. M., Kennedy, K. M., Section, J., Park, D. C., et al. (2013). Age-related differences in memory-encoding fMRI responses after accounting for decline in vascular reactivity. Neuroimage 78, 415–425. doi: 10.1016/j.neuroimage.2013.04.053

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, T. T., and Wong, E. C. (2005). A signal processing model for arterial spin labeling functional MRI. Neuroimage 24, 207–215. doi: 10.1016/j.neuroimage.2004.09.047

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, H., Xu, F., Rodrigue, K. M., Kennedy, K. M., Cheng, Y., Flicker, B., et al. (2011). Alterations in cerebral metabolic rate and blood supply across the adult lifespan. Cereb. Cortex 21, 1426–1434. doi: 10.1093/cercor/bhQ234

PubMed Abstract | CrossRef Full Text | Google Scholar

Mayhew, S. D., Mullinger, K. J., Bagshaw, A. P., Bowtell, R., and Francis, S. T. (2014). Investigating intrinsic connectivity networks using simultaneous BOLD and CBF measurements. Neuroimage 99, 111–121. doi: 10.1016/j.neuroimage.2014.05.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Mohtasib, R. S., Lumley, G., Goodwin, J. A., Emsley, H. C. A., Sluming, V., and Parkes, L. M. (2012). Calibrated fMRI during a cognitive Stroop task reveals reduced metabolic response with increasing age. Neuroimage 59, 1143–1151. doi: 10.1016/j.neuroimage.2011.07.092

PubMed Abstract | CrossRef Full Text | Google Scholar

Moradi, F., Buračas, G. T., and Buxton, R. B. (2012). Attention strongly increases oxygen metabolic response to stimulus in primary visual cortex. Neuroimage 59, 601–607. doi: 10.1016/j.neuroimage.2011.07.078

PubMed Abstract | CrossRef Full Text | Google Scholar

Mullinger, K. J., Cherukara, M. T., Buxton, R. B., Francis, S. T., and Mayhew, S. D. (2017). Post-stimulus fMRI and EEG responses: evidence for a neuronal origin hypothesised to be inhibitory. Neuroimage 157, 388–399. doi: 10.1016/j.neuroimage.2017.06.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Mullinger, K. J., Mayhew, S. D., Bagshaw, A. P., Bowtell, R., and Francis, S. T. (2014). Evidence that the negative BOLD response is neuronal in origin: a simultaneous EEG-BOLD-CBF study in humans. Neuroimage 94, 263–274. doi: 10.1016/j.neuroimage.2014.02.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Murphy, K., Birn, R. M., and Bandettini, P. A. (2013). Resting-state fMRI confounds and cleanup. Neuroimage 80, 349–359. doi: 10.1016/j.neuroimage.2013.04.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Nichols, T. E., and Holmes, A. P. (2002). Nonparametric permutation tests for functional neuroimaging: a primer with examples. Hum. Brain Mapp. 15, 1–25. doi: 10.1002/hbm.1058

PubMed Abstract | CrossRef Full Text | Google Scholar

O’Rourke, M. F., and Hashimoto, J. (2007). Mechanical factors in arterial aging: a clinical perspective. J. Am. Coll. Cardiol. 50, 1–13. doi: 10.1016/j.jacc.2006.12.050

PubMed Abstract | CrossRef Full Text | Google Scholar

Obata, T., Liu, T. T., Miller, K. L., Luh, W. M., Wong, E. C., Frank, L. R., et al. (2004). Discrepancies between BOLD and flow dynamics in primary and supplementary motor areas: application of the balloon model to the interpretation of BOLD transients. Neuroimage 21, 144–153. doi: 10.1016/j.neuroimage.2003.08.040

PubMed Abstract | CrossRef Full Text | Google Scholar

Ogawa, S., Tank, D. W., Menon, R., Ellermann, J. M., Kim, S. G., Merkle, H., et al. (1992). Intrinsic signal changes accompanying sensory stimulation: functional brain mapping with magnetic resonance imaging. Proc. Natl. Acad. Sci. U S A 89, 5951–5955. doi: 10.1073/pnas.89.13.5951

PubMed Abstract | CrossRef Full Text | Google Scholar

Parkes, L. M., Rashid, W., Chard, D. T., and Tofts, P. S. (2004). Normal cerebral perfusion measurements using arterial spin labeling: reproducibility, stability, and age and gender effects. Magn. Reson. Med. 51, 736–743. doi: 10.1002/mrm.20023

PubMed Abstract | CrossRef Full Text | Google Scholar

Podlutsky, A., Ballabh, P., and Csiszar, A. (2010). Oxidative stress and endothelial dysfunction in pulmonary arteries of aged rats. Am. J. Physiol. Heart Circ. Physiol. 298, H346–H351. doi: 10.1152/ajpheart.00972.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

Power, J. D., Barnes, K. A., Snyder, A. Z., Schlaggar, B. L., and Petersen, S. E. (2012). Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. Neuroimage 59, 2142–2154. doi: 10.1016/j.neuroimage.2011.10.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Power, J. D., Mitra, A., Laumann, T. O., Snyder, A. Z., Schlaggar, B. L., and Petersen, S. E. (2014). Methods to detect, characterize and remove motion artifact in resting state fMRI. Neuroimage 84, 320–341. doi: 10.1016/j.neuroimage.2013.08.048

PubMed Abstract | CrossRef Full Text | Google Scholar

Restom, K., Behzadi, Y., and Liu, T. T. (2006). Physiological noise reduction for arterial spin labeling functional MRI. Neuroimage 31, 1104–1115. doi: 10.1016/j.neuroimage.2006.01.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Romero, L., Walsh, V., and Papagno, C. (2006). The neural correlates of phonological short-term memory: a repetitive transcranial magnetic stimulation study. J. Cogn. Neurosci. 18, 1147–1155. doi: 10.1162/jocn.2006.18.7.1147

PubMed Abstract | CrossRef Full Text | Google Scholar

Sadaghiani, S., Uğurbil, K., and Uludağ, K. (2009). Neural activity-induced modulation of BOLD poststimulus undershoot independent of the positive signal. Magn. Reson. Imaging 27, 1030–1038. doi: 10.1016/j.mri.2009.04.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Samanez-Larkin, G. R., and D’Esposito, M. (2008). Group comparisons: imaging the aging brain. Soc. Cogn. Affect. Neurosci. 3, 290–297. doi: 10.1093/scan/nsn029

PubMed Abstract | CrossRef Full Text | Google Scholar

Simon, A. B., and Buxton, R. B. (2015). Understanding the dynamic relationship between cerebral blood flow and the BOLD signal: implications for quantitative functional MRI. Neuroimage 116, 158–167. doi: 10.1016/j.neuroimage.2015.03.080

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, S. M. (2002). Fast robust automated brain extraction. Hum. Brain Mapp. 17, 143–155. doi: 10.1002/hbm.10062

CrossRef Full Text | Google Scholar

Storti, S. F., Boscolo Galazzo, I., Montemezzi, S., Menegaz, G., and Pizzini, F. B. (2017). Dual-echo ASL contributes to decrypting the link between functional connectivity and cerebral blow flow. Hum. Brain Mapp. 38, 5831–5844. doi: 10.1002/hbm.23804

PubMed Abstract | CrossRef Full Text | Google Scholar

Studholme, C., Hawkes, D. J., and Hill, D. L. (1998). “Normalized entropy measure for multimodality image alignment,” in SPIE The International Society for Optical Engineering, ed. K. M. Hanson (San Diego, CA: Medical Imaging), 132–143.

Google Scholar

Tak, S., Polimeni, J. R., Wang, D. J. J., Yan, L., and Chen, J. J. (2015). Associations of resting-state fmri functional connectivity with flow-BOLD coupling and regional vasculature. Brain Connect. 5, 137–146. doi: 10.1089/brain.2014.0299

PubMed Abstract | CrossRef Full Text | Google Scholar

Tak, S., Wang, D. J. J., Polimeni, J. R., Yan, L., and Chen, J. J. (2014). Dynamic and static contributions of the cerebrovasculature to the resting-state BOLD signal. Neuroimage 84, 672–680. doi: 10.1016/j.neuroimage.2013.09.057

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, C. H., Low, K. A., Kong, T., Fletcher, M. A., Zimmerman, B., Maclin, E. L., et al. (2017). Mapping cerebral pulse pressure and arterial compliance over the adult lifespan with optical imaging. PLoS One 12:e0171305. doi: 10.1371/journal.pone.0171305

PubMed Abstract | CrossRef Full Text | Google Scholar

Tarantini, S., Tran, C. H. T., Gordon, G. R., Ungvari, Z., and Csiszar, A. (2017). Impaired neurovascular coupling in aging and Alzheimer’s disease: contribution of astrocyte dysfunction and endothelial impairment to cognitive decline. Exp. Gerontol. 94, 52–58. doi: 10.1016/j.exger.2016.11.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Thakral, P. P., Madore, K. P., and Schacter, D. L. (2017). A role for the left angular gyrus in episodic simulation and memory. J. Neurosci. 37, 8142–8149. doi: 10.1523/JNEUROSCI.1319-17.2017

PubMed Abstract | CrossRef Full Text | Google Scholar

Toth, P., Tarantini, S., Csiszar, A., and Ungvari, Z. (2017). Functional vascular contributions to cognitive impairment and dementia: mechanisms and consequences of cerebral autoregulatory dysfunction, endothelial impairment and neurovascular uncoupling in aging. Am. J. Physiol. Heart Circ. Physiol. 312, H1–H20. doi: 10.1152/ajpheart.00581.2016

PubMed Abstract | CrossRef Full Text | Google Scholar

Trott, D. W., Seawright, J. W., Luttrell, M. J., and Woodman, C. R. (2011). NAD(P)H oxidase-derived reactive oxygen species contribute to age-related impairments of endothelium-dependent dilation in rat soleus feed arteries. J. Appl. Physiol. 110, 1171–1180. doi: 10.1152/japplphysiol.01037.2010

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsvetanov, K. A., Henson, R. N. A., Tyler, L. K., Davis, S. W., Shafto, M. A., Taylor, J. R., et al. (2015). The effect of ageing on fMRI: correction for the confounding effects of vascular reactivity evaluated by joint fMRI and MEG in 335 adults. Hum. Brain Mapp. 36, 2248–2269. doi: 10.1002/hbm.22768

PubMed Abstract | CrossRef Full Text | Google Scholar

Tustison, N. J., Avants, B. B., Cook, P. A., Zheng, Y., Egan, A., Yushkevich, P. A., et al. (2010). N4ITK: improved N3 bias correction. IEEE Trans. Med. Imaging 29, 1310–1320. doi: 10.1109/TMI.2010.2046908

PubMed Abstract | CrossRef Full Text | Google Scholar

van Osch, M. J. P., Teeuwisse, W. M., van Walderveen, M. A. A., Hendrikse, J., Kies, D. A., and van Buchem, M. A.. (2009). Can arterial spin labeling detect white matter perfusion signal? Magn. Reson. Med. 62, 165–173. doi: 10.1002/mrm.22002

CrossRef Full Text | Google Scholar

Vecchio, F., Miraglia, F., Curcio, G., Altavilla, R., Scrascia, F., Giambattistelli, F., et al. (2015). Cortical brain connectivity evaluated by graph theory in dementia: a correlation study between functional and structural data. J. Alzheimers Dis. 45, 745–756. doi: 10.3233/JAD-142484

PubMed Abstract | CrossRef Full Text | Google Scholar

Viviani, R., Messina, I., and Walter, M. (2011). Resting state functional connectivity in perfusion imaging: correlation maps with BOLD connectivity and resting state perfusion. PLoS One 6:e27050. doi: 10.1371/journal.pone.0027050

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Z., Aguirre, G. K., Rao, H., Wang, J., Fernández-Seara, M. A., Childress, A. R., et al. (2008). Empirical optimization of ASL data analysis using an ASL data processing toolbox: ASLtbx. Magn. Reson. Imaging 26, 261–269. doi: 10.1016/j.mri.2007.07.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Wise, R. G., Harris, A. D., Stone, A. J., and Murphy, K. (2013). Measurement of OEF and absolute CMRO₂: MRI-based methods using interleaved and combined hypercapnia and hyperoxia. Neuroimage 83, 135–147. doi: 10.1016/j.neuroimage.2013.06.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Wise, R. G., Ide, K., Poulin, M. J., and Tracey, I. (2004). Resting fluctuations in arterial carbon dioxide induce significant low frequency variations in BOLD signal. Neuroimage 21, 1652–1664. doi: 10.1016/j.neuroimage.2003.11.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Woods, R. P., Grafton, S. T., Holmes, C. J., Cherry, S. R., and Mazziotta, J. C. (1998). Automated image registration: I. general methods and intrasubject, intramodality validation. J. Comput. Assist. Tomogr. 22, 139–152. doi: 10.1097/00004728-199801000-00027

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, C. W., Gu, H., Lu, H., Stein, E. A., Chen, J.-H., and Yang, Y. (2009). Mapping functional connectivity based on synchronized CMRO₂ fluctuations during the resting state. Neuroimage 45, 694–701. doi: 10.1016/j.neuroimage.2008.12.066

PubMed Abstract | CrossRef Full Text | Google Scholar

Zappe, A. C., Pfeuffer, J., Merkle, H., Logothetis, N. K., and Goense, J. B. M. (2008). The effect of labeling parameters on perfusion-based fMRI in nonhuman primates. J. Cereb. Blood Flow Metab. 28, 640–652. doi: 10.1038/sj.jcbfm.9600564

PubMed Abstract | CrossRef Full Text | Google Scholar

Zou, Q., Wu, C. W., Stein, E. A., Zang, Y., and Yang, Y. (2009). Static and dynamic characteristics of cerebral blood flow during the resting state. Neuroimage 48, 515–524. doi: 10.1016/j.neuroimage.2009.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar